Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	457,918	C→T	L165L (CTG→TTG)	clpX →	ATP‑dependent Clp protease ATP‑binding subunit ClpX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	457,918	0	C	T	96.8%	105.2 / NA	31	L165L (CTG→TTG)	clpX	ATP‑dependent Clp protease ATP‑binding subunit ClpX
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (13/17);  minor base A (1/0);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGCTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTT  >  NC_000913/457671‑458164                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       aCAGGCGAAAAAAGGTCTGGCGGTCGCGGTATACAACCATTACAAACGTCTTCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGTCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGt                                                                                                                                                                                                                                                     <  2:160541/251‑1 (MQ=255)     gCGAAAAAAGTGCTGGCGGTCGCGGTTAACAACCATTAAAAACGTCTGCGAAACGGGGGTACCAGCAAAGGCGTCGAGTTGGGCAAAAGAAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAAAGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGGCCGAAGCCGGTTAAGTGGGGGAAGAAGGTGAAAAAATCATTAAGAAGATGTTGc                                                                                                                                                                                                                                                  >  1:110948/1‑250 (MQ=255)                gCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACt                                                                                                                                                                                                                                      <  1:275576/251‑1 (MQ=255)                       gTCGCGGTATACAACCATTACAAACGTCTGCGAAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGAAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGt                                                                                                                                                                                                                               >  1:222424/1‑251 (MQ=255)                             gTATACAACAATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGa                                                                                                                                                                                                                          >  1:258642/1‑250 (MQ=255)                                 aCAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGc                                                                                                                                                                                                                      >  1:323541/1‑250 (MQ=255)                                            aaaCGTCTGCGCAACGGCGATACCAGCAATGGCGGCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCTGAAATGCGTCGACGATGTCCAGAAAGCACAGCGTGGTAt                                                                                                                                                                                                          <  1:330035/251‑1 (MQ=255)                                                       cAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATct                                                                                                                                                                                               <  2:206242/251‑2 (MQ=255)                                                             cGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACAGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGaa                                                                                                                                                                                          >  1:29645/1‑250 (MQ=255)                                                                 aCCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACATCGTGGTATTGTCTACATCGATGAAATCg                                                                                                                                                                                      <  2:29464/250‑1 (MQ=255)                                                                     gCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAg                                                                                                                                                                                 <  1:101396/251‑1 (MQ=255)                                                                      cAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACa                                                                                                                                                                                   >  1:161395/1‑248 (MQ=255)                                                                               cGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGt                                                                                                                                                                        >  1:162749/1‑250 (MQ=255)                                                                               cGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGCTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGt                                                                                                                                                                        >  1:175454/1‑250 (MQ=255)                                                                                                       gctgATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTAc                                                                                                                                                 <  1:147199/249‑1 (MQ=255)                                                                                                                                        gctggctgAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCAc                                                                                                              <  2:84471/251‑1 (MQ=255)                                                                                                                                                   cGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTg                                                                                                   <  1:90419/251‑1 (MQ=255)                                                                                                                                                       ggCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCg                                                                                               >  2:27076/1‑251 (MQ=255)                                                                                                                                                           cgcCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAgg                                                                                           >  1:134780/1‑251 (MQ=255)                                                                                                                                                                         tCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAgctgct                                                                              <  2:238180/250‑1 (MQ=255)                                                                                                                                                                                         gACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTgg                                                             <  1:305270/251‑1 (MQ=255)                                                                                                                                                                                          aCGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTggg                                                            <  2:121044/251‑1 (MQ=255)                                                                                                                                                                                            gcgACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCg                                                          >  1:239142/1‑251 (MQ=255)                                                                                                                                                                                                  acacTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAAc                                                     <  1:375216/250‑1 (MQ=255)                                                                                                                                                                                                              gCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGa                                        <  1:327274/251‑1 (MQ=255)                                                                                                                                                                                                                      tgtgGGTGAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGc                                >  1:47200/1‑251 (MQ=255)                                                                                                                                                                                                                             gAAGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGttt                          <  2:7974/250‑2 (MQ=255)                                                                                                                                                                                                                              aaGACGTTGAAAACATCATTCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGAt                        >  2:203812/1‑251 (MQ=255)                                                                                                                                                                                                                                                 tCAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGttt      <  1:58094/250‑1 (MQ=255)                                                                                                                                                                                                                                                  cAGAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTAt    <  2:2857/251‑1 (MQ=255)                                                                                                                                                                                                                                                    gAAGTTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTAttt  >  1:81343/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       ACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAACGCTGCTGGCTGAAACGCTGGCGCGCCTGCTGGATGTTCCGTTCACCATGGCCGACGCGACTACACTGACCGAAGCCGGTTATGTGGGTGAAGACGTTGAAAACATCATTCAGAAGCTGTTGCAGAAATGCGACTACGATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTT  >  NC_000913/457671‑458164

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