bcftools mpileup -Ou -f 02428.genome.fa I662_sorted_dedup_reads_realigned.bam.recal.bam| bcftools call -c| vcfutils.pl vcf2fq -d 10 -D 100 | gzip > 3dilpoid.fq.gz
做PSMC的准备文件vcfutils.pl vcf2fq总是报错Use of uninitialized value $l in numeric lt (<) at vcfutils.pl line 566.