ABI Paired End Reads

[sparks]

Hi,
Does anyone have information on the new ABI Paired end reads, especially file formats and orientations of the reads.

We are developing a colour space version of Novoalign that's performing pretty well but we would like to support paired end as well as mate pair for first release and having some sample data would help a lo, let us know if you have some reads to share, we don't need many.

If anyone would like to try beta of NovoalignCS just ask.

Colin

[epigen]

NovoalignCS sounds great. I'm currently comparing the performance of BWA and BioScope on SOLiD transcriptome data (single end, unfortunately). Having Novoalign as a third option would be nice, especially since other people in my group are happy using it on their non-colorspace reads. I always stumble across all kinds of bugs anyway, so volunteering as a beta tester seems obvious.

[sparks]

hi epigen,

That would be great if you could help. Just email me at colin at novocraft dot com I'll send you necessary info.

Colin

[nilshomer]

I have found novoalignCS very sensitive and specific when the alignments are evaluated after variant calling (SAMtools). I would definitely recommend testing this alignment tool, although it may be a little slow (like the regular novoalign) if you have 10+ sequencers (talking to you Broad/Baylor/WashU/etc.). Otherwise the multi-threaded version is the way to go, with MPI if you have it installed.

One issue is that the # of reads from a SOLiD slide is in the hundreds of millions. If you don't have MPI installed, and you still want to align the many reads from a SOLiD run in a parallel fashion, I have created a script to convert and split the reads after using the "solid2fastq" utility in BFAST. It's a perl script, but is renamed with a .txt extension to fool the HTTP server:

http://genome.ucla.edu/~nhomer/novoalignCS/bfast2novo.txt

You can use this in a piped fashion:

solid2fastq [options] <csfastas> <quals> | bfast2novo.pl - <out.prefix> <split num>

Hopefully this will help those with idle machines fill them up with novoalignCS.

[Rao]

data set is available on solidsoftwarecommunity.com

[KevinLam]

Hmmm why not approach ABI for the sample data sets?
with SOLiD4 being the new standard you might get different results if u used old datasets

[sparks]

Thanks for help and suggestions. e managed to get a few files of Colour space paired end reads and latest versions of NovoalignCS now fully support paired end and mate pair reads.

Colin