Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,233,680	T→C	intergenic (‑447/‑2227)	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_pgi	4,233,680	0	T	C	100.0%	154.2 / NA	47	intergenic (‑447/‑2227)	lysC/yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/18);  total (29/18)

ATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCG  >  NC_000913_3_bme_pgi/4233540‑4233823                                                                                                                                             |                                                                                                                                                aTTTATAAATAAAGGGTGAGCGGGGCGGGTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCAt                                                                                                                                         <  1:747647/149‑1 (MQ=255)         ataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACAtttt                                                                                                                                 >  1:52018/1‑149 (MQ=255)          taaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACCCTCAACATTACGCTAACGGCACCAAAACCATCACAttttt                                                                                                                                >  2:12536/1‑149 (MQ=255)          taaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAACCCATCACAttttt                                                                                                                                >  2:12384/1‑149 (MQ=255)               ggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCtgtg                                                                                                                           >  1:627545/1‑149 (MQ=255)                         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                                                                                                           <  2:648145/83‑1 (MQ=255)                                                                                                 ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCt                                                                                                           >  1:648145/1‑83 (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)                                                                                                                        tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACt                                                                                                                        >  2:98569/1‑47 (MQ=255)                                                                                                                        tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACt                                                                                                                        <  1:98569/47‑1 (MQ=255)                                                                                                                        tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                  >  1:483466/1‑149 (MQ=255)                                                                                                                        tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                  >  1:499433/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)                                                                                                                                             |                                                                                                                                                ATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCG  >  NC_000913_3_bme_pgi/4233540‑4233823

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: ‑