假如双端测序得到这两个文件:
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
[email protected]
FAFU
时间: 2024-10-24 09:31:37