ファイルの2番目の列が以下の数値のいずれかに等しい場合:
65,81,83,97,99,113,129,145,147,161,163,177 #To be read 65 OR 81 OR 83 OR 97 OR 99 etc..
次に、行全体を出力ファイルに出力する必要があるため、次のように入力します。
samtools view myfile.bam | awk '{for (i=65,81,83,97,99,113,129,145,147,161,163,177) if ($2==i) ; print$0} > output.bam
それはうまくいくでしょうか?