gpertea / gffread

GFF/GTF utility providing format conversions, region filtering, FASTA sequence extraction and more

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

gffread -w works but not -x and -y

Ruiqi-CUB opened this issue · comments

Hi there, thanks for developing this very useful tool!

I tried to use it to extract CDS and protein from a genome and corresponding gff file.

It works with -w flag and gives me the output. However, it doesn't work with -x or -y flag. Do you happy to know what could cause the issue?

Here is the code I used:
gffread -w(or -x or -y) output.fasta -g genome.fna genome_annot.gff

Thanks a lot

Hello, may I ask if this problem you said has been solved? I'm having the same issue.

I am checking this on the last version but I cannot reproduce the issue it works as expected on the small example files:
gffread -x transcripts_CDS.fa -g genome.fa annotation.gff

gffread -y transcripts_proteins.fa -g genome.fa annotation.gff

Note that these options require the presence of the CDS features in the user provided annotation gff/gtf,