pjotrp / bioruby-gff3-plugin Goto Github PK

Hi Pjotr,
thanks for making this happen.
I may have encountered a bug. In the following gff, my genes of interest are on the reverse strand (my genes of interest are those with "Vg" in the name)
http://fourmidable.unil.ch/temp/Si_gnF%2Escaffold10535-splittingyw.gff
(the GFF was first built by MAKER, and then manually edited using Apollo)
(my connection is slow here so the upload may take another 10 mins or so)

My sequences come out looking like this:

Si_Vg3-RA Sequence:Si_gnF.scaffold10535_2342963:2350091 (2342963:2350091)
TACACCAAGGGACAGCAGGAAGAGAATGAACCACTCATGTGTAGTCGAAACGCTTATCATTGTCAAAAATGCGCTTTTAAAGTTCATCAGAAGTTTATTAGAACATGAAATTAATTTGAAAGCTTTATTTTAACAAACTGAACCATATTTTAGCTAATTATATGTAGATAATATTAACAACGTTTTAGACTTAAAGTGTAAAACGTGTATGAATAATTTAATTTAGTTCGGTCCGAAAAACTCTGTTATATTATGCACTTTATTAAAAGTGAATTCCCTTTAAAGTTAAGTTAGACTGTTTTTTATGTTCAAGGTTAATAATATTAAACTAGTAAAGAAACAGACCGTCAGCCGCACCGGCATCGGTACGGACTGGTGCTCGTACGGACCCTTGGCGTCTTGCTCATAGTTATGAGACATAAGCAAGCCTGTGACTGACCACACCTGTGCAACTTTGTCGTTATATGACCTTAAGTTAATTTCCCACAAGAGCAGTAAGTTCATTTTAGTCTCCTCAACAACGTTCGCTTCATATAGTTAGGTGCTATACGAGTATATGTGGTTCTTAATAGCTTGCCAGGCATAAGGTTCTAAGGACTTCTCTTAGAGCTTATAGCGCTATATGGGTACAGCCCTTTCGGTAAACTCTAGTTCAACTTCGTGCCTCACTAGGCCCTAAATAATAAGCTAGCATTGCATGGATGAACCCTCCACTTATACGAGTTCCCATAGCATCCAGTCGACGTCTAGCTGTGAGTCCCGCTTTTGCGCTATCTATCGGTCTCATGAGTCTAAGGGAGGTTGAGTCTCGGAAGAAGGCGATGCAAATTTCGGTACCTCCTAAGGCAGCCACCGTTTACGCTCCAAGATATACTCTAATGCGGCAACGGGGTTGTACATCGGGTTTGCTCCGGTCTATCTCATGGATACAGCAGACATGGGTCGTTTCCAGTAGTAATACTTCAATTCTTCAACTTCTTAATACTCTTCACGGTCCTCGCTGTCGAGATGGTAATGCCATACCTGCAATTTTACTGCTTCCTTTTATACTACTTTGCTTTATTTCAACAAAGCCATTCAATATATATAAAATTAACCATAGATATAAAATGACTATTAATGAAAGAATAAGATATTAGATCAATTTGTTACACTACTGACATTAACAACTTTTAATACTTTGTATGTAAATACTTTCTAATACTTTCTTTTCTTTGCCCTTTGATAAAATAAAATAAAATTTATAGCCATGAAACAGTATATGCCCCAACAGGTTTAATGTAACGTTGGTCTCTTACTATGCTATATATCTGCTTTGAATAAACTTTAACATACTAACAATGTCCTTAATAGGTGCTGAGTGTATCAATAGTGCCCATGGAACTTTTCAAAGTGGTAAGTTAGATGACTCTACTTCTTGCTCTATTGGCAAGTTGGACTTAGTAGGCTAAGAGGGTAACCATGACAGATGTCATATTGCTTTAATTGAAACCGGTTTTATTTGTTTTAGAGGTTGAGAACCAAACCTGGTGAGCTTAATTCGTTGCAGCTCAGTTGACCTTTAGACCACATGTATAAATTATTAGGAAAGAGACTAAGGCTCGTCGCTTCTCATCCAGTCGGCTCGTAATCAGCTTTAAGGCTCGTTCTTTTGAGAAATCTCTGGTTTTTTTCTAGAAAGGTGTCAGTGTCAAGGAGTAGGTCATCGAGATCATCAAGGTCGTCGCTCCTTCTTTTACTTAGACAGTACGTTAGATTCCGCAGTAATGCTTTGTAAAAGTACCGAGGCTTGCAGGGCGACAATGGAATGAAATAGCCAAAGTTTCCGTTCTGCTAGTACTTTAGATTACTCGTGTTGCAATACGTTGAACGGTTCCTGAATGAAGAGGTTTATCGATTTCTTTATGTCTTAGGAAGACTTCCTATACTCTTATGCGACCTCTTTATACATTTGAATTTCTTAGAGTAAGCGTGGTACCTAGCGTTCGTCATGTGACTCAATCTCGTTATACATAGGGTTAAATTATTCCGTTGTCACCTTCCGCTTTTACGAACCATGTGGAATGCGCTACGACAGCATGTACGACCTTGACCTGGACGAAAACAGTGATAGCTTTTAACCGATTTCTCACCTGTTCAATTTCCGCTTCTCCGCCGTCTTGAGGAAAGATTTTAAGGGTTTTCGCAAGTGGTTGGTTGTGGCCTAATATAGTTTCTTAAGAAACATTCATAGGTAAATACGTGTAAGTATTAATATCAAGGTCTCAGCTAACGTTGTAAATGAATTTGTAGTTTTGAAAATCTCTTCTCATATACTTTAGCCGGTTTTAACCGTAGTTCGCCTTTAGATTTTCTAATTAAGGTGTCACAAAAGAACTACACATACTTTCTCTTTTTATAATTAATATTTTAATTCAAACATATATAAAATTATATATATAAATGCGAATAAAAATCTAATGTAAAATGTATGTCTAAATATGTGTTGCTGTTCAAAATTATTTTCCTGTTTAAATTTTTTTTTGCGTTATAGAAAACAATTTTGAATATGACAACACATATTCACATACGTTTTATATAAACTTTTATTCGTATTTACATATATAATTTTATGTTCAGTCTATTATAATGTTTTTATTTTTTTTCAACTTGTATTATAGTTCATAATGCATAAAACATATTGTAAAAAAACAATAGTTGTTAAAAGACGCGTGAAAGGCCATAATGTTTATAATAAAAAGAGAAAGTATATGTAGTTCTTTGTGACACCTTAATTAAAAAATCCAAAGACGAACGCCGATTAAAATTGAACTCATAAGTAAATATGTCTTATTGGAATTCTGAGTTCGTAAAAAATTAAGAAAAATTGGTTAAGTTTATAAATAATAATTTTTAATATATTTTTTTAATAAAAATTAAAATTTTAAACGTATTAAATTAGACACGTCTTTGACTAATTCTCAAGTCACCATTGAGTTGTTCTTATGCATTTACACAGGCGTGGAGATCGTAAACGACTTGACAATGCGTTATTAATACAGCAGGGTAGAATGATAGGCCATGTGTCAAAGCCGGCATACTGAGAATTCCCTTTACTTCTCTATCTGTTGATATAAAGAATAAACCGATTGGTTAACGTTGTTCCGATGGAGCTTTTGTTATGAGTCTTTTAAGTCTGTAAGTAGAAACGTGAACCACAGTGACGGGTAGGCTTCTAATAGAGTCAGAAACTCGGTATGAATCTTCCGTTCGACGGCTGCTGCTTTATAGTTGCATACGAGTACCACCGGCGAGACATACTAGATAGGTCTCTGTATGGTTTTAATCAACCTGGTTAAAAGATATTCGAGATGTACTTACTTTTGCTTCGAGTGCTTCAAGCAACGTACCGTCATGTCGTTAAATAGGACTGTCTGGGCGGTTACTAATGCAACGTTGCGCATCGCTTTATGTGATTAATGCTAGTCTCACTAGTCCACTTGAGACGGCACTTCTCATGCGATTTGTCGTAGTAATTATGCTTCGCTGGCCTTACCGCTTTAGAGCGATTGTTCCGCGCGTCACAGTCCATAGATCACTTAGGTTTCTTGATACTGTGGACCATGAGCTTCCCGATGATATATCTAAAACTTTTGACCCAAAAGTTTCCTGAATTACACTTTTACCAACGATCATTACTACGGCATGATGGGTCTATACATATGCAACCGGAACTATCATAAAAGTTGAAGGAGGCATTCGGGTGAAAGCTTCATCCTATACGCCATAGCTCGATGTCCGTCCAAATGCTAAACTAGTTACTCAACACCTTGAGGATAGTTAAACTTCTTTACTCTCTCTTTAGTGTTCCTAGTGCGCAACTCTTTGATCGCGTTCTTGAATTCTAATTTAGGCCTGTCTTCTTATTAAACCTTCCTGTACAGGACAAATTGAGCCATATGCCAAGGTACCAGATAGGAATGCTGTTCGTATCTTAATCTCTCCGGCAACGACCACTCATAGAACTTGAAATGAAGTTTTAGTATAATGAATGTTTCTATAATTAATAATAATGTAACGTCGCGACTTTTTTGAAGACTGCTCACTGAGGTTTGACTTTTGTTGACGTAAATTGTTAAAGCTTTTCTATCACTCGAAGGGCTACCTTTACCCACACGGCAAGCAGATAAGAAAGCTCGACGGTCATAAACATTTTTCACTTTAATTGAAATTTTTTCCTCTTGGGTAGTGGAGCTCCTCACCGCAAATGCTTTGCAATAAGACGTTATCTCATGTCTTCGCCAAACCAAAGTATCGTGGAAAGCTCATAGTCTTAATATAACGACCATAACTGTTTTTACCTTATTACGCTCATGGAGATTTTATGCTTTGGTTATAACTATATTTTGTCTTTTTGAAACGTAACTTCTAAGTAGGCTTGTATGGCGTTAGACCTTGCTCATGACCTAATTGTGTAATATCACAGCAAGGGAAGTGGTGCGCTGTTTTATAGAAATTAGAAGTTGGTCATAGATTACTCCCATTATGCGCAGGACACTATTGGAGTCTTTATGTATTTTACTGTTTTCTTTTTCCAGGTAAAAGATAGTTTTATCTTAGACTATGGTGCTTTTTTCTTAGTCAGAACCTTCTATAGCAATGCCCTTAGAGCTTTAGTAGATTAAGATTATTGCTCGCAATATACTTTTATCTGTGTTGTAAACTCAGATTCGTTCATCGTTTCACACTCTATGTTTAGCTGTACTGCAAGCTACGTCACTGTTAAGTGCCATTCTTAGTCGTTGGTAGCGTATTTCTCTACGTTGTGGTGTCATTTGACCTAACCTTTGGTTTGTCATTCCTTTCTTCTCTTCTTTAACATTTACAGGAGTCGCGTCCAGAGTTCAGCCCATGTCAGAAACATCGCCTGCATTCAAAGTCAGAGGGCTCCAATGTTCTGTTGTGAATACAAAAATGACAGCCGTCGCACTCTTCATTATATCTGGTCTTTAATTCTGTAATGAAAATACAATTATGATTACGTCGCGTTCTTCATTTTATACTTGAAACAATGAGTGTTCTTCATGTTATACGTATAGGATGAGGAGAGTTAAAGCTTATGCGGTAATTATTGCTTGGTTTTCTATTTAATTTCCCACATAACGCTATACCTTCTTGGACGTTATGTCCTTTACTTTAACAATAGTGGCCTTCGAGGAGTGTTTCAGGTGTTGACTCTCTATACTATCTCTTGAGGTCGTAATGGTTTGTTACATACCTTCTCTAAGTTTTCCCTTTCTTTAGACACGCTTGAACGTTATTCCGATGACTGCAACGGGTTCACTCCCTAGTCGATTTGAAAGTATAACTACGTAGGGTCGAAAGACTTTATGCGGTCTTTATACTGGTCCACTAGCCAGAGGAGTTAATATGCTTAAATAGTGTTATGTTACAAGTCGTTTTGAGACTTTGGTTGTGGTAGCAACATGTTTTGGGTACCCACTACCATGGTTGGCATGTTCTCGGCACCATGGCAACCCGGTAATTTGGAAGTCTTAGTGTTTCCGTTTCACTTTAACTACAGGACAATCTGCTTCATAGTGTTGCATTCAATGTTTTAAAACATTTTTGTTATTAATTATAAAAAATTAACAAATTCATAATTTTCTGCATTACCATTGATATAGATATACATTACTATTATTAAATACACAATTAATTCTTAGTTTAATTACGAGTTTATTACAGCGTCTTCGAGATTAGAAGATTTGAGATCTGAAACATCGTAAAGTATTAATATAAAATTTATTAAATTTATTAACATCAAGTTTACTTCATTTATTAACAAGATAGTAATCAAAAGCTACGTAATATAGATTGACTCATTTGAGAGAATACTTTACATTACATCTTATAGATAGTTTCAGAGTATATTAAATTGTAGTGACTTAAGAAAAAATGTCGAAGCACGTGAGAGCTATTGCTGTTCTAAGAATGAAAACTATTGGTCGAGATATTGCACGTTGATCCTTTTACATTCGTGCACGAAAATTGGTGGATAGGTGTTCTAAGGGTATTAGCATCTTTAATATAAGGCCTTTCGAGCTTTCACCGACAGAATCGATTTCTATGCCTGTTACTGTCGTCTTTACAAATACATACCGACCCGTTAGATCTTTAACTTGACTTCTTTCACCCATTACTAGATTTTCAACGTTATTTGCCTGTCTTGCAGCTCTATGGCCTCTTTCCGGTAGTTCTTTCGTTACCTTTACTTTAGTAAAAGCTTTAGCACGTTAATGGTCTGCCTAGGGAGAGACAGTAAAGTCTCTTCATACCTTATTGGCAGTATAAGCTGCCTTTTGTGCAAGCTAATATACATTCATTTTTAGATTTTATTGAACTATTATAAAATTCATTAATTAATATTTGTCATGCTACCATGTTCTTTCACATTATTTTTAGCTTTATTTAAACAATGATCCGTTTACCACGTTGTATAGCATTACGGTATGCACCAGAGACGCCGTTGATGCTGAGGTCCGCACTATTGCTAAAGGAGTGAGGATTCTTGACAGAGAACTGCTTTGGCCTTCTTAAACGGCGATGGATGCGATACTGTTTGCTCTTAACGGTTCCTGGACGAGGCCTTTTATTTGCTTTTCGGCTTCCTCGGTACACGTAACTTCACGGTCTCGTCGTTGTCTACTTGCAATAGTCGCTGTCTCTTCGCCCATCTTACTACTGTCTCCCGCCTTTAACCCCGATGGTAGTCAGATTAGCGTTCTTTCTTGTACCTGTTCTATCATTTTCTCCAGTGCCAGTATTTTTTATGTTAGTCTTTCTGAGCGTCCTCCCACCTAGTTTGCTCTCAGTTATATCCTTCTTTGTGTTGTAACAGATAGCGTGCTCTCAACATCTTTACCTGCTATTTTAGACAAAGTGATGATGTGGACATGGTCCTACGGAGGTTCTATGGTCCGGGCATCTTTCTCAAGGCTTCTTTATACTAGAAATAACGGACAGCTTCTTGTTACTTAGATACCTGGAGTTTGCATCCCAGCTCCTTCCGCGATTCGGTCTAAAATGAGTTTTCGGACAGTTGTTCTATGTTTTGAAGGTTTAAGGACAGAGAACGTCGCGGCGTATT

(first three nucl are "TAC" which is the complement of ATG). But it's not the reverse-complement either. So I think there is a snafu in there. These genes of interest are on the reverse strand. Maybe your parser forgets to take the complement?

Also, are you planning on letting users extract more generic labels? (eg: this gff has the "gene" label)

Thanks!
yannick

Hi,
Are there any examples of how to change IDs in a GFF3 file?

Thank you in advance,

Michal

gff3-fetch mRNA m_hapla.WS232.genomic.fa m_hapla.WS232.annotations.gff3 > m_hapla.WS232.mRNA.fa

caused

/home/wrk/.gems/gems/bio-gff3-0.8.7/lib/bio/db/gff/digest/gffparser.rb:119:in block in each_mRNA_seq': undefined methodseqname' for #Array:0x00000000de79e8 (NoMethodError)

User request for output fasta containing direction information of CDS