FASTA格式拆分工具

FASTA格式由威廉·皮尔森(William Raymond Pearson)和戴维·吉布森(David J. Gibson)开发,最初用于序列比对软件BLAST。BLAST是用于比对两个或多个序列的软件,它是生物信息学领域最常用的软件之一。

在开发BLAST之前,序列比对通常使用手动方法进行。这是一个非常耗时和费力的过程。FASTA格式的开发简化了序列比对的过程,使其成为一种更有效的方法。

FASTA格式的组成

由以下两部分组成:

  • 序列名称行:以“>”字符开头,后跟序列的名称和注释。
  • 序列行
    • 核酸:以“A”、“C”、“G”、“T”或“N”字符表示。
    • 氨基酸:以“A”、“R”、“N”、“D”、“C”、“Q”、“E”、“G”、“H”、“I”、“L”、“K”、“M”、“F”、“P”、“S”、“T”、“W”、“Y”或“V”字符表示。

示例:

>Ctg135

ACCCGCAGGGAACAAGAATTGCGCAAAAGTGACACCTACGTCACAACTAA

核苷酸对照表:

字母含义
A腺嘌呤
C胞嘧啶
G鸟嘌呤
T胸腺嘧啶
N未知碱基

氨基酸对照表:

字母含义
A丙氨酸
R精氨酸
N缬氨酸
D天冬氨酸
C半胱氨酸
Q谷氨酸
E谷氨酰胺
G甘氨酸
H组氨酸
I异亮氨酸
L亮氨酸
K赖氨酸
M甲硫氨酸
F苯丙氨酸
P脯氨酸
S丝氨酸
T苏氨酸
W色氨酸
Y酪氨酸
V缬氨酸

FASTA格式具有以下优点:

  1. 格式简单,易于理解和使用。
  2. 可以直接用记事本打开,方便编辑。

多个FASTA格式的序列聚合在一个文件时,通常使文件变得很大,分析会占用大量的内存。因此,推荐时把大文件拆解成小文件,每个文件一条序列,即每次处理一条序列。

直接将下列代码复制到RStudio中(需要联网状态),然后根据提示选择FASTA文件即可。

source("https://dataholdcn.cn/R/GBIT/GBIT.R")   # 导入分析框架
GBIT.splitFasta()   # 将FASTA文件按照序列分开

若在使用中发现BUG,请及时告知,非常感谢!

https://www.lecent.wang/d/10-yuan-chuang

评论

发表评论

了解 数据控|突破是我们的每一步 的更多信息

立即订阅以继续阅读并访问完整档案。

继续阅读