fastq文件在经过fastqc文件质检后,一般都会生成一个网页版的文件,我们可以根据文件来分析我们的测序结果的好坏,前提是我们能够读懂这个文件中显示的内容,接下来我们主要解读一下每一张图所代表的信息。
1 首先我们看一下左边的summary:绿色代表PASS;黄色代表WARN;红色代表FAIL。当出现黄色时说明需要查看结果。
2 Basic Statistics
Basic statistics是该fastq一些基本信息,主要
Filename:文件名
File type: 文件类型
Encoding:测序平台的版本和相应的编码版本号,用于计算Phred反推error P时用
Total Sequences: 输入文本的reads的数量
Sequence length: 测序长度
%GC: GC含量,表示整体序列的GC含量,由于二代测序GC偏好性高,且深度越高,GC含量会越高。
3.Per base sequence quality
横轴为read长度,纵轴为质量得分,Q = -10*log10(error P)。柱状表示该位置所有序列的测序质量的统计,柱状是25%~75%区间质量分布,error bar是10%~90%区间质量分布,蓝线表示平均数。一般要求所有位置的10%分位数大于20,即大于最多允许该位置10%的序列低于Q20。当任何碱基质量低于10,或者任何中位数低于25报WARN,需注意;当任何碱基质量低于5或者任何中位数低于20报FAIL。
4.Per base sequence content
统计在序列中的每一个位置,四种不同碱基占总碱基数的比例,检测有无AT、GC分离的现象。横轴为位置,纵轴为百分比。正常情况下四种碱基出现的频率应是接近的,且没有位置差异,因此好的样品中四条线应该是平行且接近的,由于刚开始测序仪状态不稳定,造成前几个碱基有波动。在reads 开头出现碱基组成偏离往往是我们的建库操作造成的,比如建 GBS 文库时在 reads 开头加了 barcode;barcode的碱基组成不是均一的,酶切位点的碱基组成是固定不变的,这样会造成明显的碱基组成偏离;在 reads结尾出现的碱基组成偏离,往往是测序接头的污染造成的。当所有位置的碱基比例一致现出偏差时,即四条线平行且分开,代表文库有偏差,或测序中的系统误差;当部分位置碱基的比例出现偏差时,即四条线在某些位置纷乱交织,则有overrepresented sequence的污染。当任一位置的A/T比例与G/C比例相差超过10%,报"WARN";当任一位置的A/T比例与G/C比例相差超过20%,报"FAIL",我这里的数据就不是很好。
5.Per sequence GC content
横轴表示GC含量,纵轴表示不同GC含量对应的read数,蓝线是理论分布(正态分布,通过从所测数据计算并构建理论分布),红色是实际情况,两个比较接近判为好的。曲线形状的偏差往往是由于文库的污染或是部分reads构成的子集有偏差(overrepresentedreads);形状接近正态分布但偏离理论分布的情况提示我们可能有系统偏差;如果出现两个或多个峰值,表明测序数据里可能有其他来源的DNA序列污染,或者有接头序列的二聚体污染。偏离理论分布的reads超过15%时,报"WARN";偏离理论分布的reads超过30%时,报"FAIL"。
6.Per base N content
当出现测序仪不能分辨的碱基时会产生N,横轴为碱基分布,纵轴为N比率,当任一位置N的比率超过5%报WARN,超过20%报FAIL。我这里几乎没有。
7.Sequence Length Distribution
理论上每次测序仪测出的read长度是一致的,但是由于建库等因素通常会导致一些小片段,如果报FAIL,表明此次测序过程中产生的数据不可信。
8.Sequence Duplication Levels
统计序列完全一致的reads的频率,横轴表示重复的次数,纵轴表示重复的reads的数目。一般测序深度越高,越容易产生一定程度的重复序列。
9.Overrepresented sequences
当有某个序列大量出现时,超过总reads数的0.1%时报WARN,超过1%时报FAIL。
10.Adapter Content
横轴表示碱基位置,纵轴表示百分比。当fastqc分析时没有选择参数-a adapter list时,默认使用图例中的4种通用adapter序列进行统计。若有adapter残留,后续必须去接头。