Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,256	C→T	G218G (GGC→GGT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,256	0	C	T	96.8%	81.8 / ‑5.1	32	G218G (GGC→GGT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (15/15);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAA‑GAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACC  >  NC_000913/4620016‑4620490                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           gCCGGTTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCg                                                                                                                                                                                                                                  <  1:143519/251‑1 (MQ=255)           cTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTcc                                                                                                                                                                                                                         >  1:314005/1‑250 (MQ=255)                                     ggTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTgg                                                                                                                                                                                              >  1:92808/1‑251 (MQ=255)                                     ggTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTAGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCCGCGTCCCGGGAACCGTCACGACCTgg                                                                                                                                                                                              >  1:285591/1‑251 (MQ=255)                                      gTCATTCTGGAGCAGCTGTGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGACCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGCAGAGCTGACCAACGGCGGCTACACCATCGGTCGTGTTATCGCTAGGCCGTTTATCGGTGACCAAGCACGTAACTTCCAGCGTAACGGTAACCGTCaagcaccggc                                                                                                                                                                                             >  2:67843/1‑242 (MQ=255)                                                           gAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTAGCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCaccgac                                                                                                                                                                        <  1:333415/251‑1 (MQ=255)                                                                  acaTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTTAGCCGCCAGCCCCGACCGtgctg                                                                                                                                                                 >  1:350041/1‑251 (MQ=255)                                                                         aaCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCCACAACATCGGTCGTGGTATCGCCCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTCACGGTAACCGTCACGACCTGGCTGTTGGGCCACCAGCACCAGCCGGGCTGCaaaacc                                                                                                                                                          >  1:98013/1‑246 (MQ=255)                                                                            cGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGATCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTgg                                                                                                                                                       <  1:350669/251‑1 (MQ=255)                                                                             ggCAAGCCTATTTTCTATCCCTCCGCTGTCTCCGTGTTCCAGATTGCCTGCCACGAATGATCCTTT‑GGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACATCACCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCTGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGt                                                                                                                                                      <  1:256993/251‑1 (MQ=255)                                                                                        tttCTATACCTCCGCTGCCTCCGTGTACCAGATTGCCTGCCACGAC‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCCGAAACTGGTTGATGAAAAac                                                                                                                                           >  1:303974/1‑251 (MQ=255)                                                                                         ttCTATACCTCCGCTGACTCCGTGTTCAAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAac                                                                                                                                           <  1:151379/250‑1 (MQ=255)                                                                                             atcccTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCGCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAacac                                                                                                                                         <  1:360013/245‑1 (MQ=255)                                                                                                     gCTGACTCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtggt                                                                                                                              <  1:158973/251‑1 (MQ=255)                                                                                                           tCCGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGAGGCTACAAAAACGGTCGTGTTATCGCTCGTCAGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAAACGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGGGCTGCAGAAACTGGTTGATGTAAAACACGGTCAGGGGGTTTATGt                                                                                                                        >  1:22104/1‑251 (MQ=255)                                                                                                            ccGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTc                                                                                                                       >  1:306294/1‑251 (MQ=255)                                                                                                            ccGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTc                                                                                                                       <  1:319598/251‑1 (MQ=255)                                                                                                             cGTGTTCCAGATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCCCCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCg                                                                                                                      >  1:313362/1‑251 (MQ=255)                                                                                                                      gATTGCCTGCCACGAA‑GGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCAGCTACAACATCTGTCGCGTTATCGCACCTCCGTTTATCGGTGACCAAGCCGGTGACTTCCAGCGTACCGGGAACCGTCACGCCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGGAACTGGTTGATGCAAAACACGGTTAGGTGGGTTCTTGCGGTAAAAttt                                                                                                             >  1:132679/1‑250 (MQ=255)                                                                                                                                     aaggaaccttaggccaggATAAACCCTACGAACTGTGCGAAATCGCCCGTGAAGAGCGGACCACCGGCGGCGACAACATCGGTCGGGTTATCGCTCGTCCGTTTAGCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGcc                                                                                               <  2:100362/235‑1 (MQ=255)                                                                                                                                                gtgccggATAAACTCTACGAACTGCTCGAAATCGCTCTTTGAGTGCTGACCCACGGCGGCTACAATTTCGGTCGTGTTATCGCTCGTCCGTTTATGGGCGACAAAGCCGGTAACTCCCATCGTACCGGTAACCGTCATGACCTGGCTGTTGAGCCGCGAGCACCGACCGTGCGGCAGAAGCTGGTTGATGAAAAACACGGGCAGGTGGTTTCTGTCGGGAAAATTGCGGACATCTACGCCAACTGCGGTAt                                                                                    <  2:165880/246‑1 (MQ=255)                                                                                                                                                 gccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTCAAATTGCAGACATCTACGCCAACTGCGGTATc                                                                                   >  1:311224/3‑251 (MQ=255)                                                                                                                                                     ggATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCAcc                                                                                >  1:339207/1‑250 (MQ=255)                                                                                                                                                        tAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGCCCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGGGGTTTCTGTCGGTAAAATTTCCGACATCTACGCCCACTGCGGTATCAccc                                                                               >  1:159809/1‑247 (MQ=255)                                                                                                                                                          aaccctACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaagaaa                                                                          <  1:374797/247‑1 (MQ=255)                                                                                                                                                           aCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAg                                                                         <  1:161927/251‑1 (MQ=255)                                                                                                                                                                aCGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGaaa                                                                    <  1:76168/251‑1 (MQ=255)                                                                                                                                                                                                  cAACGGCGGCTACAACATCTGTCGTGTTATCGCTCGTCCGTTTTTCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGTCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGtcacca                                  <  1:69963/251‑1 (MQ=255)                                                                                                                                                                                                         ggcTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCaa                              >  1:53314/1‑248 (MQ=255)                                                                                                                                                                                                                            ttATCGCTCGTCCGTTTAGCGGTGACAAAGCCGTTAACTTCCAGCGTACCGGTACCCGCCACGCCCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGa         <  2:314005/250‑1 (MQ=255)                                                                                                                                                                                                                                 gCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTAACGACCTGGCTGTTGAGCCGCCAGCACCGACAGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTATGTCGGTAAAATTGCAGAAATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAGGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGGGATAacac   >  1:170132/1‑251 (MQ=255)                                                                                                                                                                                                                                    cGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACAcc  <  1:190136/249‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           GCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAA‑GAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACC  >  NC_000913/4620016‑4620490

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