[Bioperl-l] [BioSQL-l] How is is_circular recorded in BioSQL (by BioPerl)?
roy.chaudhuri at gmail.com
Mon Jul 25 16:25:28 UTC 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