我有多个包含信息的文本文件。这里我展示了两个文本文件,如下所示:
Sample1.txt
Status /documents/Sample1.sorted.bam
Assigned 50945040
Unassigned_Unmapped 947866
Unassigned_MappingQuality 0
Unassigned_Chimera 0
Unassigned_FragmentLength 0
Unassigned_Duplicate 0
Unassigned_MultiMapping 49013681
Unassigned_Secondary 0
Unassigned_Nonjunction 0
Unassigned_NoFeatures 21189312
Unassigned_Overlapping_Length 0
Unassigned_Ambiguity 4430011Sample2.txt
Status /documents/Sample2.sorted.bam
Assigned 36335614
Unassigned_Unmapped 870456
Unassigned_MappingQuality 0
Unassigned_Chimera 0
Unassigned_FragmentLength 0
Unassigned_Duplicate 0
Unassigned_MultiMapping 68688141
Unassigned_Secondary 0
Unassigned_Nonjunction 0
Unassigned_NoFeatures 23746485
Unassigned_Overlapping_Length 0
Unassigned_Ambiguity 3734593对于单个文本文件,我使用grep:
grep "Assigned\|Unmapped\|MultiMapping\|NoFeatures\|Ambiguity" Sample1.txt > output.txt但我希望输出如下所示,我可以使用一个小脚本对所有文本文件和制表:
Sample1 Sample2
Assigned 50945040 36335614
Unassigned_Unmapped 947866 870456
Unassigned_MultiMapping 49013681 68688141
Unassigned_NoFeatures 21189312 23746485
Unassigned_Ambiguity 4430011 3734593发布于 2018-05-27 23:33:08
$ cat tst.awk
$2 != 0 {
printf "%s%s", (NR>1 ? $1 : "Name"), OFS
for (i=2; i<=NF; i+=2) {
gsub(/^.*\/|\..*$/,"",$i)
printf "%s%s", $i, (i<NF ? OFS : ORS)
}
}
$ paste Sample1.txt Sample2.txt | awk -f tst.awk | column -t
Name Sample1 Sample2
Assigned 50945040 36335614
Unassigned_Unmapped 947866 870456
Unassigned_MultiMapping 49013681 68688141
Unassigned_NoFeatures 21189312 23746485
Unassigned_Ambiguity 4430011 3734593要获得Excel可以理解的输出,而不是问题中显示的输出,请执行以下操作:
$ cat tst.awk
BEGIN { OFS="," }
$2 != 0 {
printf "%s%s", (NR>1 ? $1 : "Name"), OFS
for (i=2; i<=NF; i+=2) {
gsub(/^.*\/|\..*$/,"",$i)
printf "%s%s", $i, (i<NF ? OFS : ORS)
}
}
$ paste Sample1.txt Sample2.txt | awk -f tst.awk > output.csv然后双击output.csv将其打开。
https://stackoverflow.com/questions/50552072
复制相似问题