Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.We appreciate your help!

Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
Attention:
We will be out of the office on November 11th and 13th 2019, due to the U.S. holiday(Veteran's day) and due to a team event(Nov 13th). We will return to monitoring the GATK forum on November 12th and 14th respectively. Thank you for your patience.

Getting very less entries in vcf file??

Dear Developers,

I Run Unified genotyper for identify the indels and deletion.
I used Realigner target creater got forIndelRealigner.intervals. used this output for Indel Realigner and get the output bam file which i gave as input toUnified Genotyper.
I got result as

CHROM POS ID REF ALT QUAL FILTER INFO FORMAT BRCA1

gi|262359905|ref|NG_005905.2| 48155 . G A 387.77 . AC=1;AF=0.500;AN=2;BaseQRankSum=-0.466;DP=47;Dels=0.00;FS=2.721;HaplotypeScore=0.0000;MLEAC=1;MLEAF=0.500;MQ=29.31;MQ0=0;MQRankSum=-4.812;QD=8.25;ReadPosRankSum=0.535 GT:AD:DP:GQ:PL 0/1:25,22:45:99:416,0,676

But if i run without these steps (Realigner target creater, Indel Realigner) i am getting the same output.
Also i used the same dataset and run it in samtools and got more SNP.

I am looking for a deletion in the dataset.
I would like to get the result from GATK Could you please suggest on the same.

Thanks
Sridhar

Best Answer

Answers

  • sridhar28sridhar28 Member

    Hello Geraldine,

    Thank you for the detailed reply.
    Actually in the analysis we tried with a software(not samtools) for variant calling and compared the results with GATK here we dint get the deletion we observed.
    I tried Haplotypecaller earlier i got few entries more but the deletion which i expect i dint get it..

    Thanks
    Sridhar

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    I see -- you should definitely try the genotype given alleles mode, then. That will at least give you some metrics on what the caller thinks of the deletion you observe.

  • sridhar28sridhar28 Member

    Hello Geraldine,

    Cold you please tell me how to use the option GENOTYPE_GIVEN_ALLELES ??
    Because it shows error wehn i use UnifiedGenotyper and haplotypecaller..

    Thanks
    Sridhar

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    Can you tell me what command line you tried and what error you got?

  • sridhar28sridhar28 Member

    java -jar galaxy/apps/GenomeAnalysisTK-2.5-2-gf57256b/GenomeAnalysisTK.jar -R brca.fasta -T UnifiedGenotyper -I BRCA_S21_picard.bam --GENOTYPE_GIVEN_ALLELES

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    Oh, I see. You need to provide a VCF file containing the alleles you are looking for, with the --alleles argument as described here. You can also provide the same VCF to the -L argument if you want to restrict your analysis to just those sites, but that is optional.

Sign In or Register to comment.