Understanding paired-end SOLiD

[nilshomer]

BWA does not support paired end reads with SOLiD.

[AmitL]

~bump~
Does anyone work with BWA on SOLiD paired-end data?

[gringer]

Bowtie can do color-space, Bowtie2 cannot (and is not likely to in the future).

[AmitL]

But... Bowtie does not handle color-space... or does it?

[gringer]

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.

[AmitL]

Thank you very much!

Do you run such analyses with BWA?

[gringer]

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).

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.

[AmitL]

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

=]