Question: Calling variants help
1
gravatar for rquigl01
18 months ago by
rquigl0110
rquigl0110 wrote:

Hello Galaxy Biostar, For the past three weeks I've been trying to call SNPs and INDELs on my NGS samples for a summer research project. I'm a college student and all of this is new to me, but I've learned so much in these past three weeks. So far I have put my fastq files through bowtie2 on galaxy which created a sorted BAM file. I've ran QC on the fastq files and the BAM files and they have all came up in good results. I can open up the BAM files for my samples in IGV and scroll through the genome looking for mutations (I'm studying Arabidopsis Thaliana). Now for the part that I'm stuck at.... I can't find a program (that works) that will call all of the SNPs and INDELs for me. Either they don't call any INDELs at all or they miss many obvious SNPs that I can see on IGV. It's just getting very frustrating. I don't know how to run code, so doing anything outside of galaxy really isn't an option. Any help would be nice! Point me in the right direction! Is there someone who could do this step for me? Thank you galaxy community!

indels snps galaxy mpileup calling • 493 views
ADD COMMENTlink modified 18 months ago by bryantl90 • written 18 months ago by rquigl0110

What programs have you tried?

ADD REPLYlink written 18 months ago by bryantl90

mpileup and unified genotyper. Freebayes works but it doesnt open in IGV

ADD REPLYlink written 18 months ago by rquigl0110
1
gravatar for bryantl
18 months ago by
bryantl90
United States
bryantl90 wrote:

You will not find a program that will call ALL of the SNPs and INDELs. Each program has different biases and is prone to different types of errors. Here is a link to a paper that specifically looks at different pipelines for variant calling: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4671096/

Personally, I have found that freebayes works for my purposes, but the program that works best for you will be dependent on your dataset and what exactly you are looking for.

ADD COMMENTlink written 18 months ago by bryantl90

freebayes doesnt open in IGV. All it says is malformed header

ADD REPLYlink written 18 months ago by rquigl0110

mpileup and unified genotyper. Freebayes works but it doesnt open in IGV

ADD REPLYlink written 18 months ago by rquigl0110
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 16.09
Traffic: 64 users visited in the last hour