gff文件_GFF文件格式说明

gff格式是Sanger研究所定义,是一种简单的、方便的对于DNA、RNA以及蛋白质序列的特征进行描述的一种数据格式,比如序列的那里到那里是基因,已经成为序列注释的通用格式,比如基因组的基因预测,许多软件都支持输入或者输出gff格式。目前格式定义的最新版本是版本3。原始定义见SONG website

gff是纯文本文件,由tab键隔开的9列组成,以下是各列的说明:

Column 1: seqid

序列的编号,编号的有效字符[a-zA-Z0-9.:^*$@!=_?-|]

Column 2: source

注释信息的来源,比如"Genescan"、"Genbank"等,可以为空,为空用"."点号代替

Column 3: type

注释信息的类型,比如Gene、cDNA、mRNA等,或者是SO对应的编号

Columns 4 & 5: start and end

开始与结束的位置,注意计数是从1开始的。结束位置不能大于序列的长度

Column 6: score

得分,数字,是注释信息可能性的说明,可以是序列相似性比对时的E-values值或者基因预测是的P-values值。"."表示为空。

Column 7: strand

序列的方向, +表示正义链, -反义链 , ? 表示未知.

Column 8: phase

仅对注释类型为 "CDS"有效,表示起始编码的位置,有效值为0、1、2。

Column 9: attributes

以多个键值对组成的注释信息描述,键与值之间用"=",不同的键值用";"隔开,一个键可以有多个值,不同值用","分割。

注意如果描述中包括tab键以及",=;",要用URL转义规则进行转义,如tab键用 %09代替。键是区分大小写的,

以大写字母开头的键是预先定义好的,在后面可能被其他注释信息所调用。

预先定义的键包括:

ID 注释信息的编号,在一个GFF文件中必须唯一;

Name 注释信息的名称,可以重复;

Alias 别名

Parent Indicates 该注释所属的注释,值为注释信息的编号,比如外显子所属的转录组编号,转录组所属的基因的编号。值可以为多个。

Target Indicates: the target of a nucleotide-to-nucleotide or protein-to-nucleotide alignment.

Gap:The alignment of the feature to the target if the two are not collinear (e.g. contain gaps).

Derives_from:Used to disambiguate the relationship between one feature and another when the relationship

is a temporal one rather than a purely structural “part of” one. This is needed for polycistronic genes.

Note 备注

Dbxref 数据库索引

Ontology_term: A cross reference to an ontology term.

例子

### 使用 `gff3_to_tbl` 工具的方法 #### 安装依赖工具 为了能够顺利运行 `gff3_to_tbl` 脚本,通常需要安装一些必要的软件包。这些可能包括 Perl 及其模块、BioPerl 库以及其他辅助程序。 ```bash sudo apt-get update && sudo apt-get install -y perl bioperl ``` #### 下载并准备脚本 可以从 NCBI 或其他资源下载最新的 `gff3_to_tbl.pl` 脚本文件,并将其放置在一个合适的目录下以便执行: ```bash wget https://ftp.ncbi.nlm.nih.gov/genomes/TOOLS/gff3_to_tbl_v2.8.tar.gz tar zxvf gff3_to_tbl_v2.8.tar.gz cd gff3_to_tbl_v2.8/ chmod +x *.pl export PATH=$PATH:/path/to/gff3_to_tbl_directory ``` #### 准备输入文件 该工具主要处理 GFF3 文件作为输入数据源。确保拥有一个格式良好的 GFF3 文件用于转换操作。此外还需要 FASTA 格式的序列文件来补充基因组位置信息[^1]。 #### 运行命令示例 下面是一个简单的例子展示如何调用此工具来进行转换工作: ```bash perl gff3_to_tbl.pl input.gff3 output.tbl <fasta_sequence.fa> ``` 这里: - `input.gff3`: 是待解析的 GFF3 注释文件路径; - `output.tbl`: 将要生成的目标表格文件名; - `<fasta_sequence.fa>`: 提供给定特征的具体坐标范围内的碱基序列信息; 请注意,在实际应用过程中可能会遇到不同版本间的差异以及特定参数设置的需求,因此建议查阅官方文档获取最准确的帮助指南。 #### 验证输出结果 一旦完成了上述步骤之后,可以打开产生的 `.tbl` 文件查看结构化后的表单内容是否符合预期标准。如果一切正常,则说明已经成功掌握了基本的操作流程。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值