このファイルを使用しようとしています...
MT 200610-1 100 2757
MT 200610-10 100 6753
MT 200610-100 100 15173
列 1、2、4 を使用して、次のような新しいファイルを作成したいです...
chrMT 2756 2757 200610-1
chrMT 6752 6753 200610-10
chrMT 15172 15173 200610-100
最初の 3 つだけをこのように実行すればうまくいきます。
awk 'BEGIN{FS=="\t";OFS=="\t"}{print "chr"$1,$4-1,$4}' test_ace_geno1.map
chrMT 2756 2757
chrMT 6752 6753
chrMT 15172 15173
最後に 2 列目を追加しようとすると、奇妙なことが起こります...
awk 'BEGIN{FS=="\t";OFS=="\t"}{print "chr"$1,$4-1,$4,$2}' test_ace_geno1.map
200610-16 2757
200610-10 6753
200610-100 15173
セットリストを使用して vim を実行すると、通常のタブとして表示されるため、全体を台無しにしているその列で何が起こっているのかわかりません。
MT^I200610-1^I100^I2757$
MT^I200610-10^I100^I6753$
MT^I200610-100^I100^I15173$
答え1
dos
スタイル行末付き:
$ cat -e ip.txt
MT 200610-1 100 2757^M$
MT 200610-10 100 6753 ^M$
MT 200610-100 100 15173^M$
$ awk 'BEGIN{FS=="\t";OFS=="\t"}{print "chr"$1,$4-1,$4}' ip.txt
chrMT 2756 2757
chrMT 6752 6753
chrMT 15172 15173
$ awk 'BEGIN{FS=="\t";OFS=="\t"}{print "chr"$1,$4-1,$4,$2}' ip.txt
200610-16 2757
chrMT 6752 6753 200610-10
200610-100 15173
スタイルの行末と必要な出力のperl
両方を処理するには、次のようにします。dos
$ perl -lane 'print "chr$F[0]\t", $F[3]-1, "\t$F[3]\t$F[1]"' ip.txt
chrMT 2756 2757 200610-1
chrMT 6752 6753 200610-10
chrMT 15172 15173 200610-100
では、まずawk
行末をスタイルに変更し、可能な場合は で変更するか、unix
dos2unix
perl
$ perl -i -pe 's|\r\n|\n|' ip.txt
$ cat -e ip.txt
MT 200610-1 100 2757$
MT 200610-10 100 6753 $
MT 200610-100 100 15173$
その後
$ awk -v OFS='\t' '{print "chr"$1,$4-1,$4,$2}' ip.txt
chrMT 2756 2757 200610-1
chrMT 6752 6753 200610-10
chrMT 15172 15173 200610-100
または、 を使ってGNU awk
、入力レコードをDOS形式の行末処理するように設定する
gawk -v RS='\r\n' -v OFS='\t' '{print "chr"$1,$4-1,$4,$2}' ip.txt
答え2
ステートメント内のすべてをprint
希望の形式で入力し、減算などの必要な操作$NF - 1
も実行します。
awk '{print "chr" $1, ($NF - 1), $NF, $2}' file.txt
$4
ここではの代わりにを使用することもできます$NF
が、4 番目のフィールドが最後のフィールドであるため、両方とも同じフィールドを指すことになります。
例:
% awk '{print "chr" $1, ($NF - 1), $NF, $2}' file.txt
chrMT 2756 2757 200610-1
chrMT 6752 6753 200610-10
chrMT 15172 15173 200610-100