高考考试网
当前位置: 首页 高考资讯

生信分析(生信分析必须了解的4种文件格式)

时间:2023-05-24 作者: 小编 阅读量: 8 栏目名: 高考资讯

fasta文件格式FASTA文件主要用于存储生物的序列文件,例如基因组,基因的核酸序列以及氨基酸等,是最常见的生物序列格式,一般以扩展名fa,fasta,fna等。从第二行开始为序列本身,只允许使用既定的核苷酸或氨基酸编码符号。例如基因组中的单碱基突变,SNP,插入/缺失INDEL,拷贝数变异CNV,和结构变异SV等,都是利用VCF格式来存储的。

在做生物信息的过程中,经常需要进行各种文件格式。每一种生物软件都有固定的文件格式要求。因此,需要非常每一种数据的文件格式,从某种意义上来说,生物信息分析的过程就是进行各种文件格式的转换过程。例如当前很多分析都可以概括为从fastq到bam,从bam到vcf的过程。

fasta文件格式

FASTA文件主要用于存储生物的序列文件,例如基因组,基因的核酸序列以及氨基酸等,是最常见的生物序列格式,一般以扩展名fa,fasta,fna等。fasta文件中,第一行是由大于号">"开头的任意文字说明,用于序列标记,为了保证后续分析软件能够区分每条序列,单个序列的标识必须是唯一的,序列ID部分可以包含注释信息。从第二行开始为序列本身,只允许使用既定的核苷酸或氨基酸编码符号。序列部分可以在一行,也可以分成多行。

>gi|556503834|ref|NC_000913.3|:190-255 Escherichia coli str. K-12 substr. MG1655, complete genome
ATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGA
>gi|556503834|ref|NC_000913.3|:337-2799 Escherichia coli str. K-12 substr. MG1655, complete genome
ATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTC
TGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCT
GGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATT
TTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGGTTCCCGCTGGCGCAATTGAAAACTTTCG
TCGATCAGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAG

fastq文件格式

fastq文件格式是用来存储测序文件的,它是含有quality的fasta文件。

@DJB775P1:248:D0MDGACXX:7:1202:12362:49613
TGCTTACTCTGCGTTGATACCACTGCTTAGATCGGAAGAGCACACGTCTGAA

JJJJJIIJJJJJJHIHHHGHFFFFFFCEEEEEDBD?DDDDDDBDDDABDDCA
@DJB775P1:248:D0MDGACXX:7:1202:12782:49716
CTCTGCGTTGATACCACTGCTTACTCTGCGTTGATACCACTGCTTAGATCGG

IIIIIIIIIIIIIIIHHHHHHFFFFFFEECCCCBCECCCCCCCCCCCCCCCC
  • 第一行:以‘@’开头,是这一条read的名字,这个字符串是根据测序时的状态信息转换过来的,中间不会有空格,它是每一条read的唯一标识符,同一份FASTQ文件中不会重复出现,甚至不同的FASTQ文件里也不会有重复;

  • 第二行:测序read的序列,由A,C,G,T和N这五种字母构成,这也是我们真正关心的DNA序列,N代表的是测序时那些无法被识别出来的碱基;

  • 第三行:以‘ ’开头,在旧版的FASTQ文件中会直接重复第一行的信息,但现在一般什么也不加(节省存储空间);
    第四行:测序read的质量值,这个和第二行的碱基信息一样重要,它描述的是每个测序碱基的可靠程度,用ASCII码表示。

质量值体系

从表中可以看到下限有33和64两个值,我们把加33的的质量值体系称之为Phred33,加64的称之为Phred64(Solexa的除外,它叫Selexa64)。不过,现在一般都是使用Phred33这个体系,而且33也恰好是ASCII的第一个可见字符('!')

sam格式介绍

sam文件主要用来存储短序列比对的结尾,即将测序数据定位到基因组上的表示形式。

第一列:是reads ID
第二列:是flag标记的总和
第三列:比对到参考序列上的染色体号。
第四列:为在参考序列上的位置
第五列:比对的质量值,MAPQ
第六列:代表比对结果的CIGAR字符串
第七列:mate比对到的染色体号,若是没有mate,则是*
第八列:比对到参考序列上的第一个碱基位置
第九列:Template的长度,
第十列:为read的序列
第十一列:为ASCII码格式的序列质量;

VCF文件格式介绍

VCF是Variant Call Format的简称,是一种定义的专门用于存储基因序列突变信息的文本格式。在生物信息分析中会大量用到VCF格式。例如基因组中的单碱基突变,SNP, 插入/缺失INDEL, 拷贝数变异CNV,和结构变异SV等,都是利用VCF格式来存储的。将其存储为二进制格式就是BCF。

1.CHROM [chromosome]: 染色体名称,
2.POS [position]: 参考基因组突变碱基位置,如果是INDEL,位置是INDEL的第一个碱基位置。
3.ID [identifier]: 突变的名称,
4.REF [reference base(s)]:参考染色体的碱基
5.ALT [alternate base(s)]: 与参考序列比较,发生突变的碱基,
6.QUAL [quality]: Phred标准下的质量值
7.FILTER [filter status]:使用其它的方法进行过滤后得到的过滤结果
8.INFO

文件格式介绍:https://genome.ucsc.edu/FAQ/FAQformat.html#format1



,
    推荐阅读
  • 橘红好还是陈皮好(陈皮橘红与化橘红)

    化橘红化橘红与上述两种重要同科不同种。明末清初化州橘红因其疗效确实,已得到进一步的肯定,受到当时医家及官家的赏识。当前无论调剂配方或配制成药,一律配付柚类皮橘红。事实上,解放以后,化州柚产量一度萎缩,致使化州柚来源的化橘红难以为继。毛橘红总黄酮含量、柚皮苷含量均高于光橘红。这意味着,毛橘红在祛痰、镇咳、消炎作用均优于光橘红。由此可证,道地药材毛橘红质量显著优于光橘红。

  • 城市人套路深(城市套路深我要回农村)

    后辞去教职,加入外企,为公司高管。《村庄》中读历史,是在读人性,读伦理,那些记忆中有关家乡的一切,都浓缩在近百年来几代人经历的生存过程中。每个人都有一个故乡,在外界纷扰的“城市”中,找到心中的“村庄”,心静则安宁。

  • 产品迭代是怎么进行的(产品版本号是如何确定的)

    产品的版本号是如何确定的?常识类概念常见的版本号命名规则在这里先简单的和大家聊一下关于版本的常见命名规则和几种类型。在大部分情况下常见的版本号是三段式的,即X.Y.Z我们用X来表示大版本号,一般当产品出现重大更新、重写、不再向后兼容的情况时我们会在X上加1,当X是0的时候我们默认为开发或测试阶段。

  • 弘扬伟大精神凝聚奋斗伟力(伟大时代呼唤伟大精神)

    庆祝中华人民共和国成立70周年大会盛大举行。会议明确提出建立健全党和国家功勋荣誉表彰制度,决定成立党和国家功勋荣誉表彰工作委员会。作出的一系列重要指示批示为做好这项工作提供了根本遵循。2021年6月29日,人民大会堂金色大厅,党中央首次颁授“七一勋章”。对“七一勋章”获得者的高度赞誉,传遍会场内外。

  • 赵本山压不住小沈阳场面(我不敢用小沈阳)

    2009年,小沈阳凭借小品《不差钱》一炮而红,名声响彻大江南北。赵本山生气了,开始限制小沈阳的演出,并冷落他。小沈阳当时的表演舞台主要在洗浴中心、夜场和迪厅。与此同时,宋小宝也在为了生活出卖尊严。同年5月8日,“赵本山二人转大奖赛”正式启动。小沈阳参加了第一届大赛并获得了铜奖,冠军是吉林省“四大名旦”之一的闫学晶。

  • 红楼梦的简介(关于红楼梦的简介)

    下面希望有你要的答案,我们一起来看看吧!红楼梦的简介《红楼梦》是中国古代的章回体长篇小说,又名叫做《石头记》,被列为中国的古典四大名著之首,一般认为是清代作家曹雪芹所写。《红楼梦》的版本可以分为120回程本和80回脂本两个系统,程本指的是程伟元的印刷本,脂本指的是脂砚斋在不同时期所抄录的早起手稿。

  • 免费的lol专用加速器(分享好用的lol手游免费加速器)

    英雄联盟加速器推荐。《英雄联盟》游戏是一款非常火爆的英雄对战MOBA竞技网游。《英雄联盟》有很多服务地区,其中,国服、韩服、美服都是非常火爆的,这时候就需要游戏加速器来辅助加速游戏,不然会因为晚间人数众多导致游戏的延迟、卡顿、连接不上等问题的出现。

  • 培训师资格证怎么考(如何考培训是资格证)

    下面更多详细答案一起来看看吧!考试时间安排:每年5月和11月。费用:预计在1500元-3000元之间。

  • 猫能活多大(猫能活多久)

    以下内容希望对你有帮助!猫能活多大有调查统计显示,科学饲养的家猫寿命一般在15岁左右,猫咪到了八岁的时候已经进入老年期。当然这个数值并不是固定的,猫咪的寿命长短跟猫咪自己的遗传基因、主人的饲养方法、生活环境等多方面因素有关。在英国有一只名字叫做露西的猫咪,活到了39岁,好比一个人活到了172岁,是世界上最长寿的猫咪。

  • 酸辣粉皮汤怎么做(做粉皮的步骤)

    以下内容希望对你有帮助!酸辣粉皮汤怎么做鲜粉皮3张,香菜5克,竹笋片10克,米醋15克,胡椒粉3克,酱油8克,精制盐1克,味精2克,熟猪油12克。将粉皮切成1厘米宽,5厘米长的条;取大汤碗,放入胡椒粉、米醋和香菜。炒锅置旺火上,下肉汤、笋片、粉皮、酱油、盐、味精烧沸,加入熟油,冲入大汤碗里即可。