gunzip -c my.vcf.gz | awk -F ' ' '/^#/{print; next;} {n=0; for(i=10;i<=NF;++i) { if($i!="." && index($i,"0/0")==0) n++;} if(n==1) print; } '
没有评论:
发表评论