如何解决在UNIX中从另一个文件中查找一个文件的内容并保持相同的顺序
我有一个文件,例如:
Colcht.WP_006104309.1
Moopro.WP_070396948.1
Mastes.WP_027845098.1
Phowil.WP_068791039.1
Cyaapo.WP_015218744.1
Gemher.WP_017295578.1
Lyncon.WP_039726659.1
Glokil.WP_023171831.1
Noscya.WP_087539356.1
Photen.WP_073607805.1
Hydriv.WP_073598454.1
Lepoha.WP_088893428.1
Nodnod.WP_017300904.1
Noscya.WP_087540001.1
Spisub.WP_017307136.1
Scy0HK.WP_073635112.1
PlaSR0.WP_054467905.1
我使用grep -A 1 -F -f
从另一个文件进行grep操作,有关上述文本文件中每一行的信息。结果的一部分如下:
>Cyaapo.WP_015218744.1
MNIVVVGLSHKTAAVEVREKLSIPEAKIEDSIRHLLTYPHIEEVAIISTCNRLEIYAVVKETEQGVKEITQFLAEIGNLALLELRRHLFILLHQDAIRHLMRVAAGLESLVLGEGQILAQVKTTHKLGMKYNGMSRLLDRLFKQAISAGKRVRTETNIGTGAVSISSAAVELVDTKIEDLSSQKVTIIGAGKMSRLLVQHLLAKGVEDIIIVNRSHNRSQELAKQFPQANLKLNLLEDMMTMVAQSDIVFTSTGATQPILDKNNLSSLSINHSLMLVDISVPRNVASDVTELEFIKSYNVDDLKAVVAQNHASRREMAREAENLLEEEIEAFELWWQSLETVPTISCLRSKIEQIREQELEKALSRLGSEFAEKHQEVIEALTRGIVNKILHDPMVQLRTQKDIEARRKALHTLQTLFDLEVSEQLI
>Gemher.WP_017295578.1
MNIIVVGLSHKTAAVEVREKLSIPEAKIEESIKHLLSYPHIEEVGIISTCNRLEIYAVVKETEQGVKEITQFLAEIGHLSLHSLRRHLFILLHQDAIRHLMRVSAGLESLVLGEGQILAQVKNTHKLSTKYQGMGRLLDRLFKQAMSAGKRVRTETNIGTGAVSISSAAVELVDMKLDDLSRQKVSIIGAGKMSRLLVQHLISKGVSDITIVNRSVSRSKELAKQFPQIELKLNLLEEMMEIVRDSDIVFTSTGATQPILDKNNLCSIECHHSIMLVDISVPRNVASDVEELDFIVAYNVDDLKAVVAQNQASRREMAREAELLLEEEIEAFELWWQSLETVPTISCLRSKIEEIREQELEKALSRLGSEFAEKHQEVIEALTRGIVNKILHDPMVQLRSQQDIEARRKALQTLQTLFNLEISEQFG
>Glokil.WP_023171831.1
MQIAVIGLSHRTAPVEIREKVSIPEQQVAEYVSRLRSCSQIAECAILSTCNRLEIYAVLRDSEHGLREVTQFLAESKGVALPMLQRHLFTLLHQDAVMHLMRVAAGLDSIVLGEGQILAQIKVTHKLAQQGKGVDRILNQLFKAAITGGKRVREETDIGKGAISVSSAAVEMAMRKKNRRSLQDQRCLVVGAGKMGELVLRHLISKGARQIIVLNRSLEKAAQMVEQFGLMLPVATIDELGNHLGAADLIFTCTSASEPLINYERLSQVRREQPLMIFDIAVPRNVAVDVEELSNVHLFNVDHLKQVVEENRAYRQLMVQQCEDILLQQLDEFLDWWRNLEAVPTINSLRQKVETIREQELEKALSRLGTEFGEKHQGIIDSLTRAIVNKILHDPMVQLRAQRDVEARRRALQTLQTLFNLEPLGSNPEPPVL
>Hydriv.WP_073598454.1
MNIAVVGLSHKTAPVEIREKLSIQEAKLESALAHLRSYPHIIEVAIISTCNRLEIYAIATETDQGVREISQFLSEIGHIPLDRLRRYLFILLHQDAVRHLMRVAAGLESLVLGEGQILAQVKNTHKLAQKYQSLGQILDRLFKQAMTAGKRVRSETNIGTGAVSISSAAVELAHMKAENLAARRVCIIGAGKMSRLLVQHLLAKGTQQICIVNRSHRRAEELASQFPEVQLKLYPLTEMMSAVAASDIVFTSTAATEPIINRSQLEASLTRDRELMLFDISVPRNVHADVGGMESVQSYNVDDLKAVVAQNYESRRKMAQEAEALLEEEIAAFELWWRSLETVPTISCLRSKVETIREQELEKALSRLGTEFAEKHQEVIEALTRGIVNKILHEPMVQLRAQQDIEARRRCLQSLQMLFNLEIEKQVI
>Lepoha.WP_088893428.1
MNIVVVGLSHRTAPVEVREKLSIPTPQMEAAIAHLRSFPHIEEATILSTCNRLEVYVVTSETEQGVREVTQFLSEYGKISVSQLRPYLFILLHDDAVMHLMRVSAGLDSLVLGEGQILAQVKHTHKVGQQYNGIGRILNRLFNQAITAGKRVRTETSIGTGAVSISSAAVELAQLKVQHLPACRVAILGAGKMSRLVVQHLISKGATQICIVNRSLDSARELAQQFKEAEIRLHLLDEMMHVICNSDLVFTATAATEPLIDRAKLESTIDPLHSLKLFDISVPRNVHADVNELDHVQLFNVDDLKAVVAQNQESRRQMALEAENILDEEVAAFDLWWRSLETVATISELRDKVEAIRAQELEKALSRLGSEFAEKHQEVIEALTRGIVNKILHDPMVQLRAQQDIEARRRAMQTLRSLFNLEEPASNKA
但是,如您所见,第一个文件(1st: Colcht.WP_006104309.1,2nd: Moopro.WP_070396948.1,3rd: Mastes.WP_027845098.1
)中的顺序不再受尊重。
预期输出的子集为:
>Colcht.WP_006104309.1
MENYNTSNIDNVLLLKGDDIINLFKNREQDILDLVKLTYKIHGRGDSTLPHSSFLRFPDKNKERIIALPAYLGGEINTAGIKWIASFPGNLARGMERASAILIINSTETGRPQAIMEGSVISAKRTAASAALAAHFLRDRQSLVTVGLIGCGLINFETVRFLLKVRPEIETLFLYDVSLEKSDQFKRKCQQLSQNRELVILDNPDDVFKHSSVIALATTASQPHIVDISACQSDSIILHTSLRDLSPEIILSVDNIVDDIDHVCRAQTSIHLAEQKTGNRDFIRCPLSDILNGVAAPRQNNSQIAVFSPFGLGVLDLALGQLAYQLADETNVGTRLTSFFPVSWLQREDE
>Moopro.WP_070396948.1
METAYQGFAQQQPGDVIVLSASDILSLLAGREKELIEVVRQTYIAHARGESALPPSPFLRFANHPKNRIIAKPAYLGESFETAGIKWISSFPDNYQFGLLRASAVIILNSVKTGFAEAILEGSVISAKRTAASAALAARLLQSETQPESIGIIACGVINFEITRFLLAEFPTVKNLVIFDIDHERAVQYKSRCETNFETPNITIANDINTVLSSTSIISIATTETTPHIFEISACQPGSNILHISLRDFSPEVILSCDNIVDDVEHICSAQTSVHLAEQKINHRHFIRGSIGDILCGKILAKPTPSAITIFSPFGLGILDLAVAKLVHEWGIARNLGTVIPSFGCLPHE
>Mastes.WP_027845098.1
MSNKHHLSFTYLSQEDLLDAGCFDIRMVMDIAEKAMLEFERHHVIFPEKIVQIFNQATQERINCLPATLLDEKVCGVKWVSVFPMNPIEHDQQNLSAIFILSEIETGFPICVMEGTLASNMRVAAIGGLAAKYLARQDSEVIGFIGAGEQAKMHLIAMKAVCPSLKQCRVAAHVVKHEEQFIAELSRLYPEMEFVSTNTNLQKAIEDADILVTATSAQAPLLKATWVKPGTFYSHIGGWEDEFEVALQADKIVCDDWETVTHRTQTLSRMYQEGLINANNIHADLHELVSGKKAGRESQTERIYFNAVGLAYIDIAIAMAMFNRAREKQKGTQLDLQQSMVFEHLGLKSKVKL
>Phowil.WP_068791039.1
MRVISAAEVQAALDFETLVGRLRDIFRRGGEAPARQQYDIAITGEPAQTLLLAPAWQAGRHVGVQIATVTPGNGARGLPAGMGAYLLLDGRSGAPAALIDGPMLTLRRTAAASALASAYLSRPDSARLLMVGTGALAPHLIAAHAAVRPIREVLVWGRTPAKAARLAKAVKLPRVRLAWTEDLEGAVRGADIVACATLSQQPLLRGAWLRPGQHLDLVGAYRPEMRESDGEVFRRARVFVDTRAGALAEAGDLIQALAEGALSAADVAADLFELARGEKAGRRFYDQITLFKATGSALEDLAAAQLTVERA
>Cyaapo.WP_015218744.1
MNIVVVGLSHKTAAVEVREKLSIPEAKIEDSIRHLLTYPHIEEVAIISTCNRLEIYAVVKETEQGVKEITQFLAEIGNLALLELRRHLFILLHQDAIRHLMRVAAGLESLVLGEGQILAQVKTTHKLGMKYNGMSRLLDRLFKQAISAGKRVRTETNIGTGAVSISSAAVELVDTKIEDLSSQKVTIIGAGKMSRLLVQHLLAKGVEDIIIVNRSHNRSQELAKQFPQANLKLNLLEDMMTMVAQSDIVFTSTGATQPILDKNNLSSLSINHSLMLVDISVPRNVASDVTELEFIKSYNVDDLKAVVAQNHASRREMAREAENLLEEEIEAFELWWQSLETVPTISCLRSKIEQIREQELEKALSRLGSEFAEKHQEVIEALTRGIVNKILHDPMVQLRTQKDIEARRKALHTLQTLFDLEVSEQLI
有人知道我如何像上面一样grep保持第一个文件的顺序吗?
非常感谢您的帮助:)
解决方法
如果要保留第一个文件的顺序,最好解析第二个文件,存储下一行,然后解析第一个文件。
> cat tst.awk
FNR==NR && p {
a[prev]=$0
p=0
next
}
FNR==NR && $0~/^>/ {
prev=substr($0,2)
p=1
next
}
$0 in a {
print ">" $0 RS a[$0]
}
用法:
awk -f tst.awk file2 file1
如果file2
很大并且没有足够的内存,则可以将file2
替换为grep
命令的输出(仅{{1}的一部分) }。
file2
否则,您仍然可以传递awk -f tst.awk <(grep -A 1 -f file1 file2) file1
,但必须保存行的顺序,并在file1 file2
部分中完成工作。
END
用法:
> cat tst.awk
FNR==NR {
row[NR]=$0
a[$0]
next
}
p {
next_row[x]=$0
p=0
next
}
substr($0,2) in a {
x=substr($0,2)
p=1
}
END {
for (i=1;i in row;i++)
if (next_row[row[i]])
print ">" row[i] RS next_row[row[i]]
}
,
$ awk '
NR==FNR { if (/^>/) tag=$1; else a[substr(tag,2)]=tag ORS $0; next }
$1 in a { print a[$1] }
' file2 file1
>Colcht.WP_006104309.1
MENYNTSNIDNVLLLKGDDIINLFKNREQDILDLVKLTYKIHGRGDSTLPHSSFLRFPDKNKERIIALPAYLGGEINTAGIKWIASFPGNLARGMERASAILIINSTETGRPQAIMEGSVISAKRTAASAALAAHFLRDRQSLVTVGLIGCGLINFETVRFLLKVRPEIETLFLYDVSLEKSDQFKRKCQQLSQNRELVILDNPDDVFKHSSVIALATTASQPHIVDISACQSDSIILHTSLRDLSPEIILSVDNIVDDIDHVCRAQTSIHLAEQKTGNRDFIRCPLSDILNGVAAPRQNNSQIAVFSPFGLGVLDLALGQLAYQLADETNVGTRLTSFFPVSWLQREDE
>Moopro.WP_070396948.1
METAYQGFAQQQPGDVIVLSASDILSLLAGREKELIEVVRQTYIAHARGESALPPSPFLRFANHPKNRIIAKPAYLGESFETAGIKWISSFPDNYQFGLLRASAVIILNSVKTGFAEAILEGSVISAKRTAASAALAARLLQSETQPESIGIIACGVINFEITRFLLAEFPTVKNLVIFDIDHERAVQYKSRCETNFETPNITIANDINTVLSSTSIISIATTETTPHIFEISACQPGSNILHISLRDFSPEVILSCDNIVDDVEHICSAQTSVHLAEQKINHRHFIRGSIGDILCGKILAKPTPSAITIFSPFGLGILDLAVAKLVHEWGIARNLGTVIPSFGCLPHE
>Mastes.WP_027845098.1
MSNKHHLSFTYLSQEDLLDAGCFDIRMVMDIAEKAMLEFERHHVIFPEKIVQIFNQATQERINCLPATLLDEKVCGVKWVSVFPMNPIEHDQQNLSAIFILSEIETGFPICVMEGTLASNMRVAAIGGLAAKYLARQDSEVIGFIGAGEQAKMHLIAMKAVCPSLKQCRVAAHVVKHEEQFIAELSRLYPEMEFVSTNTNLQKAIEDADILVTATSAQAPLLKATWVKPGTFYSHIGGWEDEFEVALQADKIVCDDWETVTHRTQTLSRMYQEGLINANNIHADLHELVSGKKAGRESQTERIYFNAVGLAYIDIAIAMAMFNRAREKQKGTQLDLQQSMVFEHLGLKSKVKL
>Phowil.WP_068791039.1
MRVISAAEVQAALDFETLVGRLRDIFRRGGEAPARQQYDIAITGEPAQTLLLAPAWQAGRHVGVQIATVTPGNGARGLPAGMGAYLLLDGRSGAPAALIDGPMLTLRRTAAASALASAYLSRPDSARLLMVGTGALAPHLIAAHAAVRPIREVLVWGRTPAKAARLAKAVKLPRVRLAWTEDLEGAVRGADIVACATLSQQPLLRGAWLRPGQHLDLVGAYRPEMRESDGEVFRRARVFVDTRAGALAEAGDLIQALAEGALSAADVAADLFELARGEKAGRRFYDQITLFKATGSALEDLAAAQLTVERA
>Cyaapo.WP_015218744.1
MNIVVVGLSHKTAAVEVREKLSIPEAKIEDSIRHLLTYPHIEEVAIISTCNRLEIYAVVKETEQGVKEITQFLAEIGNLALLELRRHLFILLHQDAIRHLMRVAAGLESLVLGEGQILAQVKTTHKLGMKYNGMSRLLDRLFKQAISAGKRVRTETNIGTGAVSISSAAVELVDTKIEDLSSQKVTIIGAGKMSRLLVQHLLAKGVEDIIIVNRSHNRSQELAKQFPQANLKLNLLEDMMTMVAQSDIVFTSTGATQPILDKNNLSSLSINHSLMLVDISVPRNVASDVTELEFIKSYNVDDLKAVVAQNHASRREMAREAENLLEEEIEAFELWWQSLETVPTISCLRSKIEQIREQELEKALSRLGSEFAEKHQEVIEALTRGIVNKILHDPMVQLRTQKDIEARRKALHTLQTLFDLEVSEQLI
>Gemher.WP_017295578.1
MNIIVVGLSHKTAAVEVREKLSIPEAKIEESIKHLLSYPHIEEVGIISTCNRLEIYAVVKETEQGVKEITQFLAEIGHLSLHSLRRHLFILLHQDAIRHLMRVSAGLESLVLGEGQILAQVKNTHKLSTKYQGMGRLLDRLFKQAMSAGKRVRTETNIGTGAVSISSAAVELVDMKLDDLSRQKVSIIGAGKMSRLLVQHLISKGVSDITIVNRSVSRSKELAKQFPQIELKLNLLEEMMEIVRDSDIVFTSTGATQPILDKNNLCSIECHHSIMLVDISVPRNVASDVEELDFIVAYNVDDLKAVVAQNQASRREMAREAELLLEEEIEAFELWWQSLETVPTISCLRSKIEEIREQELEKALSRLGSEFAEKHQEVIEALTRGIVNKILHDPMVQLRSQQDIEARRKALQTLQTLFNLEISEQFG
>Glokil.WP_023171831.1
MQIAVIGLSHRTAPVEIREKVSIPEQQVAEYVSRLRSCSQIAECAILSTCNRLEIYAVLRDSEHGLREVTQFLAESKGVALPMLQRHLFTLLHQDAVMHLMRVAAGLDSIVLGEGQILAQIKVTHKLAQQGKGVDRILNQLFKAAITGGKRVREETDIGKGAISVSSAAVEMAMRKKNRRSLQDQRCLVVGAGKMGELVLRHLISKGARQIIVLNRSLEKAAQMVEQFGLMLPVATIDELGNHLGAADLIFTCTSASEPLINYERLSQVRREQPLMIFDIAVPRNVAVDVEELSNVHLFNVDHLKQVVEENRAYRQLMVQQCEDILLQQLDEFLDWWRNLEAVPTINSLRQKVETIREQELEKALSRLGTEFGEKHQGIIDSLTRAIVNKILHDPMVQLRAQRDVEARRRALQTLQTLFNLEPLGSNPEPPVL
>Hydriv.WP_073598454.1
MNIAVVGLSHKTAPVEIREKLSIQEAKLESALAHLRSYPHIIEVAIISTCNRLEIYAIATETDQGVREISQFLSEIGHIPLDRLRRYLFILLHQDAVRHLMRVAAGLESLVLGEGQILAQVKNTHKLAQKYQSLGQILDRLFKQAMTAGKRVRSETNIGTGAVSISSAAVELAHMKAENLAARRVCIIGAGKMSRLLVQHLLAKGTQQICIVNRSHRRAEELASQFPEVQLKLYPLTEMMSAVAASDIVFTSTAATEPIINRSQLEASLTRDRELMLFDISVPRNVHADVGGMESVQSYNVDDLKAVVAQNYESRRKMAQEAEALLEEEIAAFELWWRSLETVPTISCLRSKVETIREQELEKALSRLGTEFAEKHQEVIEALTRGIVNKILHEPMVQLRAQQDIEARRRCLQSLQMLFNLEIEKQVI
>Lepoha.WP_088893428.1
MNIVVVGLSHRTAPVEVREKLSIPTPQMEAAIAHLRSFPHIEEATILSTCNRLEVYVVTSETEQGVREVTQFLSEYGKISVSQLRPYLFILLHDDAVMHLMRVSAGLDSLVLGEGQILAQVKHTHKVGQQYNGIGRILNRLFNQAITAGKRVRTETSIGTGAVSISSAAVELAQLKVQHLPACRVAILGAGKMSRLVVQHLISKGATQICIVNRSLDSARELAQQFKEAEIRLHLLDEMMHVICNSDLVFTATAATEPLIDRAKLESTIDPLHSLKLFDISVPRNVHADVNELDHVQLFNVDDLKAVVAQNQESRRQMALEAENILDEEVAAFDLWWRSLETVATISELRDKVEAIRAQELEKALSRLGSEFAEKHQEVIEALTRGIVNKILHDPMVQLRAQQDIEARRRAMQTLRSLFNLEEPASNKA
以上内容在这些输入文件上运行:
$ head -n 50 file*
==> file1 <==
Colcht.WP_006104309.1
Moopro.WP_070396948.1
Mastes.WP_027845098.1
Phowil.WP_068791039.1
Cyaapo.WP_015218744.1
Gemher.WP_017295578.1
Lyncon.WP_039726659.1
Glokil.WP_023171831.1
Noscya.WP_087539356.1
Photen.WP_073607805.1
Hydriv.WP_073598454.1
Lepoha.WP_088893428.1
Nodnod.WP_017300904.1
Noscya.WP_087540001.1
Spisub.WP_017307136.1
Scy0HK.WP_073635112.1
PlaSR0.WP_054467905.1
==> file2 <==
>Cyaapo.WP_015218744.1
MNIVVVGLSHKTAAVEVREKLSIPEAKIEDSIRHLLTYPHIEEVAIISTCNRLEIYAVVKETEQGVKEITQFLAEIGNLALLELRRHLFILLHQDAIRHLMRVAAGLESLVLGEGQILAQVKTTHKLGMKYNGMSRLLDRLFKQAISAGKRVRTETNIGTGAVSISSAAVELVDTKIEDLSSQKVTIIGAGKMSRLLVQHLLAKGVEDIIIVNRSHNRSQELAKQFPQANLKLNLLEDMMTMVAQSDIVFTSTGATQPILDKNNLSSLSINHSLMLVDISVPRNVASDVTELEFIKSYNVDDLKAVVAQNHASRREMAREAENLLEEEIEAFELWWQSLETVPTISCLRSKIEQIREQELEKALSRLGSEFAEKHQEVIEALTRGIVNKILHDPMVQLRTQKDIEARRKALHTLQTLFDLEVSEQLI
>Gemher.WP_017295578.1
MNIIVVGLSHKTAAVEVREKLSIPEAKIEESIKHLLSYPHIEEVGIISTCNRLEIYAVVKETEQGVKEITQFLAEIGHLSLHSLRRHLFILLHQDAIRHLMRVSAGLESLVLGEGQILAQVKNTHKLSTKYQGMGRLLDRLFKQAMSAGKRVRTETNIGTGAVSISSAAVELVDMKLDDLSRQKVSIIGAGKMSRLLVQHLISKGVSDITIVNRSVSRSKELAKQFPQIELKLNLLEEMMEIVRDSDIVFTSTGATQPILDKNNLCSIECHHSIMLVDISVPRNVASDVEELDFIVAYNVDDLKAVVAQNQASRREMAREAELLLEEEIEAFELWWQSLETVPTISCLRSKIEEIREQELEKALSRLGSEFAEKHQEVIEALTRGIVNKILHDPMVQLRSQQDIEARRKALQTLQTLFNLEISEQFG
>Glokil.WP_023171831.1
MQIAVIGLSHRTAPVEIREKVSIPEQQVAEYVSRLRSCSQIAECAILSTCNRLEIYAVLRDSEHGLREVTQFLAESKGVALPMLQRHLFTLLHQDAVMHLMRVAAGLDSIVLGEGQILAQIKVTHKLAQQGKGVDRILNQLFKAAITGGKRVREETDIGKGAISVSSAAVEMAMRKKNRRSLQDQRCLVVGAGKMGELVLRHLISKGARQIIVLNRSLEKAAQMVEQFGLMLPVATIDELGNHLGAADLIFTCTSASEPLINYERLSQVRREQPLMIFDIAVPRNVAVDVEELSNVHLFNVDHLKQVVEENRAYRQLMVQQCEDILLQQLDEFLDWWRNLEAVPTINSLRQKVETIREQELEKALSRLGTEFGEKHQGIIDSLTRAIVNKILHDPMVQLRAQRDVEARRRALQTLQTLFNLEPLGSNPEPPVL
>Hydriv.WP_073598454.1
MNIAVVGLSHKTAPVEIREKLSIQEAKLESALAHLRSYPHIIEVAIISTCNRLEIYAIATETDQGVREISQFLSEIGHIPLDRLRRYLFILLHQDAVRHLMRVAAGLESLVLGEGQILAQVKNTHKLAQKYQSLGQILDRLFKQAMTAGKRVRSETNIGTGAVSISSAAVELAHMKAENLAARRVCIIGAGKMSRLLVQHLLAKGTQQICIVNRSHRRAEELASQFPEVQLKLYPLTEMMSAVAASDIVFTSTAATEPIINRSQLEASLTRDRELMLFDISVPRNVHADVGGMESVQSYNVDDLKAVVAQNYESRRKMAQEAEALLEEEIAAFELWWRSLETVPTISCLRSKVETIREQELEKALSRLGTEFAEKHQEVIEALTRGIVNKILHEPMVQLRAQQDIEARRRCLQSLQMLFNLEIEKQVI
>Lepoha.WP_088893428.1
MNIVVVGLSHRTAPVEVREKLSIPTPQMEAAIAHLRSFPHIEEATILSTCNRLEVYVVTSETEQGVREVTQFLSEYGKISVSQLRPYLFILLHDDAVMHLMRVSAGLDSLVLGEGQILAQVKHTHKVGQQYNGIGRILNRLFNQAITAGKRVRTETSIGTGAVSISSAAVELAQLKVQHLPACRVAILGAGKMSRLVVQHLISKGATQICIVNRSLDSARELAQQFKEAEIRLHLLDEMMHVICNSDLVFTATAATEPLIDRAKLESTIDPLHSLKLFDISVPRNVHADVNELDHVQLFNVDDLKAVVAQNQESRRQMALEAENILDEEVAAFDLWWRSLETVATISELRDKVEAIRAQELEKALSRLGSEFAEKHQEVIEALTRGIVNKILHDPMVQLRAQQDIEARRRAMQTLRSLFNLEEPASNKA
>Colcht.WP_006104309.1
MENYNTSNIDNVLLLKGDDIINLFKNREQDILDLVKLTYKIHGRGDSTLPHSSFLRFPDKNKERIIALPAYLGGEINTAGIKWIASFPGNLARGMERASAILIINSTETGRPQAIMEGSVISAKRTAASAALAAHFLRDRQSLVTVGLIGCGLINFETVRFLLKVRPEIETLFLYDVSLEKSDQFKRKCQQLSQNRELVILDNPDDVFKHSSVIALATTASQPHIVDISACQSDSIILHTSLRDLSPEIILSVDNIVDDIDHVCRAQTSIHLAEQKTGNRDFIRCPLSDILNGVAAPRQNNSQIAVFSPFGLGVLDLALGQLAYQLADETNVGTRLTSFFPVSWLQREDE
>Moopro.WP_070396948.1
METAYQGFAQQQPGDVIVLSASDILSLLAGREKELIEVVRQTYIAHARGESALPPSPFLRFANHPKNRIIAKPAYLGESFETAGIKWISSFPDNYQFGLLRASAVIILNSVKTGFAEAILEGSVISAKRTAASAALAARLLQSETQPESIGIIACGVINFEITRFLLAEFPTVKNLVIFDIDHERAVQYKSRCETNFETPNITIANDINTVLSSTSIISIATTETTPHIFEISACQPGSNILHISLRDFSPEVILSCDNIVDDVEHICSAQTSVHLAEQKINHRHFIRGSIGDILCGKILAKPTPSAITIFSPFGLGILDLAVAKLVHEWGIARNLGTVIPSFGCLPHE
>Mastes.WP_027845098.1
MSNKHHLSFTYLSQEDLLDAGCFDIRMVMDIAEKAMLEFERHHVIFPEKIVQIFNQATQERINCLPATLLDEKVCGVKWVSVFPMNPIEHDQQNLSAIFILSEIETGFPICVMEGTLASNMRVAAIGGLAAKYLARQDSEVIGFIGAGEQAKMHLIAMKAVCPSLKQCRVAAHVVKHEEQFIAELSRLYPEMEFVSTNTNLQKAIEDADILVTATSAQAPLLKATWVKPGTFYSHIGGWEDEFEVALQADKIVCDDWETVTHRTQTLSRMYQEGLINANNIHADLHELVSGKKAGRESQTERIYFNAVGLAYIDIAIAMAMFNRAREKQKGTQLDLQQSMVFEHLGLKSKVKL
>Phowil.WP_068791039.1
MRVISAAEVQAALDFETLVGRLRDIFRRGGEAPARQQYDIAITGEPAQTLLLAPAWQAGRHVGVQIATVTPGNGARGLPAGMGAYLLLDGRSGAPAALIDGPMLTLRRTAAASALASAYLSRPDSARLLMVGTGALAPHLIAAHAAVRPIREVLVWGRTPAKAARLAKAVKLPRVRLAWTEDLEGAVRGADIVACATLSQQPLLRGAWLRPGQHLDLVGAYRPEMRESDGEVFRRARVFVDTRAGALAEAGDLIQALAEGALSAADVAADLFELARGEKAGRRFYDQITLFKATGSALEDLAAAQLTVERA
>Cyaapo.WP_015218744.1
MNIVVVGLSHKTAAVEVREKLSIPEAKIEDSIRHLLTYPHIEEVAIISTCNRLEIYAVVKETEQGVKEITQFLAEIGNLALLELRRHLFILLHQDAIRHLMRVAAGLESLVLGEGQILAQVKTTHKLGMKYNGMSRLLDRLFKQAISAGKRVRTETNIGTGAVSISSAAVELVDTKIEDLSSQKVTIIGAGKMSRLLVQHLLAKGVEDIIIVNRSHNRSQELAKQFPQANLKLNLLEDMMTMVAQSDIVFTSTGATQPILDKNNLSSLSINHSLMLVDISVPRNVASDVTELEFIKSYNVDDLKAVVAQNHASRREMAREAENLLEEEIEAFELWWQSLETVPTISCLRSKIEQIREQELEKALSRLGSEFAEKHQEVIEALTRGIVNKILHDPMVQLRTQKDIEARRKALHTLQTLFDLEVSEQLI
,
您可以尝试仅用显示的示例编写的以下内容吗?其中Input_file2是其中一行从>
开始的文件。我尚未测试它是否可以工作。
awk '
FNR==NR{
a[$0]
next
}
/^>/ && (substr($0,2) in a){
found=1
val=$0
next
}
found!=""{
print val ORS $0
found=""
}' Input_file1 Input_file2
说明: 添加以上详细说明。
awk ' ##starting awk program from here.
FNR==NR{ ##Checking condition if FNR==NR which is true when Input_file1 is being read.
a[$0] ##Creating array a with index of current line here.
next ##next will skip all further statements from here.
} ##Closing FNR==NR block here.
/^>/ && (substr($0,2) in a){ ##Checking condition if line starts from > and sub-string from 2nd character to till all line is present in array a then do following.
found=1 ##Setting found to 1 here.
val=$0 ##Creating val which has current line value here.
next ##next will skip all further statements from here.
}
found!=""{ ##Checking if found is SET then do following.
print val ORS $0 ##Printing variable val ORS and current Line here.
found="" ##Nullifying found here.
}' Input_file1 Input_file2 ##Mentioning Input_file names here.
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。