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

Hello,

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 • 238 views
ADD COMMENTlink modified 9 months ago by colindaven2.5k • written 9 months ago by gsc5290
1

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 9 months ago by _r_am31k
0
gravatar for colindaven
9 months ago by
colindaven2.5k
Hannover Medical School
colindaven2.5k 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>
Options:
ADD COMMENTlink written 9 months ago by colindaven2.5k

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 9 months ago by genomax92k

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 9 months ago by colindaven2.5k
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: 1299 users visited in the last hour