The program gffread can be used to validate, filter, convert and perform various other operations on GFF files. Because the program shares the same GFF parser code with Cufflinks, Stringtie, and gffcompare, it could be used to verify that a GFF file from a certain annotation source is correctly "understood" by these programs. Thus the gffread utility can be used to simply read the transcripts from the file, and optionally print these transcripts back, in either GFF3 (default) or GTF2 format (with the -T option), while discarding any non-essential attributes, optionally fixing some potential issues with the input file(s)
gffread is invoked using the following:
- Input (s)
Full path to a multi-fasta file with the genomic sequences for all input mappings (reference genome - can either be a custom or a reference genome from DE).
- The annotation file (gtf file)
Fasta file with spliced exons for each GFF transcript (Default - output.fa)
Please work through the documentation and add your comments on the bottom of this page, or email comments to firstname.lastname@example.org. Thank you.
The test data for gffread-2.2.1 consists of a sample gtf file and 1st chromosome from Arabidopsis thaliana
Run SplAdder as following:
- Input (s)
- Sample gtf file - Sample_cuffcompare_out.gtf
- Chromosome 1 fasta file - TAIR10_chr1.fasta
- output.fa fasta file with spliced exons for each GFF transcript