vqslod

[shajoezhu]

for i in $(seq 1 1 14); do bcftools view -f "PASS" --include 'VQSLOD > 5' PG0389-C.wg.vcf.gz | grep -v "#" | grep Pf3D7_$(printf %02d $i)_v3 | wc -l; done