0

だから私は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 のすべての行と比較するにはどうすればよいですか?

ベスト、セルグ

4

1 に答える 1