技巧从GFF⽂件提取CDS并翻译成蛋⽩Usage:
gffread <input_gff> [-g <genomic_qs_fasta> | <dir>][-s <q_info.fsize>]
[-o <outfile.gff>] [-t <tname>] [-r [[<strand>]<chr>:]<start>..<end> [-R]]
[-CTVNJMKQAFGUBHZWTOLE] [-w <exons.fa>] [-x <cds.fa>] [-y <tr_cds.fa>]
[-i <maxintron>]
i can do itFilters and/or converts GFF3/GTF2 records.
<input_gff> is a GFF file, u '-' if the GFF records will be given at stdin
Options:
-g full path to a multi-fasta file with the genomic quences
美国怪谈剧情
for all input mappings, OR a directory with single-fasta filesopt
(one per genomic quence, with file names matching quence names)
-
lesspressurebetterlifes <q_info.fsize> is a tab-delimited file providing this info
for each of the mapped quences:
<q-name> <q-length> <q-description>
曼彻斯特大学(uful for -A option with mRNA/EST/protein mappings)
-i discard transcripts having an intron larger than <maxintron>
-r only show transcripts overlapping coordinate range <start>..<end>
(on chromosome/contig <chr>, strand <strand> if provided)
-R for -r option, discard all transcripts that are not fully
长裤的英文contained within the given range
-U discard single-exon transcripts
绅士风度英文
-C coding only: discard mRNAs that have no CDS feature
-
F full GFF attribute prervation (all attributes are shown)
-G only par additional exon attributes from the first exon
and move them to the mRNA level (uful for GTF input)
-A u the description field from <q_info.fsize> and add it
as the value for a 'descr' attribute to the GFF record
-O process also non-transcript GFF records (by default non-transcript
records are ignored)
-V discard any mRNAs with CDS having in-frame stop codons
-H for -V option, check and adjust the starting CDS pha
if the original pha leads to a translation with an
in-frame stop codon
华东理工大学排名
-
B for -V option, single-exon transcripts are also checked on the
opposite strand
-N discard multi-exon mRNAs that have any intron with a non-canonical
splice site connsus (i.e. not GT-AG, GC-AG or AT-AC)
-J discard any mRNAs that either lack initial START codon
or the terminal STOP codon, or have an in-frame stop codon
(only print mRNAs with a fulll, valid CDS)
--no-pudo: filter out records matching the 'pudo' keyword
-M/--merge : cluster the input transcripts into loci, collapsing matching
transcripts (tho with the same exact introns and fully contained)
-d <dupinfo> : for -M option, write collapsing info to file <dupinfo>
-
-cluster-only: same as --merge but without collapsing matching transcripts
-K for -M option: also collap shorter, fully contained transcripts
with fewer introns than the container
-Q for -M option, remove the containment restriction:
(multi-exon transcripts will be collapd if just their introns match,
while single-exon transcripts can partially overlap (80%))
--force-exons: make sure that the lowest level GFF features are printed as
"exon" features
-E expo (warn about) duplicate transcript IDs and other potential
problems with the given GFF/GTF records
-D decode url encoded characters within attributes
-
Z merge clo exons into a single exon (for intron size<4)
-w write a fasta file with spliced exons for each GFF transcript
-x write a fasta file with spliced CDS for each GFF transcript
-W for -w and -x options, also write for each fasta record the exon
runboxcoordinates projected onto the spliced quence
-y write a protein fasta file with the translation of CDS for each record
-L Enmbl GTF to GFF3 conversion (implies -F; should be ud with -m)vsc
-m <chr_replace> is a reference (genomic) quence replacement table with
this format:
<original_ref_ID> <new_ref_ID>
