[Bioperl-l] [BioSQL-l] How is is_circular recorded in BioSQL (by BioPerl)?

Roy Chaudhuri roy.chaudhuri at gmail.com
Mon Jul 25 12:25:28 EDT 2011


>> Actually, when round-tripping a GenBank file through BioSQL, the LOCUS line
>> molecule type ends up in lower case, which makes me wonder if it is coming
>> from alphabet in the biosequence table.
>
> If so, that may break for viral GenBank files where the LOCUS line may say
> RNA, but the sequence is given using acgt (i.e. the DNA alphabet).

Just tried it and that does seem to be the case. It's not the only thing 
that breaks on round tripping, for example circular genomes become linear.



More information about the Bioperl-l mailing list