Question: GFF3 to GTF
0
gravatar for silas008
3.6 years ago by
silas008100
Brazil
silas008100 wrote:

I'm tryning to convert a c elegans mirBase gff3 file to gtf file with gffread (by cufflinks) but my gtf output is blank.

The code:

gffread miRNA_WBcel215.gff3 -T -o miRNA_WBcel215.gtf

Who knows why?

Thank you

gff gtf • 4.0k views
ADD COMMENTlink modified 13 months ago by RamRS24k • written 3.6 years ago by silas008100

Hmm.. the command seems fine. Could you add some lines of WBcel215.gff3 in order to see how it looks like? For other hand, maybe you'd like to try the other solutions (apart from the gffread utility) proposed in this post.

ADD REPLYlink modified 13 months ago by RamRS24k • written 3.6 years ago by iraun3.6k
0
gravatar for Juke-34
3.6 years ago by
Juke-342.8k
Sweden
Juke-342.8k wrote:

You could try the simple Perl approach to check if it works:

use Bio::Tools::GFF;

# $gff corresponds to your gff3 input file
my $gff_in = Bio::Tools::GFF->new(-file => $gff, -gff_version => 3);

# $gtf corresponds to your gtf output file
my $gtf_out = Bio::Tools::GFF->new(-file => $gtf, -gff_version => 2.5);

while( my $feature = $gff_in->next_feature) {
    $gtf_out->write_feature($feature);
}
ADD COMMENTlink modified 13 months ago by RamRS24k • written 3.6 years ago by Juke-342.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1440 users visited in the last hour