This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Process Soapalign Results Before Call Snps And Indels By Samtools?

Hi. I want to use SOAP to align my data to references and then call SNPs and Indels by using Samtools. Problems is that there are two output files related to pair aligned reads and single aligned reads produced by soap. So how to ideal with this two files before I use samtools? I tried to just combine two files by cat command but I am not sure if it will influences the final results. So is there any suggestions? Thanks.

samtools
• 3,215 views • link
updated by lh3 lh3 • written by Haiping • Haiping

1 answer

There is a soap2sam.pl script in samtools, but I would suggest you change to another mapper that native supports SAM output. Recommended: bwa, novoalign or stampy.

• 0 views • link
written by lh3 lh3

Thanks. Actually I do hava novoalign resutls. I just want compare the results between soap and novoalign.

• 0 views • link
written by Haiping • Haiping

Novoalign is almost certainly better. No need to compare in fact.

• 0 views • link
written by lh3 lh3

Log in to answer this question.

More posts like this