Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,233,794	G→A	intergenic (‑561/‑2113)	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,794	0	G	A	100.0%	88.2 / NA	28	intergenic (‑561/‑2113)	lysC/yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/13);  total (15/13)

GTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGA  >  NC_000913_3_bme_pgi/4233651‑4233936                                                                                                                                                |                                                                                                                                               gTTTACAATCAACATTCCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAg                                                                                                                                           <  2:369626/149‑1 (MQ=255)   ttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCt                                                                                                                                         >  2:34815/1‑149 (MQ=255)                 aCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTggg                                                                                                                           <  1:369956/149‑1 (MQ=255)                                aaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGATACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgc                                                                                                            >  1:485942/1‑149 (MQ=255)                                    cATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGGAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAACGCATTGGGCGGAGTCGGCcccgcccgc                                                                                                        >  2:80943/1‑148 (MQ=255)                                         caTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgc                                                                                                   <  2:485942/149‑1 (MQ=255)                                                           gcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCCTTAGGCGGAATCGGCCCCGGGCGATATGCGCGCGGCATTTCGCGCCg                                                                                 >  1:3986/1‑149 (MQ=255)                                                              aCAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAcc                                                                              >  2:511447/1‑149 (MQ=255)                                                                 aTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccg                                                                           <  1:199087/149‑1 (MQ=255)                                                                        aaaGTCACAATTCTCAAAATCAGAAGAGCATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAAAGGTAACAAAGATGAAGAACCATTGGGCGGAATAGGCCCCGCGCGATATGGGCGCGGCATTCAGCGCCGACACAGGCCGGtt                                                                    >  2:1750/1‑149 (MQ=255)                                                                              acaATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCg                                                              >  2:535950/1‑149 (MQ=255)                                                                                  ttCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGatat                                                          >  1:272307/1‑149 (MQ=255)                                                                                       aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCCATATTCGct                                                     >  2:331190/1‑149 (MQ=255)                                                                                             agaagaGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCt                                               <  2:8090/149‑1 (MQ=255)                                                                                               aagaGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGACATTCAGCGCCGACCCCGGCCGGGTCGGGCGATATTCGCTCTGCCTgg                                             >  1:418283/1‑149 (MQ=255)                                                                                                   gTATGGCTAATGGCAGGAGACGCGACGAAGCTTAGAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGgatga                                         <  2:424548/149‑1 (MQ=255)                                                                                                    tATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGgatgat                                        >  1:226147/1‑149 (MQ=255)                                                                                                        gCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtcttc                                    >  2:353978/1‑149 (MQ=255)                                                                                                         cTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtcttct                                   <  1:353978/149‑1 (MQ=255)                                                                                                          tAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTg                                  >  2:487061/1‑149 (MQ=255)                                                                                                           aaTGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGt                                 <  1:298037/149‑1 (MQ=255)                                                                                                                    agaCGGGAGGAAGCTTGAAGCAATGGTCACAAAGCTGAAAAAGCATTGGGCGGAATAGGCCTCGCGTGATATGCGCGCGGCATTCAGTGCCGACCTCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATTATCTTCTGTTCGACtttt                        <  2:143131/149‑4 (MQ=255)                                                                                                                        gcgACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGaa                    <  2:272307/149‑1 (MQ=255)                                                                                                                          gacgaAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGt                  >  2:221864/1‑149 (MQ=255)                                                                                                                             gaAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTAGATGATCTTATGTTCGACTGGTCGAAGTGcc               >  1:113406/1‑149 (MQ=255)                                                                                                                                gCTGGAAGCAACGGTCACAAAGCGGAAAAAGCGTTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGCCAGGTAGAAGTGTCggg            <  2:65579/149‑1 (MQ=255)                                                                                                                                          ccGGTCACAAAGGTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATGCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAAcgacga  <  2:335086/148‑1 (MQ=255)                                                                                                                                          aCGGTCACAAAGCTGAAAAAGCATTGGGCAGAATCGGCCCCGCGCGATATGCGCTCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAAcgacga  <  1:469499/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                               GTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGA  >  NC_000913_3_bme_pgi/4233651‑4233936

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