Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,773,173 Δ1 bp coding (893/1914 nt) mtlA → mannitol‑specific PTS enzyme IICBA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,773,1730G.97.2% 470.6 / 0.4 110coding (893/1914 nt)mtlAmannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base . (55/51);  total (59/51)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GATTGAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGT  >  NC_000913/3772928‑3773414
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
gATTGAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGGGGTCTGGTTTCTCCGGCATCTCCGGGtttat                                                                                                                                                                                                                                              >  2:112165/1‑248 (MQ=255)
   tGAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGATAAACCGTCTGCGGGCGGTGCGGCACTCATCCACTTCCTGGGGGGCATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGATCACGCTGACTATCCTGGGCGGAGGACTGGTTTCTCCGGCATCTCCGGTTTTATCCt                                                                                                                                                                                                                                          >  1:53717/1‑251 (MQ=255)
    gAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCtt                                                                                                                                                                                                                                         >  1:214451/1‑251 (MQ=255)
     aaGCTAACCCATGTCCAGGTATGGGCGTGCTGCTGGCGTACATGGTCTTTGGTCGGGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACTAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGTGGTATGACTGGCGTGTTCACGCTGACTATCCTTGGCGGTGGGCGGGTTTCGCCGGCAACTCCGGTTTTATCCTTg                                                                                                                                                                                                                                        <  2:89489/251‑1 (MQ=255)
     aaGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCGTTGGGCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGAGCCTCGCAGTCATCCTCGGCGGTATGACTGTCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTg                                                                                                                                                                                                                                        <  2:368229/251‑1 (MQ=255)
         tAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTTTCCTTGCTGt                                                                                                                                                                                                                                    >  2:472566/1‑251 (MQ=255)
             ccaggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGCCTGGCGTGTTCAAGATGACTATCCTGGGCGGTGGGCTGGTTTCTCCGGCATCTCCGGTTTTTTCCTTGCTGCACTg                                                                                                                                                                                                                                >  1:294162/1‑251 (MQ=255)
                ggtccaggtATGGGCGTGCTGCTGGCGTCCATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCCATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTTTCTTTGCCGTACTGGCg                                                                                                                                                                                                                             >  2:48808/1‑251 (MQ=255)
                ggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGGGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTCCTGGCg                                                                                                                                                                                                                             >  2:364053/1‑251 (MQ=255)
                    ccggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCCCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTCCTGGCgctg                                                                                                                                                                                                                          >  2:449394/3‑247 (MQ=255)
                    caggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGa                                                                                                                                                                                                                         <  1:153434/251‑1 (MQ=255)
                      gttATGGGCGTGCTGCTGGCGTACATGTTCTTGGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGaca                                                                                                                                                                                                                       <  1:414924/249‑1 (MQ=255)
                       gtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGTTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGaccc                                                                                                                                                                                                                      >  2:160130/1‑249 (MQ=255)
                               gTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTATCGCTACACAGTCTGCGGGCGGAGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGCGTTCCCGCTGACTATCCTGGGCGGGGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAACAggg                                                                                                                                                                                                              >  2:245777/1‑250 (MQ=255)
                                 gctgctGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGGCTGGTTTCTCCGGCCTCTCCGGGTTTATCCTTGCTGTACTGGCGATGACCCCaaacggagc                                                                                                                                                                                                            >  2:448222/1‑245 (MQ=255)
                                                     ttGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTGTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTACGCTAACATcg                                                                                                                                                                                           <  2:414212/248‑1 (MQ=255)
                                                     tgggTCGTGGTAGCGCTAATCAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTAGCCACGAAATCTTCTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCGGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTGTATCCTTGCTGTACTGGCGATGACCCCAAATTGTGCTTACTTCGCTGACATcg                                                                                                                                                                                           <  2:113411/246‑1 (MQ=255)
                                                      tGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTTTCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCggg                                                                                                                                                                                       >  1:325919/1‑251 (MQ=255)
                                                        gTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGgtg                                                                                                                                                                                     >  1:583769/1‑251 (MQ=255)
                                                           gTGGTAGCGCTAAACAGTCTGCTGGCGTTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGgtgtgt                                                                                                                                                                                  <  2:325919/251‑1 (MQ=255)
                                                             ggTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGgtgtgtgg                                                                                                                                                                                >  1:216466/1‑250 (MQ=255)
                                                                   gcTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGGCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgcggcg                                                                                                                                                                          >  1:47322/1‑251 (MQ=255)
                                                                   cggcgaacaactatgactggagCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGAGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATACTGGGCGGGGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGAAGACAACAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgcggcg                                                                                                                                                                          >  2:627261/22‑251 (MQ=255)
                                                                       aaCAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATgg                                                                                                                                                                      <  1:112165/251‑1 (MQ=255)
                                                                            gcggcgggcggTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCGGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGtc                                                                                                                                                                 <  2:114588/248‑1 (MQ=255)
                                                                             cTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCCTCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGGGTGTGTGCGGCGATGGCTggct                                                                                                                                                                >  2:614140/1‑248 (MQ=255)
                                                                                 gggcggTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCtt                                                                                                                                                            >  1:592922/1‑251 (MQ=255)
                                                                                        gCGGCAATCATCCCCTTCCTGGGGGGTATCCACGAAATCTCCTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGt                                                                                                                                                      <  2:158598/250‑1 (MQ=255)
                                                                                              atcatcCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCt                                                                                                                                               >  1:172944/1‑251 (MQ=255)
                                                                                                    cacgtCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTg                                                                                                                                         <  1:505035/247‑1 (MQ=255)
                                                                                                    cACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTTGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGGGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTg                                                                                                                                         >  1:454606/1‑251 (MQ=255)
                                                                                                    cACTTCATGGGGGGCATCCACGCGATCTACTTCCCGTATCTGCTGATGAATCCGCGTCTGACCCTCACAGTCATCCTCGGCGGTATAACTGGCGTGTTCCCACTGACTATCCTGGGCGGTGGTCTGGTTTCCCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCCAAAGGTGCTTACTTCGCTTACATCGCGGGTGTGGGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTCTTTTg                                                                                                                                         >  2:256096/1‑251 (MQ=255)
                                                                                                     aCTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGc                                                                                                                                        >  1:563570/1‑251 (MQ=255)
                                                                                                          cTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGaa                                                                                                                                    >  1:216253/1‑250 (MQ=255)
                                                                                                            ggggggTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCGCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGaaaa                                                                                                                                  <  2:287203/250‑1 (MQ=255)
                                                                                                             gggggTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGAGATCACGCTGAATATCCTGGGCGGTGGTCTGGGTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTAGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAAcc                                                                                                                                >  1:269886/1‑251 (MQ=255)
                                                                                                                 gtctCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCTCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCGGGGCTGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCGTACTGCGCTAACATCGCGGGTGTGTGTTCGGCGATGGCTGTCTCCGTCTTTTTCTCTGCTATTTTGCTGAAATCCAGCa                                                                                                                            <  2:418213/248‑1 (MQ=255)
                                                                                                                 gTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGGCTGGGTTCTCCGGCATCTCCGGTTTTATCCTTGGTTTACTGGCGATGACACCCAAAGGGGCTTACTTCGCTAACATCGCGGGTGTGGGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGc                                                                                                                             >  1:443257/1‑250 (MQ=255)
                                                                                                                   aTCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCTCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGTGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATGTCCTGTTGTATCCTTGCGGTACTTGCTATGACATCAAAAGGTGCTTGCGGCGCTAACATCGCGGGTGTGTGGGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTAGTTTGCTGAAAACCAGCaa                                                                                                                           <  1:98402/250‑1 (MQ=255)
                                                                                                                    tCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTATGGTTTCTTCGGCCTCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGCTGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGcca                                                                                                                           >  1:516120/1‑247 (MQ=255)
                                                                                                                     ccACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAACCCAGCAAAGt                                                                                                                        <  1:119893/251‑1 (MQ=255)
                                                                                                                       aCGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGa                                                                                                                      <  1:324170/251‑1 (MQ=255)
                                                                                                                         gAAATCTACTTCCCGTATGCGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTATGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGCTGACACCAAAAGGGGCTTACTTCGCTAAACTCGCGGGGGTGGGTGCGGCGATGGCAGTCTCCTTCGTTGTCTATGCTATTTTGCTGTAAACCAGCAAAGTGaa                                                                                                                     >  2:399655/1‑250 (MQ=255)
                                                                                                                          aaaTCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACCATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCACAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGaa                                                                                                                     >  1:598853/1‑249 (MQ=255)
                                                                                                                           aaTCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTTGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAaaga                                                                                                                  <  1:399655/251‑1 (MQ=255)
                                                                                                                            aTCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAaaga                                                                                                                  >  1:407247/1‑250 (MQ=255)
                                                                                                                              cTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGGGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAaagacg                                                                                                                >  2:490320/1‑248 (MQ=255)
                                                                                                                               tACTTACCGTCTGTGCTGATGCATCCGCGTCTCATCCTCGCAGTCATCCCAGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGCCCGCGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAACAGGTGCTTACTTCCCTAACATCGCGGGGGTGTGGGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAAACAGAACAGTGAaagaaga                                                                                                               >  2:592275/1‑250 (MQ=255)
                                                                                                                                  ttCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGGCACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCACAGTGAaagaaga                                                                                                               >  1:385178/1‑247 (MQ=255)
                                                                                                                                     ccGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCATGGGCGGTGGTCTGGGTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTGCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAgatgat                                                                                                        <  1:612993/251‑1 (MQ=255)
                                                                                                                                         aTGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATAttt                                                                                                    >  2:456842/1‑250 (MQ=255)
                                                                                                                                           gtgCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAATAAGATGATATTGaa                                                                                                  <  1:153532/251‑1 (MQ=255)
                                                                                                                                             gcttatgaATCCGCGTCTGATCCTCGCAGTCTTCCTCGGCGGTATGTCTGGCGTGTTCGCGCTGACGATCCTGGGCGGTGGTCTGGGTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGaat                                                                                                 <  2:518448/246‑2 (MQ=255)
                                                                                                                                             gCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAagc                                                                                                >  1:41531/1‑251 (MQ=255)
                                                                                                                                             gCTGATGAATCCGCGTCTGATCCTAGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTTACTATCCTGGGCGGAGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGGTGTCTCTGCTATTTTGCTGAAAACCCGCAAAGTGAAAGAAGAAGATGATATTGGagc                                                                                                >  2:55269/1‑251 (MQ=255)
                                                                                                                                               tgatgaATCCGCGTCTGCTCCTCGCAGTCCTCCTCGGCGGTATGTCTGGCGTGTTCACGCTGACCATCCTGGGCGGTCGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGGGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTTTTTTGCTGAAAAcaa                                                                                                                               >  2:428242/1‑216 (MQ=255)
                                                                                                                                               tgatgaATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCGCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCa                                                                                                                               <  1:428242/218‑1 (MQ=255)
                                                                                                                                                   gaATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAAAAGATGATATTGAAGCAGCAAc                                                                                          >  1:278344/1‑251 (MQ=255)
                                                                                                                                                   gaATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGGGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGAGGAAATTGAAGCAGCaa                                                                                           >  2:23834/1‑250 (MQ=255)
                                                                                                                                                   gaATCCGAGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAAc                                                                                          <  1:304691/251‑1 (MQ=255)
                                                                                                                                                    aaTCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTTGTTTCGCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACt                                                                                         <  1:621927/251‑1 (MQ=255)
                                                                                                                                                              tGATCCTCGCATTCATCCTCGGCGGAATGACTGGCGTGTTTACGCTGACTTTCCTGGGCGGGGGTCTGGTTTCTCCGGCATCTCCGGGTTTATCCGTGCTGTACTGGCGATGACACCAAAAGGTGCTTAATTCGCTAACATCCAGGGTGTGTGTGCGGGGAGAGCTGACTCCTTCGTTGGCTCTGCTATTTTGCTGAAAACCCGCAAAGTGAAAGAAGAAGATGCTATTGGAGCAGCAAATCGTCGTaggc                                                                               >  2:293392/1‑248 (MQ=255)
                                                                                                                                                              tGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCCAAGTGAAagaagaaga                                                                                                            >  1:14788/1‑222 (MQ=255)
                                                                                                                                                              tGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAagaagaaga                                                                                                            <  2:14788/222‑1 (MQ=255)
                                                                                                                                                               gATCCTCGCAGTCATCCTCGGCGGTATGACTGGCTTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGc                                                                               <  1:606539/250‑1 (MQ=255)
                                                                                                                                                               gATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTa                                                                                  >  2:384988/1‑247 (MQ=255)
                                                                                                                                                               gATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTAGTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTAt                                                                                 <  2:216253/248‑1 (MQ=255)
                                                                                                                                                                aTCCTCGCAGGCATCCGCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTGCGCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTGCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACGCGTCGTATGCAg                                                                             <  2:455057/251‑1 (MQ=255)
                                                                                                                                                                  ccTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGGGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTAATGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGGTGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGCAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAgg                                                                            >  1:66428/1‑250 (MQ=255)
                                                                                                                                                                   cTCGCAGTCATCATAGGCGGTATGACTGGCGTGTTCAAGCTGACTATCCTGGGCCGGGGTCTGGTTTCTCCGGCATCTCCGGATTTATACTTGCTGTACTGGCGATGACACCAAAAGGTGCTTAATTCGCTAAAATAGCGGGGGTGTTTGCGGCGATGGCTGTCTCCTTCGTTATCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCCACTCGTCGTATGCAggcc                                                                          >  2:593454/1‑249 (MQ=255)
                                                                                                                                                                     cGCAGTCAGCCTCGGCGGTATGCCTGGCGTGGTCACGCTGACTATCCTGGGCTGTGGCCTGGTTGCCCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTGTGTGCTATTTGGCTGAAACCCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACAt                                                                        <  1:256096/251‑1 (MQ=255)
                                                                                                                                                                       cAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCCGCAAAGTGAAAGAAGAAGATGATATTGAag                                                                                                 >  1:518137/1‑224 (MQ=255)
                                                                                                                                                                       cAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAag                                                                                                 <  2:518137/224‑1 (MQ=255)
                                                                                                                                                                       cAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGa                                                                      <  2:179518/251‑1 (MQ=255)
                                                                                                                                                                         gTCATACTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGCGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGGGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCGGCAAAGTGAAAGAAGAAGCTGATATTGAAGCAGCAACtcttccgatgcgggaaatgaa                                                                     >  2:115336/1‑231 (MQ=255)
                                                                                                                                                                         gTCAGCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGaaa                                                                    <  2:587533/251‑1 (MQ=255)
                                                                                                                                                                          tCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGaa                                                                     >  1:570745/1‑249 (MQ=255)
                                                                                                                                                                           cATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGc                                                                  <  2:242514/251‑1 (MQ=255)
                                                                                                                                                                              ccTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATTACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTg                                                                <  1:420424/250‑1 (MQ=255)
                                                                                                                                                                                 cggcggTATTACTTGCGTGGTCACTCGGACTATCCTGGGCGGTGGTCTGGGTTCTCCGGCATCTCCGGTTTTAGCCTTGCTGGACTGGCGATGACAACAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGGGGGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGTTAGTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGAGATGGAAGCAGCAACTCGTCGTATGCAAGACATGAAAGCTGAGt                                                             <  2:608468/250‑1 (MQ=255)
                                                                                                                                                                                        aTGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAgg                                                      <  2:624380/250‑1 (MQ=255)
                                                                                                                                                                                          gACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGgcgc                                                   >  2:475620/1‑251 (MQ=255)
                                                                                                                                                                                          gACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCACCGGCATCTCCGGTTTTATACTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGAGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAAACAGCAAAGTGAAAGAAGACGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGATGAGTCTAAAggagc                                                   >  1:251631/1‑248 (MQ=255)
                                                                                                                                                                                                   gtgcACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGTGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCt                                          <  1:537357/248‑1 (MQ=255)
                                                                                                                                                                                                   gtTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGCGTCTAAAGGCGCATCTCCGCt                                          >  2:18218/1‑251 (MQ=255)
                                                                                                                                                                                                      cACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCCGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTc                                       >  2:499103/1‑251 (MQ=255)
                                                                                                                                                                                                       aCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTct                                      >  1:471039/1‑251 (MQ=255)
                                                                                                                                                                                                             aCTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACGGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCGCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGAGATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGc                                <  1:245777/251‑1 (MQ=255)
                                                                                                                                                                                                                aTCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCg                               <  1:615045/249‑1 (MQ=255)
                                                                                                                                                                                                                 tCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGc                                >  2:604043/1‑247 (MQ=255)
                                                                                                                                                                                                                 tCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGc                                <  1:604043/247‑1 (MQ=255)
                                                                                                                                                                                                                          gtggtCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACg                   >  1:386502/1‑251 (MQ=255)
                                                                                                                                                                                                                              tCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCt               >  2:484378/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ggTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCt               <  1:48808/248‑1 (MQ=255)
                                                                                                                                                                                                                                  gTTTCTCCGGCATCTCCGGTTTTATACTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGc           >  2:428019/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTTTCTCCGGCATCGCCGGTTTTATCCGTGCTGTAGTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTGCGTTGTCTCTGCTATGTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTTAAGCAGCAACTCGTCGGATGCAGGTCATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTGTGCTTGCGATGTGACTAACTGCCTGAGc           <  2:328044/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tttCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGcc          >  1:40287/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tcccGGTATCTCCGGTTTTATCCGTGCTGTACTGGCGATGACAACAAAAGGTGCTTACTTCGCTATCAGCGCGGGTGTGTGTGCGGCGATGGCCGTCGCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGAGGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAGGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCac        <  2:341226/248‑1 (MQ=255)
                                                                                                                                                                                                                                         cGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCATCTTTTTCTCTGCTATTGGGCTTAAAACCAGCAAAGTGAAAGAAGAAGATTATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCAGCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCacgtac    <  1:268343/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCacgtacgt  <  1:109214/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
GATTGAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGT  >  NC_000913/3772928‑3773414

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