Question: How to extract fasta sequences from assembled transcripts generated by Stringtie
3
gravatar for seta
3.2 years ago by
seta1.3k
Sweden
seta1.3k wrote:

Hi all,

I used STAR and stringtie for mapping reads to reference genome and assembly. As you know, the generated assembled transcripts by stringtie are in gtf format. Now, I want to have fasta sequence of assembled transcript. I used gffread, but all sequences had the same header! maybe it's not compatible with stringtie. Could you please help me out to convert assembled transcripts by stringtie in gtf format to fasta format?

Thanks

ADD COMMENTlink modified 11 months ago by Juke344.5k • written 3.2 years ago by seta1.3k

use gffread, you can find it in cufflink package

ADD REPLYlink written 11 months ago by dukecomeback40
2
gravatar for zzqr
2.9 years ago by
zzqr40
zzqr40 wrote:

The stringtie_merged.gtf file have seqname, start, end strand info. So, you can use R GRanges object and getSeq function from GenomicRanges and BSgenome packages to retrive sequences.

ADD COMMENTlink written 2.9 years ago by zzqr40
1
gravatar for lakhujanivijay
2.9 years ago by
lakhujanivijay5.1k
India
lakhujanivijay5.1k wrote:

You can also use bedtools getfasta to fetch sequences from GTF or BED files.


UPDATE

Here is the perfect solution

ADD COMMENTlink modified 2.7 years ago • written 2.9 years ago by lakhujanivijay5.1k

I used this, but I run into the following error

"Error (GFaSeqGet): subsequence cannot be larger than 465 Error getting subseq for gene1 (465..1503)!"

Did you had any issues using gffread?

Thanks

ADD REPLYlink written 22 months ago by spriyar10
0
gravatar for Juke34
11 months ago by
Juke344.5k
Sweden
Juke344.5k wrote:

I use agat_sp_extract_sequences.pl from AGAT.

agat_sp_extract_sequences.pl --cdna --gff input.gtf --fasta genome.fa -o output.fa

ADD COMMENTlink modified 8 months ago • written 11 months ago by Juke344.5k
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: 1464 users visited in the last hour