vcf中染色体号列前添加chr并跳过header:
awk '{
if($0 !~ /^#/)
print "chr"$0;
else if(match($0,/(##contig=<ID=)(.*)/,m))
print m[1]"chr"m[2];
else print $0
}' test_nochr.vcf >test_addchr.vcf
vcf中染色体号列前添加chr并跳过header:
awk '{
if($0 !~ /^#/)
print "chr"$0;
else if(match($0,/(##contig=<ID=)(.*)/,m))
print m[1]"chr"m[2];
else print $0
}' test_nochr.vcf >test_addchr.vcf