SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > SOLiD



Similar Threads
Thread Thread Starter Forum Replies Last Post
Error using tophat mapping AB SOLiD color space data Luyi Tian Bioinformatics 7 03-24-2016 06:36 AM
The optimal options for TopHat in mapping 35bp SOLiD reads holywoool Bioinformatics 0 09-22-2011 07:10 PM
Problems with mapping SOLID color space sequences to hg19 using TopHat davidehs Bioinformatics 1 06-24-2011 02:21 PM
Mapping of SOLiD reads JohanS SOLiD 2 05-26-2011 03:04 AM
Poor R3 (35bp) mapping boscoliveira SOLiD 8 03-04-2011 08:09 PM

Reply
 
Thread Tools
Old 11-29-2016, 02:35 AM   #1
Ilarius
Junior Member
 
Location: Barcelona

Join Date: Nov 2016
Posts: 2
Default mapping 35bp solid reads with tophat unsuccessful

Hi, I recently downloaded some solid short reads (35 bp), single end, and convert them in cfasta with ahi-dump.

I found that there are very few mappers that support solid read mapping in color space and finally i decided to use tophat using bowtie1 (which supported color space mapping).
Unfortunately I less less than 30 % of the reads mapping.
My previous experience has always been with 100 bp paired end reads, so I am not completely sure if I am using the right parameters for short read mapping:

tophat -G $annotation --segment-length 17 --segment-mismatches 1 -g 1 --bowtie1 -p 12 -o $OUTPUT_FOLDER --coverage-search --color --quals $genome $READ $QUAL

maybe is the segment-mismatch and maximum multi-hits too restrictive? (-g is 20 by default). What for your experience are the best parameters for short read mapping?

Thanks
Ilario
Ilarius is offline   Reply With Quote
Old 11-30-2016, 03:43 AM   #2
colindaven
Senior Member
 
Location: Germany

Join Date: Oct 2008
Posts: 391
Default

In general, in comparison to Illumina PE reads don't expect high mapping percentages. From memory I think 50-70% was pretty good, as the data were not extensively filtered (on the 5500xl machine at least).

You probably don't want to use Tophat for such short reads as I don't think it will be able to find anchors/seeds successfully.

My best experiences with SOLiD were with the commercial program NovoalignCS. You might be able to get a free trial.

Otherwise, Shrimp2 was also not too bad.

Why not use just bowtie1 instead of tophat2 with bowtie1 ?
colindaven is offline   Reply With Quote
Old 12-06-2016, 03:54 AM   #3
Ilarius
Junior Member
 
Location: Barcelona

Join Date: Nov 2016
Posts: 2
Default

I did use tophat2 with bowtie1, anyway if I let the default -g option (i think it is 20) it takes ages to finish.

As regards shrimp I didn't find any working link to download it!
Ilarius is offline   Reply With Quote
Old 08-04-2017, 03:47 AM   #4
colindaven
Senior Member
 
Location: Germany

Join Date: Oct 2008
Posts: 391
Default

I recently found another non-commercial aligner for SOLiD data - it is called subread. I have only tried for Illumina data but it's very good there.

Again, I wouldn't use tophat for such short read data.
colindaven is offline   Reply With Quote
Old 08-06-2017, 08:44 PM   #5
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 794
Default

Quote:
What for your experience are the best parameters for short read mapping?
Anything that gets you quicker to re-analysing your target sample with a non-SOLiD system. You shouldn't be mapping SOLiD reads at all, especially not 35bp ones; interpreting the results is just too confusing (even for the people writing code to do the mapping).
gringer is offline   Reply With Quote
Reply

Tags
parameter tuning, short read alignment, solid

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:04 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO