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

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

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

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



,
    推荐阅读
  • 四川仁寿的人口有多少(仁寿人说湖广填四川)

    当今仁寿的人口祖籍,绝大部分都是外省各地,其中以湖广为最多,湖广中又以麻城和孝感乡为第一,这与整个四川移民的状况几乎完全吻合。总之,麻城孝感乡人最多,这是不争的事实。这富、荣二邑,在明朝大致包括川南地区的泸州北部、荣昌、隆昌、富顺等地。在光绪《井研县志·氏族》中,有34个姓氏家族都是明初以来迁入该地的。据《仁寿县志》乾隆、嘉庆版载:“旧额人丁八千三百七十八户”。

  • 关于随礼这件事你怎么看(话说随礼你随礼了吗)

    婚姻礼俗,馈赠风俗。送礼时,集体合伙,分摊钱财。是在男女结婚前,亲友们按约定的数额出钱,集中之后做贺礼,交付将要结婚的男女。旧时随礼只限于男家的亲友,现代女家的亲友也从此俗。好像不久后,也会迎来自己的婚礼,也会收到礼钱。当时每人送了一百块,参加了她专门为同事准备的婚宴。婚宴上,算是认识了。再过了不久,发现被拉黑了,消息发送不出去了。随着时代的不同,送礼也是水涨船高。

  • 腰椎轻度骨质增生该怎么治疗(得了腰椎骨质增生应该怎么治疗)

    由于腰椎的特殊生理位置的限制,因此手术破坏了正常的腰椎骨生理结构,造成手术损伤大,易造成腰椎术后失稳,术后瘢痕组织粘连,术中误伤神经根等一系列不良反应。凭借着这些优点,椎间孔镜技术已成为治疗腰椎骨质增生的最终选择。

  • 飞机刹车主要靠发动机反推吗(飞机减速主要靠反推还是刹车)

    如果跑道是干的,开不开反推,对飞机着陆后的制动距离几乎没有影响。第二会影响发动机短舱的设计,后期还要考虑反推的维护成本。最重要的,发动机反推并不是适航规章要求的,但是机轮刹车是。而且这个清单是需要国家的民航局和适航局批准的。但是这个飞机可是英国女王的专机,很排面的。BAE146从1983年生产到了2002年,一共产了三四百架,是英国最成功的的支线客机了。

  • 中学时代的重要性(中学时代对我们有什么意义)

    中学时代的重要性中学时代对我们有什么意义,中学阶段是人生发展的一个新阶段,可以为我们的一生奠定重要基础。中学时代见证着一个人从少年到青年的生命进阶。这将为我们的人生长卷打上更加丰富而厚实的底色。中学时代是孩子成长过程中很重要的一环,是孩子青春期认识社会的一段关键时期,对孩子性格的养成有十分重要的意义。

  • 关于春节的介绍(关于春节的简介)

    关于春节的介绍春节,即中国农历新年,俗称新春、新岁、岁旦等,口头上又称过年、过大年。春节历史悠久,由上古时代岁首祈岁祭祀演变而来。在春节期间,全国各地均有举行各种庆贺新春活动,带有浓郁的各地域特色。百节年为首,春节是中华民族最隆重的传统佳节。受到中华文化的影响,世界上一些国家和地区也有庆贺新春的习俗。据不完全统计,已有近20个国家和地区把中国春节定为整体或者所辖部分城市的法定节假日。

  • 让人绝望的末日类电影(在5部雪地生存类电影中)

    凶手假借摄影师寻找模特儿或搭便车等名义,在长达12年的期间,前后连续绑架并残杀了17名年轻的应召女郎,并用小飞机将尸体弃置在人迹罕至的荒原里。带来补给物资的罗西前来告诉完成研究报告的法利,不久后将带他离开荒原,并计划引入外来投资者在这里建立渡假村,法利对此不予置评。事实上,不只是对于狼族,人类的种种行为已经严重破坏了此地区的生态平衡。

  • 怎样鼓励高三孩子自信(怎么缓解高三孩子的压力)

    下面更多详细答案一起来看看吧!怎样鼓励高三孩子自信父母的信任和鼓励,对孩子来说,就是心情的稳定剂,是保证孩子能够良好发挥的基石。所以,只要家长和孩子一起努力,订立合理目标,明确前进方向,确实执行措施,让孩子看到自己的进步和实力,信心就会自然而然的产生。坚持每天锻炼身体1小时,但要减少剧烈运动。

  • 唯美好听的微信昵称女(唯美好听的微信昵称女有哪些)

    温暖人心小太阳冰纨,今天小编就来聊一聊关于唯美好听的微信昵称女?接下来我们就一起去研究一下吧!唯美好听的微信昵称女温暖人心小太阳冰纨倒影年华我爱你,别介意玉婳瀚海阑干深海爱人情太深岛未晴,心未眠勾勒出你的微笑一米阳光似水流年冰冷艳城听风的倾诉心已成伤却未碎水墨江山

  • 乘法的初步认识学习内容解读(乘法初步认识的学习指导)

    整理了前几年给家长写的关于乘法初步认识这个单元的学习指导,希望对家长有用。),相同加数有4个,就在乘号后面写上4,2×4=8,2表示相同加数是2,4表示有相同的4个2,整道算式就表示4个2相加是8。乘法算式的各部分名称。加法中是“加数、加数、和”,乘法中是“乘数、乘数、积”,积就比和字多两笔。第二个乘数是5,就摆5组,一共是5个4。摆学具是帮助理解巩固乘法意义的很好的方法。