[gt-users] gff3 parser
Gordon Gremme
gremme at gmail.com
Fri Feb 13 13:14:06 CET 2009
> ok. i think i've resolved this:
> http://gist.github.com/62770
>
> for a.gff of:
> ##gff-version 3
> 1 ucb gene 2 2 . - . ID=geneA
> 1 ucb gene 2 2 . + . ID=geneB
>
> and b.gff of:
> gff-version 3
> 1 ucb gene 4 4 . - . ID=geneA
>
> it gives:
>
> $ bin/gt gff3 -retainids a.gff b.gff 2>/dev/null
> ##gff-version 3
> ##sequence-region 1 2 2
> 1 ucb gene 2 2 . - . ID=geneA
> ###
> 1 ucb gene 2 2 . + . ID=geneB
> ###
> ##sequence-region 1 4 4
> 1 ucb gene 4 4 . - . ID=geneA.1
> ###
>
>
>
> and if i add another c.gff that's same as b.gff, it creates geneA.2
Ok, great! Could you add a test to testsuite/gt_gff3_include.rb to
test the new behaviour and make sure your code doesn't break the old
behaviour? It currently breaks some tests, just call testsuite.rb
-keywords gt_gff3 to check for yourself.
Gordon
More information about the gt-users
mailing list