File: split1.maf

package info (click to toggle)
last-align 1645-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 14,684 kB
  • sloc: cpp: 44,403; python: 5,212; ansic: 1,938; sh: 710; makefile: 457
file content (201 lines) | stat: -rw-r--r-- 16,701 bytes parent folder | download | duplicates (3)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
#
# a=31 b=1 A=21 B=3 e=190 d=94 x=189 y=48 z=189 D=1e+09 E=0.516058
# R=01 u=0 s=2 S=1 M=0 T=0 m=100 l=1 n=100 k=1 w=1000 t=4.75495 j=3 Q=0
# /home/mcfrith/data/genome/calMil1/last/calMil1-uMAM8
# Reference sequences=21204 normal letters=936953443
# lambda=0.214243 K=0.246122
#
#     A   C   G   T   M   S   K   W   R   Y   B   D   H   V
# A   5  -7  -3  -6   2  -5  -4   2   3  -6  -5   1   1   1
# C  -8   5  -7  -2   2   2  -4  -4  -7   3   1  -5   1   0
# G  -2  -7   5  -8  -4   2   2  -4   3  -7   0   1  -5   1
# T  -6  -3  -7   5  -4  -5   2   2  -6   3   1   1   1  -5
# M   2   2  -5  -4   2   0  -4   0   0   0  -1  -1   1   1
# S  -4   2   2  -4   0   2   0  -4   0   0   1  -1  -1   1
# K  -4  -5   2   2  -4   0   2   0   0   0   1   1  -1  -1
# W   2  -5  -5   2   0  -5   0   2   0   0  -1   1   1  -1
# R   3  -7   3  -7   0   0   0   0   3  -7  -2   1  -1   1
# Y  -7   3  -7   3   0   0   0   0  -7   3   1  -1   1  -2
# B  -5   1   1   1  -1   1   1  -1  -1   1   1   0   0   0
# D   1  -5   1   0  -1  -1   1   1   1  -2  -1   1  -1   0
# H   0   1  -5   1   1  -1  -1   1  -2   1   0  -1   1  -1
# V   1   1   1  -5   1   1  -1  -1   1  -1   0   0   0   1
#
# Coordinates are 0-based.  For - strand matches, coordinates
# in the reverse complement of the 2nd sequence are used.
#
# name start alnSize strand seqSize alignment
#
# m=1 s=190
#
# Query sequences=195
a score=211 mismap=0.00255
s KI635889   4462115 58 +   6429418 TACCTATAGATTGGTTGACAGCTTTCATGCACAAATTTTCCACCAGAATCTTTCCTGG
s chr1     248943367 58 - 248956422 TGCCTATGGACTGGTTGACAGCCTTCATGCACAGGTTCTCCACCAGAGCCTTTCTTGG
p                                   #$)-1356899::::::::::::::::::::::::::::::::::::9888640,+(%

a score=331 mismap=1e-10
s KI635917 1140607 159 +   4339253 ACCTTGTTcttgttccttctgcttcttcttcttctccatcttaatcAGCTTTGCATTGCGCAGTTTGGCCTTTCCAATGCCTCCGGCCTGACGAATGGATTCCAGCAGACTAGCCCTTCCATCCGAGGGTTTCACAACCTCCTTGGGTGCCCCCTGGAC
s chr1       15793 155 + 248956422 ACCTTGCTCctgctccttctgctgctg---cttctccaGCTTTCGCTCCTTCATGCTGCGCAGCTTGGCCTTGCCGATGCCCCCAGCTTGGCGGATGGACTCTAGCAGAGTGG-CCAGCCACCGGAGGGGTCAACCACTTCCCTGGGAGCTCCCTGGAC
p                                  %*.2578<>?BDFJMOPPQQQQQQQQQQQQUY\^___________```bbccccdejnsw{~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~{zwtpkgbb``^^^^^^^^^^^^^^^^]]\\\\[[[[ZZXUTROKGFDBA>;73.*%

a score=520 mismap=1e-10
s KI635917 1143406 188 +   4339253 TCTGTGCTGAACGTGGGCAGTTCAGGCATGTTGCTTCCTGGCACGGAGGGCGCGATGCCAGGTCCCAGGTCAGCACTGTACATGAGGTCATCTGCAATCCCAGGCAGATCCGGTAGGTAGGATGGGACGTCAATCTCTGGGACCTGACCCAGGTCTGGTACATAGAAGTAATTCTCAGTAGTCTGAAA
s chr1       16876 185 + 248956422 TCAGTGTGGAAGGTGGGCAGTTCTGGAATGGTGC---CAGGGGCAGAGGGGGCAATGCCGGGGCCCAGGTCGGCAATGTACATGAGGTCGTTGGCAATGCCGGGCAGGTCAGGCAGGTAGGATGGAACATCAATCTCAGGCACCTGGCCCAGGTCTGGCACATAGAAGTAGTTCTCTGGGACCTGCAA
p                                  #$%()**+.1237<@EILOPQQQQQQQQQQQQQQQQQRRSTTUXY]begiimoquz}~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~}xsojigc_[VRMID?;:740+'&$$$$$$#"""

a score=340 mismap=1e-10
s KI635917 1144102 140 +   4339253 ACCCGTCTTTCTAGTTGCTCTCGCGTGGTGATTGAGAGCGGGGCATCAAACGGCTTCTCTTCCTTCTCAGTCTCCAGTGCCATGTGTGTCTTTGTCACAGCTCCTGCTAACGGATCCAAGAGAACGTACTTCTTGTAGCT
s chr1       17230 140 + 248956422 ACCTGCTGTTCCAGCTGCTCTCTCTTGCTGATGGACAAGGGGGCATCAAACAGCTTCTCCTCTGTCTCTGCCCCCAGCATCACATGGGTCTTTGTTACAGCACCAGCCAGGGGGTCCAGGAAGACATACTTCTTCTACCT
p                                  #%%&&&''+.12589=BEHJKLLLMOPQTWXYYZZ[\\\afjosx|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~|xsonljihfecbbba`_]ZVURONMLIHFC?;62-((&$$#

a score=334 mismap=1e-10
s KI635917 1144725 145 +   4339253 TACGGATTCTCAGTAGTGTTGAACAGCAGCAGGGAGCTGATGGAGCTAATGTTCCGAGGAAGGCTTCCCAACCCCTCCTCCGTCTCATCCTCCTCTCGAGACTTGGTGTTCAGACACACGGGGTAATGCTTTGCCTTTTCCTGAG
s chr1       17602 145 + 248956422 TACAGGTTCTCGGTGGTGTTGAAGAGCAGCAAGGAGCTGACAGAGCTGATGTTGCTGGGAAGACCCCCAAGTCCCTCTTCTGCATCGTCCTCGGGCTCCGGCTTGGTGCTCACGCACACAGGAAAGTCCTTCAGCTTCTCCTGAG
p                                  $(*+-/37;?ABEHINRW[_cegglpuy}~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~}yxvuttssrplhdcccccccccccccccba`^[ZZZXVRNLJFEEDCBBBA>>===<::862.*%

a score=340 mismap=1e-10
s KI635917 1146887 159 +   4339253 CACACCTGCACAGCCTTATCATCCAGAGGCTTCAGCTTGCTGTGGATCTTGTAGCGTGGGCGCTTCTGAGAAGCAGGATCATCCGCGCCCGAGAATATGGAAGTGTAATCCTGTAGCCGGTCAGGTGCAGGGTACTTAGCACTGGAAAAGACCTGCAGA
s chr1       17909 159 + 248956422 CAGACCTGCAGGGCCCGCTCGTCCAGGGGGCGGTGCTTGCTCTGGATCCTGTGGCGGGGGCGTCTCTGCAGGCCAGGGTCCTGGGCGCCCGTGAAGATGGAGCCATATTCCTGCAGGCGCCCTGGAGCAGGGTACTTGGCACTGGAGAACACCTGTGGA
p                                  #$%).1456667888889;=>@AAAAAAAAAAABGKPTWZ[\`eimqstx{~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~zvrmhd_^\XTPKGB><:6652/+'%$#

a score=240 mismap=7.38e-06
s KI635917 1147176 117 +   4339253 TACTTTTGTGGCTTTCTTGCTGCCTTTGATCCTCTCGATCTTGGCCTGTGCCAGATTGATGCGGTCACTGATGGACTGCAGGCGAGTCTGATTCCTCTCCACGTTTTGGGAAACCCT
s chr1       18264 117 + 248956422 TACCTTGATGGCCTTCTTGCTGCCCTTGATCTTCTCAATCTTGGCCTGGGCCAAGGAGACCTTCTCTCCAATGGCCTGCACCTGGCTCCGGCTCTGCTCTACCTGCTGGGAGATCCT
p                                  $'()*++-159;=AEJNQRSTTTTTTTTTTTTTTTTTTTTTTTTTTSSQQPOLHGGFFFFFFFFFFFFFFFFFFFFEEECA@@@@@@@@@@??????>>=<;987777642.-+*($

a score=252 mismap=4.27e-07
s KI635917 1161884 119 +   4339253 ACCTCGTGAAGACATCATTGGAGATTTTTTGGAGATACTGAAGAGCATCAGTGATCTGGTGAATGGACTCTTCCCGGCGCAGGTCCGGCTGGATTATTGGCACCGTGTACACCTGTCCT
s chr1       24735 119 + 248956422 ACCTGCTGAAGATGTCTCCAGAGACCTTCTGCAGGTACTGCAGGGCATCCGCCATCTGCTGGACGGCCTCCTCTCGCCGCAGGTCTGGCTGGATGAAGGGCACGGCATAGGTCTGACCT
p                                  #%&&''+/35667778889:>ACDEFIKLNOOSUVZ]^___``````````````````````````````````````````_]ZYWSOKFB=987666542/.-,,+*))))('&%#

a score=192 mismap=0.109
s KI635862   2418195 96 +  12728579 CACAGGCCGCTGCAGCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCCGA-GCAGCC-GGATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248924753 98 - 248956422 CATAGCTCACTGCAGCCTCGATCTCCCAAGCTCAAGTGATCCTCCCGCCTCAGCCTCCCATGTAGCTGGGACTACAGGTGTGTGACACCATGACTGGC
p                                   #$%%&&'(()*****************************************************************************))(''&&%%$#

a score=201 mismap=0.031
s KI635862   2418214 77 +  12728579 GAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC--GAGCAGCCGGATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248920976 79 - 248956422 GAACTCCTGGCCTCAAGCGATCCTCCCACCTCAGCCTCCCAAAGTGTTGGGATTATAGACATGAGCCACTGCACCTGGC
p                                   %),.//00000000000000000000000000000000000000000000000000000000000///////..-*)($

a score=214 mismap=0.00137
s KI635862   2418214 80 +  12728579 GAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC--GAGCAGCCGGATAACAGGCGCGCGCCACCGCGCCCGGCGAG
s chr1     248916719 82 - 248956422 GAACTCCTGACCTCAGGTGATCCACCCGCCTCAGCCTCCCAAAGTGTTGGGATTACAGGCACGAGCCACTGCATCCAGCCAG
p                                   %*.269;<==========================================================<<:87410/.+*($$#

a score=214 mismap=0.00209
s KI635862   2418214 77 +  12728579 GAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCCGAGCAGCC--GGATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248879450 79 - 248956422 GAATTCCTGGCCTCAAGTGATCTGCCCGCCTCAGCCTCCCAAAATGCTGGGATCACAGGCATGAGCCACCACGCCCGGC
p                                   %(+,158:;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;::862.*%

a score=198 mismap=0.0234
s KI635862   2418209 85 +  12728579 GCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC--GAGCAGCCGGATAACAGGCGCGCGCCACCGCGCCCGGCGAG
s chr1     248845257 87 - 248956422 GTCCCGAACTCCTGACCTTAGGTGATCCACCAGCATCGGCCTCCCAAAGTGCTGGGATTACAGGTGTGAGCCACTGTGTCCAGCCAG
p                                   #$&'(+.001111111111111111111111111111111111111111111111111111111111111111100//..-*)'$$#

a score=190 mismap=0.228
s KI635862   2418217 74 +  12728579 CTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCCG-AGCAGCCGG-ATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248845129 76 - 248956422 CTCCTGGGTTCAAGCGATTCTCCTGCCTCAGCCTCCCTAGTAGCTGGTATTACAGGTGCCTGCCACCACACCCAGC
p                                   $%&''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''&&%$#

a score=194 mismap=0.127
s KI635862   2418218 73 +  12728579 TCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC-GAGCAGCC-GGATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248826215 75 - 248956422 TCCTGGGTTCAAGCGATTCTCCTGCCTCAGCCTCCCAAGTAGCTGGGACTACAGGTGTGTGCCACCACACCCGGC
p                                   $'()))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))('$

a score=192 mismap=0.12
s KI635862   2418202 89 +  12728579 CGCTGCAGCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC-GAGCAGCC-GGATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248811964 91 - 248956422 CACTGCAGCCTCCACCTGCCAGGCTCAAGCAATTCTCCTGCCTCAGCCTCCCGAGTAGCTGGGATTACAGGTAAGTGCCACCACACCCAGC
p                                   #$'())))))))))))))*****************************************************************))))(&%#

a score=191 mismap=0.124
s KI635862   2418206 85 +  12728579 GCAGCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC-GAGCAGCC-GGATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248811595 87 - 248956422 GCAACCTCCACCTTCTGGGTTCAAGCAATTCTCCTGCCTCAGCCTCCCGAGTAGTTGGGATTACAGGCATGTGCCACCATGCCCAGC
p                                   "##$$$$$$%%'(())))))*********************************************************))))))(&%#

a score=207 mismap=0.00555
s KI635862   2418202 89 +  12728579 CGCTGCAGCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC-GAGCAGCC-GGATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248806800 91 - 248956422 CACTGCAACCTTCACCTTCTGGGTTCAAGCAATTCTCCTGCCTCAGCCTCCCGAGTAGCTGGGATTACAGGTGCCCGCCACCATACCCAGC
p                                   #%(,.///000001145567777777777777777777777777777777777777777777777777777777777776530/..-+(&$

a score=192 mismap=0.117
s KI635862   2418195 96 +  12728579 CACAGGCCGCTGCAGCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTC-CGAGCAGCCGG-ATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248803921 98 - 248956422 CATAGCTCACTGCAGCCTTGACCTCCTAGGCTAAAGCAATCCTCCCACCTTAGCCTCTCCAGTAGCTGGAACTACAGGCATGCATCACCATGTCCAGC
p                                   #$$%&&'(()*******************************************************************************)))(('&%#

a score=330 mismap=0.75
s KI636063 409336 105 +    789923 AAAATATGGAACGCTTCACGAATTTGCGTGTCATCCTTGCGCAGGGGCCATGCTAATCTTCTCTGTATCGTTCCAATTTTAGTATATGTGCTGCCGAAGCGAGCA
s chr1     157784 102 + 248956422 aaaataTGGAATGCTTCACAAATTTGCATGTCATTCTTTCACAGAGGCCGTGCCAA---TCTCTCTATTGTTCCAACTTAAGTATGTGTGCTACTGAGGCAAGCA
p                                 !"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!

a score=203 mismap=0.0127
s KI635862   2418202 85 +  12728579 CGCTGCAGCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTC-CGAGCAGCC-GGATAACAGGCGCGCGCCACCGCGCC
s chr1     248798099 87 - 248956422 CACTGCAGCCTCAACCTCCTGGGCTCAAGCAATCCTCCTGCCTCAGCCTCACAAGTAGCTGGGACTACAGGTGCTTGTCACCACACC
p                                   #%)-01233333333333333333333333333333333333333333333333333333333333333333322211110.*)'&$

a score=201 mismap=0.0167
s KI635862   2418211 79 +  12728579 CCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC--GAGCAGCCGGATAACAGGCGCGCGCCACCGCGCCCGG
s chr1     248796986 81 - 248956422 CTAGAACTTTTGACCTCAGGTGATCCGCCCGCCTCAGCCTCCCAAAGTGCTGGGATTACAGGCGTGATCCACCGCGCCCAG
p                                   #%),/11222222222222222222222222222222222222222222222222222222222222222222221/-)%#

a score=209 mismap=0.00211
s KI635862   2418209 87 +  12728579 GCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC--GAGCAGCCGGATAACAGGCGCGCGCCACCGCGCCCGGCGAGGG
s chr1     248792490 89 - 248956422 GTCTAGAACTCCTGACCTTGGGTGATCCACCCGCCTCGGCCTCCCAAAGTGCTGGGATTACAGGCATGAGCCACTGCGCCCAGCCTGGG
p                                   #$'(-158:;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;:9740,+($$$#"

a score=271 mismap=6.21e-09
s KI635901   1311390 127 +   5048387 TCTTTTTTTCCTTCTGTCAACACTTTTAGCTGAAGCCATAAAAACGACAAGTCCCGTCATAATCCCATACCCACTCTTTCAGTCCAATGCTGAAGATCTTAGTGTTGAAGGTTTACCCGAAGGTATT
s chr1     248791606 123 - 248956422 TCCCTttttctttttttt----CTTATAGCTCGAGCTGTAAAAGCCAAAGGTCCGGTGATGATCCCATACCCTTTTTTCCAGTCTCATGTTGAAGATTTTTATGTAGAAGGCCTtcccaaaggaatt
p                                    "##$&'''''''''''''''''+./048:;<<=ACEFGKNQRSSSSTTWX\_abbdefikloqrrrssssssssssssssssssrrrrrrrrqolid`^\XXWVTQQOMIEA@?>>=;8752.)('$

a score=202 mismap=0.0217
s KI635862   2418204 87 +  12728579 CTGCAGCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC-GAGCAGCC-GGATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248790144 89 - 248956422 CTGCAACCTCCACCTCCTGGGTTCAGGCAATTCTCCTGCCTCAGCCTCCCGAGTAGCTGGGATTACAGGTGCCTGCCACCATGCCCAGC
p                                   $')*+++++++,-/011111111111111111111111111111111111111111111111111111111111111111000/.+(&$

a score=191 mismap=0.145
s KI635862 2418202 87 +  12728579 CGCTGCAGCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCT-CCGAGCAGCC-GGATAACAGGCGCGCGCCACCGCGCCCG
s chr1      167211 89 + 248956422 CACTGCAACCTCTGCCTCCTGGGTTCAAGCGATTCTCCTGCCTTGGCCTCCCGAATAGCTGGGATTACAGACATGCGCCACCACACCCG
p                                 "#%&&&'''''''''()))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))((('&$

a score=192 mismap=0.162
s KI635862   2418214 77 +  12728579 GAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC--GAGCAGCCGGATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248786969 79 - 248956422 GAACTCCTGACCTCAAGTGATCCACCCATCTCGGCCTCCCAAAGTGCTGGGATTACAGGCATGAGCCACTGAGCCCAGC
p                                   $&((((((((((((((((((((((((((((((((((((((((((((((((((((((((((((((((((((((((('&%#

a score=194 mismap=0.0422
s KI635862   2418197  97 +  12728579 CAGGCCGCTGCAGCC-CAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC-GAGCAGCCG-GATAACAGGCGCGCGCCACCGCGCCCGGCGAG
s chr1     248786816 100 - 248956422 CAGTTCACTGCAACCTCTGCCCCCCTGGGTTCAAACAATTCTCCTGTCTCAGCCTCCTGAGTAGCTGAGATTACAGGCATGTGCCACCATGCCCGGCTAG
p                                    "###$%%&&&''''''''''')*,-.....................................................................-+($$#

a score=199 mismap=0.0162
s KI635862   2418217 87 +  12728579 CTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC-GAGCAGCCGG-ATAACAGGCGCGCGCCACCGCGCCCGGCGAGGGTTGCGGCT
s chr1     248781694 89 - 248956422 CTTTTGGGTTCAAGGGATTCTCCTGCCTCAGCCTCCCGAGTAGCTGGGATTACAGGCATGTGCCACCACACCCGGCTAATGTTGTAGTT
p                                   #$%&)*++,/12222222222222222222222222222222222222222222222222222222222222221/-,+++**(&%$#"

a score=226 mismap=0.000106
s KI635862   2418202 89 +  12728579 CGCTGCAGCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC-GAGCAGCCGG-ATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1     248781029 91 - 248956422 CACTGCAACCTCCACCTCCCGGGTTCAAGCGATTCTCCTGCCTCAGGCTCCCGAGTAGCTGGGATTACAGGCGCGCGCCACCATGCCCAGC
p                                   #%),./0000000226:<>>@@ABEGHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHGFC@;75530,('$

a score=194 mismap=0.0709
s KI635862 2418197 94 +  12728579 CAGGCCGCTGCAGCCCAGAACTCCTGGCCTCAAGCGATCCACCAGCCTCAGCCTCC-GAGCAGCCGG-ATAACAGGCGCGCGCCACCGCGCCCGGC
s chr1      175660 96 + 248956422 CAGCTCACTGCAATCTCTGCCCCCTGGGTTCAAGCAATTCTCCTTCCTCAGCCTCCTGAGTAGCTGGGATTACAGGCATGCACCACCACGCCTGGC
p                                 "###$%%&''''''''''''))++,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,+))'$

a score=211 mismap=0.00255
s KI635889   4462115 58 +   6429418 TACCTATAGATTGGTTGACAGCTTTCATGCACAAATTTTCCACCAGAATCTTTCCTGG
s chr1     248772848 58 - 248956422 TGCCTATGGACTGGTTGACAGCCTTCATGCACAGGTTCTCCACCAGAGCCTTTCTTGG
p                                   #$)-1356899::::::::::::::::::::::::::::::::::::9888640,+(%

a score=392 mismap=1e-10
s KI635917 1140607 159 +   4339253 ACCTTGTTcttgttccttctgcttcttcttcttctccatcttaatcAGCTTTGCATTGCGCAGTTTGGCCTTTCCAATGCCTCCGGCCTGACGAATGGATTCCAGCAGACTAGCCCTTCCATCCGAGGGTTTCACAACCTCCTTGGGTGCCCCCTGGAC
s chr1      186314 156 + 248956422 ACCTTGCTCctgctccttctg---cttcttcttctccAGCTTTCGCTCCTTCATGCTGCGCAGCTTGGCCTTGCCGATGCCCCCAGCTTGGCGGATGGACTCTAGCAGAGTGGCCCGGCCACCGGAGGGGTCGACCACTTCCCTGGGAGCTCCCTGGAC
p                                  %*.2578;>?BDEHJKKKKKKKKKPUY^bgkortuuuuuuuuuuuuvvwxxyyyz{~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~zvvuutsrpponkgcbaa``___^]\YVUROKGFDBA>;73.*%

a score=513 mismap=1e-10
s KI635917 1143406 188 +   4339253 TCTGTGCTGAACGTGGGCAGTTCAGGCATGTTGCTTCCTGGCACGGAGGGCGCGATGCCAGGTCCCAGGTCAGCACTGTACATGAGGTCATCTGCAATCCCAGGCAGATCCGGTAGGTAGGATGGGACGTCAATCTCTGGGACCTGACCCAGGTCTGGTACATAGAAGTAATTCTCAGTAGTCTGAAA
s chr1      187398 185 + 248956422 TCAGTGTGGAAGGTGGGCAGTTCTGGAATGGTGC---CAGGGGCAGAGGGGGCAATGCCGGGGCCCAGGTCGGCAATGTACATGAGGTCGTTGGCAATGCCGGACAGGTCAGGCAGGTAGGATGGAACATCAATCTCAGGCACCTGGCCCAGGTCTGGCACATAGAAGTAGTTCTCTGGGACCTGCAA
p                                  #$%()**+.1237<@EILOPQQQQQQQQQQQQQQQQQRRSTTUXY]begiimoquz}~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~}xsojigc_[VRMID?;:740+'&$$$$$$#"""