Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,380,807	G→A	*220* (TAG→TAA)	ycjU →	beta‑phosphoglucomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,380,807	0	G	A	100.0%	144.3 / NA	45	*220* (TAG→TAA)	ycjU	beta‑phosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/25);  total (20/25)

ACGCCATTAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAGCAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTCG  >  NC_000913/1380680‑1380929                                                                                                                                |                                                                                                                           aCGCCATTAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCa                                                                                                                 >  1:1379783/1‑139 (MQ=255)   gCCATTAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAAc                                                                                                               <  2:1324036/139‑1 (MQ=255)        tAACGCCAGCGGTATGTGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGc                                                                                                          >  2:1565319/1‑139 (MQ=255)             ccAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGc                                                                                                     <  2:1379783/139‑1 (MQ=255)              cAGCGGTATGTGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                    <  1:1565319/139‑1 (MQ=255)               aGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCtt                                                                                                   >  1:79352/1‑139 (MQ=255)               aGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                    <  2:265786/138‑1 (MQ=255)               aGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                    >  1:265786/1‑138 (MQ=255)                gCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCttt                                                                                                  <  1:265396/139‑1 (MQ=255)                 cGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTc                                                                                                 >  2:883022/1‑139 (MQ=255)                    tATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGtt                                                                                              <  1:1322954/139‑1 (MQ=255)                          cTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                    >  1:1313638/1‑127 (MQ=255)                          cTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                    <  2:1313638/127‑1 (MQ=255)                                      cGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAAt                                                                            <  2:79352/139‑1 (MQ=255)                                       ggCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATc                                                                           >  2:206180/1‑139 (MQ=255)                                         cgcgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTa                                                                         >  2:253526/1‑139 (MQ=255)                                             ggCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCttt                                                                                                  >  1:1166265/1‑110 (MQ=255)                                             ggCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCttt                                                                                                  <  2:1166265/110‑1 (MQ=255)                                             ggCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGAt                                                                     >  1:471640/1‑139 (MQ=255)                                             ggCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGAt                                                                     <  2:471640/139‑1 (MQ=255)                                                  aaCCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCa                                                                <  1:1496826/139‑1 (MQ=255)                                                   aCCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAg                                                               <  1:1408786/139‑1 (MQ=255)                                                     cGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGt                                                             <  1:253526/139‑1 (MQ=255)                                                             aaTTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATaa                                                                   >  2:436027/1‑125 (MQ=255)                                                             aaTTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATaa                                                                   <  1:436027/125‑1 (MQ=255)                                                                      tGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACt                                            <  2:452168/139‑1 (MQ=255)                                                                      tGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACt                                            <  1:179144/139‑1 (MQ=255)                                                                       gCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACaa                                                                                                                         >  2:1362851/1‑61 (MQ=255)                                                                       gCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACaa                                                                                                                         <  1:1362851/61‑1 (MQ=255)                                                                            cAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCt                                                                                                         <  2:123248/72‑1 (MQ=255)                                                                            cAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCt                                                                                                         >  1:123248/1‑72 (MQ=255)                                                                                          acCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGa                                                                                                                     >  2:105683/1‑46 (MQ=255)                                                                                          acCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGa                                                                                                                     <  1:105683/46‑1 (MQ=255)                                                                                            cTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACa                                                                                                                          <  1:353281/39‑1 (MQ=255)                                                                                            cTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                    >  2:601039/1‑61 (MQ=255)                                                                                            cTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                    <  1:601039/61‑1 (MQ=255)                                                                                              ggCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGtt                                                                                              >  2:935874/1‑65 (MQ=255)                                                                                              ggCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGtt                                                                                              <  1:935874/65‑1 (MQ=255)                                                                                               gCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTACCGATAAaga                   >  1:1533555/1‑139 (MQ=255)                                                                                                         tCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCgtg         >  2:744149/1‑139 (MQ=255)                                                                                                                tCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGt    >  1:869727/1‑137 (MQ=255)                                                                                                                tCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGt    <  2:869727/137‑1 (MQ=255)                                                                                                                tCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTCg  <  1:883022/139‑1 (MQ=255)                                                                                                                       aaaCGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAgg                                                              >  2:1029226/1‑72 (MQ=255)                                                                                                                       aaaCGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAgg                                                              <  1:1029226/72‑1 (MQ=255)                                                                                                                                |                                                                                                                           ACGCCATTAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAGCAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTCG  >  NC_000913/1380680‑1380929

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑