Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,962,959 68 bp→GCGGCCGC coding (131‑198/408 nt) ygdB ← DUF2509 domain‑containing protein YgdB

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2962959 2963026 68 86 [0] [0] 86 ygdB DUF2509 domain‑containing protein YgdB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 29629580 (0.000)81 (1.220)
+GCGGCCGC
74/478 0.0 100% coding (199/408 nt) ygdB DUF2509 domain‑containing protein YgdB
?NC_000913 2963027 = 0 (0.000)coding (130/408 nt) ygdB DUF2509 domain‑containing protein YgdB

GCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGAATCCTGTCACGCTATAGCTGGCATTCACCACGGTTTGCGGTTC  >  NC_000913/2962709‑2963206
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
gCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTAACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGg                                                                                                                                                                                                                                                         <  1:285983‑M1/251‑2 (MQ=255)
  cAACATTAACTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCACCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCAAGTAAACGCAAACAAACCTGGGCATAGGTTTCAGCGTACTGCGAGCCCTGAACTGCTGGgcg                                                                                                                                                                                                                                                       >  2:486022‑M1/1‑248 (MQ=255)
   aaCATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcgg                                                                                                                                                                                                                                                      >  2:235746‑M1/1‑247 (MQ=255)
    aCATTAGCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGATGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGATTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTCTCAGCGTACTGCCAGCACTGAACTGCTGGgcggc                                                                                                                                                                                                                                                     >  2:362674‑M1/1‑246 (MQ=255)
      aTTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGCAGCTGACATAACGCCCTCTCTTTCAGCGGACAACAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACACCGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAACCAAACCTGGGCATCGTTTTCAGCGTACTGCGAGCACTGAACTGCTGGtcggccg                                                                                                                                                                                                                                                   >  2:387987‑M1/1‑244 (MQ=255)
         aCCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcac                                                                                                                                                                                                                                                >  1:201821‑M1/1‑241 (MQ=255)
           cTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAACCAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCACTAAGGCTTCATTATCTCACAGCAACCGCAAACAAACCTGGGCATCGGTTTCAACGTACTGCGAGCACTGAACTGCTGGgcggccgcactg                                                                                                                                                                                                                                              >  2:324071‑M1/1‑239 (MQ=255)
              cGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgact                                                                                                                                                                                                                                           >  2:223775‑M1/1‑236 (MQ=255)
                    gCTAAAGCCTTGTTGATTCTTCAGGTAAGCGGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCTTGGCGAAAAAACAATGTTTCCATCAATTACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCGCTGAACTGCTGGgcggccgcactgtctttccat                                                                                                                                                                                                                                     <  2:328064‑M1/251‑22 (MQ=255)
                       aaaGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgct                                                                                                                                                                                                                                  >  1:312769‑M1/1‑227 (MQ=255)
                       aaaGACTTGTTGATTCGTCAGGGAAGATGACATAACGCCCTCGCTTTCAGCGGACAAAAATCGCTCCGGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCCCAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCTCTGAACTGCTGGgcggccgcactgcctttccatg                                                                                                                                                                                                                                    <  2:349937‑M1/249‑23 (MQ=255)
                            cttaacgattcttCAGGCCAGCTGACATAACGTCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTACATTATCTACCAGTAAACGCAAACACACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAAATGCCAGgcggccgccctgcctttccatgctcacgc                                                                                                                                                                                                                             >  2:62723‑M1/7‑222 (MQ=255)
                                 tgATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgaga                                                                                                                                                                                                                         <  1:26437‑M1/250‑34 (MQ=255)
                                  gATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCTTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagagg                                                                                                                                                                                                                       <  1:186289‑M1/251‑36 (MQ=255)
                                      cttcAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAACATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGGCTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaa                                                                                                                                                                                                                    >  2:358263‑M1/1‑212 (MQ=255)
                                        tcAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaac                                                                                                                                                                                                                 >  1:47803‑M1/1‑210 (MQ=255)
                                         cAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaac                                                                                                                                                                                                                 <  1:516851‑M1/250‑42 (MQ=255)
                                            ggAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcg                                                                                                                                                                                                             <  1:324071‑M1/251‑46 (MQ=255)
                                                    aCATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgtt                                                                                                                                                                                                     >  1:389267‑M1/1‑198 (MQ=255)
                                                      aTAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttgg                                                                                                                                                                                                   >  2:348187‑M1/1‑196 (MQ=255)
                                                          cGCCCTCTCTTTCAGCGGACAAAAATTGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAGACGCCTTCATAACAGACAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcattgcctttccttgctcacgcaagaggcaaaactgcgatcctgttggctc                                                                                                                                                                                                >  1:426204‑M1/1‑192 (MQ=255)
                                                                  ctTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgt                                                                                                                                                                                       >  1:292339‑M1/1‑184 (MQ=255)
                                                                   tttCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgt                                                                                                                                                                                       >  1:281325‑M1/1‑183 (MQ=255)
                                                                       aGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaat                                                                                                                                                                                    >  1:166728‑M1/1‑179 (MQ=255)
                                                                        gCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagc                                                                                                                                                                                 <  1:189228‑M1/251‑74 (MQ=255)
                                                                         cgtccaaaaATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGAAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACGGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaacttcgatcctgttggctcattccttgtaatagca                                                                                                                                                                                <  2:471482‑M1/247‑75 (MQ=255)
                                                                               aaaaTCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactac                                                                                                                                                                          >  2:15796‑M1/1‑171 (MQ=255)
                                                                               aaaaTCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactac                                                                                                                                                                          <  2:147436‑M1/251‑81 (MQ=255)
                                                                                               cgTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcat                                                                                                                                                          >  1:286602‑M1/1‑155 (MQ=255)
                                                                                               cgTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcat                                                                                                                                                          >  1:41501‑M1/1‑155 (MQ=255)
                                                                                                gTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatc                                                                                                                                                         >  2:507956‑M1/1‑154 (MQ=255)
                                                                                                     gAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGCACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggac                                                                                                                                                    >  1:316112‑M1/1‑149 (MQ=255)
                                                                                                      aaaaaaCAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATGATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggacc                                                                                                                                                   <  1:62723‑M1/251‑104 (MQ=255)
                                                                                                       aaaaaCAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTTATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggacc                                                                                                                                                   >  1:173828‑M1/1‑147 (MQ=255)
                                                                                                        aaaaCAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccag                                                                                                                                                 <  2:88501‑M1/251‑106 (MQ=255)
                                                                                                        aaaaCAATGTGTCCAGCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTAGCTGCCAGTAATCGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGggggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctaatgccttgtaatagcaaactacccaaaaccagcagcatcaggaccag                                                                                                                                                 <  2:101488‑M1/251‑106 (MQ=255)
                                                                                                         aaaCAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacacaaaaccagccgcatcaggacc                                                                                                                                                   >  2:18089‑M1/1‑145 (MQ=255)
                                                                                                          aaCAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccag                                                                                                                                                 >  1:545031‑M1/1‑144 (MQ=255)
                                                                                                              aagttTCCAGCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCAAAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctggtggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccag                                                                                                                                           <  2:481240‑M1/249‑112 (MQ=255)
                                                                                                                 tttCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtg                                                                                                                                         <  1:251798‑M1/250‑114 (MQ=255)
                                                                                                                   tCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgac                                                                                                                                       <  1:513735‑M1/250‑116 (MQ=255)
                                                                                                                      aTCAATGACGTCGCCTGTTCGGCACAACGAAACGCCTTCATAACCAGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAGACAAACCTGGGCATCGGTTTCATGGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgaggggcaaaactgcgatcctgttggctcagtccttgtaatatcaaactacccaaaacgagcatcatccggaccagagccagtgacgaa                                                                                                                                    >  2:95227‑M1/1‑132 (MQ=255)
                                                                                                                        cAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaac                                                                                                                                  >  1:59229‑M1/1‑130 (MQ=255)
                                                                                                                              cTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtactagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaaaacccttt                                                                                                                           >  2:168466‑M1/1‑124 (MQ=255)
                                                                                                                                   ccTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcg                                                                                                                      >  1:9757‑M1/1‑119 (MQ=255)
                                                                                                                                   ccTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcg                                                                                                                      >  1:536183‑M1/1‑119 (MQ=255)
                                                                                                                                       ttCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttc                                                                                                                  >  1:267957‑M1/1‑115 (MQ=255)
                                                                                                                                          gCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcac                                                                                                                >  2:127513‑M1/1‑112 (MQ=255)
                                                                                                                                             acaACGAAACGCCTTCACAACCGTCAATGAATAAGGCTTGATTATCTGCCAGTAAACGCAACCAAACCTGTGGATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgacgttccatgctcacgcgagaggcaaaactgcgatcctgtttgctcattccttgtactagcaaattacccaaaaccagcagcatcaggcccagagccagtgacgaaacaccctgtgcgcggttcacaggt                                                                                                            <  2:52968‑M1/251‑143 (MQ=255)
                                                                                                                                               aaCGAAACGCGTTCATAACCGGCAATCAATAAGTCTTCATTAGCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCATCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttg                                                                                                          <  2:351950‑M1/251‑145 (MQ=255)
                                                                                                                                                aCGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCCTCGGTTTCAGCGTACTGCGAGCACTGCACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaacatgcgatcctgttggctcattccttgtcatagcaaactaaccaaaaccagcagcctcaggaccagagccagtgacgcaacaccctttgcgcggctcacagggtg                                                                                                          >  2:238011‑M1/1‑106 (MQ=255)
                                                                                                                                                   aatcgccTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatc                                                                                                      >  2:499261‑M1/4‑103 (MQ=255)
                                                                                                                                                   aaaCGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatc                                                                                                      >  1:378780‑M1/1‑103 (MQ=255)
                                                                                                                                                        ccTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtc                                                                                                 <  1:10805‑M1/251‑154 (MQ=255)
                                                                                                                                                         cTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtc                                                                                                 >  2:176812‑M1/1‑97 (MQ=255)
                                                                                                                                                               aaCCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctat                                                                                          >  1:358509‑M1/1‑91 (MQ=255)
                                                                                                                                                               aaCCGGCAAGCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatccggttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctat                                                                                          <  2:536183‑M1/251‑161 (MQ=255)
                                                                                                                                                                   ggCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagct                                                                                      >  1:438706‑M1/1‑87 (MQ=255)
                                                                                                                                                                    gCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctg                                                                                     >  1:202900‑M1/1‑86 (MQ=255)
                                                                                                                                                                             aaGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcacc                                                                            >  1:193965‑M1/1‑77 (MQ=255)
                                                                                                                                                                              aGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccc                                                                           >  2:340619‑M1/1‑76 (MQ=255)
                                                                                                                                                                                  ttCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggt                                                                       >  1:44946‑M1/1‑72 (MQ=255)
                                                                                                                                                                                    cagtatcTGCAATTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGATCTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcgtgtcacgctatagctggcattcaccacggttt                                                                     <  1:438785‑M1/248‑182 (MQ=255)
                                                                                                                                                                                      ttATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgc                                                                   >  1:132430‑M1/1‑68 (MQ=255)
                                                                                                                                                                                          cTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggtt                                                               >  1:230521‑M1/1‑64 (MQ=255)
                                                                                                                                                                                            gCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttca                                                             <  1:358263‑M1/251‑190 (MQ=255)
                                                                                                                                                                                                 tAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagactt                                                        >  2:390231‑M1/1‑57 (MQ=255)
                                                                                                                                                                                                     cGCAACCAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactg                                                    >  1:94040‑M1/1‑53 (MQ=255)
                                                                                                                                                                                                      gCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactgg                                                   >  1:438498‑M1/1‑52 (MQ=255)
                                                                                                                                                                                                        aaacaaacCTGGGCATCGGTTTCAGCGTACTGAGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgcaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacatccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagactcactggc                                                  >  2:555181‑M1/1‑50 (MQ=255)
                                                                                                                                                                                                          acaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagc                                               >  1:176360‑M1/1‑48 (MQ=255)
                                                                                                                                                                                                            aaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcac                                             <  2:230521‑M1/251‑206 (MQ=255)
                                                                                                                                                                                                            aaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcac                                             <  2:41501‑M1/251‑206 (MQ=255)
                                                                                                                                                                                                               ccTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgc                                           >  2:432361‑M1/1‑43 (MQ=255)
                                                                                                                                                                                                                cTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcat                                         >  1:174729‑M1/1‑42 (MQ=255)
                                                                                                                                                                                                                  gggCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcata                                        <  1:387987‑M1/250‑211 (MQ=255)
                                                                                                                                                                                                                      aTCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcatattaac                                   >  1:169403‑M1/1‑36 (MQ=255)
                                                                                                                                                                                                                      aTCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcatattaac                                   >  1:396396‑M1/1‑36 (MQ=255)
                                                                                                                                                                                                                      aTCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcatattaac                                   <  2:552046‑M1/251‑216 (MQ=255)
                                                                                                                                                                                                                                  gTACTGCGAGCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctccttccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcatattaaccgtcaacaccgg                       >  1:497423‑M1/1‑24 (MQ=255)
                                                                                                                                                                                                                                           gCACTGAACTGCTGGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcatattaaccgtcaacaccggcgagaagcc              <  2:312769‑M1/251‑237 (MQ=255)
                                                                                                                                                                                                                                                     gctgGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcatattaaccgtcaacaccggcgagaagccgctgacatcc    >  1:132348‑M1/1‑5 (MQ=255)
                                                                                                                                                                                                                                                     gctgGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcatattaaccgtcaacaccggcgagaagccgctgacatcc    <  1:15796‑M1/251‑247 (MQ=255)
                                                                                                                                                                                                                                                     gctgGgcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcatattaaccgtcaacaccggcgagaagccgctgacat      <  1:141646‑M1/249‑245 (MQ=255)
                                                                                                                                                                                                                                                        gggcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcatattaaccgtcaacaccggcgagaagccgctgacatcctg  <  1:473070‑M1/250‑249 (MQ=255)
                                                                                                                                                                                                                                                        gggcggccgcactgactttccatgctcacgcgagaggcaaaactgcgatcctgttggctcattccttgtaatagcaaactacccaaaaccagcagcatcaggaccagagccagtgacgaaacacccttttcgcggttcacaggttgaatcctgtcacgctatagctggcattcaccacggtttgcggttcagacttactggcagcacgcatattaaccgtcaacaccggcgagaagccgctgacatcctg  >  2:36914‑M1/1‑2 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
GCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGAATCCTGTCACGCTATAGCTGGCATTCACCACGGTTTGCGGTTC  >  NC_000913/2962709‑2963206

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: