Bug Bulletin: The recent 3.2 release fixes many issues. If you run into a problem, please try the latest version before posting a bug report, as your problem may already have been solved.

Indel calling in pooled samples

sashasasha Posts: 8Member

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 Posts: 8Member

    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 Posts: 8Member
    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

    Post edited by sasha on
Sign In or Register to comment.