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

GGGCATTGACGCCATTAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAGCAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTCGGCCCGTCGGG  >  NC_000913/1380672‑1380939                                                                                                                                        |                                                                                                                                     gggCATTGACGCCATTAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACaa                                                                                                                                   <  1:1203170/139‑1 (MQ=255)  ggCATTGACGCCATTAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACaaa                                                                                                                                  >  2:1412104/1‑139 (MQ=255)                tAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGc                                                                                                                    <  2:1103642/139‑1 (MQ=255)                      cAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                              <  2:65528/139‑1 (MQ=255)                                     ggTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACatat                                                                                               >  2:1236270/1‑139 (MQ=255)                                       tGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTc                                                                                             >  1:1098108/1‑139 (MQ=255)                                               ggCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATc                                                                                     >  1:1103290/1‑139 (MQ=255)                                                   cgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACg                                                                                 <  2:1098108/139‑1 (MQ=255)                                                           aCCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAg                                                                         <  2:164613/139‑1 (MQ=255)                                                              ggggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTg                                                                      <  2:1103290/139‑1 (MQ=255)                                                                ggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGa                                                                                                                               >  1:55879/1‑80 (MQ=255)                                                                ggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGa                                                                                                                               <  2:55879/80‑1 (MQ=255)                                                                ggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTa                                                                                   >  1:85192/1‑124 (MQ=255)                                                                ggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCa                                                                    >  2:538472/1‑139 (MQ=255)                                                                ggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCACAGCTTTCGTTACAACATATTCAAAAAATCTa                                                                                   <  2:85192/124‑1 (MQ=255)                                                                       ttACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACa                                                                                                                                    >  2:520129/1‑68 (MQ=255)                                                                       ttACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATgtggtg                                                             <  1:1236270/139‑1 (MQ=255)                                                                        tACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGa                                                            <  1:538472/139‑1 (MQ=255)                                                                         aCTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGaa                                                           >  2:1064987/1‑139 (MQ=255)                                                                            gTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                              <  2:952667/85‑1 (MQ=255)                                                                            gTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                              >  1:952667/1‑85 (MQ=255)                                                                                          aaTCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGa                                                            <  2:144499/121‑1 (MQ=255)                                                                                          aaTCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGa                                                            >  1:144499/1‑121 (MQ=255)                                                                                               ctcacCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAACTTTCGTTACAACATATTCAAAAAATCTACGATCACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGcc                                     >  1:965641/1‑139 (MQ=255)                                                                                                 cacCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGa                                                                                                                               <  2:1208202/47‑1 (MQ=255)                                                                                                 cacCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGa                                                                                                                               >  1:1208202/1‑47 (MQ=255)                                                                                                    cTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATaa                                                                             >  2:42232/1‑94 (MQ=255)                                                                                                    cTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATaa                                                                             <  1:42232/94‑1 (MQ=255)                                                                                                                aTCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                              >  1:603235/1‑49 (MQ=255)                                                                                                                aTCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                              <  2:603235/49‑1 (MQ=255)                                                                                                                 tCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGc                                                                     <  2:549852/89‑1 (MQ=255)                                                                                                                 tCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGc                                                                     >  1:549852/1‑89 (MQ=255)                                                                                                                        tCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATCACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAAtt                          <  1:934177/125‑1 (MQ=255)                                                                                                                        tCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGtt                          >  2:934177/1‑125 (MQ=255)                                                                                                                                  cGTATAACAAAGGAATCAACATGGCTCATCTTTCGTTACAACATATTCAAAAAAGCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTCGGCCCGTCggg  <  2:965641/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     GGGCATTGACGCCATTAACGCCAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAGCAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTCGGCCCGTCGGG  >  NC_000913/1380672‑1380939

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