map.txt文件有两列:
$ head map.txt
NC_013669.1 NW_024465699.1
NC_013669.1 NW_024464918.1
NC_013669.1 NW_024465023.1
NC_013669.1 NW_024465553.1
NC_013669.1 NW_024465573.1
NC_013669.1 NW_024467106.1
NC_013669.1 NW_024467122.1
NC_013669.1 NW_024466199.1
NC_013669.1 NW_024465561.1
NC_013669.1 NW_024465035.1
需求:
根据第二列匹配到文件GCF_017591425.1_NIBS_Ocur_1.0_genomic.fna中的行,然后用第一列进行替换(前面加>)
#!/bin/bashcat map.txt | while read a1 a2;dorep=`grep -w $a2 GCF_017591425.1_NIBS_Ocur_1.0_genomic.fna`sed -i "s/$rep/>${a1}/g" GCF_017591425.1_NIBS_Ocur_1.0_genomic.fnadone
