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_pgi	4,233,708	0	C	T	100.0%	63.6 / NA	21	intergenic (‑475/‑2199)	lysC/yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/10);  total (11/10)

GCGGTTGTCAACGATGGGGTCATGCGGATTTTTCA‑TCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGC  >  NC_000913_3_bme_pgi/4233564‑4233839                                                                                                                                                  |                                                                                                                                    gtggTTGTCAACGATGGGGTCATCCGGATTTTTCACACCATACCCACTGGT‑AGGAGTTCAGCTAATAAGTGCTTCACGGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTa                                                                                                                                 <  1:124320/147‑1 (MQ=255)    gTTGTCAACGATGGGGTCATGCGGATTTTTCA‑TCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTCAGTAAGCATGTCTAAAACCATCACATTTTTTTGTGACTGGTGCTACaa                                                                                                                              >  1:423742/1‑149 (MQ=255)             gATGGGGTCATGCGGATTTTTCA‑TCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCaaa                                                                                                                     <  1:391430/149‑1 (MQ=255)                           gactgttCA‑TCCACTCCTGGCCGTCAGTAGTCCAGCTAATAAATGCTTCACTGCGCTAAGGGTTCACACTCAACATTACGCCAACGGCACTAAAACCATCACATTTTTCCGTTACCGGTGCTACAATCTTCCAAAGTCACAATCCTCaa                                                                                                       <  2:327598/144‑1 (MQ=255)                              ttttCA‑TCCGCTCCTGGCGTTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTCGTGCTACAATCTTCCAAAGTCACAATTCTCAatat                                                                                                    <  2:306035/149‑1 (MQ=255)                                     tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaag                                                                                              >  1:8090/1‑149 (MQ=255)                                        aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGt                                                                                           >  2:405728/1‑149 (MQ=255)                                                  gTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATg                                                                                 >  2:199087/1‑149 (MQ=255)                                                    cAGTAGTTCAGTTAATAAATGCTTCACGGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGc                                                                               <  2:192740/149‑1 (MQ=255)                                                    cAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCCCTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGc                                                                               >  2:469499/1‑149 (MQ=255)                                                                     aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAAAGACGCGACGAAGc                                                              >  2:396928/1‑149 (MQ=255)                                                                          ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAAGCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGaa                                                         <  2:214655/149‑1 (MQ=255)                                                                                  gcTGAGTGTTTACACGCAACATTACGCTAACGGCACTAAAACCATAACATTTTTCTGTGACTGATGCTAAAATATTCCAAAGTCGCAATTCTCAAAATCAGAAGATAGTTGCTCATGGCGAGAGACGCGAAGACACTTAAAACAACgat                                                 >  2:446810/1‑147 (MQ=255)                                                                                    tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCCAAATTCTCAAAATCAGAAGAGTATTGATAATGCCACGAGACGGGACGAAGCTTGACGCAACGGTca                                               >  2:481535/1‑148 (MQ=255)                                                                                    tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTca                                               <  1:405728/149‑2 (MQ=255)                                                                                         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)                                                                                                                                                  |                                                                                                                                    GCGGTTGTCAACGATGGGGTCATGCGGATTTTTCA‑TCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGC  >  NC_000913_3_bme_pgi/4233564‑4233839

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