我想修改这个命令并创建一个命令来过滤带有" value "标志和超过2个"pass"。任何建议吗?
这个命令只能使用一个PASS:
awk '{if(($5=="val") && ($0 ~ /PASS/ )) {print $0}}' sample.vcf
假设$5是标志字段:
awk '{if(($5=="val") && ($0 ~ /^.*PASS.*PASS.*$/ )) {print $0}}' sample.vcf
BTW{if(($5=="val") && ($0 ~ /PASS/ )) {print $0}}
将不匹配任何行,因为if $5 == "val" $0
永远不会匹配~ /PASS/