[Bioperl-l] trans-spliced genes & gff2 & chado
Charles Hauser
chauser at duke.edu
Mon Jul 21 14:00:16 EDT 2003
I am generating gff2 to run GBrowse and have a trans-spliced gene to
represent. I can write code to deal with this particular instance, but
was wondering if there is a generic solution?
- is there a 'standard/accepted format' to display trans-spliced genes?
- in chado I believe one would generate a feature for each segment -
correct?
Charles
gene join(32737..32824,complement(174205..174384),
complement(69520..71506))
/gene="psaA"
/note="trans splicing"
CDS join(32737..32825,complement(174205..174384),
complement(69520..71506))
while( my $seq = $seqio->next_seq ) {
foreach my $f ( $seq->top_SeqFeatures() ) {
<snip>
$out->write_feature($f);
}
}
More information about the Bioperl-l
mailing list