Question: SOLiD paired end reads mapping problems
1
gravatar for MiguelMorard
4.2 years ago by
MiguelMorard20
Spain,Valencia
MiguelMorard20 wrote:

Hello,

I have Solid paired end reads (F3, F5) from DNA-Seq I want to map on a reference genome (yeast). I used Bowtie with the -C parameter and then I tried different values for the -v but I can't map more than 51% of the reads. 

I would like to know if somebody uses other alignment tools with colorspace data that could perform better than bowtie.

Thanks a lot

dna bowtie sequence alignment solid • 1.6k views
ADD COMMENTlink modified 4.2 years ago by brentp22k • written 4.2 years ago by MiguelMorard20
2
gravatar for brentp
4.2 years ago by
brentp22k
Salt Lake City, UT
brentp22k wrote:

That's solid data for you. Though 50% is pretty low. How long are the reads?

Back when I was subject to such horrible things, I did a comparison of various aligners on colorspace data https://github.com/brentp/bowfast/tree/master/aligner-compare

 

Bowtie should work well with the options here: https://github.com/brentp/bowfast/blob/master/aligner-compare/run-bowtie.sh

You can see that other aligners do better but they are more of a pain to deal with.

  • bfast has huuuuge indexes
  • shrimp requires splitting the data (but you can follow my run-shrimp.sh script
  • novoalign isnt' free

sympathies.

 

ADD COMMENTlink written 4.2 years ago by brentp22k

The F3 are around 75bp and the F5 are around 35bp. I will try your parameters. I also found this http://cushaw3.sourceforge.net/homepage.htm#latest yesterday but I haven't had time to test it.

 

ADD REPLYlink written 4.2 years ago by MiguelMorard20
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: 2411 users visited in the last hour