Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Viral Genome SNP and Indel Detection
0
gravatar for anjali.gopal91
3.6 years ago by
anjali.gopal91 • 50
United States
anjali.gopal91 • 50 wrote:

Hi all,

I'm trying to call snps and indels for HIV, and I'm wondering if there are suggestions for good mutation detection software for viruses. I expect to have fairly large indels (20-30bp in a 9kb genome), as well as the regular snps, across multiple  samples that are pooled together in one run.

I'm probably going to try SNVer (http://snver.sourceforge.net/), LoFreq (http://sourceforge.net/p/lofreq/wiki/Home/), or Breseq (http://barricklab.org/twiki/bin/view/Lab/ToolsBacterialGenomeResequencing), but I'm also wondering if there are any other software packages that would specifically be good at detecting indels. 

Are there any recommendations for software? 

Thanks!
sequencing snp alignment • 1.3k views
ADD COMMENT • link •
modified 3.6 years ago by Brian Bushnell ♦ 16k • written 3.6 years ago by anjali.gopal91 • 50
0
gravatar for Brian Bushnell
3.6 years ago by
Brian Bushnell ♦ 16k
Walnut Creek, USA
Brian Bushnell ♦ 16k wrote:

If you are interested in detecting long indels, I highly recommend using BBMap for alignment, as it is the most accurate for detecting such mutations in my tests.  Samtools mpileup and GATK are both capable of calling variants from sam files, and are widely used.  I have also heard good things about FreeBayes, but note that it does not support the current SAM specification (which is already several years old), so you will need to add the "sam=1.3" flag to BBMap's command line to generate old-style cigar strings that FreeBayes can process.

 
ADD COMMENT • link written 3.6 years ago by Brian Bushnell ♦ 16k
Please log in to add an answer.

Similar posts • Search »

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