GFF read utility (included as a part of cufflinks)

[syintel87]

From reference genome, if I want to extract sequences based on annotation file in order to create new fasta file,
does gffread utility exclude intron regions?

(reference.fasta + annotation.gtf => gffread utility => annotation + exon sequence)

Thank you in advance.

[yzzhang]

Hello, have you done with gffread successfully? I used tophat- cufflinks to analyze my 12 samples, and got a merged.gtf using cuffmerge. I wanted to discard those single exon transcripts and extract the coding sequence and protein sequence. My merged.gtf is 180 Mb and my genome is 280 Mb. I just ran gffread command and I got the gff file and exon sequence using -w argument, but I just got empty cds file and protein file when I used -x and -y argument.
Did you face this condition?
Any help is appreciated.
Best Regards,
yun

[syintel87]

I am sorry, but it did not seem to work properly in my case, so that I wrote my own code using python in order to extract CDS from reference genome file.