您的位置:首页 > 其它

如何计算测序深度

2014-11-27 00:05 218 查看
假如双端测序得到这两个文件:

L2Y-0_TGACCA_L002_R1.fastq.gz

L2Y-0_TGACCA_L002_R2.fastq.gz

如果用wc -l 看gz文件的行数, 返回的结果是不对的

一般用fastqc看下质量怎么样:

/share/hiseq/SampleProcess/FastQC/fastqc 20_GTGGCC_L000_R1.fastq.gz

后面只需要放一个,自动两个都会跑出来

然后查看结果, 有总序列数目:



这个是单端的。

那么覆盖深度就为:

(total_sequences * 2) * 测序的长度 / 基因组的大小

freemao

free_mao@qq.com

FAFU
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: