[Bioperl-l] parsing ACE file

Florent Angly florent.angly at gmail.com
Thu Jan 24 22:06:29 UTC 2008


That would be the module Bio::Assembly::IO::ace
It works fine as far as I know.
To parse an assembly, use Bio::Assembly::IO: 
http://doc.bioperl.org/bioperl-live/Bio/Assembly/IO.html
Regards,
Florent

Phillip San Miguel wrote:
> Stefano Ghignone wrote:
>> Dear All,
>>     dealing with an assembly .ace file and a list of contigs (from 
>> that assembly), how can I extract from the .ace file the read names 
>> forming each listed contig? Is there any module doing this job?
>>
>> Any suggestion about how to start is welcome...
>> Cheers
>>
>> Stefano
>>
>>   
> perl -ne 'next unless (/^(?:CO)|(?:RD)/);print' acefile.ace
>
> will give you a list of each the contigs followed by the reads in each 
> contig, if "acefile.ace" is a phrap ace file.
>
> There is a bioperl module for handling phrap ace file, but I'm not 
> sure what its current status is. Last time I looked (probably a couple 
> of years ago) it seemed to have been abandoned half-finished.
>




More information about the Bioperl-l mailing list