Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,689 2 bp→AG intergenic (‑456/‑2217) 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,6890CA100.0% 133.1 / NA 38intergenic (‑456/‑2218)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (22/16);  total (22/16)
*NC_000913_3_bme_pgi4,233,6900AG100.0% 130.3 / NA 38intergenic (‑457/‑2217)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/16);  total (22/16)

AAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCC  >  NC_000913_3_bme_pgi/4233546‑4233827
                                                                                                                                               ||                                                                                                                                         
aaataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCAttt                                                                                                                                      >  4:147244/1‑150 (MQ=255)
   taaaGGGTGAGCGGGGCGGTTGTCANCGATGGGGNCATNCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTc                                                                                                                                   >  3:338320/1‑150 (MQ=255)
                        gTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCtt                                                                                                              <  2:144797/150‑1 (MQ=255)
                                 ggggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCtg                                                                                                                                 <  1:458626/122‑1 (MQ=255)
                                 ggggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCtg                                                                                                                                 >  2:458626/1‑122 (MQ=255)
                                         gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                         <  4:140214/138‑1 (MQ=255)
                                         gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                         >  3:140214/1‑138 (MQ=255)
                                          cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTGCAATCTTCCAAAGTCACAATTCTCa                                                                                            >  3:307733/1‑150 (MQ=255)
                                           ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaa                                                                                           >  1:83347/1‑150 (MQ=255)
                                                tttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                      <  3:31786/150‑1 (MQ=255)
                                                 ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGCTTACCCTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTATGACTGGCGCTACAATCTTCCAAAGTCACACTTCTCAAAATCag                                                                                     >  3:332984/1‑150 (MQ=255)
                                                       cACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGt                                                                               <  3:4766/150‑1 (MQ=255)
                                                        aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                                                              >  2:445031/1‑150 (MQ=255)
                                                                     agtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTGCAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                 <  4:307733/150‑1 (MQ=255)
                                                                      gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaa                                                                                   >  4:62992/1‑131 (MQ=255)
                                                                      gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaa                                                                                   <  3:62992/131‑1 (MQ=255)
                                                                       tagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAagag                                                                                <  4:104049/133‑1 (MQ=255)
                                                                       tagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAagag                                                                                >  3:104049/1‑133 (MQ=255)
                                                                         gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAgaga                                                             <  3:290736/150‑1 (MQ=255)
                                                                              gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGcga                                                        >  4:34157/1‑150 (MQ=255)
                                                                                          ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAg                                            >  3:112638/1‑150 (MQ=255)
                                                                                               tGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACg                                       >  2:472961/1‑150 (MQ=255)
                                                                                                       gggTTTACACTCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCaaa                               >  4:163835/1‑150 (MQ=255)
                                                                                                                 tCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCtt                                                                                                              >  2:561439/1‑61 (MQ=255)
                                                                                                                 tCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCtt                                                                                                              <  1:561439/61‑1 (MQ=255)
                                                                                                                 tCAACATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                     >  1:477651/1‑150 (MQ=255)
                                                                                                                   aaCATTACGCTAACGGCACTAAAACCATAGCATTTTTTTTTTACTTTCACTACAATCTTCCAAAATCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGCCGCGACAAAAATTGAATCAACGGTCGCAAAGCTGAAAAAGCa                   >  2:203966/1‑150 (MQ=255)
                                                                                                                   aaCATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTg                                               <  3:90762/122‑1 (MQ=255)
                                                                                                                   aaCATTACGCTAACGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTg                                               >  4:90762/1‑122 (MQ=255)
                                                                                                                    ataTTACGCTAAAGGCACTAAAACAATAGCATTTTTCTGTGAGTGGTGATACAATCTTCCAAAGTAACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAt                  <  4:338320/148‑1 (MQ=255)
                                                                                                                           gCTAACGGCACTAAAACCATAGCATTTTTCTGTGAc                                                                                                                             <  3:229790/36‑1 (MQ=255)
                                                                                                                           gCTAACGGCACTAAAACCATAGCATTTTTCTGTGAc                                                                                                                             >  4:229790/1‑36 (MQ=255)
                                                                                                                              aaCGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAAt        >  4:87522/1‑150 (MQ=255)
                                                                                                                               aCGGCACTAGAAGCATAGCAGTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAATATCAGAAGGGTAGTGCTAATGGCAGGAGACCCGGCGAAGCTTCAATCAACGCTCTTATTGTTTAAAAAGCATTGGGCGGAATc       <  4:175651/150‑1 (MQ=255)
                                                                                                                                cGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTca                                                                                                     <  3:90311/55‑1 (MQ=255)
                                                                                                                                cGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTca                                                                                                     >  4:90311/1‑55 (MQ=255)
                                                                                                                                cGGCACTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGACAAAAAGCTAAAAAAGCATTGGGCGGAATCg      <  2:83347/150‑1 (MQ=255)
                                                                                                                                    aCTAAAACCATAGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGccc  >  1:553493/1‑150 (MQ=255)
                                                                                                                                               ||                                                                                                                                         
AAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCC  >  NC_000913_3_bme_pgi/4233546‑4233827

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