Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,096,754	G→A	G25D (GGT→GAT)	rdgB →	dITP/XTP pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,096,754	0	G	A	96.9%	93.8 / NA	32	G25D (GGT→GAT)	rdgB	dITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/11);  minor base C (1/0);  total (21/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGGTTGACGGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAA  >  NC_000913/3096507‑3096982                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                     cgGTTGACGGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGTTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGGGATGTAAAAAGGCAAACGTGGCCGCCACAAACAAATTAAAAGCATTAATCCGCAACAAAGCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCt                                                                                                                                                                                                                                   <  2:227454/251‑1 (MQ=255)        cGGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATc                                                                                                                                                                                                                            >  1:263600/1‑251 (MQ=255)              ggccaACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGcc                                                                                                                                                                                                                      >  1:69453/1‑251 (MQ=255)                        cATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGTTTGCCAAAAGCCAGTTGGTGATTAAAAAAGGCGACCTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGAcc                                                                                                                                                                                                            <  2:263600/251‑1 (MQ=255)                                             aagcctttcttgtttGCCAAAAGCCAGGTGGTGATTGAAAAAGGAGAAGAGGGCCGCCACAAACAAATAAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGAACACAGTTGTCCTCGCAACAGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACCGACCTCGGCGATGATTCCGCTga                                                                                                                                                                                         >  1:227950/14‑249 (MQ=255)                                                 aaTTCAGGGTTGCCAAAAGACAGGTGGTGATTGACAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCACGCCAGAAATCGCGGCGTTAATTAATTAGGTCTCCTATGCAAACAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAAc                                                                                                                                                                                   >  2:38048/1‑251 (MQ=255)                                                        ggTTGCAAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTg                                                                                                                                                                            >  1:51423/1‑251 (MQ=255)                                                                       ggtggtGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTGCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAAc                                                                                                                                                             <  1:26655/251‑1 (MQ=255)                                                                       ggtggtGATTGAAAAAGGCGAACTGGCCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAAc                                                                                                                                                             <  1:31688/251‑1 (MQ=255)                                                                                                                      tCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGc                                                                                                              >  2:74292/1‑251 (MQ=255)                                                                                                                            aTCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGc                                                                                                        >  1:226324/1‑251 (MQ=255)                                                                                                                             tCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATACTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAACGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGcc                                                                                                       >  1:29043/1‑251 (MQ=255)                                                                                                                              ccgaaacaaATCACGCCAGAAATCGAGGCGTTAATTAATTAGGAATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTAACGGCAATTGCcg                                                                                                      >  1:193133/5‑251 (MQ=255)                                                                                                                              ccGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCACAAGTGACCGCTTTACCGGCAATTGCcg                                                                                                      >  1:193141/1‑251 (MQ=255)                                                                                                                                gCAACAAATCCCGCCAGCAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGCCCTTTATCGAAAACGCGCTTCTGCAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCcgcc                                                                                                    >  2:223880/1‑249 (MQ=255)                                                                                                                                              cAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGGGGCCCAAACAGACCGCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGc                                                                                      <  2:168980/251‑1 (MQ=255)                                                                                                                                                                     taGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGgcgcgc                                                               >  1:243615/1‑251 (MQ=255)                                                                                                                                                                      aGGTATCCTATTCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGcc                                                              <  2:24154/251‑1 (MQ=255)                                                                                                                                                                                  cAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTTGCGGCGCGCCGGGCCTTTACTc                                                  >  1:110857/1‑251 (MQ=255)                                                                                                                                                                                     aaaGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATCCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcgc                                               >  1:174122/1‑251 (MQ=255)                                                                                                                                                                                                    aCCGGCAATGTAGGTACTGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTCGATGGGCTTGGGGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTg                                 >  2:16294/1‑250 (MQ=255)                                                                                                                                                                                                       ggCAATGTCGGTCAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCACAAACAGACCTCGGCGTTGCTTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTg                                                                                                                            >  2:21080/1‑156 (MQ=255)                                                                                                                                                                                                       ggCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTg                                                                                                                            <  1:21080/156‑1 (MQ=255)                                                                                                                                                                                                             gTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAc                       >  2:53713/1‑251 (MQ=255)                                                                                                                                                                                                                   aaaGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCa                 <  1:53713/251‑1 (MQ=255)                                                                                                                                                                                                                            gAGCTTGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCt        <  1:249703/251‑1 (MQ=255)                                                                                                                                                                                                                               cTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCACTTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGCTCAAAAAAATATGc      >  2:178469/1‑250 (MQ=255)                                                                                                                                                                                                                               cTGGCGTCGCTGCTTAGCGAATTCGATCTTGATATCGTGGCCCAAACCGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGCGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGCTCAAAAGAATCTGc      >  2:171985/1‑250 (MQ=255)                                                                                                                                                                                                                                 ggCGTCGCTGCTTAGCGACTTCGCTCTTGATATCGTGGCCCACACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGCGACCGCTTTACCGGCCATTGCCGCCGCCTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGCTATTCCGGTGCAGACGCGCCCGACCACAAGAACCTGccc    >  2:150287/1‑248 (MQ=255)                                                                                                                                                                                                                                  gCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAATCCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCaaaa  <  1:7770/251‑1 (MQ=255)                                                                                                                                                                                                                                      cGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCa     <  2:198207/244‑1 (MQ=255)                                                                                                                                                                                                                                      cGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCa     >  1:198207/1‑244 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                     CGGTTGACGGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAA  >  NC_000913/3096507‑3096982

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