[Bioperl-l] Bio::SeqIO::scf traces scrambled?

Chris Fields cjfields at illinois.edu
Thu Jun 18 20:27:02 UTC 2009


On Jun 18, 2009, at 3:02 PM, Charles Tilford wrote:

> Cook, Malcolm wrote:
>> Charles,
>>
>> Another possible stopgap that might work for you, if you're working  
>> with AB1 chromatograms and have ABIs kb-basecaller turned on, is to  
>> use Bio::Trace::ABIF
>>
>> 	http://search.cpan.org/dist/Bio-Trace-ABIF/lib/Bio/Trace/ABIF.pm
>>
>> It works great and includes implementation of ABIs algorithm  
>> allowing to (re)compute trace clear ranges using kc-basecallers  
>> quality scores and any windowing/quality parameters.
>>
>> Its not in the bioperl project but it is an easy install from CPAN.
>>
> Thanks - we installed that a few weeks ago, and it was on my list of  
> things to try, but I had not gotten to it yet since I was getting  
> data out of the SCF SeqIO module. Even though the SeqIO::scf data  
> looks ok, the fact that I need to unscramble it makes me nervous...  
> Thanks, too, for the example code. I'll try out the Bio::Trace::ABIF  
> module and see if it works with our files.
>
> Thanks,
> CAT

You definitely shouldn't need to unscramble it; my guess is this is a  
legit bug that just has gone unnoticed.  I see that you have filed a  
ticket on it so we can at least track it.  Thanks!

chris



More information about the Bioperl-l mailing list