How To Process Soapalign Results Before Call Snps And Indels By Samtools?
1
0
Entering edit mode
13.0 years ago
Haiping ▴ 110

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 • 2.6k views
ADD COMMENT
2
Entering edit mode
13.0 years ago
lh3 33k

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.

ADD COMMENT
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

Haha.OK. Thanks.

ADD REPLY

Login before adding your answer.

Traffic: 3633 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6