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
•
written
by
gsc529 •
1 answer
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:
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.
Log in to answer this question.
-
is there any way to normalize the counts and compare it for the splicing quantification?
written by PK •
-
input file in rmats
written by Lambodarswain316 •
-
How robust is rMATS when swapping samples?
written by Apex92 •
-
how to visualize correctly significant splicing events- inconsistencies between rMATS results and I…
written by Reuth •
-
novel AS events detection from RNA seq
written by Sara
-
alternative splicing: very different outputs from rMATS vs altanalyze
written by smurph50 •
-
Alternative Splicing Analysis
written by noahhelton98 •
-
Event-based software for alternative splicing with three conditions?
written by htyu1lmzta •
-
Alternatives to LifeScope for realigning SOLiD BAM files
written by roysomak4 •
-
How Can I Know The Length Of Mapped Reads From Bam File?
written by Jordan
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.