使用awk
cat test.sam |awk '{if ($0~/^@/){print} else if($0~/NH:i:1\>/){print}}' >test_filter.sam
本文标题：过滤sam文件，留下唯一比对reads
本文链接：https://www.haomeiwen.com/subject/qykfzktx.html
网友评论