fasta中如何根据相同的序列号将序列合并

小小脑

>super821_319450_333026
LMDELRGAQFLSKIDLRSGYNQIRVREQDIPKTAFRCHYGHFEFLVMPFGLTNAPATFQS
CMNHILRRHLRKFLLVFFDNILIYSRTWEEHFQHIEAVLHILEEQQFYAKLSKCEFGLTE
KLYLGHI
>super8297_1747273_1757914
RIDDLFDQMKGATVFSKIDLRFGYHQLRIKEDDIPKTAFQTXFGHYEFVVVPFGLTNAPI
VFTSLMNGVFKKYLDQFVQVFLDGILVYSKNKKEHEEHLRVVLTCLRKNQLYEKLSKRSF
FQNKIHYLGHII
>super821_319450_333026
NSAWEIIPRSEGKTVVGSRWIYKVKQAADESVEKYKARFVAXGFSXIEGIDYEETFAPVA
RYSSIRTILSLSAQMGWHXXXXXXXXXXXQMDVKTAFLNGVIEEEVYIEQPEGFEIFSSE
SHVCRLKIALYGLKQAPHAWYTRIDSYFTGLGFSKSEEDLNLYQIVVEGKLLIIVLYVDD
LILTGDEQLILSCKADLAKEFEMKDLGLLHYFLGLEIWQHSDGLFVSQGKYAREILEKFK
MHGCKPIDTPL
>super8277_107828_122351
MEQPLGFEVQDRDTYVCRLKKALYGLKQAPRAWYERIDSYLMKLGFTRSDADPNLYFKVD
NDKPLILVLYVDDLFLTCANPLIHQCKTELPSEFEIKDLGLMHYFLGLEVRQKPGEIFLS
QGKYVVKLLERFGMVDCKSVTTPM
>super8319_7037_15604
CMDYRQLNKLTVKKKYPLPRIDDLFNQVRGAKVFSKIDLRFRYHQIXIKDEDICKTTFHT
SIYGHYEFVVLPFGLTNAPATFMRLMHGIFHQFLDNFVLILYFILIFIDDILIYSKSEEE
HKXHLRIVLEMLRKHQLYAKFSKCDFFKTKIQYLGHV
>super8277_107828_122351
DDLFDQMKGATVFLKIDLRSGYHQFKISESDVSKTAFRTRFRHYEFIVVPFGLTNAPAVF
MNLMNNVLRKYLDKFVQVFLDDILIYSKNENEHEEHLRIVLSCLREHKLYGKLSKCAFFQ
如上序列，第一条和第三条，第四条和第六条的序列号时一致的，
问题：如何根据相同的序列号将序列合并？

zxy877298415

回复 1# 小小脑

合并后的效果是什么样子的？

小小脑

回复 2# zxy877298415

不好意思没说清楚。本来是想直接序列合并放一起就行了。但是发现这样不合理，所以我想是把重复给编个号。如下
[url=][/url]        [url=][/url] >super821_319450_333026
LMDELRGAQFLSKIDLRSGYNQIRVREQDIPKTAFRCHYGHFEFLVMPFGLTNAPATFQS
CMNHILRRHLRKFLLVFFDNILIYSRTWEEHFQHIEAVLHILEEQQFYAKLSKCEFGLTE
KLYLGHI
>super8297_1747273_1757914
RIDDLFDQMKGATVFSKIDLRFGYHQLRIKEDDIPKTAFQTXFGHYEFVVVPFGLTNAPI
VFTSLMNGVFKKYLDQFVQVFLDGILVYSKNKKEHEEHLRVVLTCLRKNQLYEKLSKRSF
FQNKIHYLGHII
>super821_319450_333026_1
NSAWEIIPRSEGKTVVGSRWIYKVKQAADESVEKYKARFVAXGFSXIEGIDYEETFAPVA
RYSSIRTILSLSAQMGWHXXXXXXXXXXXQMDVKTAFLNGVIEEEVYIEQPEGFEIFSSE
SHVCRLKIALYGLKQAPHAWYTRIDSYFTGLGFSKSEEDLNLYQIVVEGKLLIIVLYVDD
LILTGDEQLILSCKADLAKEFEMKDLGLLHYFLGLEIWQHSDGLFVSQGKYAREILEKFK
MHGCKPIDTPL
>super8277_107828_122351
MEQPLGFEVQDRDTYVCRLKKALYGLKQAPRAWYERIDSYLMKLGFTRSDADPNLYFKVD
NDKPLILVLYVDDLFLTCANPLIHQCKTELPSEFEIKDLGLMHYFLGLEVRQKPGEIFLS
QGKYVVKLLERFGMVDCKSVTTPM
>super8319_7037_15604
CMDYRQLNKLTVKKKYPLPRIDDLFNQVRGAKVFSKIDLRFRYHQIXIKDEDICKTTFHT
SIYGHYEFVVLPFGLTNAPATFMRLMHGIFHQFLDNFVLILYFILIFIDDILIYSKSEEE
HKXHLRIVLEMLRKHQLYAKFSKCDFFKTKIQYLGHV
>super8277_107828_122351_1
DDLFDQMKGATVFLKIDLRSGYHQFKISESDVSKTAFRTRFRHYEFIVVPFGLTNAPAVF
MNLMNNVLRKYLDKFVQVFLDDILIYSKNENEHEEHLRIVLSCLREHKLYGKLSKCAFFQ
>super821_319450_333026_2
NSAWEIIPRSEGKTVVGSRWIYKVKQAADESVEKYKARFVAXGFSXIEGIDYEETFAPVA
RYSSIRTILSLSAQMGWHXXXXXXXXXXXQMDVKTAFLNGVIEEEVYIEQPEGFEIFSSE
SHVCRLKIALYGLKQAPHAWYTRIDSYFTGLGFSKSEEDLNLYQIVVEGKLLIIVLYVDD
LILTGDEQLILSCKADLAKEFEMKDLGLLHYFLGLEIWQHSDGLFVSQGKYAREILEKFK
MHGCKPIDTPL
这样如何操作呢？请指教

zxy877298415

回复 3# 小小脑

awk -vRS=">" -F "\n" 'a[$1]++{$1=$1"_"a[$1]-1}1'  OFS="\n" ORS=">" file

小小脑

回复 4# zxy877298415

谢谢大神指点，但是文件的最后会有个'>'.这个是为什么呢？

zxy877298415

回复 5# 小小脑

awk -vRS=">" -F "\n" 'a[$1]++{$1=$1"_"a[$1]-1}NF{printf ">"$0}'  OFS="\n"  file

cd_ly

awk '/>/&&++a[$1]>1{$0=$0"_"++b[$0]}1' file

wh7211

回复 3# 小小脑

1. awk '/>/{a[$0]++;$0=a[$0]==1?$0:$0"_"a[$0]-1}1' file

复制代码