Header Leaderboard Ad

Collapse

How do I convert Trinity output to GTF/GFF?

Collapse

Announcement

Collapse

SEQanswers June Challenge Has Begun!

The competition has begun! We're giving away a $50 Amazon gift card to the member who answers the most questions on our site during the month. We want to encourage our community members to share their knowledge and help each other out by answering questions related to sequencing technologies, genomics, and bioinformatics. The competition is open to all members of the site, and the winner will be announced at the beginning of July. Best of luck!

For a list of the official rules, visit (https://www.seqanswers.com/forum/sit...wledge-and-win)
See more
See less
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • How do I convert Trinity output to GTF/GFF?

    Hello everyone and thank you for your time,
    This is a reiteration of an unanswered question:
    TRINITY:Input (Fastq) and Output files (GFF)

    How can I convert (FASTA) transcripts produced from Trinity into a GTF/GFF file? I am not interested in alternative workflows for particular things, I just want a GTF/GFF annotation file. Also, I'm not familiar with any long sequence alignment algorithms (would BWA be best??) and a simple conversion to GTF/GFF. Can you please suggest a simple workflow?

  • #2
    Answered my own question

    Well, since this problem went unanswered, I figured out a way to resolve this issue. It involves all of the things I thought it might, and there may be a simpler way. If you have a better method, please reply so that your solution will be easy to find.

    Solution:
    bwa-mem of Trinity.fasta to reference genome
    samtools -Sbh | bam2bed | bed12ToBed6 | bedToGenePred stdin stdout | genePredToGtf file stdin assembly.gtf

    Comment


    • #3
      Hi
      I have the same problem.
      How can I replace Trinity gene numbers to real gene names ?
      I have blast hits file from TransDecoder of Trinity.
      My question is that how to run DEG to get real gene names, not the Trinity names (Tr1...Trn..)?

      Comment

      Latest Articles

      Collapse

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 06-01-2023, 08:56 PM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 06-01-2023, 07:33 AM
      0 responses
      9 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 05-31-2023, 07:50 AM
      0 responses
      6 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 05-26-2023, 09:22 AM
      0 responses
      11 views
      0 likes
      Last Post seqadmin  
      Working...
      X