Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,235,408	Δ1 bp	intergenic (‑2175/‑499)	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,235,407	0	A	.	100.0%	120.1 / NA	28	intergenic (‑2174/‑500)	lysC/yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (16/12);  total (16/12)

TGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAT  >  NC_000913_3_bme_pgi/4235263‑4235549                                                                                                                                                  |                                                                                                                                               tGCCTTCGATCAATGGGGCATGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAATGAGCGAGTGGTCGCGACACTTCCACGCAAGGGCTTGTCGCGCAT‑GTGTCTGCGCGAAGGAAAGCTTGatcac                                                                                                                                            >  2:454479/1‑145 (MQ=255)    ttcccagcAATGTGTCGGGGAACGGGGTAAGTGACTGGGGACGGAACTTTTGCCGGTTGTATCGGGCAAGGACGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgt                                                                                                                                        <  2:152911/142‑1 (MQ=255)     ttCGATCAATGGGCCGGGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGGTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCCTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgtc                                                                                                                                       <  1:288820/149‑1 (MQ=255)                gggCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCGGGTGCAGGTCTCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGc                                                                                                                            >  1:123505/1‑149 (MQ=255)                  gCGGGGAACTGGGTAAGGGACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCg                                                                                                                          <  1:264251/149‑1 (MQ=255)                        agcTGGGTAAGGAATTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCGTCCACGCAAGGGCTTCTAGCGCTT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                    <  2:136984/147‑1 (MQ=255)                               tAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                                                                             >  1:362338/1‑149 (MQ=255)                                  ggAATTGGCGACGGGACTGTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                          <  2:257033/149‑1 (MQ=255)                                   gatcgggCGACGGAACTGGTGGCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                         <  2:251568/144‑1 (MQ=255)                                       tGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccgcct                                                                                                     >  2:483434/1‑147 (MQ=255)                                                aaCTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCAAcc                                                                                            >  2:259307/1‑149 (MQ=255)                                                   ttttACCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                         >  1:403855/1‑149 (MQ=255)                                                       gCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAGAGCTTG‑ATCATCGTCGATATGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                     <  2:488719/149‑1 (MQ=255)                                                                aTCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                                                                            >  1:265788/1‑149 (MQ=255)                                                                 tCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGAAAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                           >  2:483457/1‑149 (MQ=255)                                                                  cGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgcg                                                                          <  2:362338/149‑1 (MQ=255)                                                                    gggaaggaaggGGCGCGGGGTCGCGGAGCTTCCACGCAAGGGCGTGTCGCGCAT‑CTGCAGGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAATCGATGCCTGATGCGACGCGGt                                                                        <  1:442602/146‑1 (MQ=37)                                                                    ggCAAGGAAGGGGCGAGAGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAAGCGTTCAAGACGCATCCGGCACCCGATGCATGATGCGAcgggat                                                                        >  2:103130/1‑145 (MQ=255)                                                                              gggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                              >  1:144045/1‑149 (MQ=255)                                                                                     tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAG‑CTGCATGCGCGAAAGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCCCGCCGCATACGGCAACCGATGCCGCATGCGACGGAGTCGCGACGTCTCAGGCCt                                                       >  2:99547/1‑149 (MQ=37)                                                                                        tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                    >  1:313647/1‑149 (MQ=255)                                                                                        tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCCCGCCCCATCAGGCCTACa                                                    >  2:293850/1‑149 (MQ=38)                                                                                                cTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATg                                            >  1:311084/1‑149 (MQ=255)                                                                                                      cGCAAGGGCTTGTCGCGCAT‑CTGCAGGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCCTCTTATCAGGCCTACAGGTCGATGCCGata                                      <  2:519848/149‑1 (MQ=255)                                                                                                             gCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACGGGTCGATGCCGATATGTACAt                               <  2:473707/149‑1 (MQ=255)                                                                                                                 ttCGCGCATGCAGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGt                            >  1:302269/2‑149 (MQ=25)                                                                                                                     cgcAT‑CTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATGAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                       >  1:352709/1‑149 (MQ=255)                                                                                                                           cTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCAGTTAAGTCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGCCGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTCCATCGTATTCGGCatt                  >  2:358601/1‑147 (MQ=37)                                                                                                                                           aaGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCTTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat  <  1:17900/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                               TGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCAT‑CTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAT  >  NC_000913_3_bme_pgi/4235263‑4235549

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