Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,786,5360AG28.8% 141.8 / 12.1 67E426G (GAA→GGA) gpmM2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (7/40);  new base G (19/0);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 3.80e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCC  >  NC_000913/3786292‑3786785
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                           
aaCGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCACCTCACAGACAGAGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGAACGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTCGGTCCCACCGGGGTGATGGAAGCGGCGGGTAAAGCGGGGGGAGCGCTGGGTCACAGCGTGGgagaagg                                                                                                                                                                                                                                                      >  2:421953/1‑250 (MQ=255)
    gTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTcgcg                                                                                                                                                                                                                                                  <  1:248685/251‑1 (MQ=255)
      ggCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCCCCGAAAGAGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGGTCACTGCGTGGGAGAAGGCGCGa                                                                                                                                                                                                                                                 >  1:68435/1‑250 (MQ=255)
            gaagaGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGt                                                                                                                                                                                                                                          <  2:95577/251‑1 (MQ=255)
                agTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGATCACTGCGTGGGAGCAAGTGCGAAAgagggttgg                                                                                                                                                                                                                                      >  1:296983/1‑243 (MQ=255)
                      tCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCgt                                                                                                                                                                                                                                 <  1:29620/250‑1 (MQ=255)
                             cGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAACCGGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCTGTTAAAGCTGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGAc                                                                                                                                                                                                                         <  1:8202/251‑1 (MQ=255)
                                    cGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTTGATCCCTGCGTGGGAGAAGTCGCGCAAACGGGTGAGGCCgtgggggggccacggct                                                                                                                                                                                                                  >  1:2895/1‑239 (MQ=255)
                                     gCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGATGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAAAGGTTGAAGCGGTGGAGTACTGCGTGGGAGAAGACGCGGAAGCGGGTGAGGACGGGGGTGGACAAatgctg                                                                                                                                                                                                                 >  1:182424/1‑251 (MQ=255)
                                          cTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGATCACATCGTGGGAGAAGACGCGAAAACGGTTGAGTCCCTGGGGGGGCAAATGCTGATCAc                                                                                                                                                                                                            >  1:238682/1‑251 (MQ=255)
                                           tGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGACAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGCACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTTAAGCCGATGAAGCGGTGGGTAACTGCGGGGAAAGAGCCCCGCAAGCGGGTGAGGCCGTGGGGGGGCAActggcggtcacc                                                                                                                                                                                                           >  1:233845/1‑242 (MQ=255)
                                            gATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCAcc                                                                                                                                                                                                           <  2:434559/250‑1 (MQ=255)
                                            gATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCg                                                                                                                                                                                                          <  2:107571/251‑1 (MQ=255)
                                             aTCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGAGCACTGCGGGGGAGAAGGCGCGAAAGCGGGTGAGTCCGTGGGTGGGCAAACGATGATCACCGc                                                                                                                                                                                                         >  1:315492/1‑251 (MQ=255)
                                                  cTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGAGGTTGAAGGGCTGGGTCACTGCGTGGGAGAAGTCGCGAAAGAGGGTGAGTCCCTGGGGGGGCAACCGCTGGTAACCCCTGAcc                                                                                                                                                                                                    >  1:357194/1‑251 (MQ=255)
                                                      ccGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTg                                                                                                                                                                                                       <  1:129993/244‑1 (MQ=255)
                                                            gTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGCGCTGACCGCAAAACTGGTTGCGGCCATCAAGCGCGGCAAATACGACACCATCCTCTGTAACTCTCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTTAAGCGGCTGACGCGCCGGGACCCTGCGCGGGAGAAGGCGCCAAAGCGGTTGAGTCCGTTGGGGGGCCACTGCTCATCCCCCCTGACCACGGTAACGc                                                                                                                                                                                          >  1:272777/1‑251 (MQ=255)
                                                               gCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGa                                                                                                                                                                                       <  1:275309/251‑1 (MQ=255)
                                                                  acctacGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCa                                                                                                                                                                                    <  1:147574/251‑1 (MQ=255)
                                                                              cAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGGAGCGCGGGGTAACTGCGGGGGAGAAGGCGCGAAAGCGGTTGAGGCCGTGGGGGGGAAAATGCTGAGCACGCCGGACCCCGGTAAAGCTGAGCAGATGGGCGGTcc                                                                                                                                                                        >  2:239224/1‑251 (MQ=255)
                                                                              cAACCGGAAAGGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATcc                                                                                                                                                                        <  2:69009/251‑1 (MQ=255)
                                                                                 cctgAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCTGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCTAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGc                                                                                                                                                                     <  1:438167/248‑1 (MQ=255)
                                                                                      aaTGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGGGGACAACTGCTGATCACCGCGGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACgg                                                                                                                                                                <  2:272777/251‑1 (MQ=255)
                                                                                          aGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCa                                                                                                                                                            <  1:269693/251‑1 (MQ=255)
                                                                                                  agagCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGGTCACTGCGGGGGAGAAGGCGCGAAAACGGTTGAGGCCGGGGGGGGGCAACTGCTGATCACCGCTGACCACCGGAAAGCTGAGCGGAGGCGCGATTCGGCCACGGGGCCggcacac                                                                                                                                                     >  2:66531/1‑250 (MQ=255)
                                                                                                   gagCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAggcacacac                                                                                                                                                   <  2:197054/251‑1 (MQ=255)
                                                                                                        gACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGGGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGGTGCGCGATCCGGCAACGGGTCAggcacacacggc                                                                                                                                                <  1:434423/249‑1 (MQ=255)
                                                                                                             aaaaaCTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTTCGTGGGAGAAGGCGCGAAAGCGGGTGAGTCCGGGGGGGGACAACTGCTGATCACCGCTGACCCCGGGAACGCTGAGCAGATGCGCGGTCCGGCCACGGGGCCGGCACCCACGGCCCACAcc                                                                                                                                        >  2:110999/1‑251 (MQ=255)
                                                                                                              aaaaCTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCa                                                                                                                                       <  2:233845/251‑1 (MQ=255)
                                                                                                               aaaCTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTGGGAGCGCTGGGTCACTGCGTGGGAGGAAGCGCGAAAACGGGTGAGGCCGGGGGGGGGGAACTGCTGATCACCCCTGACCCCGGGAAAGCCGTGCAGGTGCGCGGTCCGGCCACGGGGCAGGGACACCCCGcacccacccc                                                                                                                                      >  2:271469/1‑249 (MQ=255)
                                                                                                                     gTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAACGGGTGAAGGGCTGGGTCACTGCGGGGGAGAAGGCGCGAAAGCGGGTGGATCCGGGGGGGGACAAATGCTGATCACCCCCGGCCACGGGAACGCTGGGCAGAAGCGGCATACGGCAACGGGGGCGGCACACAAGGCCCAcaaccaacatcc                                                                                                                                 >  2:355681/1‑247 (MQ=255)
                                                                                                                           gCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTcc                                                                                                                          <  1:452732/251‑1 (MQ=255)
                                                                                                                            ccATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGGTAAAGAGGGTGAAGCGCTGGATCACAGCGGGGGAGAAGGGGCGCAAGCGGGTGGGGCCGGGGGGGGGCAAATGCTGATAACCCCCGACCCAGGGAACGCGGAGCCGAAGCGCGATCCGGCCACGGGGCAGGCACACACGGGACACCCAAACCTGCCAGTTccc                                                                                                                         >  2:58109/1‑250 (MQ=255)
                                                                                                                                cAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTg                                                                                                                      <  1:327076/250‑1 (MQ=255)
                                                                                                                                cAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGGTAAAACGGTTGAAGCGCGGGAACAACGCGAGGGAGAAGGCGCGGAAGCGGTTGAGTCCGCGGGGGGGGAACTGCTGATTCCCGCTGACCAACGTAACGGCGCGCCGATGCGCGATCCGGCAACGGGGCCGGCACACACGGCACACACCCACCTGCCAGTTCCGCTgg                                                                                                                     >  1:466622/1‑250 (MQ=255)
                                                                                                                                 aaGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGAt                                                                                                                    <  1:81471/251‑1 (MQ=255)
                                                                                                                                 aaGAGCGGCAAATACGACACCATCATCTGTAACAATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGAt                                                                                                                    <  1:107575/251‑1 (MQ=255)
                                                                                                                                   gagCGGCAAATACGACACCATCATCAGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGAttt                                                                                                                  <  1:387799/251‑1 (MQ=255)
                                                                                                                                   gagCGGCAAATACGAAACGAACATCTGGAACTAGCCGAACGGCGACATGGGAGGTCACACAGGGGGGATGGAAGCGGCGGTTAAAGCGGGGGAAGCGCTGGTTCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAGCTGCTGATCACCGCGGACCACGGTAACGCTGAGCATATGCGCGATCCGGCATCGGGTCAAGCATACACGGTACACGCCAACCTGGCAGGTGCGCGGAttt                                                                                                                  <  2:460163/251‑1 (MQ=255)
                                                                                                                                     gCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTAc                                                                                                                <  2:106688/251‑1 (MQ=255)
                                                                                                                                          aaaTACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAACGGTTGGAGCGCTGGGTCACTGCGTGGGAGAAGGCGCGAAAGAGGGTGAGTCCGTGGGGGGGCAACTGCCGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGCACCGGCAAAGGGGCCGGCACACACGGCACACACCCACCTGCCCGGTCCGCTGAttttcagtgg                                                                                                           >  2:69156/1‑244 (MQ=255)
                                                                                                                                            aTACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGCTAAAGCCGTTGAAGCGCAGGATCACTGCGTGGGAGAAGTCGCGAAAGCGGGTGAGTCCGTGGGGGGGCAACTGCTGATCACCGCTGACCCCGGGAAAGCCGAGCAGATGCGCGATCCGGCAACGGGGCAGGCACACACGGCACACACCAAACTGCCAGGTCCACTGAttttcggtgggg                                                                                                         >  1:422791/1‑242 (MQ=255)
                                                                                                                                              ccGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGAt                                                                                                       <  2:72772/250‑1 (MQ=255)
                                                                                                                                                gACACCATCATCTGTACCTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGAGGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATaa                                                                                                     <  2:2895/251‑1 (MQ=255)
                                                                                                                                                     catcatCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGAGGGTGGAAAGCTGGGTCACTGCGTGGAAGAAGTCGCGGAAACGGGTGAGTCCGGGGGGGGGCAAATGCTGATCACCGCTGACCACGGGAACGCCGGGGAGAGTCGCGATCCGGCAACGGGGCCGGCACACAAGGCACACACCAACCCGCCCGGTCCGCTGATTTACGTTGGTGGTAAGaaac                                                                                                >  1:170117/1‑249 (MQ=255)
                                                                                                                                                             gTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAg                                                                                          <  1:368254/249‑1 (MQ=255)
                                                                                                                                                               aaCTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGGAAGCGGTTTAGGCCGTGGGGGGGCAACTGCTGGTCACCCCTGGCCACGGGAACGCTGGGCAGATGCGCGATTCGGCCACGGGTCAGGGACACACGGCACACACCAACCTGCCCGTTCCGCTGATTTTCGTTGGTGATAAGAACGTGAAAGCGGt                                                                                      >  1:324235/1‑251 (MQ=255)
                                                                                                                                                                 cTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGGGGCAGAAGTCGCGAAAGCGGGTGAGTCCGGGGGGGGACAAATGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGGCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGGTAAGAACGTGAAAGCgggtg                                                                                    >  1:39940/1‑248 (MQ=255)
                                                                                                                                                                   aTCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGAACACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGaa                                                                                  <  2:315492/251‑1 (MQ=255)
                                                                                                                                                                            gCGAAATGGTAGGTCACACCGGGGTGAGGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCaaa                                                                         <  2:182424/251‑1 (MQ=255)
                                                                                                                                                                             cGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGCGCTGGGTCACTGCGGGGGAGAAGGCGCGAAAGCGGGTGAGTCCGTGGGGGGACCACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAAAGGGGCAGGCCCAAACGGCACCCACCAACCCGCCCGTTCCGCTGATTTACGCTGGTGATAAGAACGGGAAAGCGGGTGAAGGCGGCaaa                                                                         >  2:359206/1‑250 (MQ=255)
                                                                                                                                                                                 aTGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGGAGAAGTCGCGAAAGCGGTTTAGTCCGTGGGGGGGCAACTGCTGGTCACCGCCGACCACGGGAACGCTGAGCAGATGCGCGATCCGGCAACGGGGCAGGCACCCACGGCACACACCAACCCGCCAGTTCCCCTGCTTTTCGTTGGGGGTAAGAACGGGAAAGCGGTTGAAGGAGGCAAACTTTc                                                                    >  1:395561/1‑251 (MQ=255)
                                                                                                                                                                                  gggTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCt                                                                   <  2:415296/250‑1 (MQ=255)
                                                                                                                                                                                   ggTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTg                                                                  <  1:372115/251‑1 (MQ=255)
                                                                                                                                                                                     tAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAACGGTTGAAGCGCTGGGTCACTGCGTGGAAGAAGTCGCGAAAGCGGGTGAGTCCGTGGGGGGACAACTGCTGGTCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCCACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGGTGAAGGCGGCAAACTTTCTGAc                                                                >  1:300308/1‑251 (MQ=255)
                                                                                                                                                                                      aGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGTTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCGGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCTCACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGAGAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACa                                                               <  1:453911/251‑1 (MQ=255)
                                                                                                                                                                                        gTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGGTCACTGCGTGGGAGAAGTCGCGAAAGAGGGTGAGTCCGTGGGGGGACAACTGCTGATCACCGCCGACCCCGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCCAACCCCAACCTGCCAGTTCCCCTGATTTACGTTGGTGATAAGAACGCGAAAACGGGTGAAGGCGGCAAAATTTCTGACATc                                                             >  1:81820/1‑251 (MQ=255)
                                                                                                                                                                                           ccacCAGGGTGAGGGAAGCGGCGGTTAAAGCGGTTGACGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATcgcg                                                          <  2:347742/250‑1 (MQ=255)
                                                                                                                                                                                                     gATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCAtg                                                 <  1:293835/250‑1 (MQ=255)
                                                                                                                                                                                                          aaGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCg                                           <  1:359206/251‑1 (MQ=255)
                                                                                                                                                                                                              ggcggTTAAAGCGGTTGAAGCGCTGGAGCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGa                                       <  2:70383/251‑1 (MQ=255)
                                                                                                                                                                                                                cggTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCGGACATCGCGCCGACCATGTTGTCGCTGATg                                     <  1:176122/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACTTCGCGCCGACCATGTTGTCGCTGATgg                                    >  2:370644/1‑251 (MQ=255)
                                                                                                                                                                                                                    tAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCCCGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGt                                  <  2:177724/250‑1 (MQ=255)
                                                                                                                                                                                                                              gAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTGGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATcc                        <  1:365342/250‑1 (MQ=255)
                                                                                                                                                                                                                                  cgcTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGAACGAGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGAAAAGGGTCAAGCAAACACGGAACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTATGACATCGCGCCGACAATGTTGTCGATCATGGGAATGGAACTCACGCa                    >  1:427352/1‑250 (MQ=255)
                                                                                                                                                                                                                                       gATCACTGCGTGGAAGAAGTCTCGAAAGCGGTTGAGGCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGAt              <  2:128725/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGt       <  1:316945/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   gaagaaGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGcc  <  1:150566/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaagaagTCGCTAAAGCGGTTGTGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGTTAACGCTGAGCAGATGCGCGATCCTGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTTCGTTGGTGATAAGAACGTGAAAGCTGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATTTTGGCGCTGATGTGGATGGAAATCCCGCAAGAGATGACTGGTAAGcc  <  1:278355/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                           
AACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCC  >  NC_000913/3786292‑3786785

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: