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%	107.8 / NA	34	*220* (TAG→TAA)	ycjU	beta‑phosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/16);  total (18/16)

GGCATTGACGCCATTAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAGCAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTG  >  NC_000913/1380673‑1380922                                                                                                                                       |                                                                                                                    ggCATTGACGCCATTAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACaaa                                                                                                                 <  1:361044/139‑1 (MQ=255)                                 cTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTacaaca                                                                                 >  1:1360149/1‑139 (MQ=255)                                  tCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACat                                                                                >  1:1360150/1‑139 (MQ=255)                                                cgcgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACaa                                                                                                                  >  2:334110/1‑91 (MQ=255)                                                cgcgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACaa                                                                                                                  <  1:334110/91‑1 (MQ=255)                                                cgcgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCa                                                                                                <  2:374940/109‑1 (MQ=255)                                                cgcgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCa                                                                                                >  1:374940/1‑109 (MQ=255)                                                cgcgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                             <  1:100319/112‑1 (MQ=255)                                                cgcgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                             >  2:100319/1‑112 (MQ=255)                                                   gggCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGa                                                               <  2:1414312/139‑1 (MQ=255)                                                    ggCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAg                                                                                                                >  1:626880/1‑89 (MQ=255)                                                    ggCTTAACCGGGGCGCAATTACTGCTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAg                                                                                                                <  2:626880/89‑1 (MQ=255)                                                                 cgcAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCAtg                                                 >  2:1509238/1‑139 (MQ=255)                                                                                cTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCa                                                         >  1:1338456/1‑116 (MQ=255)                                                                                cTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCa                                                         <  2:1338456/116‑1 (MQ=255)                                                                                   cAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAAcc                               >  1:1182144/1‑139 (MQ=255)                                                                                   cAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAAcc                               >  1:1386815/1‑139 (MQ=255)                                                                                         aaTCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGaaa                         >  2:100629/1‑139 (MQ=255)                                                                                         aaTCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGaaa                         >  2:451388/1‑139 (MQ=255)                                                                                           tcactcacCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCttt                                                                                           >  2:1241645/1‑71 (MQ=255)                                                                                           tcactcacCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCttt                                                                                           <  1:1241645/71‑1 (MQ=255)                                                                                               tcacCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCg                   <  2:1360149/139‑1 (MQ=255)                                                                                               tcacCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCg                   <  2:1360150/139‑1 (MQ=255)                                                                                                 acCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCttt                                                                                           >  1:1223360/1‑65 (MQ=255)                                                                                                 acCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCttt                                                                                           <  2:1223360/65‑1 (MQ=255)                                                                                                        cgcgGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACg                                                                >  2:136159/1‑85 (MQ=255)                                                                                                        cgcgGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACg                                                                <  1:136159/85‑1 (MQ=255)                                                                                                           ggTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCa       <  2:44363/139‑1 (MQ=255)                                                                                                           ggTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCa       <  1:451388/139‑1 (MQ=255)                                                                                                               aTCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCgt   <  1:100629/139‑1 (MQ=255)                                                                                                                tCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCgtg  >  1:870659/1‑139 (MQ=255)                                                                                                                tCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCgtg  >  1:1440160/1‑139 (MQ=255)                                                                                                                     cTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACtt                                    >  2:1486094/1‑100 (MQ=255)                                                                                                                     cTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACtt                                    <  1:1486094/100‑1 (MQ=255)                                                                                                                                       |                                                                                                                    GGCATTGACGCCATTAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAGCAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTG  >  NC_000913/1380673‑1380922

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