如何調(diào)整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


    我只能表示無盡的膜拜。。。。。。。。。。。一用就用三種方法。謝謝大神。