だから私は2つのファイルを持っています.1つのVCFは次のようになります
88 Chr1 25 C - 3 2 1 1
88 Chr1 88 A T 7 2 1 1
88 Chr1 92 A C 16 4 1 1
次のような遺伝子を持つ別の
GENEID Start END
GENE_ID 11 155
GENE_ID 165 999
2番目のファイルの2番目と3番目の位置の範囲内に遺伝子の位置(VCFファイルの3列目)があるかどうかを見て、それを印刷するスクリプトが必要です。
私がこれまでに行ったことは、ファイルに参加して行うことでした
awk '{if (3>$12 && $3< $13) print }' > out
私がしたことは、結合されたファイルの現在の行のみを比較することです (値が同じ行にある場合にのみ印刷されます)。列 3 のすべての行を列 12 および 13 のすべての行と比較するにはどうすればよいですか?
ベスト、セルグ