csfasta to fasta?

[brachysclereid]

Is anyone aware of a script that will convert a color space fasta to a fasta file?

[BAMseek]

One thing to keep in mind - you could convert the colors into bases, but an error in the color call would throw off all the subsequent bases after that position. It may be better to align the data in colorspace (converting the reference sequence into color space), and then finding the most parsimonious conversion of the color space reads into base space reads given the reference.

[brachysclereid]

csfasta to fasta

I would still like to do this. In theory those problems would result in kmers occurring only once or a few times in the data set and could be removed.


Does such a script exist?

[kevleb]

you should use encodeFasta.py script in corona from SOLiD.

[pmiguel]

I would still like to do this. In theory those problems would result in kmers occurring only once or a few times in the data set and could be removed.

Alas, no. Without knowing the base n-1, there are 4 base space sequences that may derive from a color space sequence starting at base n. One of those base space sequences is correct, the other 3 are incorrect. The result is that any error at any position sends one down the incorrect path. However, these will not create wildly divergent kmers because there are only 3 of them.

It might be possible to collapse the resulting "homoconvertamers" knowing their derivation. But it would likely be easier to just do all the analysis in color space (and take advantage of the extra error detection capabilities of color space) and do the conversion to base space as a final step.

--
Phillip

[westerman]

I definately agree with Phillip and BAMSeek as to the advisability of doing your work in colorspace instead of in basespace. As for your specific question, there have been several cs-to-bs converters posted on SeqAnswers. Gringer posted one recently (which I did not like). Other people have posted theirs. Search for forum for the keywords "basespace colorspace"