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

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

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

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



,
    推荐阅读
  • 大猩猩被人类囚禁30年(被囚禁30年的红毛猩猩抑郁离世)

    电影里,被人类当作实验品的红毛猩猩凯撒,在服用基因药物后拥有了高智商。6个月后,审批手续办好,夏特克被琳恩带回了家,从此过上了与祖辈们截然不同的生活。猩猩要和人类上学的消息,成为当地轰动一时的新闻。琳恩火速赶往校长办公室。警察赶来,将夏特克团团围住。可她的申请被驳回,理由是要保证当地居民的安全。那些日子,琳恩继续为了夏特克的自由而不断地东奔西走,终于迎来了一个好消息,此时已过去漫长的11年。

  • 小儿发烧抽搐是什么样(小儿发烧抽搐是怎么回事)

    众所周知,新生儿本身免疫力要比成年人要差一些,所以很容易发烧感冒。并且,新生儿有时候发烧时还会出现抽搐的情况。尤其是在发烧中的宝宝,夜间护理更是十分重要。

  • 老墙抹灰技巧(墙面翻新粉刷技巧分享)

    旧墙面在翻新前一定要做墙面基底处理,粉化严重的要全部铲除,铲除后抹上一层白水泥。墙面底漆的作用是抗碱防潮、封固基材、保护面漆,提高面漆的质感和遮盖力。所以墙面的修补最好是整块墙面,或是尽可能地选择较大面积进行修补,这样可以尽量避免产生色差,修补后的质量也更好。墙面色差是修补墙面时需要尽量避免的。为此,应该注意按照油漆桶的说明配比、涂刷,最好估算好用量。

  • 如何做猕猴桃水果沙拉(制作的步骤是什么)

    以下内容希望对你有帮助!如何做猕猴桃水果沙拉首先将猕猴桃用清水洗一遍,打一盆水,加入一小勺食用盐,搅匀后再把所有的水果放进去泡上10分钟,这个步骤可以有效的去除果皮上残留的农药!再洗三片大一点的生菜叶摆在果盘上。将洗好的水果捞起来,滤干水分,把苹果去皮去核,和香蕉,猕猴桃切成丁状,草莓对半切,把从超市买的沙拉酱和酸奶一起放在碗里,均匀搅拌。

  • 牛肉中哪个部位的肉最营养最嫩 牛肉的哪个部位营养价值最高

    大部分都是脂肪含量低的精肉,全是红肉,没有半点的肥油。适合煎、炒、炸、牛排,老少皆宜。菲利是切自沙朗肉中段,可切成牛排或烧烤。牛里脊肉是脊骨里面的一条瘦肉。

  • 中秋节送什么礼物给父母最合适(可以送父母这两款暖心好物)

    中秋节送什么礼物给父母最合适中秋节马上就要到了,虽然三天时间很短暂,但是能跟爸爸妈妈在一起做饭、做家务、聊天说笑应该是最幸福的事情了。而说到准备中秋礼物这件事,现在能给二老的更应该是健康、实用,而不该追求所谓的面子而买一些华而不实的礼品。

  • 抖音letmeseeyou得给得给是什么歌(letmeseeyou得给得给的歌名)

    下面内容希望能帮助到你,我们来一起看看吧!抖音letmeseeyou得给得给是什么歌“letmeseeyou得给得给”歌名是《Swagger》。

  • 传祺ga6高级版(传祺GA6是一台好开)

    而1.5T爱信6AT的组合最低油耗为6.1L,这个油耗表现甚至比本田思域还要低一些,同时在质量上也要更加可靠。经测试,在城市道路和快速高速道路各半的情况下,广汽传祺车的油耗表现都比较令人满意,符合一台中型轿车应有的水准。在性能上,传祺GA6搭载了,百公里油耗仅需6.3L,因此很多朋友都觉得,整车才是家用好车。油耗方面,百公里综合油耗为6.8L,在同级的自主紧凑型车中属于优秀的水准。

  • 肉价高企(肉价为啥这么高)

    近期猪肉价格较快上涨,引起社会关注。据有关部门最新通报,今年肉类供应总体有保障。毋庸讳言,眼下,生猪生产形势濒临困境。本轮非洲猪瘟来势凶猛,全世界尚无有效疫苗。合理范围内的猪价波动,有利于产业的优胜劣汰和转型升级。美国是生猪生产大国,近年来也经历了4次生猪价格波动,每次4至6年。我国生猪价格波动比美国还要剧烈些,缘故在于产业素质不够高,供应和需求衔接不够好。

  • 关于李白送别诗的诗句(总结了5首关于李白的诗词)

    今天给大家推荐这五首李白传唱度很高的诗歌,希望大家能够喜欢。飞流直下三千尺,疑是银河落九天。桃花潭水深千尺,不及汪伦送我情。想渡黄河,冰雪堵塞了这条大川;要登太行,莽莽的风雪早已封山。烹羊宰牛且为乐,会须一饮三百杯。岑夫子,丹丘生,将进酒,杯莫停。古来圣贤皆寂寞,惟有饮者留其名。五花马,千金裘,呼儿将出换美酒,与尔同销万古愁。