- 1、本文档共16页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
一、生物信息分析流程
获得原始测序序列(SequencedReads)后,在有相关物种参考序列或参考基因组的情
况下,通过如下流程进行生物信息分析:
二、项目结果说明
1原始序列数据
TM
高通量测序(如illuminaHiSeq2000/MiSeq等测序平台)测序得到的原始图像数据文
件经碱基识别(BaseCalling)分析转化为原始测序序列(SequencedReads),我们称之为
RawData或RawReads,结果以FASTQ(简称为fq)文件格式存储,其中包含测序序列(reads)
的序列信息以及其对应的测序质量信息。
FASTQ格式文件中每个read由四行描述,如下:
@EAS139:136:FC706VJ:2:2104:15343:1973931:Y:18:ATCACG
GCTCTTTGCCCTTCTCGTCGAAAATTGTCTCCTCATTCGAAACTTCTCTGT
+
@@CFFFDEHHHHFIJJJ@FHGIIIEHIIJBHHHIJJEGIIJJIGHIGHCCF
其中第一行以“@”开头,随后为illumina测序标识符(SequenceIdentifiers)和描
述文字(选择性部分);第二行是碱基序列;第三行以“+”开头,随后为illumina测序标
识符(选择性部分);第四行是对应序列的测序质量(Cocketal.)。
illumina测序标识符详细信息如下:
EAS139Uniqueinstrumentname
136RunID
FC706VJFlowcellID
2Flowcelllane
2104Tilenumberwithintheflowcelllane
15343x-coordinateoftheclusterwithinthetile
197393y-coordinateoftheclusterwithinthetile
1Memberofapair,1or2(paired-endormate-pairreadsonly)
YYifthereadfailsfilter(readisbad),Notherwise
180whennoneofthecontrolbitsareon,otherwiseitisanevennumber
ATCACGIndexsequence
第四行中每个字符对应的ASCII值减去33,即为对应第二行碱基的测序质量值。如果
TM
测序错误率用e表示,illuminaHiSeq2000/MiSeq的碱基质量值用Q表示,则有下列
phred
关系:
文档评论(0)