Originally posted by unique379
I am confused with your comment regarding getting fastq format in decode format? There are some serious conversion issues if one directly converts color space to base space. Much better if you leave everything in color space, and this means you will have to index the reference in colorspace (using the -C parameter in bowtie-build).
Did you know you can use bowtie to map with existing csfasta + qual data files? Here is an example of a command line we have previously used:
bowtie -C -f --integer-quals -l 21 -n 1 --nomaqround --chunkmbs 256 --best --col-cqual --col-keepends --sam --quals input.qual bowtie_ref_index input.csfasta output.sam
You (and others) might be interested in a new software tool I have developed to help analyse miRNA data. Visit miRspring.victorchang.edu.au for detailed information. The website describes the miRspring document as well as providing some basic information on preparing BAM files.
Good luck! Would be interested to hear your findings when comparing to lifescope output.
D
Leave a comment: