Question: Methods to convert SOLiD reads to nucleotide space
gravatar for gsc529
13 months ago by
gsc5290 wrote:


I would like to know if there are any methods that can be used to convert colorspace reads from BAM files generated from colorspace reads on the SOLiD system to nucleotide space?

I am looking to use analyze some data generated a number of years ago using both SOLiD 3 and SOLiD 4 systems for alternative splicing. I currently have the aligned BAM files that were generated using LifeScope software but I do not, however, have the csfastq's. I am hoping to run the data through rMATS, or a similar package, to look at alternative splicing events in two different study groups (case vs control). Are there any such things similar to rMATS that would be better suited to handle colorspace?

Any suggestions/guidance would be greatly appreciated! Thank you

sequencing rna-seq assembly • 269 views
ADD COMMENTlink modified 13 months ago by colindaven2.6k • written 13 months ago by gsc5290

Titles should be brief and clear. Titles convey what the post is about and should not contain text that goes in the body of the post. In your case, "Methods to convert SOLiD reads to nucleotide space" would have been a much better title than this wall of text.

ADD REPLYlink written 13 months ago by Ram32k
gravatar for colindaven
13 months ago by
Hannover Medical School
colindaven2.6k wrote:

Because the reads were aligned against a reference, you can just use BAM 2 fastq methods to extract the sequences. That way you get a FASTQ and avoid colourspace altogether, which I would not recommend using for RNA-seq projects in any case.

You can use a modern samtools for this:

samtools fastq

No input file specified.
Usage: samtools fastq [options...] <in.bam>
ADD COMMENTlink written 13 months ago by colindaven2.6k

Do you know if this works with SOLiD data and generate correct results?

In this day and age there is really no good reason to work with colorspace data. OP should ideally find a new/more current dataset.

ADD REPLYlink written 13 months ago by GenoMax96k

Yep, it works, yep, they should use other data - but that wasn't the question.

I'd probably trust it for alignments, but not for SNV calls as you're not using the colour space information.

ADD REPLYlink written 13 months ago by colindaven2.6k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1937 users visited in the last hour