Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • mattanswers
    Member
    • Oct 2009
    • 65

    #16
    I would also be very interested in generating a gtf file for arabidopsis. Would it be possible to post the code ?

    Comment

    • celeste8
      Junior Member
      • Aug 2010
      • 2

      #17
      Download TAIR10_GFF3_genes.gff inTAIR, and then use the Perl script below to convert this gff file to gtf which is used as conference annotation Cufflinks.

      #!/usr/bin/perl

      use warnings;
      use strict;
      while (<>) {
      chomp;
      my @parts = split /\t/;
      if ( !(($parts[2] eq 'three_prime_UTR')|($parts[2] eq 'five_prime_UTR')|($parts[2] eq 'exon')|($parts[2] eq 'CDS')) ) {
      next;
      }
      elsif ( $parts[2] eq 'three_prime_UTR') {
      $parts[2] = '3UTR';
      }
      elsif ($parts[2] eq 'five_prime_UTR') {
      $parts[2] = '5UTR';
      }
      $parts[8] =~ s/^Parent=((\w+)\.\w*).*/gene_id \"$2\"; transcript_id \"$1\";/sg;
      $_ = join "\t" , @parts;
      print "$_\n";
      }

      Comment

      • mattanswers
        Member
        • Oct 2009
        • 65

        #18
        celeste8, thank you very much !

        Comment

        • samhokin
          Member
          • Nov 2013
          • 20

          #19
          Thanks, celeste8!!
          Sam Hokin
          Computational Scientist, Carnegie and NCGR

          Comment

          Latest Articles

          Collapse

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, 06-09-2026, 11:58 AM
          0 responses
          17 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 06-05-2026, 10:09 AM
          0 responses
          27 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 06-04-2026, 08:59 AM
          0 responses
          38 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 06-02-2026, 12:03 PM
          0 responses
          61 views
          0 reactions
          Last Post SEQadmin2  
          Working...