[BioC] [Bioc-sig-seq] SOLiD color space data

Steve Lianoglou mailinglist.honeypot at gmail.com
Sun Dec 12 16:57:25 CET 2010


Hi Neel,,

On Sat, Dec 11, 2010 at 9:26 PM, Neel Aluru <naluru at whoi.edu> wrote:
> Thanks, Martin. I also have read and heard about doing everything in color space. But, all the short read aligners convert the csfasta files to fastq before mapping to the genome.

I don't actually think this is true. I've just had to deal with some
SOLiD data, and bowtie does everything in colorspace (you can give it
the csfasta and qual files separately -- look for their -Q flag). Note
that you have to align to a special color space index.

I'd also be very surprised if BWA doesn't handle color space reads,
since (i) you can specify building a color space specific indiex (in
the bwa index ... command); and (ii) their paper says they handle
color space reads. I think the fastq file you would use as input to
BWA would just have the sequences in color space.

The SHRiMP aligner also does colorspace ...

-- 
Steve Lianoglou
Graduate Student: Computational Systems Biology
 | Memorial Sloan-Kettering Cancer Center
 | Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact



More information about the Bioconductor mailing list