Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,637	C→T	T345T (ACC→ACT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,637	0	C	T	97.1%	93.4 / NA	35	T345T (ACC→ACT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (11/23);  total (11/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAG  >  NC_000913/4620393‑4620873                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             aCGCCAACTGCGGTATCACCAATAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCt                                                                                                                                                                                                                                           <  1:379881/248‑1 (MQ=255)    ccAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCAc                                                                                                                                                                                                                                     <  1:11820/251‑1 (MQ=255)    ccAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCAc                                                                                                                                                                                                                                     >  1:83880/1‑251 (MQ=255)       aCTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGGGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTTCTGGCCACGGt                                                                                                                                                                                                                                  >  2:109560/1‑251 (MQ=255)          gCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTgt                                                                                                                                                                                                                               <  1:142992/251‑2 (MQ=255)          gCGGTATCACCAAGAAAGTGAAAGCAACGGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTgt                                                                                                                                                                                                                               <  1:13019/251‑2 (MQ=255)            ggTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGCGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGCTGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTgtg                                                                                                                                                                                                                              >  1:22146/1‑248 (MQ=255)              tATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATc                                                                                                                                                                                                                           >  1:322595/1‑251 (MQ=255)                              aaaGCCACTTGCCTGGCAGCGCGGTTTGACGCCAACATCAAAGAGATGAAAGAAGCGGGTGATCACACCACCGTCTTCCCCAACTTCGATGATTTCGACTCTTCCTGGTGCCCCCGTCGCGACGTCTCCGGTTATGCTGCGGGTTTGGAGCTGTTCGACCTCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGCCGACATCCTGATCCTCACTGCTGAACACGGTTGTGATCCGACCTGGACCGGTAc                                                                                                                                                                                                           <  2:313362/251‑1 (MQ=255)                                   aaCTGGTCTGGACGCGCTGTGTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGTTAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGTTACTGAcc                                                                                                                                                                                                      <  2:53314/250‑1 (MQ=255)                                       ggCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACcaca                                                                                                                                                                                                   <  1:342911/250‑1 (MQ=255)                                                           cGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTAc                                                                                                                                                                              >  1:365885/1‑251 (MQ=255)                                                                 caTCAAAGAGATGAAAGAATCGGTTGATAACCCCAGCGTCGACTCCAACTTCGTTGACTTCGAATCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGGCCTCCGTCTGCCGGATCTGATGGCTCTGCTTCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACAGTCCGGTACTGGttt                                                                                                                                                                        <  2:322595/251‑3 (MQ=255)                                                                  aTCAAAGAGAAGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCAGCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGTCATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGttta                                                                                                                                                                       <  1:47358/251‑4 (MQ=255)                                                                   acaataagaTGAAAGAATGGGTTGATAACGCCATCGTCTTCACCAACTTCGTTGTCTTCGACTCTTCCTGGGGCCCCCGTCGCGACGTCGCCGGTTATGCTTCGGGTCTGGAGCTGTTCGGCCGCCGTCTGCCGGCGCTGATGTCTCTGCTGCGCGATGACGACATCCTTATCCTCACTGCTGACCACGGTTGTGATCCGTCCTGGACCGGTACTGACCACACTCGTGAACACATTCCGGTAcggggttac                                                                                                                                                                      <  2:159809/245‑9 (MQ=255)                                                                        tcgataaaaGAAACGGGAGATAACACCATCGTCTTCCCCAACTTCGTTGACTTCGACTCTTCCTGGGACCACCGTCGCGACGTCGCCGGTTATGCTGCTGGACTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGGAGACGACATCCTGATCCTCAGTGAGGACCACGGGTGTGATCCTGCCTGGACCGTGACTGACCACACGCGTGAACACATTCCGGTACTgggttacggacc                                                                                                                                                                 >  2:71523/7‑241 (MQ=255)                                                                        gagaTGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTGCGTTGACTTCGACTCTTCCTGGGGCCACCGGCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGccc                                                                                                                                                                 <  1:57351/251‑1 (MQ=255)                                                                                      gggTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCgg                                                                                                                                                   >  1:137656/1‑251 (MQ=255)                                                                                      gggTGATAACACCAGCGTCTTCACCACCTTCGTTGACTTCGACCCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCTGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGATCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCgg                                                                                                                                                   <  1:71523/251‑1 (MQ=255)                                                                                           aTAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCa                                                                                                                                              <  1:38140/251‑1 (MQ=255)                                                                                               ccccATCGTCTTCCCCAACGTCGTTGATTTCGACTCCTCCTGGGGCCACCGGCGCGACGTCGCCGGTTATTCTGCGGGTCTGGTGCGTTTCGACCGCCGTCTGCCGGAGCTGATGGCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTgg                                                                                                                                          <  2:291890/249‑1 (MQ=255)                                                                                                      gTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACAGGGCTCACTGGGTCaccg                                                                                                                                   >  1:216077/1‑248 (MQ=255)                                                                                                          ccACCAACTGCTTTGATTTCGACCATGCCTGGGGCCCCCGTCGCGACGTCGCCGGTTAGTCTGCGGGTCTTGAGCTGTTCGACCTCCGTCTGCCGGAGCGGAGGTCTCTGCTGCGCGATGACGTCATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACAGTCCGGTACTGGTTTGCGGCCCGAAAGTAAAACCGGGCGCACTTGGTCACCGTg                                                                                                                                 <  2:350041/248‑1 (MQ=255)                                                                                                               cacttcggtgtgttCGGCTCTTCCCGGGGCCACCTACGGGACGTCGCAGGTTATGCTGCTGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAGCCTGATCCTCACTGCTGTACACGTTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTTACCGTGAAACCtt                                                                                                                          <  2:55212/239‑1 (MQ=255)                                                                                                                aCTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTc                                                                                                                         <  1:46653/251‑1 (MQ=255)                                                                                                                   tcggtgtcttggACTCTGCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCGGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCGGCTGCGCGATGACGACATCCTGATCCTCACTGCGGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTc                                                                                                                         <  2:83880/237‑1 (MQ=255)                                                                                                                    cGTTGACTTCGACTCGTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGATCGACCGCCGTCTGCCGGAGCTGCTGGCTCTGCTGCGAGATGACGACAGCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACAACCCGCGTGCACACTTTCAGGTACTGGCTGACAGCCCGAAAGGAAAACAGGGCTCACTGGGGCAACATGAAACCTACGCgg                                                                                                                     >  2:194992/1‑251 (MQ=255)                                                                                                                       tGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGATACCTTCGCggca                                                                                                                   >  1:35403/1‑248 (MQ=255)                                                                                                                                   tgccTGTGGCCACCTTCGCGGCGTGATCGGTTATGCTGCGGGGCTGTTGCTGGTCGACCGCAGTCTGCCGGAGTGGATGTCTCTGCTGCGCGATGACGACATCCTGCTCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGAct                                                                                                       <  2:306294/248‑1 (MQ=255)                                                                                                                                                        cggtggctgTTATGCGTCGGGTCTGGTGCTGTTCTACCGCCGGCTGCCGGAGCTGATGTCTCTGTTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACtt                                                                                 <  2:311224/243‑1 (MQ=255)                                                                                                                                                                                ggatcTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCCCCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAAGATTTTGGTACTTCTGATATGGGATATGGCAAAGc                                                           <  1:87469/245‑1 (MQ=255)                                                                                                                                                                                      tttCGACCGCCGTCGGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACACCCTGATCCTCACTGCTGACCACTGTTGTGATCCGACCTGGTCCGGTACTGACCTCACTCGTGAACTCATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTc                                                    <  1:384064/249‑1 (MQ=255)                                                                                                                                                                                        tcgcacgccgTCTGCCGTATCTGATGTCTCCGCAGCGCGATGACGTTCTCCTGATCCTCACTGCTGTCCACGTTTGTGTTACGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCTGTTCTGGTTTACGGCCCGAATGTAAAACCGGGCTCACTTGGTCACCTTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAACTATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTg                                                  <  2:12333/245‑1 (MQ=255)                                                                                                                                                                                                                               aTGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccgggcccctcc          >  1:265939‑M1/1‑213 (MQ=255)                                                                                                                                                                                                                                       aTCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacaggga  <  1:41907‑M1/251‑47 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             ACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAG  >  NC_000913/4620393‑4620873

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