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%	152.9 / NA	51	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 (32/19);  total (32/19)

GGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGC  >  NC_000913_3_bme_pgi/4233566‑4233853                                                                                                                                               |                                                                                                                                                  ggTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTGCa                                                                                                                                             >  1:255841/1‑149 (MQ=255)      tCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCtt                                                                                                                                        >  2:775416/1‑149 (MQ=255)       cAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTc                                                                                                                                       <  2:537910/149‑1 (MQ=255)           gATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGAGCTAAGGGTTTATACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCaaa                                                                                                                                   <  2:721904/149‑1 (MQ=255)             tGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGt                                                                                                                                 >  2:62383/1‑149 (MQ=255)               gggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTca                                                                                                                               <  1:450412/149‑1 (MQ=255)                 gTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAACGTcaca                                                                                                                             >  1:128341/1‑149 (MQ=255)                 gTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTcaca                                                                                                                             >  1:131229/1‑149 (MQ=255)                    aTGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGATCGGGGTGTACACTCAACATTACGCTAACGGCACCAAAACCAGCACATTTTTCTGTGACTGGTGCTACAAGCTTCCAAAGTCACAAtg                                                                                                                          >  1:491768/1‑148 (MQ=255)                       cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                                                       >  2:793215/1‑149 (MQ=255)                           tttttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTAACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTCCAATCTTGCAAAGTCACAATTCTCaaaa                                                                                                                   >  1:764409/1‑149 (MQ=255)                            ttttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAAt                                                                                                                  <  1:62383/149‑1 (MQ=255)                              ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                                                <  1:501090/149‑1 (MQ=255)                                  tCCACTACTGGCGGTCACTAGTTCAGCTAATAAATGCTTCACTTCGCTAAGGGTTTACACTCAACATGATGCTAACGGCACCAAAACCATCACGTTTTTCTGTGACTGGTGCTAGAATCTTCCAAAGTGACAAGTCTCAAAAGCagaag                                                                                                            >  2:686454/1‑149 (MQ=255)                                   ccACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                                           <  1:392832/149‑1 (MQ=255)                                   ccACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                                           <  1:219055/149‑1 (MQ=255)                                     aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCAAAATTCTCAAAATCAGAAGAGt                                                                                                         >  1:472361/1‑149 (MQ=255)                                         cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTg                                                                                                     >  2:145453/1‑149 (MQ=255)                                          tGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGAATTGc                                                                                                    >  2:527178/1‑149 (MQ=255)                                            gCGGTCAGTAGTTCAGATAATAAATGCTTCACTGCGTTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTa                                                                                                  <  1:670667/149‑1 (MQ=255)                                                         cAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcg                                                                                     >  2:563935/1‑149 (MQ=255)                                                          cgCTCATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgc                                                                                    <  1:11790/148‑1 (MQ=255)                                                            cTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGcga                                                                                  >  1:437988/1‑149 (MQ=255)                                                                    tGCTTCACTGCGATAAGGATTTACAATCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTCTAATCTTCCAAAGTCACAATTCTCAAAACTGGAGTAGTATTACTAATGGCAAGAGACGCGACGAAGCtt                                                                          <  2:142712/149‑1 (MQ=255)                                                                       ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGaa                                                                       <  1:810203/149‑1 (MQ=255)                                                                               gcTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAaga                                                                                         <  1:287976/123‑1 (MQ=255)                                                                               gcTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAaga                                                                                         >  2:287976/1‑123 (MQ=255)                                                                               gcTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGt                                                               >  2:364191/1‑149 (MQ=255)                                                                                 tAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGtag                                                             >  1:336518/1‑147 (MQ=255)                                                                                    gggTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCACAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACAGTAGCaa                                                          >  2:232030/1‑149 (MQ=255)                                                                                      gTTTACACTCAACATTACGCTAACGGCCCCCACTCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTACTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAg                                                        >  2:391846/1‑149 (MQ=255)                                                                                      gTTTACACTCAACATGACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAg                                                        <  2:524610/149‑1 (MQ=255)                                                                                       tttACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGACGCTTGAAGCAACGGTAGCaagcc                                                       >  2:205867/1‑146 (MQ=255)                                                                                           cacTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGaa                                                   <  2:540508/149‑1 (MQ=255)                                                                                              tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGcga                                                                                  <  2:526336/115‑1 (MQ=255)                                                                                              tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGcga                                                                                  >  1:526336/1‑115 (MQ=255)                                                                                                        gCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCGTGAAGCAACGATAGCAAAGCTGAAAAAGCATTGGGCg                                      >  2:36610/1‑149 (MQ=255)                                                                                                            aCGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAAt                                  >  2:785882/1‑149 (MQ=255)                                                                                                            aCGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAAt                                  >  1:82942/1‑149 (MQ=255)                                                                                                            aCGGCACCAAAACCATCACATTTTTCTGAGACTGGTGCCACCATCGTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGTCGAAGCTTGAAGCAACGGTAGCATAGCTGAACCAGCATTGGGCGGagt                                  >  1:32544/1‑147 (MQ=255)                                                                                                              ggCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCCACGCAGCTTGAAGCAACGATCGCAACGCTGAAAAAGCATGTGCCGGAATCg                                >  2:227051/1‑149 (MQ=255)                                                                                                               gCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                               >  1:149871/1‑149 (MQ=255)                                                                                                                 accaaaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCCCAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcc                             <  1:224308/146‑1 (MQ=255)                                                                                                                 accaaaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGACGAGTATTGCTAATGGCAACAGACGCGACGAAGCTTGAACCAACGGTAGCACACCCGAAAACGCATTGGCCGGAATCGGcc                             >  2:369109/4‑149 (MQ=255)                                                                                                                   caaaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCAAAATTGTCTAAATCAGAAGAGTCTTGCTAATGGCACGAGACGCGACGAAGCTTGAAGCAACGGTAGCAGAGCTGAAAAAGCCTTGGTCGGAATCGGgccc                           >  2:151183/2‑149 (MQ=255)                                                                                                                                tttttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGc                                       <  1:16907/124‑1 (MQ=255)                                                                                                                                tttttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGc                                       >  2:16907/1‑124 (MQ=255)                                                                                                                                  tttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGACCATTAATGCGACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgcgc            >  1:405241/1‑149 (MQ=255)                                                                                                                                           cTGGTGCTACAATCTTCCAAAGGCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAg   <  1:733702/149‑1 (MQ=255)                                                                                                                                           cTGGTGCCACAATCTTCCAAAGTCACAATTTTCAAAATAAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAg   <  2:469033/149‑1 (MQ=255)                                                                                                                                            tGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgc  >  2:691302/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                  GGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGC  >  NC_000913_3_bme_pgi/4233566‑4233853

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