如何调整fastaID的位置

321wangke321

有这样一个fasta文件

>MYLAATYLSNGSELLLEILGPGIIGGLFLPILGALPDAMLILVSGLSGTTETAQSQVSVGMGLLAGSTVMLLTVIWGSCILVGKCDIENSIAIDNKDTKGLSLTGSGVSTDIWTSYAARIMAISVIPFLIVQLPQLLNSTSGRHLAVLISLIVSVVLLVSYCLYQVFQPWIQRRKIAYAKHKHVISGLLQHLTNQERLLTVDGEPDEEIIKKLFKTIDQDSDGYLSFSELRALIVGIRFDEIELDKNDAVNKVMKDFDTSHDSRISVDEFWVGVSKWLKEAKRAGNANSDPGHRTMKFLSDFHSRTKEEHALLGAGDQSDEVVEGAENPKWTSLKAGLMLLVGTLIAAAFADPLVDAVDNFSDATSIPTFFISFIALPLATNSSEAVSAIIFASRKKIRTASLTFSELYGAVTMNNILCLSVFLALVYIRGLSWDFSAEVLVILIVCVVMGILGSLRTNYPLWTCSIAFLLYPFSLALVYVLDYVLGWS*
ppa004833m.g

>MGTSTTISVLFLMLIALHSPIHANPIKTIVVIVMENRSFDHMLGWMKRINPEINGVDGSEANPLSTTDPNSKRFFFKNESHYVDPDPGHSFQAIREQIFGSDNTSADPPPMNGFAQQAFSMDNTTNMSQDVMNGFNPDMVAVYKTLVSEFAVCDRWFASVPSSTQPNRLYVHSGTSAGATSNIPALLAKGYPQRTIFENVDDAGISFGIYYQNIPATLFYRNLRKLKYISKFHMYDTFKSHAKQGKLPGYTVVEQRYFDVKIAPANDDHPSHDVFQGQMFVKEVYETLRASPQWNDTLLIITYDEHGGFYDHVATPVHGVPSPDGIVGPEPFLFQFDRLGVRVPTIMVSPWIEKGTVLHGPNGSPFPTSEFEHSSIPATVRKIFNLSAPLTKRDEWAGTFEGIVQTRTEPRTDCPEQLPTPVKIRTSEPKEDAKLTEFQQEMTQLAAVLKGDNIFTSYPKKTGKEMTVKDGKEYADDAVKRFFEAGLYAKRMGVSEEQIVQMRPSLTTRPSPIVNQTP*
ppa004282m.g

如何变成
>ppa004833m.g
MYLAATYLSNGSELLLEILGPGIIGGLFLPILGALPDAMLILVSGLSGTTETAQSQVSVGMGLLAGSTVMLLTVIWGSCILVGKCDIENSIAIDNKDTKGLSLTGSGVSTDIWTSYAARIMAISVIPFLIVQLPQLLNSTSGRHLAVLISLIVSVVLLVSYCLYQVFQPWIQRRKIAYAKHKHVISGLLQHLTNQERLLTVDGEPDEEIIKKLFKTIDQDSDGYLSFSELRALIVGIRFDEIELDKNDAVNKVMKDFDTSHDSRISVDEFWVGVSKWLKEAKRAGNANSDPGHRTMKFLSDFHSRTKEEHALLGAGDQSDEVVEGAENPKWTSLKAGLMLLVGTLIAAAFADPLVDAVDNFSDATSIPTFFISFIALPLATNSSEAVSAIIFASRKKIRTASLTFSELYGAVTMNNILCLSVFLALVYIRGLSWDFSAEVLVILIVCVVMGILGSLRTNYPLWTCSIAFLLYPFSLALVYVLDYVLGWS*


>ppa004282m.g
MGTSTTISVLFLMLIALHSPIHANPIKTIVVIVMENRSFDHMLGWMKRINPEINGVDGSEANPLSTTDPNSKRFFFKNESHYVDPDPGHSFQAIREQIFGSDNTSADPPPMNGFAQQAFSMDNTTNMSQDVMNGFNPDMVAVYKTLVSEFAVCDRWFASVPSSTQPNRLYVHSGTSAGATSNIPALLAKGYPQRTIFENVDDAGISFGIYYQNIPATLFYRNLRKLKYISKFHMYDTFKSHAKQGKLPGYTVVEQRYFDVKIAPANDDHPSHDVFQGQMFVKEVYETLRASPQWNDTLLIITYDEHGGFYDHVATPVHGVPSPDGIVGPEPFLFQFDRLGVRVPTIMVSPWIEKGTVLHGPNGSPFPTSEFEHSSIPATVRKIFNLSAPLTKRDEWAGTFEGIVQTRTEPRTDCPEQLPTPVKIRTSEPKEDAKLTEFQQEMTQLAAVLKGDNIFTSYPKKTGKEMTVKDGKEYADDAVKRFFEAGLYAKRMGVSEEQIVQMRPSLTTRPSPIVNQTP*

劳烦大神出手，万分感谢。

haooooaaa

1. sed -n '/>/{s/>//;x;n;s/^/>/;G;p}' b

复制代码

1. sed '/^>/{N;s/>\(.*\)\n\(.*\)/>\2\n\1/}' b
   
2. 
   
3. awk 'a{print ">"$0"\n"a;a=""}/>/{a=substr($0,2)}' b

复制代码

1. t = ""
   
2. for i in open("b"):
   
3.     if i[0] == ">":
   
4.         t = i[1:].strip()
   
5.     elif i[0] != ">" and i.strip():
   
6.         print ">%s\n%s\n" % (i.strip(), t)

复制代码

321wangke321

回复 2# haooooaaa


    我只能表示无尽的膜拜。。。。。。。。。。。一用就用三种方法。谢谢大神。