Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,141,192 Δ1 bp coding (1030/2502 nt) frwA ← putative PTS enzyme: Hpr, enzyme I and II components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,141,1890C.93.5% 818.2 / 49.2 184coding (1033/2502 nt)frwAputative PTS enzyme: Hpr, enzyme I and II components
Reads supporting (aligned to +/- strand):  ref base C (6/6);  new base . (87/85);  total (93/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 194 total aligned reads displayed.

CTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGC  >  NC_000913_3_bme_pgi/4141040‑4141332
                                                                                                                                                     |                                                                                                                                               
  gcgcTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTccc                                                                                                                                               >  1:2030347/1‑150 (MQ=255)
   cgcTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTccct                                                                                                                                              <  2:2030347/150‑2 (MQ=255)
   cgcTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTccct                                                                                                                                              <  2:751326/150‑2 (MQ=255)
          aCGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCaa                                                                                                                                      >  2:2465691/1‑150 (MQ=255)
              aGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACcac                                                                                                                                  >  2:2050148/1‑150 (MQ=255)
                  aGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTTTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAAc                                                                                                                                     <  4:1495988/143‑1 (MQ=255)
                   gCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCg                                                                                                                             >  4:1260721/1‑150 (MQ=255)
                     gATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGc                                                                                                                            >  4:918964/1‑150 (MQ=255)
                         ttAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCgg                                                                                                                       <  1:2439499/150‑1 (MQ=255)
                              ggCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTg                                                                                                                  <  1:885701/150‑1 (MQ=255)
                                  aTTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGt                                                                                                              <  4:787445/150‑1 (MQ=255)
                                  aTTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGt                                                                                                              >  3:1292574/1‑150 (MQ=255)
                                   tttCCATGCGGATACCGTCAGCGGTACGGGCTTGTTTGGTCAGCCAGAGACGCTGTTGCTCACGCAGGGCGGCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTc                                                                                                             >  4:1116869/1‑150 (MQ=255)
                                   tttCAATGCGGATACCGTCAGCGGTNCGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTc                                                                                                             >  1:2231383/1‑150 (MQ=255)
                                   tttCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTc                                                                                                             >  2:1471134/1‑150 (MQ=255)
                                        aTGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGatat                                                                                                        >  1:1774269/1‑150 (MQ=255)
                                         tGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGAtata                                                                                                       <  1:1572127/150‑1 (MQ=255)
                                              aTACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCg                                                                                                  <  2:2270106/150‑1 (MQ=255)
                                               tACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGt                                                                                                 <  3:61179/150‑1 (MQ=255)
                                                 ccGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAAGCGCCCCGGCGTTGCCGTCGATATAAATCGttt                                                                                               <  2:1439649/150‑1 (MQ=255)
                                                    tCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgtt                                                                                            >  1:812220/1‑150 (MQ=255)
                                                     cAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgttg                                                                                           <  1:1033357/150‑1 (MQ=255)
                                                       gCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTtgct                                                                                         <  4:302943/150‑1 (MQ=255)
                                                         ggTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTtgctgc                                                                                       <  2:1711196/150‑1 (MQ=255)
                                                         ggTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCCTCGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTtgctgc                                                                                       <  3:434329/150‑1 (MQ=255)
                                                           tACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTc                                                                                                                                                 >  3:1401475/1‑91 (MQ=255)
                                                           tACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTc                                                                                                                                                 <  4:1401475/91‑1 (MQ=255)
                                                            aCGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCAc                                                                                    <  1:382244/150‑1 (MQ=255)
                                                             cGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACg                                                                                   >  2:2543304/1‑150 (MQ=255)
                                                              gggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATc                                                                                                                              >  1:1449332/1‑106 (MQ=255)
                                                              gggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACgg                                                                                  >  1:2065915/1‑150 (MQ=255)
                                                               ggCTTGTTGGGTCAGCCAGCCACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTGGCGGATAATAACGAGCTACGGCTTCCCTGGCTCAACCTCTAGCGCCCCGGCGTTGCCGTCGAGATAAATCGTTTGTTGCTGCCACGGa                                                                                 <  2:790963/150‑1 (MQ=255)
                                                               ggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGa                                                                                 >  1:333071/1‑150 (MQ=255)
                                                                gCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGa                                                                                                           >  4:1500783/1‑123 (MQ=255)
                                                                gCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAg                                                                                >  3:1089902/1‑150 (MQ=255)
                                                                gCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAg                                                                                >  2:216653/1‑150 (MQ=255)
                                                                    gttgGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAg                                                                            >  2:1595372/1‑150 (MQ=255)
                                                                        ggTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGc                                                                         >  3:1199574/1‑150 (MQ=255)
                                                                         gTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGt                                                                                                              >  1:1678289/1‑111 (MQ=255)
                                                                          tCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAg                                                                            <  1:1189410/144‑1 (MQ=255)
                                                                          tCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATc                                                                      <  1:1471134/150‑1 (MQ=255)
                                                                          tCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATc                                                                      >  1:613743/1‑150 (MQ=255)
                                                                            aGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGcc                                                                                                                <  2:309603/106‑1 (MQ=255)
                                                                            aGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCaa                                                                    <  1:1305583/150‑1 (MQ=255)
                                                                             gCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAat                                                                   >  2:1401651/1‑150 (MQ=255)
                                                                              ccAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAata                                                                  >  2:118670/1‑150 (MQ=255)
                                                                                aGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATc                                                                >  2:444830/1‑150 (MQ=255)
                                                                                  acacGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCca                                                              >  1:898431/1‑150 (MQ=255)
                                                                                  acacGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCca                                                              <  2:442322/150‑1 (MQ=255)
                                                                                  acacGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCAACAATATCca                                                              >  3:423922/1‑150 (MQ=255)
                                                                                    acGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACgg                                                                                  <  2:889509/128‑1 (MQ=255)
                                                                                        tGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCGCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaa                                                        >  3:1188188/1‑150 (MQ=255)
                                                                                            gCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCt                                                                                                                                         <  3:860397/65‑1 (MQ=255)
                                                                                             cTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGc                                                   <  2:2216426/150‑1 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGttt                                                                                               <  2:114602/105‑1 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGttt                                                                                               >  1:114602/1‑105 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                  >  1:1949615/1‑150 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                  >  1:1406237/1‑150 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCAc                                                                                    <  2:2496032/114‑1 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                >  1:2000356/1‑150 (MQ=255)
                                                                                                 cGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTg                                               >  2:2370914/1‑150 (MQ=255)
                                                                                                   cAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACGATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGa                                             <  1:1726547/150‑1 (MQ=255)
                                                                                                    aGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAgg                                                                           <  3:224395/119‑1 (MQ=255)
                                                                                                    aGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAgg                                                                           >  4:224395/1‑119 (MQ=255)
                                                                                                        cGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCAt                                                                       >  3:1145954/1‑119 (MQ=255)
                                                                                                        cGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCAt                                                                       <  4:1145954/119‑1 (MQ=255)
                                                                                                        cGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGtt                                        >  1:2182864/1‑150 (MQ=255)
                                                                                                          tCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTga                                      >  2:1510345/1‑150 (MQ=255)
                                                                                                          tCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTga                                      >  4:631605/1‑150 (MQ=255)
                                                                                                           ccTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaa                                     <  3:1338734/150‑1 (MQ=255)
                                                                                                             tGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaac                                    <  3:1584995/150‑1 (MQ=255)
                                                                                                             tGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaac                                                       >  1:2504362/1‑130 (MQ=255)
                                                                                                             tGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaac                                                       <  2:2504362/130‑1 (MQ=255)
                                                                                                              gTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaacga                                  >  2:548923/1‑150 (MQ=255)
                                                                                                               tACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGtt                                        <  1:1052340/143‑1 (MQ=255)
                                                                                                               tACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGtt                                        >  2:1052340/1‑143 (MQ=255)
                                                                                                               tACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGt                                         >  1:552305/1‑142 (MQ=255)
                                                                                                               tACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaacgaa                                 <  1:2465691/150‑1 (MQ=255)
                                                                                                                  gcgGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTaa                                                                             >  2:1980174/1‑103 (MQ=255)
                                                                                                                    gggCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGc                            <  2:333071/150‑1 (MQ=255)
                                                                                                                        ttcttGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaa                                                        >  2:1828578/1‑119 (MQ=255)
                                                                                                                        ttcttGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaa                                                        <  1:1828578/119‑1 (MQ=255)
                                                                                                                        ttcttGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAgg                        <  1:2041790/150‑1 (MQ=255)
                                                                                                                          cttGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGAt                      >  2:27757/1‑150 (MQ=255)
                                                                                                                             gCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCAc                   <  1:1493353/150‑1 (MQ=255)
                                                                                                                              cTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGTGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCAcc                  >  1:293422/1‑150 (MQ=255)
                                                                                                                              atGATAATCACGCGCTCCGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGATGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCACCCAGCGTTGGAATGTTGAACGAACGGGCAGGGATCAcc                  <  2:1315536/149‑1 (MQ=255)
                                                                                                                                 ataataACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGc                                                   <  3:1496866/115‑1 (MQ=255)
                                                                                                                                  taataaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCaa                          <  4:132021/138‑1 (MQ=255)
                                                                                                                                  taataaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATc                     >  2:1660146/1‑143 (MQ=255)
                                                                                                                                  taataaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATc                     <  1:1660146/143‑1 (MQ=255)
                                                                                                                                     taaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGccc                                                                                                                          <  1:1367447/39‑1 (MQ=255)
                                                                                                                                      aaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTa               >  1:1753238/1‑145 (MQ=255)
                                                                                                                                      aaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTa               <  2:1753238/145‑1 (MQ=255)
                                                                                                                                        cGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgt                                                                                             <  3:994668/65‑1 (MQ=255)
                                                                                                                                            cTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTAt              >  2:1203042/1‑140 (MQ=255)
                                                                                                                                            cTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAggtggt    <  1:444830/150‑1 (MQ=255)
                                                                                                                                              aCGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTgc  >  2:687607/1‑150 (MQ=255)
                                                                                                                                                ggCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGt                >  1:1847399/1‑135 (MQ=255)
                                                                                                                                                ggCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGt                <  2:1847399/135‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGC  >  NC_000913_3_bme_pgi/4141040‑4141332

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: