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

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

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

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



,
    推荐阅读
  • 做生意欠工人工资(打工3个月倒欠公司6万)

    只差半秒,刘攀将葬身在边境的断崖之下。上午10点左右,刘攀所在的D组开始上班,通常一直要工作到第二天凌晨两点,这还取决于当月业绩完成情况。办案民警向受害人取证单纯的她不知道,屏幕外的刘攀正照着公司制作的话术本一步步将她引入陷阱。刘攀只能在心里祈祷女孩发现端倪,趁早收手,哪怕对方狠狠骂他一顿。

  • 华为matepadpro什么时候更新鸿蒙(华为MatePad11樱语粉上手)

    于大学生而言,平板因其便携性,逐步取代PC成校园主力设备,但大学生群体需求多元、使用场景多变对平板适用性提出新的挑战。屏幕方面,华为MatePad11采用了一块10.95英寸16:10的全面屏,支持120Hz高刷新率、DCI-P3电影级色域,屏占比达到86%。华为MatePad11搭载了7250mAh电池、支持22.5W快充、4喇叭大振幅哈曼卡顿联合调音的扬声器,影音兼具续航持久。鸿蒙生态使华为MatePad11实现了自身价值最大化,也打出了华为平板在市场中的差异化表现。

  • 圆润身材怎样变成扁平身材(从有肚腩的假人模特开始)

    无论如何,至少对品牌而言,这个尊重身材多样性的举动引来了关注。此外,服装的功能性如收小腹、保护胸部、减少背部勒痕等也被最大程度呈现。文章发表后引发激烈批评,大量网友前往其社交账号下要求道歉,目前这则新闻已被删除。英国调研机构Statista的数据显示,大码女性占据了美国女性人口的近70%,潜在消费金额为240亿美元。

  • 冰晶可以一直放在冰箱里冷冻吗(冰晶的作用)

    因为冰晶属于耐冻的物品,冰箱冷冻室内的低温不会对冰晶本身造成破坏,一直放在冰箱内冷冻一般也不会对冰箱产生危害。冰晶是空调扇的必备组成部分,其是空调扇制冷的核心,如果空调扇无冰晶,那么其降温功能会受到阻碍,且冰晶是无可代替的,冰块虽同样具有降温功效,但是其效果不够明显,且使用不长久,又不能够反复使用。

  • 月子餐吃什么既下奶又有营养(月子期间吃这3道菜可以补钙和促进钙吸收)

    另外,还要多吃水果,可促进钙质被有效吸收和利用。此外,口蘑中的微量元素也可提高机体的免疫力。原料:口蘑、豌豆各100克,高汤、盐、水淀粉、植物油各适量;制作:1、口蘑洗净,切成小丁;豌豆洗净。

  • 恶魔公主撞到腹黑少爷(公主殿下复仇记:初遇司徒依然)

    公主殿下复仇记:初遇司徒依然欢迎大家来到范小稀的公主殿下复仇记漫画,这一集讲述的是冷梦幽初遇同父异母的妹妹还是在任性酒吧不远处,冷梦幽看完信息刚想回去,手机却被人撞掉在地上,而罪魁祸首还没有一点歉意冷梦幽冷冷地说道:两人没有丝。

  • 滚石历年封面(滚石10年)

    去年,“滚石唱片”推出了“滚石40”的品牌,之后像“滚石撞乐队40团拼经典”,以及最近曹格推出的翻唱20首广东歌经典的《我为歌狂》,都是这个品牌在录音作品层面的出品。《滚动1990》这张合辑有三首歌,但也可以说有四十九首歌,前者说的是音轨,后者说的是出现的作品总数量,因为这也是将四十九首歌分成三首串烧联唱的合辑。在唱片的封面左上角,有“滚石10年”特别设计的Logo,而后面则有“纪念盘3”的字样。

  • 滹沱河读音(滹沱河怎么读)

    下面内容希望能帮助到你,我们来一起看看吧!

  • 一千零一夜故事100篇(一千零一夜)

    学校教务主任询问初当老师的麻小娟。2015年7月,麻小娟从湖南第一师范学院音乐学专业毕业,考进常德市鼎城区周家店中学。第二天麻小娟恢复了活力。《食堂里的那盏灯》2016年下学期,麻小娟调入蔡家岗中学。星期一,蔡家岗中学举行升旗仪式。麻小娟发现学生的红领巾有的脏了,有的戴法不对。第二天课间操,麻老师发现同学们队伍整齐,动作标准,鲜艳的红领巾飘扬在学生们胸前。

  • NT结果能看出胎儿性别吗(孕nt能看出胎儿性别)

    NT检查是通过B超扫描胎儿颈项透明带的检查,根据其数值来预估是否有患唐氏综合症的风险。NT检查结果主要作用是用来评估胎儿畸形的风险高低的,对与鉴别胎儿是没有任何作用的。生男生女健康最重要目前对于孕妇的检查中,虽然可以通过B超或羊水穿刺检查鉴别出胎儿的性别,但是国家规定下,医生是不会直接或者间接向孕妇及家人透露胎儿的性别信息的。