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	97.9%	142.1 / ‑4.1	47	*220* (TAG→TAA)	ycjU	beta‑phosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (25/21);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

CAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAGCAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTCGGCCCGTCGGGCT  >  NC_000913/1380693‑1380941                                                                                                                   |                                                                                                                                       cAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                                <  1:238590/139‑1 (MQ=255)    cGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTc                                                                                                             >  2:327166/1‑139 (MQ=255)    cGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTc                                                                                                             >  1:888608/1‑139 (MQ=255)         tGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTac                                                                                                        >  2:62125/1‑139 (MQ=255)         tGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTac                                                                                                        <  1:62125/139‑1 (MQ=255)             cTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAAc                                                                                                                           <  2:382737/116‑1 (MQ=255)             cTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAAc                                                                                                                           >  1:382737/1‑116 (MQ=255)                 gTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATAtt                                                                                                >  2:1043426/1‑139 (MQ=255)                     ggATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                                <  1:1027713/119‑1 (MQ=255)                     ggATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                                >  2:1027713/1‑119 (MQ=255)                         cGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAAt                                                                                        <  1:814761/139‑1 (MQ=255)                            cgcgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                                <  2:1299704/112‑1 (MQ=255)                            cgcgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCt                                                                                                                >  1:1299704/1‑112 (MQ=255)                            cgcgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTa                                                                                     >  1:982626/1‑139 (MQ=255)                            cgcgGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTa                                                                                     >  2:322102/1‑139 (MQ=255)                                      aCCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACa                                                                                                                                      >  1:92038/1‑80 (MQ=255)                                          gggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGc                                                                       <  2:888608/139‑1 (MQ=255)                                           ggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACaaa                                                                                                                                    >  1:108291/1‑77 (MQ=255)                                           ggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACaaa                                                                                                                                    <  2:108291/77‑1 (MQ=255)                                           ggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTaca                                                                                                       >  1:549097/1‑106 (MQ=255)                                           ggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTaca                                                                                                       <  2:549097/106‑1 (MQ=255)                                           ggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATAGCTCAGCTTTCGTTacaac                                                                                                     >  2:310436/1‑108 (MQ=255)                                           ggCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATAGCTCAGCTTTCGTTacaac                                                                                                     <  1:310436/108‑1 (MQ=255)                                             cgcAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGtt                                                                                                          <  2:671144/101‑1 (MQ=255)                                             cgcAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGtt                                                                                                          >  1:671144/1‑101 (MQ=255)                                             cgcAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCAtg                                                                    <  1:327166/139‑1 (MQ=255)                                             cgcAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCAtg                                                                    <  1:1265873/139‑1 (MQ=255)                                                   tACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGa                                                              <  1:322102/139‑1 (MQ=255)                                                         tGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATCACAAAGGAATCAAc                                                                                                                           <  1:869815/72‑1 (MQ=255)                                                         tGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAAc                                                                                                                           >  2:869815/1‑72 (MQ=255)                                                          gCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGGTAACCAGGTGCATGTGGTGAAGGACtt                                                       <  1:1043426/139‑1 (MQ=255)                                                          gCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACtt                                                       >  1:778271/1‑139 (MQ=255)                                                                          ctcacCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATaa                                                                               >  2:687855/1‑99 (MQ=255)                                                                          ctcacCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATaa                                                                               <  1:687855/99‑1 (MQ=255)                                                                               cTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAGCAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATaa                                  >  2:122174/1‑139 (MQ=255)                                                                                       ggTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAg                                                                           >  1:78570/1‑90 (MQ=255)                                                                                       ggTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAg                                                                           <  2:78570/90‑1 (MQ=255)                                                                                       ggTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACtt                                                       <  2:516132/110‑1 (MQ=255)                                                                                       ggTTATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACtt                                                       >  1:516132/1‑110 (MQ=255)                                                                                         ttATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAAt                                           >  2:1240811/1‑120 (MQ=255)                                                                                         ttATCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAAt                                           <  1:1240811/120‑1 (MQ=255)                                                                                            tCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCa                          <  1:666542/134‑1 (MQ=255)                                                                                            tCGGCCTTCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCa                          >  2:666542/1‑134 (MQ=255)                                                                                                  ttCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTc               >  2:607512/1‑139 (MQ=255)                                                                                                   tCTGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTCg              >  1:201146/1‑139 (MQ=255)                                                                                                     tGGCAAAACGTATAACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTCGGc            >  2:1181741/1‑139 (MQ=255)                                                                                                               tataACAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTCGGCCCGTCGGGCt  >  2:724334/1‑139 (MQ=255)                                                                                                                   |                                                                                                                                       CAGCGGTATGCGCTCGGTGGGGATCGGCGCGGGCTTAACCGGGGCGCAATTACTGTTGCCTTCAACGGAATCACTCACCTGGCCGCGGTTATCGGCCTTCTGGCAAAACGTATAGCAAAGGAATCAACATGGCTCAGCTTTCGTTACAACATATTCAAAAAATCTACGATAACCAGGTGCATGTGGTGAAGGACTTCAACCTGGAAATTGCCGATAAAGAGTTCATCGTGTTTGTCGGCCCGTCGGGCT  >  NC_000913/1380693‑1380941

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