Question: Standard Practice of making VCF file from SAM file
0
gravatar for hafiz.talhamalik
18 days ago by
Pakistan
hafiz.talhamalik30 wrote:

Can any one guide me to about what standard GATK tools to use when making a vcf file from allignment sam file ? I'm kinda stuck. confused which one to use and when. ???

snp gatk • 147 views
ADD COMMENTlink modified 11 days ago by Biostar ♦♦ 20 • written 18 days ago by hafiz.talhamalik30

Where are you stuck? What confuses you? Please tell us how far you've gotten in your quest to solve this.

ADD REPLYlink written 18 days ago by RamRS19k

I have used PICARD for SortSam, and MarkDuplicates. Now which one to use next ?

ADD REPLYlink modified 18 days ago • written 18 days ago by hafiz.talhamalik30

GATK Best Practices says BaseRecalibration is next. Have you looked at that resource?

ADD REPLYlink written 18 days ago by RamRS19k

I am trying to follow this one. But BaseRecalibration requires input of -knownSites latest_dbsnp.vcf. and i have whole genome of bacteria so what should i do with this parameter ? and what about next step ?

ADD REPLYlink written 18 days ago by hafiz.talhamalik30
1

Ah, bacterial genome - I think you should have mentioned that in your original question, I'd have warned you against GATK - I think GATK might not be great for non-human genomes. Maybe try samtools/bcftools?

ADD REPLYlink written 18 days ago by RamRS19k

oh okz. I used samtools/bcftools pipeline. I made a vcf file too but i wanted to counter check that with GATK pipeline. one more thing I have a genome of around mb, and samtools pipeline gave me a vcf file with around 2300 variants. Is it normal ? or did I made a mistake ?

ADD REPLYlink written 18 days ago by hafiz.talhamalik30

If you were sequencing what was supposed to be the exact same strain of bacteria as your reference, that's kind of high. If it's not, that number might be appropriate. Eyeball some of the putative SNPs in IGV, see if you can verify visually that the reads show the SNPs in the vcf.

ADD REPLYlink written 18 days ago by swbarnes24.5k

okz. that helped. Thanks alot

ADD REPLYlink written 17 days ago by hafiz.talhamalik30
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1736 users visited in the last hour