Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,708 C→T intergenic (‑475/‑2199) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,233,7080CT97.7% 299.4 / ‑2.2 87intergenic (‑475/‑2199)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (46/38);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCA  >  NC_000913_3_bme_pgi/4233562‑4233851
                                                                                                                                                  |                                                                                                                                                
gggCGGTTGTCAACGAAGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCt                                                                                                                                               >  2:555333/1‑150 (MQ=255)
   cGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACa                                                                                                                                            <  2:117284/150‑1 (MQ=255)
      ttGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATc                                                                                                                                         <  1:555333/150‑1 (MQ=255)
         tCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTc                                                                                                                                      <  2:344134/150‑1 (MQ=255)
                 ggggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTca                                                                                                                              <  4:331619/150‑1 (MQ=255)
                 ggggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGATCTACAATCGACCAACGTca                                                                                                                              >  2:553395/1‑150 (MQ=255)
                        tGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtct                                                                                                                       <  1:448755/150‑1 (MQ=255)
                          cGGATTTTTCATCCACTCCTGGCGGTCAGGAGGTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTACCGGCACTAAAACCATCACCTTTTTCTGTGACGTGTGCTAAAATCCTCGCAAGTCATAATTCTCa                                                                                                                     >  4:197680/1‑150 (MQ=255)
                           ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaa                                                                                                                    <  3:233844/150‑1 (MQ=255)
                            gATTTTTCATCCACTCCTGGCGGTCAGTAGATCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACGGGTGCTACAATCTTCCAAAGTCACAATTCTCaaa                                                                                                                   <  2:519082/150‑1 (MQ=255)
                                 ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCag                                                                                                              <  1:75921/150‑1 (MQ=255)
                                    aTCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaag                                                                                                           >  2:627334/1‑150 (MQ=255)
                                    aTCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaag                                                                                                           >  4:661084/1‑150 (MQ=255)
                                      ccACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAagag                                                                                                         >  4:128663/1‑150 (MQ=255)
                                            cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                   >  4:441949/1‑150 (MQ=255)
                                                  gTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                                               <  1:487104/132‑1 (MQ=255)
                                                  gTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                                               >  2:487104/1‑132 (MQ=255)
                                                   tCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGc                                                                                            >  1:415126/1‑150 (MQ=255)
                                                     agtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                                          >  4:308218/1‑150 (MQ=255)
                                                      gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATc                                                                                                                                         <  2:80970/102‑1 (MQ=255)
                                                      gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATc                                                                                                                                         >  1:80970/1‑102 (MQ=255)
                                                      gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                   <  4:81030/140‑1 (MQ=255)
                                                      gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                   >  3:81030/1‑140 (MQ=255)
                                                       tagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAaga                                                                                        <  2:141043/150‑1 (MQ=255)
                                                         gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAgaga                                                                                      >  3:473550/1‑150 (MQ=255)
                                                          ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                   <  1:112517/136‑1 (MQ=255)
                                                          ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                   >  2:112517/1‑136 (MQ=255)
                                                             aGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                                                      >  1:399903/1‑114 (MQ=255)
                                                             aGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                                                      <  2:399903/114‑1 (MQ=255)
                                                             aGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgct                                                                                  >  1:74735/1‑149 (MQ=255)
                                                              gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgcta                                                                                 <  1:627334/150‑3 (MQ=255)
                                                                taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacg                                                                               >  1:770648/1‑150 (MQ=255)
                                                                     aaTGCTTCACTGCGCTAATGGTTTACGCTCAACGTGACGCTAACTTCACTATTTTCAGCACATTTGTCTGGGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCt                                                                          <  3:228732/150‑1 (MQ=255)
                                                                     aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCt                                                                          >  2:696275/1‑150 (MQ=255)
                                                                     aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCt                                                                          >  4:192720/1‑150 (MQ=255)
                                                                        gCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGa                                                                       <  3:661084/150‑1 (MQ=255)
                                                                         cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                      <  3:209731/150‑1 (MQ=255)
                                                                         cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                      >  4:209731/1‑150 (MQ=255)
                                                                          ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                     <  4:393303/150‑1 (MQ=255)
                                                                             aCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAc                                                                                     >  2:913814/1‑131 (MQ=255)
                                                                             aCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAc                                                                                     <  1:913814/131‑1 (MQ=255)
                                                                              cTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAAc                                                                 <  4:473550/150‑1 (MQ=255)
                                                                               tGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACg                                                                >  2:340071/1‑150 (MQ=255)
                                                                                 cgcTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGt                                                              >  1:1020812/1‑150 (MQ=255)
                                                                                   cTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCg                                                            <  3:308218/150‑1 (MQ=255)
                                                                                    tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                      >  2:725474/1‑139 (MQ=255)
                                                                                    tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                      <  1:725474/139‑1 (MQ=255)
                                                                                    tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                                           >  2:718151/1‑150 (MQ=255)
                                                                                     aaGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCa                                                          >  2:60033/1‑150 (MQ=255)
                                                                                        ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                                       >  4:586594/1‑150 (MQ=255)
                                                                                        ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                                       >  2:733070/1‑150 (MQ=255)
                                                                                        ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                                       <  2:289449/150‑1 (MQ=255)
                                                                                               acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                                >  1:259096/1‑150 (MQ=255)
                                                                                               acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                                >  1:242688/1‑150 (MQ=255)
                                                                                               acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                                >  1:242687/1‑150 (MQ=255)
                                                                                                 tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaaa                                              >  2:439068/1‑149 (MQ=255)
                                                                                                  cAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaaac                                             >  3:643999/1‑148 (MQ=255)
                                                                                                         aCGCTNACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGc                                      <  3:197680/150‑1 (MQ=255)
                                                                                                         aCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGc                                      <  2:903838/150‑1 (MQ=255)
                                                                                                          cGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCg                                     <  2:242687/150‑1 (MQ=255)
                                                                                                            cTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGa                                   <  1:60033/150‑1 (MQ=255)
                                                                                                                cGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCg                               >  4:226365/1‑150 (MQ=255)
                                                                                                                   cACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGcc                            >  3:437061/1‑150 (MQ=255)
                                                                                                                     cTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACACTTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGcccc                          >  1:727028/1‑150 (MQ=255)
                                                                                                                      tacaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCCAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCcg                         >  2:1012133/4‑150 (MQ=255)
                                                                                                                      tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATc                                >  3:118805/1‑143 (MQ=255)
                                                                                                                      tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATc                                <  4:118805/143‑1 (MQ=255)
                                                                                                                      tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCcg                         >  3:41368/1‑150 (MQ=255)
                                                                                                                         aaCCATCACATTTTTCTGTGACTGGTGCTACAATCt                                                                                                                                        <  4:181418/36‑1 (MQ=255)
                                                                                                                         aaCCATCACATTTTTCTGTGACTGGTGCTACAATCt                                                                                                                                        >  3:181418/1‑36 (MQ=255)
                                                                                                                         aaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc                      >  4:300348/1‑150 (MQ=255)
                                                                                                                          aCCATCACATTTCTCTTTGACTGGTGCTACAATTTTCCAAAGGCACGATTCACAAAATCAGAAGAGTGTTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCCACGGCCGCAAAGCTGACAAAACATTGGGCGGAATCGGCCCcgcgcg                     <  2:267660/150‑1 (MQ=255)
                                                                                                                           ccATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGa                   <  3:192720/150‑1 (MQ=255)
                                                                                                                           ccATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGa                   <  2:770648/150‑1 (MQ=255)
                                                                                                                            cATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGat                  >  2:583125/1‑150 (MQ=255)
                                                                                                                                acaTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATgc              >  4:385796/1‑150 (MQ=255)
                                                                                                                                acaTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATAGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATgc              >  3:179406/1‑150 (MQ=255)
                                                                                                                                 caTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                                           >  1:1444/1‑105 (MQ=255)
                                                                                                                                 caTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                                           <  2:1444/105‑1 (MQ=255)
                                                                                                                                    ttttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATgcgcgc          >  3:414801/1‑150 (MQ=255)
                                                                                                                                     tttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGcgcgcg         <  3:456605/150‑1 (MQ=255)
                                                                                                                                      ttCTGTGACTGGTGCTACAAGCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCgg        <  2:615171/150‑1 (MQ=255)
                                                                                                                                        cTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                                               >  1:296942/1‑46 (MQ=255)
                                                                                                                                        cTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                                               <  2:296942/46‑1 (MQ=255)
                                                                                                                                        cTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCGCGCGCATATCGCGCggagcc    <  2:1020812/150‑5 (MQ=255)
                                                                                                                                         tgtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGCATATCGCGCggggc     <  2:504694/150‑4 (MQ=255)
                                                                                                                                          gtgACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCGCGCGCATATCGCGCggggccga  <  1:718151/150‑7 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCA  >  NC_000913_3_bme_pgi/4233562‑4233851

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