SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > SOLiD



Similar Threads
Thread Thread Starter Forum Replies Last Post
Mapping SOLiD colorspace paired end reads NestorNotabilis SOLiD 10 12-12-2012 07:14 PM
SOLiD paired-end mapping yaira04 SOLiD 6 12-12-2012 05:41 PM
SOLiD paired-end read analysis with BFAST 0.7.0a idonaldson Bioinformatics 4 03-20-2012 09:31 PM
BFAST for SOLiD paired end reads epigen Bioinformatics 31 09-03-2011 06:20 AM
BFAST indexing for 50+36 SOLiD paired end? krobison Bioinformatics 6 11-02-2010 06:25 AM

Reply
 
Thread Tools
Old 02-07-2012, 11:12 AM   #1
AmitL
Member
 
Location: Israel, Jerusalem

Join Date: Aug 2011
Posts: 35
Default Understanding paired-end SOLiD

Hello everyone,
I am working with paired-end SOLiD data, and I am not sure how to handle it.

1) What exactly is the difference between the F3 and F5 reads? Is the F5 reversed?

2) I'm using BWA for the analysis. How do I use it with paired-end? In the documentation it says 'file 1' and 'file 2', but I don't know which is which.
I have tried to run it both ways and in both times I got very few variants - 3000, 5000 while expecting 50000.
I have a color-space reference index, and I used the -c parameter in the alignment. Am I missing something?

3) Does the csfasta conversion work differently with paired-end data?

Thanks a lot,
Amit

=]
AmitL is offline   Reply With Quote
Old 02-08-2012, 12:27 AM   #2
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 838
Default

Quote:
1) What exactly is the difference between the F3 and F5 reads? Is the F5 reversed?
They are strand-specific reads. For the same RNA sequence, they will map in different directions to each other.

I think the 5' read is reverse and done first, so if you're doing an analysis taking reads in the usual order you need to specify RF orientation in the mapping program (I'm not very sure on this point).

Quote:
3) Does the csfasta conversion work differently with paired-end data?
Nope. You shouldn't do this conversion with either. Map in colour-space and let the mapping tool give you the equivalent base-space sequence.
gringer is offline   Reply With Quote
Old 02-09-2012, 05:48 AM   #3
AmitL
Member
 
Location: Israel, Jerusalem

Join Date: Aug 2011
Posts: 35
Default

Thank you very much!

Do you run such analyses with BWA?
AmitL is offline   Reply With Quote
Old 02-09-2012, 06:05 AM   #4
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 838
Default

No, I haven't [yet] bothered with looking at another program. While I might get a higher proportion of mapped reads, I've got enough problems working with error in the reads I have from Bowtie.
gringer is offline   Reply With Quote
Old 02-09-2012, 06:39 AM   #5
AmitL
Member
 
Location: Israel, Jerusalem

Join Date: Aug 2011
Posts: 35
Default

But... Bowtie does not handle color-space... or does it?
AmitL is offline   Reply With Quote
Old 02-09-2012, 07:09 AM   #6
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 838
Default

Bowtie can do color-space, Bowtie2 cannot (and is not likely to in the future).
gringer is offline   Reply With Quote
Old 02-16-2012, 06:20 AM   #7
AmitL
Member
 
Location: Israel, Jerusalem

Join Date: Aug 2011
Posts: 35
Default

~bump~
Does anyone work with BWA on SOLiD paired-end data?
AmitL is offline   Reply With Quote
Old 02-16-2012, 06:23 PM   #8
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

BWA does not support paired end reads with SOLiD.
nilshomer is offline   Reply With Quote
Old 02-17-2012, 01:48 AM   #9
AmitL
Member
 
Location: Israel, Jerusalem

Join Date: Aug 2011
Posts: 35
Default

I know the new version does not work with color-space.
Are you saying that it's true for the older versions too?
AmitL is offline   Reply With Quote
Old 02-17-2012, 08:05 AM   #10
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

The new version does not work with color space. The older version works with mate pairs and color space, but not paired end and color space.
nilshomer is offline   Reply With Quote
Old 02-20-2012, 07:30 AM   #11
AmitL
Member
 
Location: Israel, Jerusalem

Join Date: Aug 2011
Posts: 35
Default

How do you know it's meant for mate-pairs?

On the documentation it says:
Quote:
Generate alignments in the SAM format given paired-end reads.
AmitL is offline   Reply With Quote
Old 02-20-2012, 07:36 AM   #12
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Quote:
Originally Posted by AmitL View Post
How do you know it's meant for mate-pairs?

On the documentation it says:
Check out the source code "bwape.c" and search for SOLID. It uses the F3/R3 "pairing", which is mate pair.
nilshomer is offline   Reply With Quote
Old 02-20-2012, 10:55 AM   #13
AmitL
Member
 
Location: Israel, Jerusalem

Join Date: Aug 2011
Posts: 35
Default

Alright. Thanks a lot!
AmitL is offline   Reply With Quote
Old 04-18-2012, 12:07 AM   #14
yksikaksi
Member
 
Location: Scandinavia

Join Date: Dec 2009
Posts: 20
Default

Have you consider to map your csfasta reads with SHRiMP2?
yksikaksi is offline   Reply With Quote
Old 04-18-2012, 04:02 AM   #15
AmitL
Member
 
Location: Israel, Jerusalem

Join Date: Aug 2011
Posts: 35
Default

Yes.
Shrimp2 collapses because of non-existing input problems. I have checked the input manually and it was ok. It works with all other tools.

I preferred not to deal with such a buggy tool.

But thanks!
AmitL is offline   Reply With Quote
Old 04-18-2012, 04:29 AM   #16
yksikaksi
Member
 
Location: Scandinavia

Join Date: Dec 2009
Posts: 20
Default

Quote:
Originally Posted by AmitL View Post
Yes.
Shrimp2 collapses because of non-existing input problems. I have checked the input manually and it was ok. It works with all other tools.

I preferred not to deal with such a buggy tool.

But thanks!
Hum... it is strange. I managed to map my csfasta reads with SHRiMP2 but not BWA.
yksikaksi is offline   Reply With Quote
Old 04-18-2012, 12:17 PM   #17
golharam
Member
 
Location: Philadelphia, PA

Join Date: Dec 2009
Posts: 55
Default

Just out of curiosity, why are you not using LifeScope?
golharam is offline   Reply With Quote
Old 04-19-2012, 02:50 AM   #18
AmitL
Member
 
Location: Israel, Jerusalem

Join Date: Aug 2011
Posts: 35
Default

Did you try it?
Its interface is very uncomfortable and it cannot be fully operated through shell scripts.
AmitL is offline   Reply With Quote
Old 04-19-2012, 07:58 AM   #19
golharam
Member
 
Location: Philadelphia, PA

Join Date: Dec 2009
Posts: 55
Default

I agree. I only use LifeScope to map data, and I use the command-line interface. The GUI is horrible. Once I have the BAM files I do the rest of my processing outside of LifeScope.

I ask because I keep hearing LifeScope does a better job mapping reads than other tools. I haven't tried others myself though.
golharam is offline   Reply With Quote
Old 04-19-2012, 09:07 AM   #20
AmitL
Member
 
Location: Israel, Jerusalem

Join Date: Aug 2011
Posts: 35
Default

If you say so, I might give it a second chance.

Thank you
AmitL is offline   Reply With Quote
Reply

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 12:32 PM.


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