文本格式转换

mswsg

输入文件格式如下：

1. C4 Alignment:
   
2. ------------
   
3.          Query: 11117564
   
4.         Target: 3B|422103491-422104237|422103490-422104238
   
5.          Model: est2genome
   
6.      Raw score: 2560
   
7.    Query range: 0 -> 588
   
8.   Target range: 2 -> 748
   
9. 
   
10.    1 : CCGCCGCCGCCATGTATTCCTCTCCGACCTCGCCTAGAGCAGCCAGCAGGATGCCACCGCCGTG :  64
    
11.        |||||||||||||||| |||||||||||||||||| |||||||||||| |||||| ||||||||
    
12.    3 : CCGCCGCCGCCATGTAATCCTCTCCGACCTCGCCTGGAGCAGCCAGCAAGATGCCGCCGCCGTG :  66
    
13. 
    
14.   65 : GCCGCCTCTCTCCCTTCTCTTCTCTTTTTCCCTACCTTTCTTTCTCCTCTAAA-T------TGC : 121
    
15.        |||||||||||||||||||| |||||||||||||||||||| ||||||||||| |      ||
    
16.   67 : GCCGCCTCTCTCCCTTCTCTCCTCTTTTTCCCTACCTTTCTCTCTCCTCTAAAGTCTAAACTGT : 130
    
17. 
    
18. 122 : GTGCAGGTCCTCTCCGCAGGTGCTGCCATGGCC-CGTCATGGAGCTCGCCCCCGTGGGATGCAT : 184
    
19.        ||||||||||||||| |||||||||||||||||  |||||||||||||||||||||||||||||
    
20. 131 : GTGCAGGTCCTCTCCACAGGTGCTGCCATGGCCGTGTCATGGAGCTCGCCCCCGTGGGATGCAT : 194
    
21. 
    
22. 185 : GCCCCTCCAACGACGCCCTGTCCAG  >>>> Target Intron 1 >>>>  GCACTACGAG : 219
    
23.        |||||||||| ||||||||||||||++          67 bp          ++||||  ||||
    
24. 195 : GCCCCTCCAATGACGCCCTGTCCAGgt.........................agGCACCGCGAG : 296
    
25. 
    
26. 220 : CGCACTGGACGGCAGGAAGGGCTGCGTCAAGGAGCGACGCATTGCT  >>>> Target Intr : 266
    
27.        ||| ||||||||||  ||||| |||||||||||||||||||||| |++          83 bp
    
28. 297 : CGCGCTGGACGGCAAAAAGGGTTGCGTCAAGGAGCGACGCATTGTTgt................ : 345
    
29. 
    
30. 267 : on 2 >>>>  CGTGGACATCGCCTAGCGATCAGATTTGATGTTGAGAAAGACCCGGCAGGAAG : 318
    
31.                 ++||||||||||||||||| ||||||||||||||||||||||| |||||||||||
    
32. 346 : .........agCGTGGACATCGCCTAGCAATCAGATTTGATGTTGAGAAAGATCCGGCAGGAAG : 478
    
33. 
    
34. 319 : ATCTCATCCTGGAACCCAGGTGTGCTTGATTTGATATTGAGAATTCTACTTGATAGGTATGTAC : 382
    
35.        |||||||||||||||||||| ||||||||||||||||||||||||||||||| |||||||||||
    
36. 479 : ATCTCATCCTGGAACCCAGGAGTGCTTGATTTGATATTGAGAATTCTACTTGCTAGGTATGTAC : 542
    
37. 
    
38. 383 : ATGGATAGCGTACGGTGTATCCAAAGCTTTGCTTAACCTGAAATTTGTTAAAACGAATCTCTGG : 446
    
39.        || |||||||||||||||||||||||||||| |||||||||||||||||||| | |||||| |
    
40. 543 : ATAGATAGCGTACGGTGTATCCAAAGCTTTGTTTAACCTGAAATTTGTTAAAGCTAATCTCAGA : 606
    
41. 
    
42. 447 : TTTAGTAGAAGCCCACAGCGGTTGCAATACTAGAGTTATATCTATTGTGGGCAACATTGAAAAA : 510
    
43.        ||||||||||||||||||||||||||||| |||||||||||||| |||||||||||||||||||
    
44. 607 : TTTAGTAGAAGCCCACAGCGGTTGCAATATTAGAGTTATATCTAGTGTGGGCAACATTGAAAAA : 670
    
45. 
    
46. 511 : TTGCATATATCTAATCATCGTATATGTTATCACCTCTGAAACTGTGCAACATTACTACAAGATA : 574
    
47.        ||||||||||||||||||||||| ||||||||||||| ||||||||||||||||||||||||||
    
48. 671 : TTGCATATATCTAATCATCGTATGTGTTATCACCTCTAAAACTGTGCAACATTACTACAAGATA : 734
    
49. 
    
50. 575 : AATCCACGAAGCAG : 588
    
51.        ||||||||||||||
    
52. 735 : AATCCACGAAGCAG : 748
    
53. 
    
54. C4 Alignment:
    
55. ------------
    
56.          Query: 11440519
    
57.         Target: 3B|67788210-67788687|67788209-67788688 [revcomp]
    
58.          Model: est2genome
    
59.      Raw score: 2385
    
60.    Query range: 0 -> 477
    
61.   Target range: 479 -> 2
    
62. 
    
63.    1 : TAACCGCTCGAGATATTCGGTTACCGCGCGGTTACCACGTTTATCGCTCCCCCACGATAAACAC :  64
    
64.        ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
    
65. 479 : TAACCGCTCGAGATATTCGGTTACCGCGCGGTTACCACGTTTATCGCTCCCCCACGATAAACAC : 416
    
66. 
    
67.   65 : TCATACCGAGCAAAAAATCCCGATTTTTTTGAAATATTTGAATTCAAACGATCGCCTTGTAGTT : 128
    
68.        ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
    
69. 415 : TCATACCGAGCAAAAAATCCCGATTTTTTTGAAATATTTGAATTCAAACGATCGCCTTGTAGTT : 352
    
70. 
    
71. 129 : AAATAAGCTGCATCTGTGTTTTAAACAGAGAACTTATCGTGGCCCAGTGGCAAAAGGTGTTTTA : 192
    
72.        ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
    
73. 351 : AAATAAGCTGCATCTGTGTTTTAAACAGAGAACTTATCGTGGCCCAGTGGCAAAAGGTGTTTTA : 288
    
74. 
    
75. 193 : TCGCCAGACGTTGGTCGAGGGTTCGACTCGGGCTACGCACTGCTTTTTTAAACTTTTAGAGATG : 256
    
76.        ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
    
77. 287 : TCGCCAGACGTTGGTCGAGGGTTCGACTCGGGCTACGCACTGCTTTTTTAAACTTTTAGAGATG : 224
    
78. 
    
79. 257 : CAAAAGTGTTTACCTATGAAAAATTCTGGAGAAGACAGGGATCGAACTCGGGGCGCCTGAGCAA : 320
    
80.        ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
    
81. 223 : CAAAAGTGTTTACCTATGAAAAATTCTGGAGAAGACAGGGATCGAACTCGGGGCGCCTGAGCAA : 160
    
82. 
    
83. 321 : GGGATAGCTGCGGGTTGAGAGCCACTACGGTAGAACAACCTTACTGACTATTAATAGTTAAAAG : 384
    
84.        ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
    
85. 159 : GGGATAGCTGCGGGTTGAGAGCCACTACGGTAGAACAACCTTACTGACTATTAATAGTTAAAAG :  96
    
86. 
    
87. 385 : CTCATATATCTAACTCAAAAAAATTGAATTCAAAATTTGGTTTTAAATTTCGTCCGAATTTTTT : 448
    
88.        ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
    
89.   95 : CTCATATATCTAACTCAAAAAAATTGAATTCAAAATTTGGTTTTAAATTTCGTCCGAATTTTTT :  32
    
90. 
    
91. 449 : TCGGTAATTCGTGGTTACCGAGAATCTCG : 477
    
92.        | | |||||||||||||||||||||||||
    
93.   31 : TCGGTAATTCGTGGTTACCGAGAATCTCG :   3

复制代码

获取标红的结果，结果显示见最后：
C4 Alignment:
------------
         Query: 11117564
        Target: 3B|422103491-422104237|422103490-422104238
         Model: est2genome
     Raw score: 2560
   Query range: 0 -> 588
  Target range: 2 -> 748

   1 : CCGCCGCCGCCATGTATTCCTCTCCGACCTCGCCTAGAGCAGCCAGCAGGATGCCACCGCCGTG :  64
       |||||||||||||||| |||||||||||||||||| |||||||||||| |||||| ||||||||
   3 : CCGCCGCCGCCATGTAATCCTCTCCGACCTCGCCTGGAGCAGCCAGCAAGATGCCGCCGCCGTG :  66

  65 : GCCGCCTCTCTCCCTTCTCTTCTCTTTTTCCCTACCTTTCTTTCTCCTCTAAA-T------TGC : 121
       |||||||||||||||||||| |||||||||||||||||||| ||||||||||| |      ||
  67 : GCCGCCTCTCTCCCTTCTCTCCTCTTTTTCCCTACCTTTCTCTCTCCTCTAAAGTCTAAACTGT : 130

122 : GTGCAGGTCCTCTCCGCAGGTGCTGCCATGGCC-CGTCATGGAGCTCGCCCCCGTGGGATGCAT : 184
       ||||||||||||||| |||||||||||||||||  |||||||||||||||||||||||||||||
131 : GTGCAGGTCCTCTCCACAGGTGCTGCCATGGCCGTGTCATGGAGCTCGCCCCCGTGGGATGCAT : 194

185 : GCCCCTCCAACGACGCCCTGTCCAG  >>>> Target Intron 1 >>>>  GCACTACGAG : 219
       |||||||||| ||||||||||||||++          67 bp          ++||||  ||||
195 : GCCCCTCCAATGACGCCCTGTCCAGgt.........................agGCACCGCGAG : 296

220 : CGCACTGGACGGCAGGAAGGGCTGCGTCAAGGAGCGACGCATTGCT  >>>> Target Intr : 266
       ||| ||||||||||  ||||| |||||||||||||||||||||| |++          83 bp
297 : CGCGCTGGACGGCAAAAAGGGTTGCGTCAAGGAGCGACGCATTGTTgt................ : 345

267 : on 2 >>>>  CGTGGACATCGCCTAGCGATCAGATTTGATGTTGAGAAAGACCCGGCAGGAAG : 318
                ++||||||||||||||||| ||||||||||||||||||||||| |||||||||||
346 : .........agCGTGGACATCGCCTAGCAATCAGATTTGATGTTGAGAAAGATCCGGCAGGAAG : 478

319 : ATCTCATCCTGGAACCCAGGTGTGCTTGATTTGATATTGAGAATTCTACTTGATAGGTATGTAC : 382
       |||||||||||||||||||| ||||||||||||||||||||||||||||||| |||||||||||
479 : ATCTCATCCTGGAACCCAGGAGTGCTTGATTTGATATTGAGAATTCTACTTGCTAGGTATGTAC : 542

383 : ATGGATAGCGTACGGTGTATCCAAAGCTTTGCTTAACCTGAAATTTGTTAAAACGAATCTCTGG : 446
       || |||||||||||||||||||||||||||| |||||||||||||||||||| | |||||| |
543 : ATAGATAGCGTACGGTGTATCCAAAGCTTTGTTTAACCTGAAATTTGTTAAAGCTAATCTCAGA : 606

447 : TTTAGTAGAAGCCCACAGCGGTTGCAATACTAGAGTTATATCTATTGTGGGCAACATTGAAAAA : 510
       ||||||||||||||||||||||||||||| |||||||||||||| |||||||||||||||||||
607 : TTTAGTAGAAGCCCACAGCGGTTGCAATATTAGAGTTATATCTAGTGTGGGCAACATTGAAAAA : 670

511 : TTGCATATATCTAATCATCGTATATGTTATCACCTCTGAAACTGTGCAACATTACTACAAGATA : 574
       ||||||||||||||||||||||| ||||||||||||| ||||||||||||||||||||||||||
671 : TTGCATATATCTAATCATCGTATGTGTTATCACCTCTAAAACTGTGCAACATTACTACAAGATA : 734

575 : AATCCACGAAGCAG : 588
       ||||||||||||||
735 : AATCCACGAAGCAG : 748

C4 Alignment:
------------
         Query: 11440519
        Target: 3B|67788210-67788687|67788209-67788688 [revcomp]
         Model: est2genome
     Raw score: 2385
   Query range: 0 -> 477
  Target range: 479 -> 2

   1 : TAACCGCTCGAGATATTCGGTTACCGCGCGGTTACCACGTTTATCGCTCCCCCACGATAAACAC :  64
       ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
479 : TAACCGCTCGAGATATTCGGTTACCGCGCGGTTACCACGTTTATCGCTCCCCCACGATAAACAC : 416

  65 : TCATACCGAGCAAAAAATCCCGATTTTTTTGAAATATTTGAATTCAAACGATCGCCTTGTAGTT : 128
       ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
415 : TCATACCGAGCAAAAAATCCCGATTTTTTTGAAATATTTGAATTCAAACGATCGCCTTGTAGTT : 352

129 : AAATAAGCTGCATCTGTGTTTTAAACAGAGAACTTATCGTGGCCCAGTGGCAAAAGGTGTTTTA : 192
       ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
351 : AAATAAGCTGCATCTGTGTTTTAAACAGAGAACTTATCGTGGCCCAGTGGCAAAAGGTGTTTTA : 288

193 : TCGCCAGACGTTGGTCGAGGGTTCGACTCGGGCTACGCACTGCTTTTTTAAACTTTTAGAGATG : 256
       ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
287 : TCGCCAGACGTTGGTCGAGGGTTCGACTCGGGCTACGCACTGCTTTTTTAAACTTTTAGAGATG : 224

257 : CAAAAGTGTTTACCTATGAAAAATTCTGGAGAAGACAGGGATCGAACTCGGGGCGCCTGAGCAA : 320
       ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
223 : CAAAAGTGTTTACCTATGAAAAATTCTGGAGAAGACAGGGATCGAACTCGGGGCGCCTGAGCAA : 160

321 : GGGATAGCTGCGGGTTGAGAGCCACTACGGTAGAACAACCTTACTGACTATTAATAGTTAAAAG : 384
       ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
159 : GGGATAGCTGCGGGTTGAGAGCCACTACGGTAGAACAACCTTACTGACTATTAATAGTTAAAAG :  96

385 : CTCATATATCTAACTCAAAAAAATTGAATTCAAAATTTGGTTTTAAATTTCGTCCGAATTTTTT : 448
       ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
  95 : CTCATATATCTAACTCAAAAAAATTGAATTCAAAATTTGGTTTTAAATTTCGTCCGAATTTTTT :  32

449 : TCGGTAATTCGTGGTTACCGAGAATCTCG : 477
       |||||||||||||||||||||||||||||
  31 : TCGGTAATTCGTGGTTACCGAGAATCTCG :   3

python，结果展示：
>11117564 3B|422103491-422104237|422103490-422104238
CCGCCGCCGCCATGTAATCCTCTCCGACCTCGCCTGGAGCAGCCAGCAAGATGCCGCCGCCGTG
GCCGCCTCTCTCCCTTCTCTCCTCTTTTTCCCTACCTTTCTCTCTCCTCTAAAGTCTAAACTGT
GTGCAGGTCCTCTCCACAGGTGCTGCCATGGCCGTGTCATGGAGCTCGCCCCCGTGGGATGCAT
GCCCCTCCAATGACGCCCTGTCCAGGCACCGCGAG
CGCGCTGGACGGCAAAAAGGGTTGCGTCAAGGAGCGACGCATTGTT
CGTGGACATCGCCTAGCAATCAGATTTGATGTTGAGAAAGATCCGGCAGGAAG
ATCTCATCCTGGAACCCAGGAGTGCTTGATTTGATATTGAGAATTCTACTTGCTAGGTATGTAC
ATAGATAGCGTACGGTGTATCCAAAGCTTTGTTTAACCTGAAATTTGTTAAAGCTAATCTCAGA
TTTAGTAGAAGCCCACAGCGGTTGCAATATTAGAGTTATATCTAGTGTGGGCAACATTGAAAAA
TTGCATATATCTAATCATCGTATGTGTTATCACCTCTAAAACTGTGCAACATTACTACAAGATA
AATCCACGAAGCAG
>11440519 3B|67788210-67788687|67788209-67788688 [revcomp]
TAACCGCTCGAGATATTCGGTTACCGCGCGGTTACCACGTTTATCGCTCCCCCACGATAAACAC
TCATACCGAGCAAAAAATCCCGATTTTTTTGAAATATTTGAATTCAAACGATCGCCTTGTAGTT
AAATAAGCTGCATCTGTGTTTTAAACAGAGAACTTATCGTGGCCCAGTGGCAAAAGGTGTTTTA
TCGCCAGACGTTGGTCGAGGGTTCGACTCGGGCTACGCACTGCTTTTTTAAACTTTTAGAGATG
CAAAAGTGTTTACCTATGAAAAATTCTGGAGAAGACAGGGATCGAACTCGGGGCGCCTGAGCAA
GGGATAGCTGCGGGTTGAGAGCCACTACGGTAGAACAACCTTACTGACTATTAATAGTTAAAAG
CTCATATATCTAACTCAAAAAAATTGAATTCAAAATTTGGTTTTAAATTTCGTCCGAATTTTTT
TCGGTAATTCGTGGTTACCGAGAATCTCG

sunzhiguolu

回复 1# mswsg
perl 的行吗?

   

mswsg

你有时间的话也可以， 回复 2# sunzhiguolu


   

sunzhiguolu

回复 3# mswsg
这画面也太暴力了!!!

   

sunzhiguolu

回复 3# mswsg
试下,

1. #!/usr/bin/perl -w
   
2. use strict;
   
3. 
   
4. $/ = '';
   
5. my $head = '>';
   
6. while (<>){
   
7.     my @aData;
   
8.     if (@aData = /(?|Query:\s(\S+)|Target:(\s\S+))/g){
   
9.         print $head, @aData, "\n";
   
10.         next;
    
11.     }
    
12.     print @aData, "\n" if ((@aData = split /\n/) and (@aData = $aData[-1] =~ /[A-Z]+/g));
    
13. }
    

复制代码

perl abc.pl a
----------------------------------------------------------------------------
>11117564 3B|422103491-422104237|422103490-422104238
CCGCCGCCGCCATGTAATCCTCTCCGACCTCGCCTGGAGCAGCCAGCAAGATGCCGCCGCCGTG
GCCGCCTCTCTCCCTTCTCTCCTCTTTTTCCCTACCTTTCTCTCTCCTCTAAAGTCTAAACTGT
GTGCAGGTCCTCTCCACAGGTGCTGCCATGGCCGTGTCATGGAGCTCGCCCCCGTGGGATGCAT
GCCCCTCCAATGACGCCCTGTCCAGGCACCGCGAG
CGCGCTGGACGGCAAAAAGGGTTGCGTCAAGGAGCGACGCATTGTT
CGTGGACATCGCCTAGCAATCAGATTTGATGTTGAGAAAGATCCGGCAGGAAG
ATCTCATCCTGGAACCCAGGAGTGCTTGATTTGATATTGAGAATTCTACTTGCTAGGTATGTAC
ATAGATAGCGTACGGTGTATCCAAAGCTTTGTTTAACCTGAAATTTGTTAAAGCTAATCTCAGA
TTTAGTAGAAGCCCACAGCGGTTGCAATATTAGAGTTATATCTAGTGTGGGCAACATTGAAAAA
TTGCATATATCTAATCATCGTATGTGTTATCACCTCTAAAACTGTGCAACATTACTACAAGATA
AATCCACGAAGCAG
>11440519 3B|67788210-67788687|67788209-67788688
TAACCGCTCGAGATATTCGGTTACCGCGCGGTTACCACGTTTATCGCTCCCCCACGATAAACAC
TCATACCGAGCAAAAAATCCCGATTTTTTTGAAATATTTGAATTCAAACGATCGCCTTGTAGTT
AAATAAGCTGCATCTGTGTTTTAAACAGAGAACTTATCGTGGCCCAGTGGCAAAAGGTGTTTTA
TCGCCAGACGTTGGTCGAGGGTTCGACTCGGGCTACGCACTGCTTTTTTAAACTTTTAGAGATG
CAAAAGTGTTTACCTATGAAAAATTCTGGAGAAGACAGGGATCGAACTCGGGGCGCCTGAGCAA
GGGATAGCTGCGGGTTGAGAGCCACTACGGTAGAACAACCTTACTGACTATTAATAGTTAAAAG
CTCATATATCTAACTCAAAAAAATTGAATTCAAAATTTGGTTTTAAATTTCGTCCGAATTTTTT
TCGGTAATTCGTGGTTACCGAGAATCTCG


   

sunzhiguolu

大神, 我来领分了...

haooooaaa

1. import re
   
2. with open("a","rb") as f:
   
3.     content = f.read()
   
4.     for i in re.split('C4 Alignment', content)[1:]:
   
5.         print ">" + "".join(re.findall(r'Query: (.*)\n *Target: (.*)',i)[0])
   
6.         print re.sub("[^A-Z\n]","","\n".join(re.findall(r'\|+[^\n]+\n *\d+ : ([^:]+)',i)))

复制代码

1. awk '/Query:|Target:/{sub(" +"$1,"");printf m++%2?$0"\n":$1}/ +[0-9]+ :/&&i++%2{gsub("[^A-Z]","");print}'

复制代码

Linux_manne

代码比较丑.. 匆匆写了下..

1. import re
   
2. 
   
3. p1 = "Query:(.*)$"
   
4. p2 = "Target:(.*)$"
   
5. p3 = "(\|+)$"
   
6. p4 = "([A-Z]+)"
   
7. 
   
8. f = open('./sample.txt')
   
9. result = []
   
10. fl = f.readlines()
    
11. ins = []
    
12. xl = enumerate(map (lambda x: x.strip(), fl ))
    
13. for index, line in xl:
    
14.         match1 = re.search(p1,line)
    
15.         match2 = re.search(p2,line)
    
16.         match3 = re.search(p3,line)
    
17. 
    
18.         if match1:
    
19.                 if len(ins) > 1:
    
20.                         result.append(ins)
    
21.                 ins = []
    
22.                 ins.append(match1.groups()[0])
    
23.                
    
24.         if match2:
    
25.                 ins.append(match2.groups()[0])
    
26.         if match3 and (index + 1) < len(fl):
    
27.                 match4  = re.search(p4, fl[index + 1])
    
28.                 ins.append(match4.groups()[0])
    
29. result.append(ins)       
    
30. for i in result:
    
31.         print ''.join(i)

复制代码

sunzhiguolu

这回对了,

1. #!/usr/bin/perl -w
   
2. 
   
3. $/ = '';
   
4. while (<>){
   
5.     my @aData;
   
6.     if (@aData = /(?|Query:\s(\S+)|Target:(\s\S+(?:\s\[revcomp])?))/g){
   
7.         print '>', @aData, "\n";
   
8.         next;
   
9.     }
   
10.     print @aData, "\n" if ((@aData = split /\n/) and (@aData = $aData[-1] =~ /[A-Z]+/g));
    
11. }
    

复制代码

>11117564 3B|422103491-422104237|422103490-422104238
CCGCCGCCGCCATGTAATCCTCTCCGACCTCGCCTGGAGCAGCCAGCAAGATGCCGCCGCCGTG
GCCGCCTCTCTCCCTTCTCTCCTCTTTTTCCCTACCTTTCTCTCTCCTCTAAAGTCTAAACTGT
GTGCAGGTCCTCTCCACAGGTGCTGCCATGGCCGTGTCATGGAGCTCGCCCCCGTGGGATGCAT
GCCCCTCCAATGACGCCCTGTCCAGGCACCGCGAG
CGCGCTGGACGGCAAAAAGGGTTGCGTCAAGGAGCGACGCATTGTT
CGTGGACATCGCCTAGCAATCAGATTTGATGTTGAGAAAGATCCGGCAGGAAG
ATCTCATCCTGGAACCCAGGAGTGCTTGATTTGATATTGAGAATTCTACTTGCTAGGTATGTAC
ATAGATAGCGTACGGTGTATCCAAAGCTTTGTTTAACCTGAAATTTGTTAAAGCTAATCTCAGA
TTTAGTAGAAGCCCACAGCGGTTGCAATATTAGAGTTATATCTAGTGTGGGCAACATTGAAAAA
TTGCATATATCTAATCATCGTATGTGTTATCACCTCTAAAACTGTGCAACATTACTACAAGATA
AATCCACGAAGCAG
>11440519 3B|67788210-67788687|67788209-67788688 [revcomp]
TAACCGCTCGAGATATTCGGTTACCGCGCGGTTACCACGTTTATCGCTCCCCCACGATAAACAC
TCATACCGAGCAAAAAATCCCGATTTTTTTGAAATATTTGAATTCAAACGATCGCCTTGTAGTT
AAATAAGCTGCATCTGTGTTTTAAACAGAGAACTTATCGTGGCCCAGTGGCAAAAGGTGTTTTA
TCGCCAGACGTTGGTCGAGGGTTCGACTCGGGCTACGCACTGCTTTTTTAAACTTTTAGAGATG
CAAAAGTGTTTACCTATGAAAAATTCTGGAGAAGACAGGGATCGAACTCGGGGCGCCTGAGCAA
GGGATAGCTGCGGGTTGAGAGCCACTACGGTAGAACAACCTTACTGACTATTAATAGTTAAAAG
CTCATATATCTAACTCAAAAAAATTGAATTCAAAATTTGGTTTTAAATTTCGTCCGAATTTTTT
TCGGTAATTCGTGGTTACCGAGAATCTCG

mswsg

谢谢您，正在学习，和我想的思路有点不一样，正在学习。今天没分了，明天加回复 7# haooooaaa