Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,414,002	G→A	12.0%	R178H (CGT→CAT)	ackA →	acetate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,414,002	0	G	A	12.0%	63.7 / 2.3	33	R178H (CGT→CAT)	ackA	acetate kinase
Reads supporting (aligned to +/- strand):  ref base G (13/16);  new base A (2/2);  total (15/18)
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.89e-01

AAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATCCGCAACGGTAAA  >  NC_000913/2413864‑2414141                                                                                                                                           |                                                                                                                                            aagaagCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACCAAGAGCACGGCATCCGCCGTTAc                                                                                                                                         >  2:391635/1‑143 (MQ=255) aagaagCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCATTAc                                                                                                                                         >  1:330705/1‑143 (MQ=255) aagaagCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCATTAc                                                                                                                                         >  1:50512/1‑143 (MQ=255)          tGAAATCTTTCCCACAGCTGAAAGAAACAACAGGAGTTGTATTTTACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTTCAAAGAGCACGGCATCCGTCGTTACGgcgcgccc                                                                                                                                >  1:245889/1‑141 (MQ=255)                 tttCCCACCGCTGTAAGACAAAAACGTTGCTGTATTTGACTCCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACTAAGCGCACGGCATCCGTCGTTACGGCGCGCACGGCACCa                                                                                                                         <  2:206118/143‑1 (MQ=255)                   tCCCACAGCTGAAAGACAAAAACGTTGCTGTAGTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTCCGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGc                                                                                                                       <  2:325196/143‑1 (MQ=255)                           cTGAAAGACAAAAACGTTGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTa                                                                                                               <  2:101449/143‑1 (MQ=255)                              aaaGACAAAAACGTTGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCCCCAGCCACCTCTATGt                                                                                                            >  2:130357/1‑143 (MQ=255)                                  aCAAAAACGTTGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCCCGCCATCCGCCGTTCCGGCCCGCCCCGCCCCCGCCACTCCTCTGTCAcc                                                                                                        >  2:324073/1‑143 (MQ=255)                                  aCAAAAACGTTGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAAcc                                                                                                        <  1:405369/143‑1 (MQ=255)                                           ttGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCATTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCg                                                                                               <  2:50512/143‑1 (MQ=255)                                           ttGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCATTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCg                                                                                               <  2:330705/143‑1 (MQ=255)                                            tGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCgg                                                                                              <  1:115968/143‑1 (MQ=255)                                                 tATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCaaaa                                                                                         <  1:540148/143‑1 (MQ=255)                                                    tccgcgccGCGTCCACCCAGACGCTGCCGTGAGCGACTCACCTCGCCGCCCGGCCCTACCCCCTACCCAAAGAGCGCGGCATCGGTCGTTACGGCGCGCACGGCACCAGCCACTCCTATGTAACCCAGGAAGCGGCAAAAATg                                                                                      <  2:212899/137‑1 (MQ=255)                                                         tccGCGTTCCATCGGACTATGCCGGAAGAGTCTTGCCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGaa                                                                                 <  2:218980/142‑1 (MQ=255)                                                                tCCACCAGACGATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCg                                                                          <  1:4755/143‑1 (MQ=255)                                                                                     agTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACatcatc                                                     >  1:192859/1‑143 (MQ=255)                                                                                     agTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGAGGCAAAACTGCTGAACAAACCGGTAGAAGAACTGAACatcatc                                                     >  1:448253/1‑143 (MQ=255)                                                                                      gTCTTACCCCTGCTCCTTGCTTCACCACCTGTACGAAGAGCCCGGCATCCGTCGTCACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCTGCAAAAATGCTGAACAAACCGGTAGAACAACTGAtcatcatca                                                    <  2:319981/143‑1 (MQ=255)                                                                                       tCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCAc                                                   <  1:338341/143‑1 (MQ=255)                                                                                       tCTTACCTCGACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCAc                                                   <  1:220137/143‑1 (MQ=255)                                                                                         ttACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGGTACGGCGCGCACGGCACCAACCCCTCCTATGCAACCCAGCAAGCGGCAAAAATGCCGAACAAACCGGTGGAACGACACAACAGCGTCACCt                                                 >  2:424249/1‑143 (MQ=255)                                                                                             ctctACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAGCCGGCAGAAGCACTGAACATCATCACCTgacg                                             >  1:127501/1‑140 (MQ=255)                                                                                                  cGCCCTGCCTTACAACCTGTACAAAGCGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTgg                                        <  1:288505/143‑1 (MQ=255)                                                                                                    cccTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCGGCCACCTGGGc                                      >  1:303790/1‑143 (MQ=255)                                                                                                            tcccaCCTGCACACAGAGCACGGCATCCGTCGTTCCGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACggtgg                              <  2:510122/139‑1 (MQ=255)                                                                                                              cAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGtt                            <  2:151993/143‑1 (MQ=255)                                                                                                                       cAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGACACCAGGGCAACGGTGGTTCCGTTTCTg                   >  2:58886/1‑143 (MQ=255)                                                                                                                        aaaGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGc                  >  1:553355/1‑143 (MQ=255)                                                                                                                        aaaGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCCATGTAACCCAGGAAGCGGCCAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAAAGCAGCTTACGTCTCTGc                  >  2:196907/1‑143 (MQ=255)                                                                                                                             gtaCGGCATCCGTCGTTTGGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATcc             <  2:215198/141‑1 (MQ=255)                                                                                                                              cACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAAGGGGGGTTCCGGTTCTGCTATCCg            >  1:232377/1‑143 (MQ=255)                                                                                                                                        gtcgtTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGATGAACAAGCCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTGTCCGCAACGGtgac  >  1:518025/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                            AAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATCCGCAACGGTAAA  >  NC_000913/2413864‑2414141

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