The current GATK version is 3.7-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Get notifications!


You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

Got a problem?


1. Search using the upper-right search box, e.g. using the error message.
2. Try the latest version of tools.
3. Include tool and Java versions.
4. Tell us whether you are following GATK Best Practices.
5. Include relevant details, e.g. platform, DNA- or RNA-Seq, WES (+capture kit) or WGS (PCR-free or PCR+), paired- or single-end, read length, expected average coverage, somatic data, etc.
6. For tool errors, include the error stacktrace as well as the exact command.
7. For format issues, include the result of running ValidateSamFile for BAMs or ValidateVariants for VCFs.
8. For weird results, include an illustrative example, e.g. attach IGV screenshots according to Article#5484.
9. For a seeming variant that is uncalled, include results of following Article#1235.

Did we ask for a bug report?


Then follow instructions in Article#1894.

Formatting tip!


Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block as demonstrated here.

Jump to another community
Picard 2.10.4 has MAJOR CHANGES that impact throughput of pipelines. Default compression is now 1 instead of 5, and Picard now handles compressed data with the Intel Deflator/Inflator instead of JDK.
GATK version 4.beta.2 (i.e. the second beta release) is out. See the GATK4 BETA page for download and details.

Indel calling in pooled samples

Hi,
I am calling indel in pooled samples using this command:
java -jar -Xmx2g /PATH/2.1.13/GenomeAnalysisTK.jar -l INFO -T UnifiedGenotyper -I pool1.bam -I pool2.bam --out INDEL.vcf -R /reference.fa -glm INDEL

Currently i donot have any information of already known indels.
1.Do i need to first realign (RealignerTargetCreator and IndelRealigner) and then call indels even for pooled data?
2. How different will this be for calling indel on individual sample?

Looking forward for your suggesions.
with thanks sasha

Best Answers

Answers

  • sashasasha Member

    Thank you it did answer my question. In the second part i was just wondering if we have individual samples (instead of pool) then will it differ from calling indels? But i think in both cases we need to do the local realignment process (as you mentioned above).

  • sashasasha Member
    edited April 2013

    One more question related to this topic, Now that i have realigned (RealignerTargetCreator and IndelRealigner) and finally got a realigned.bam file by using this command:
    java -jar -Xmx2g /PATH/2.1.13/GenomeAnalysisTK.jar -l INFO -T IndelRealigner -I pool1.bam -I pool2.bam -I pool3.bam ... -L chr1 -targetIntervals pools_chr1.intervals -R reference.fa -o pools_realignedBam_chr1.bam

    So now i should use this pools_realignedBam_chr1.bam to call indels using:
    java -jar -Xmx2g /PATH/2.1.13/GenomeAnalysisTK.jar -l INFO -T UnifiedGenotyper -I pools_realignedBam_chr1.bam --out INDEL.vcf -R reference.fa -glm INDEL

    My question is that should i realign all pool1.bam pool2.bam … seperately? or its the right way by first realigning all pools together and then passing a single realigned file to unified genotyper??

    Looking forward for your guidance.
    Sasha

Sign In or Register to comment.