格式變換

linaiting

我這有兩種字符序列，想轉(zhuǎn)變成標(biāo)準(zhǔn)的格式，如下：
   
字符序列1：
>String1
atggaggtcaatggtcatcacaaaccaagaagagaatataatggccgagagtgcaatggtgtacaatcagtaaacaattttggcgatatcgatccatggacagcatgggcatacaaaccgcgcacagtttcattgctactgatgggaacatgctttttaatttgggcaagtggtgctcttgatccagaaagaagcttctctgttgaccgcgtttcatctgttaaaaggggtgtctttgcaatgattgctgtttttttggcttattcatttcttcaggcaccttctactgtacttattagaccacatcctgccatttggcggctggtccctaacttttttgcttttccagacccgtgatgatgctaggcaatttatgaagtatcttcaccctgatcttggtgttgaattacctgaaagatcttatggaaccgactgccgcatatatgtacctgatcatccgaaaagcaggtttaacaatgtttatgagatcctttttgatgagtttgttattgcccatatcctgctataatgataagaaacacttttatgggtattatcaattggctttgagctaatggagctcacctttcgccatatgctgccaaactttaatgagtgttggtgggatagcatcgtactagacatattaatctgcaattggtttggtatttgggctggaatgaagactgtgagatactttgatgggaggacatatgaatgggttggtttgagtcgtcaacccaatattatcagtaaggtaaaaaggacgctaggccagttcacaccagcacagtgggacaaagatgagtggtaccctctgcttggcccttggagattcatccaggtgctgagcctatgtattgtttccatgattgttgaacttaaca

想轉(zhuǎn)變成的目標(biāo)格式是這樣的：



其實(shí)就是每10個(gè)字符空一格，一行60個(gè)字符，然后在每行末尾標(biāo)明總的字符數(shù)

字符序列2：
>String2
Met Glu Val Asn Gly His His Lys Pro Arg Arg Glu Tyr Asn Gly Arg Glu Cys Asn Gly Val Gln Ser Val Asn Asn Phe Gly Asp Ile Asp Pro Trp Thr Ala Trp Ala Tyr Lys Pro Arg Thr Val Ser Leu Leu Leu Met Gly Thr Cys Phe Leu Ile Trp Ala Ser Gly Ala Leu Asp Pro Glu Arg Ser Phe Ser Val Asp Arg Val Ser Ser Val Lys Arg Gly Val Phe Ala Met Ile Ala Val Phe Leu Ala Tyr Ser Phe Leu Gln Ala Pro Ser Thr Val Leu Ile Arg Pro His Pro Ala Ile Trp Arg Leu Val His Gly Met Ala Val Val Tyr Leu Val Ala Leu Thr Phe Leu Leu Phe Gln Thr Arg Asp Asp Ala Arg Gln Phe Met Lys Tyr Leu His Pro Asp Leu Gly Val Glu Leu Pro Glu Arg Ser Tyr Gly Thr Asp Cys Arg Ile Tyr Val Pro Asp His Pro Lys Ser Arg Phe Asn Asn Val Tyr Glu Ile Leu Phe Gly Trp

想轉(zhuǎn)變成以下格式：



每三個(gè)字符看成一個(gè)單位，每隔5個(gè)單位標(biāo)上號(hào)（另起一行標(biāo)號(hào)），其實(shí)原文件每個(gè)單位之間是沒有空格的，只是為了方面看，我空了隔

這看似不難，卻把我弄得焦頭爛額，用shell和perl都有嘗試過，卻弄不出頭緒，達(dá)不到要求，希望大蝦們幫助一下小弟，非常感謝啊

liion631818

字符序列1

1. 
   
2. awk -F "" '{for(i=1;i<=NF;i++){printf($i);if(i%10==0)printf(" ");if(i%60==0)printi}}END{if(i%60)space=sprintf("%"66-(i-1)%60-((i-1)%60)/10"s", " "); print space, i-1}'
   

復(fù)制代碼

waker

dd if=urfile cbs=10 conv=unblock|xargs -n6|awk '{print $0,a+=length($0)-NF+1}'

yinyuemi

回復(fù) 1# linaiting

1. awk --re-interval '{l=length;printf substr($0,1,10)" ";a=1;while(i=sub(/^.{10}/,"",$0)){a+=i;if(a%6){printf substr($0,1,10)" "}else{printf substr($0,1,10)" "(sprintf("%8s",60*++j))"\n"}};print l}' dna
   

復(fù)制代碼

expert1

第一個(gè)簡單，設(shè)FS “”
逐個(gè)打印，每10個(gè)字符打印一個(gè)“ ”，計(jì)數(shù)60個(gè)，打印"\n",然后打印60N（N=1,2,3，，，）,
這個(gè)就是繁瑣，貌似2樓不錯(cuò)。
string 2
有點(diǎn)麻煩。

linaiting

字符序列1
liion631818 發(fā)表于 2011-09-21 14:24

    結(jié)果達(dá)到我想要的啊，很感謝大俠的幫忙��！

linaiting

回復(fù)  linaiting
yinyuemi 發(fā)表于 2011-09-21 14:37

     "dna"，呵呵，yinyuemi 可能也是學(xué)生物的，代碼好像不能計(jì)算未達(dá)到60n的字符數(shù)（最后一行無法顯示總數(shù)），不過改進(jìn)之后應(yīng)該是可以的。

   不知道各位大蝦對(duì)第二個(gè)字符序列的處理有沒有什么idea�。�？

waker

"dna"，呵呵，yinyuemi 可能也是學(xué)生物的，代碼好像不能計(jì)算未達(dá)到60n的字符數(shù)（最后一行無法顯 ...
linaiting 發(fā)表于 2011-09-21 15:06

    把第一個(gè)的方法套到第二個(gè)就行了

linaiting

dd if=urfile cbs=10 conv=unblock|xargs -n6|awk '{print $0,a+=length($0)-NF+1}'
waker 發(fā)表于 2011-09-21 14:27

    運(yùn)行之后，好像只能顯示總字符數(shù)，不過學(xué)習(xí)了，很感謝幫忙啊，不知大俠對(duì)String2有沒有什么想法啊

linaiting

把第一個(gè)的方法套到第二個(gè)就行了
waker 發(fā)表于 2011-09-21 15:10

    我試試