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,7080CT98.0% 146.9 / NA 49intergenic (‑475/‑2199)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/19);  minor base G (0/1);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCA  >  NC_000913_3_bme_pgi/4233564‑4233847
                                                                                                                                                |                                                                                                                                           
ccGGTTGTCAACGGTGGGGTCAAGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTCCACTCAACATTACGCAAACGACACCAAAACTATCACATTTTTCTGTTACTGgtgcaa                                                                                                                                         >  2:195827/2‑147 (MQ=255)
                     aTGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACAGTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAAtt                                                                                                                    <  1:73684/149‑1 (MQ=255)
                        cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTCATACATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                                                 <  1:20496/149‑1 (MQ=255)
                          gATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaa                                                                                                               >  1:390179/1‑149 (MQ=255)
                            tttttCACCCAGTCCTCGCGGTCAGGAGTTCAGCTAATAAAGGCTTCACTGCGCTAAGGGTTTACACTCAACATGACGGTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAGGTCACAATTCTCaaat                                                                                                             <  2:195576/149‑2 (MQ=255)
                                   tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaag                                                                                                      >  2:611321/1‑149 (MQ=255)
                                    ccACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                                     >  1:124880/1‑149 (MQ=255)
                                      aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGt                                                                                                   >  1:711408/1‑149 (MQ=255)
                                         ccTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAtt                                                                                                <  2:350619/149‑1 (MQ=255)
                                                ggcAGTAGTTCATCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGGCACAATTCTCAAAATCAGAAGAGTATTGCTAATg                                                                                         <  2:359051/147‑1 (MQ=255)
                                                gTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATg                                                                                         >  1:508730/1‑149 (MQ=255)
                                                 tCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTAAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATgg                                                                                        >  1:764675/1‑149 (MQ=255)
                                                  cAGTAGTTCAGCAAATAAAAGCTGCACCGCTCTAAGGGTTTACACTCAACATCACGCTAACGGCACCAAAACCAACACCTTTTACTGTGACTGGTGCTACAATCTTCCAAAGTCACAATCCTCAAAATCAGAAGAGTATTGCTAATGGc                                                                                       <  2:6656/149‑1 (MQ=255)
                                                    gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATGCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                                     <  2:599675/149‑1 (MQ=255)
                                                    ggtgtTCAGCTAATAAATGCTTCACTGCGCGAAGGGTTTACACTCTACATTACGCTAACGGCACCAAAACCATCCCAGTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                                     <  2:357181/146‑1 (MQ=255)
                                                       gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTTCTAATGGCAagag                                                                                  >  1:159310/1‑149 (MQ=255)
                                                           aGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGCCTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATACGAAAAGTATTGCTACCGGCAAAAGAcgc                                                                              >  2:679358/1‑149 (MQ=255)
                                                           aGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgc                                                                              >  1:151997/1‑149 (MQ=255)
                                                                       cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGa                                                                  <  2:52018/149‑1 (MQ=255)
                                                                        ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCt                                                                                                                                   >  1:648145/1‑83 (MQ=255)
                                                                        ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCt                                                                                                                                   <  2:648145/83‑1 (MQ=255)
                                                                               agcTACGCGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACgg                                                          <  1:346260/148‑1 (MQ=255)
                                                                                     gggTTTACACTCAACATTACGCTAACGGCCCCACCATCATCACATTTTTGTCTGCCTGGTCATCAACTCTTCAACAGTTACAACTCTGAAAATCAGAAGAGTTGTGCTAATTGCCAGAGACGCGGCTAAGCGTGATGCAACGGTCGCaa                                                    >  2:678168/1‑149 (MQ=255)
                                                                                          tACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATAGCAAGAGACGCGACGACGCTTGAAGCAACGGTCGCAAAGCTg                                               >  1:85632/1‑149 (MQ=255)
                                                                                             acTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaa                                            >  1:153308/1‑149 (MQ=255)
                                                                                              cTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATTGCAAGAGACGCGACGAAACTTGAAGCAACGGTCGCAAAGCTGaaaa                                           >  2:61111/1‑149 (MQ=255)
                                                                                              cTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAAAGGAAAGAGACGCGACGAAGCTTGAAGCAACGGGCGCAAAGCTGaaaa                                           >  2:589168/1‑149 (MQ=255)
                                                                                               tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                          >  1:499433/1‑149 (MQ=255)
                                                                                               tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                          >  1:483466/1‑149 (MQ=255)
                                                                                                cAACATAATCCTAATGGCACCAAAACCATCACATTTATCTTTGAGCAAGACTACAATCTTCCAAAGTCACAATTCGCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAATCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                         <  2:85632/149‑1 (MQ=255)
                                                                                                 aaCATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGc                                        >  1:206949/1‑149 (MQ=255)
                                                                                                 aaCATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTCAAGCAACGGTCGCAACGCAGAAAAAGc                                        >  2:464837/1‑149 (MQ=255)
                                                                                                  aCATTACGCCAACGGCACCAAAACCATCACATTTTGCTGTTACTGGTGCTAAAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAGTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCa                                       <  2:328988/149‑1 (MQ=255)
                                                                                                     ttACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAACCATTg                                    >  2:646251/1‑149 (MQ=255)
                                                                                                      taagatAACGGTACCAAAACCATGACATTTTTCTGTGAATGGTGATAAAATAATCCAAGGTGAAAATGCTCAAAACCAGTAGAGTATTGCTAATGGCAAGAGACGCGACGAGGCTTTTTGCAACGGTAAAAAAAAAAAAAAAGCATTgg                                   <  2:114674/144‑1 (MQ=255)
                                                                                                       aCGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAAGCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAGCGGTCGCAAAGCTGAAAAAGCATTggg                                  <  1:195827/149‑1 (MQ=255)
                                                                                                         gCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATATTCCAAAGTCACAAATCTAAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGCCGCAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTTGGCg                                >  1:711904/1‑149 (MQ=255)
                                                                                                            aaCGGCACCAAAACCATCACATATTAATAAAAATGGTGCTACAATCTGCCAAAGTCACAAATATCAAAACCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACAAAAAAAAAAAAAAAAAAGCATTGGGCGGaa                             <  2:711408/149‑1 (MQ=255)
                                                                                                             aCGGCACCAAAACCATCACATTTTTCTGTGACTAGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGAAAGAGACGCGACGAAGATTGAAGCACAGGTCGCAAAGCTGAAAAAGCATTGGGCGAAAt                            >  2:22291/1‑149 (MQ=255)
                                                                                                               ggCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCAAAAGTCACAATTCTAAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCAACGAAGCTTGAAGCAACGGTAGCAACACTGACAACGACTTAGGCGGAATCg                          >  2:320611/1‑149 (MQ=255)
                                                                                                                     aaaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACAGTGCGCGGAATCGGCCCcg                    >  2:808793/1‑149 (MQ=255)
                                                                                                                     aaaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTACAACAGACTTGGGCGGAATCGGCCCcg                    >  2:808165/1‑149 (MQ=255)
                                                                                                                          cATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa               <  1:454548/149‑1 (MQ=255)
                                                                                                                            tCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGata             >  1:443200/1‑149 (MQ=255)
                                                                                                                             cacaTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGatat            <  2:483466/149‑1 (MQ=255)
                                                                                                                              acaTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAACGCATTGGGCGGAATCGGCCCCGCGCGATATg           >  1:48777/1‑149 (MQ=255)
                                                                                                                               caTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgc          >  1:692342/1‑149 (MQ=255)
                                                                                                                               caTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATga          >  2:108408/1‑148 (MQ=255)
                                                                                                                                   tttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgcgc      <  1:528554/149‑1 (MQ=255)
                                                                                                                                       tgtgACTGGTGCTACAATCTACCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCa  <  1:209678/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
GCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCA  >  NC_000913_3_bme_pgi/4233564‑4233847

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