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)
          aCGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCaa                                                                                                                                      >  2:2465691/1‑150 (MQ=255)
              aGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACcac                                                                                                                                  >  2:2050148/1‑150 (MQ=255)
                  aGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTTTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAAc                                                                                                                                     <  4:1495988/143‑1 (MQ=255)
                   gCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCg                                                                                                                             >  4:1260721/1‑150 (MQ=255)
                     gATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGc                                                                                                                            >  4:918964/1‑150 (MQ=255)
                         ttAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCgg                                                                                                                       <  1:2439499/150‑1 (MQ=255)
                              ggCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTg                                                                                                                  <  1:885701/150‑1 (MQ=255)
                                  aTTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGt                                                                                                              <  4:787445/150‑1 (MQ=255)
                                  aTTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGt                                                                                                              >  3:1292574/1‑150 (MQ=255)
                                   tttCCATGCGGATACCGTCAGCGGTACGGGCTTGTTTGGTCAGCCAGAGACGCTGTTGCTCACGCAGGGCGGCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTc                                                                                                             >  4:1116869/1‑150 (MQ=255)
                                   tttCAATGCGGATACCGTCAGCGGTNCGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTc                                                                                                             >  1:2231383/1‑150 (MQ=255)
                                   tttCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTc                                                                                                             >  2:1471134/1‑150 (MQ=255)
                                        aTGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGatat                                                                                                        >  1:1774269/1‑150 (MQ=255)
                                         tGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGAtata                                                                                                       <  1:1572127/150‑1 (MQ=255)
                                              aTACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCg                                                                                                  <  2:2270106/150‑1 (MQ=255)
                                               tACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGt                                                                                                 <  3:61179/150‑1 (MQ=255)
                                                 ccGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAAGCGCCCCGGCGTTGCCGTCGATATAAATCGttt                                                                                               <  2:1439649/150‑1 (MQ=255)
                                                    tCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgtt                                                                                            >  1:812220/1‑150 (MQ=255)
                                                     cAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgttg                                                                                           <  1:1033357/150‑1 (MQ=255)
                                                       gCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTtgct                                                                                         <  4:302943/150‑1 (MQ=255)
                                                         ggTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTtgctgc                                                                                       <  2:1711196/150‑1 (MQ=255)
                                                         ggTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCCTCGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTtgctgc                                                                                       <  3:434329/150‑1 (MQ=255)
                                                           tACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTc                                                                                                                                                 >  3:1401475/1‑91 (MQ=255)
                                                           tACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTc                                                                                                                                                 <  4:1401475/91‑1 (MQ=255)
                                                            aCGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCAc                                                                                    <  1:382244/150‑1 (MQ=255)
                                                             cGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACg                                                                                   >  2:2543304/1‑150 (MQ=255)
                                                              gggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATc                                                                                                                              >  1:1449332/1‑106 (MQ=255)
                                                              gggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACgg                                                                                  >  1:2065915/1‑150 (MQ=255)
                                                               ggCTTGTTGGGTCAGCCAGCCACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTGGCGGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCTCTAGCGCCCCGGCGTTGCCGTCGAGATAAATCGTTTGTTGCTGCCACGGa                                                                                 <  2:790963/150‑1 (MQ=255)
                                                               ggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGa                                                                                 >  1:333071/1‑150 (MQ=255)
                                                                gCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGa                                                                                                           >  4:1500783/1‑123 (MQ=255)
                                                                gCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAg                                                                                >  3:1089902/1‑150 (MQ=255)
                                                                gCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAg                                                                                >  2:216653/1‑150 (MQ=255)
                                                                    gttgGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAg                                                                            >  2:1595372/1‑150 (MQ=255)
                                                                        ggTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGc                                                                         >  3:1199574/1‑150 (MQ=255)
                                                                         gTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGt                                                                                                              >  1:1678289/1‑111 (MQ=255)
                                                                          tCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAg                                                                            <  1:1189410/144‑1 (MQ=255)
                                                                          tCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATc                                                                      <  1:1471134/150‑1 (MQ=255)
                                                                          tCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATc                                                                      >  1:613743/1‑150 (MQ=255)
                                                                            aGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGcc                                                                                                                <  2:309603/106‑1 (MQ=255)
                                                                            aGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCaa                                                                    <  1:1305583/150‑1 (MQ=255)
                                                                             gCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAat                                                                   >  2:1401651/1‑150 (MQ=255)
                                                                              ccAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAata                                                                  >  2:118670/1‑150 (MQ=255)
                                                                                aGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATc                                                                >  2:444830/1‑150 (MQ=255)
                                                                                  acacGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCca                                                              >  1:898431/1‑150 (MQ=255)
                                                                                  acacGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCca                                                              <  2:442322/150‑1 (MQ=255)
                                                                                  acacGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCAACAATATCca                                                              >  3:423922/1‑150 (MQ=255)
                                                                                    acGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACgg                                                                                  <  2:889509/128‑1 (MQ=255)
                                                                                        tGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCGCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaa                                                        >  3:1188188/1‑150 (MQ=255)
                                                                                            gCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCt                                                                                                                                         <  3:860397/65‑1 (MQ=255)
                                                                                             cTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGc                                                   <  2:2216426/150‑1 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGttt                                                                                               <  2:114602/105‑1 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGttt                                                                                               >  1:114602/1‑105 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                  >  1:1949615/1‑150 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                  >  1:1406237/1‑150 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCAc                                                                                    <  2:2496032/114‑1 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                >  1:2000356/1‑150 (MQ=255)
                                                                                                 cGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTg                                               >  2:2370914/1‑150 (MQ=255)
                                                                                                   cAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACGATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGa                                             <  1:1726547/150‑1 (MQ=255)
                                                                                                    aGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAgg                                                                           <  3:224395/119‑1 (MQ=255)
                                                                                                    aGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAgg                                                                           >  4:224395/1‑119 (MQ=255)
                                                                                                        cGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCAt                                                                       >  3:1145954/1‑119 (MQ=255)
                                                                                                        cGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCAt                                                                       <  4:1145954/119‑1 (MQ=255)
                                                                                                        cGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGtt                                        >  1:2182864/1‑150 (MQ=255)
                                                                                                          tCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTga                                      >  2:1510345/1‑150 (MQ=255)
                                                                                                          tCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTga                                      >  4:631605/1‑150 (MQ=255)
                                                                                                           ccTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaa                                     <  3:1338734/150‑1 (MQ=255)
                                                                                                             tGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaac                                    <  3:1584995/150‑1 (MQ=255)
                                                                                                             tGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaac                                                       >  1:2504362/1‑130 (MQ=255)
                                                                                                             tGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaac                                                       <  2:2504362/130‑1 (MQ=255)
                                                                                                              gTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaacga                                  >  2:548923/1‑150 (MQ=255)
                                                                                                               tACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGtt                                        <  1:1052340/143‑1 (MQ=255)
                                                                                                               tACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGtt                                        >  2:1052340/1‑143 (MQ=255)
                                                                                                               tACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGt                                         >  1:552305/1‑142 (MQ=255)
                                                                                                               tACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaacgaa                                 <  1:2465691/150‑1 (MQ=255)
                                                                                                                  gcgGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTaa                                                                             >  2:1980174/1‑103 (MQ=255)
                                                                                                                    gggCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGc                            <  2:333071/150‑1 (MQ=255)
                                                                                                                        ttcttGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaa                                                        >  2:1828578/1‑119 (MQ=255)
                                                                                                                        ttcttGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaa                                                        <  1:1828578/119‑1 (MQ=255)
                                                                                                                        ttcttGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAgg                        <  1:2041790/150‑1 (MQ=255)
                                                                                                                          cttGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGAt                      >  2:27757/1‑150 (MQ=255)
                                                                                                                             gCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCAc                   <  1:1493353/150‑1 (MQ=255)
                                                                                                                              cTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGTGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCAcc                  >  1:293422/1‑150 (MQ=255)
                                                                                                                              atGATAATCACGCGCTCCGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGATGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCACCCAGCGTTGGAATGTTGAACGAACGGGCAGGGATCAcc                  <  2:1315536/149‑1 (MQ=255)
                                                                                                                                 ataataACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGc                                                   <  3:1496866/115‑1 (MQ=255)
                                                                                                                                  taataaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCaa                          <  4:132021/138‑1 (MQ=255)
                                                                                                                                  taataaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATc                     >  2:1660146/1‑143 (MQ=255)
                                                                                                                                  taataaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATc                     <  1:1660146/143‑1 (MQ=255)
                                                                                                                                     taaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGccc                                                                                                                          <  1:1367447/39‑1 (MQ=255)
                                                                                                                                      aaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTa               >  1:1753238/1‑145 (MQ=255)
                                                                                                                                      aaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTa               <  2:1753238/145‑1 (MQ=255)
                                                                                                                                        cGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgt                                                                                             <  3:994668/65‑1 (MQ=255)
                                                                                                                                            cTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTAt              >  2:1203042/1‑140 (MQ=255)
                                                                                                                                            cTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAggtggt    <  1:444830/150‑1 (MQ=255)
                                                                                                                                              aCGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTgc  >  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: ‑