SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Error while extracting fasta seqeunces from stringtie merged gtf. rajarapupriya RNA Sequencing 0 10-05-2018 08:55 AM
Extracting HLA-E allele sequences from BAM for gene editing emily202 Bioinformatics 0 03-14-2017 08:11 AM
Stringtie is breaking up an annotated gene iriet RNA Sequencing 0 04-08-2016 11:52 AM
Extracting all microbial sequences from NT kga1978 Bioinformatics 3 11-17-2012 02:32 PM
Script for extracting random sub-set of sequences silver_steve Bioinformatics 9 01-16-2012 10:07 PM

Reply
 
Thread Tools
Old 08-04-2020, 09:18 AM   #1
AlexSTorson
Junior Member
 
Location: Canada

Join Date: Jul 2020
Posts: 1
Default Extracting gene sequences from Stringtie assembly

Hi all,

Forgive me if this question has been posted before (I'm new to this forum). In the past, when doing transcript-level RNA-seq analyses, I have used the gffread utility to pull transcript sequences for annotation and downstream functional analyses.

i.e.
Code:
gffread -w transcripts.fa -g /path/to/genome.fa transcripts.gtf
However, I am currently working on a gene-level analysis and am having issues doing something similar - especially in regard to evaluating novel loci in my Stringtie assemblies. The way I see it, I have one of two options:

1) Use my RefSeq IDs for non-novel loci and annotate novel loci separately and then merge these annotation lists prior to functional analysis. (If I go this route, I still need to generate a FASTA with the novel loci anyway)

or

2) Pull all gene/loci sequences and re-annotate everything

Is there a (relatively) simple work-around for this?
AlexSTorson is offline   Reply With Quote
Reply

Tags
annotation, rnaseq, stringtie

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:27 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO