Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,832,618	G→C	G64A (GGC→GCC)	xthA →	exodeoxyribonuclease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,832,618	0	G	C	96.7%	103.5 / ‑6.7	30	G64A (GGC→GCC)	xthA	exodeoxyribonuclease III
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (17/12);  minor base T (1/0);  total (18/12)
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

TGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGGCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCG  >  NC_000913/1832541‑1832693                                                                              |                                                                            tGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTAtgccgt                                                                          >  1:305495/1‑81 (MQ=255)       tatGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCgct                                                                    <  2:305471/81‑1 (MQ=255)        atGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCgctg                                                                   >  2:305518/1‑81 (MQ=255)        atGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCgctg                                                                   >  1:305494/1‑81 (MQ=255)         tGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCgctg                                                                   <  2:305472/80‑1 (MQ=255)             tCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGAc                                                              >  2:305501/1‑81 (MQ=255)             tCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGAc                                                              >  1:305509/1‑81 (MQ=255)                    gAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCaaa                                                          >  1:305506/1‑78 (MQ=255)                           tGGCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGa                                                >  1:305497/1‑81 (MQ=255)                             gCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGAtt                                              <  1:305473/81‑1 (MQ=255)                              cGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTg                                             <  2:305474/81‑1 (MQ=255)                               gAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGc                                            >  2:305505/1‑81 (MQ=255)                               gAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGc                                            >  1:305504/1‑81 (MQ=255)                                 aGCTCGGCCACAACGCGCCTTATCACCGTCAGAAAGGCCATTACGTCGTGGCGCTGCCGCCCAAAGAGACGCCGATTGct                                           >  2:305490/1‑79 (MQ=255)                                   cTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTg                                        >  1:305499/1‑81 (MQ=255)                                      ggCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGt                                     >  2:305498/1‑81 (MQ=255)                                      ggCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGt                                     >  1:305512/1‑81 (MQ=255)                                         tACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTcgc                                  >  2:305507/1‑81 (MQ=255)                                           cAACGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCgg                                >  1:305519/1‑81 (MQ=255)                                             aCGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCt                              >  1:305510/1‑81 (MQ=255)                                              cGTGTTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCtt                             <  1:305475/81‑1 (MQ=255)                                                tgtTTTATCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTc                           <  2:305476/81‑1 (MQ=255)                                                      aTCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGt                      >  1:305508/1‑80 (MQ=255)                                                      aTCACGGGCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGt                      <  1:305477/80‑1 (MQ=255)                                                           gggCAGAAAGGCCATTATGCCGTGGCGCTGCAGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTgacgac                <  2:305478/81‑1 (MQ=255)                                                             gCAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCTCGGCTTTCCCGGTGacgacg               >  1:305522/1‑80 (MQ=255)                                                              cAGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGaa             <  1:305479/81‑1 (MQ=255)                                                               aGAAAGGCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAag            <  1:305480/81‑1 (MQ=255)                                                                     gCCATTATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGgcgc      <  1:305481/81‑1 (MQ=255)                                                                         ttATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAgcg  <  1:305482/81‑1 (MQ=255)                                                                         ttATGCCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAgcg  >  2:305511/1‑81 (MQ=255)                                                                              |                                                                            TGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATCACGGGCAGAAAGGCCATTATGGCGTGGCGCTGCTGACCAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCG  >  NC_000913/1832541‑1832693

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