Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,413,791	T→C	S108P (TCT→CCT)	ackA →	acetate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,413,791	0	T	C	97.1%	112.3 / NA	34	S108P (TCT→CCT)	ackA	acetate kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/14);  minor base G (1/0);  total (20/14)
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

GCTCAACTTTATCGTTAATACTATTCTGGCACAAAAACCAGAACTGTCTGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGTCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTGCTGTAT  >  NC_000913/2413667‑2413914                                                                                                                             |                                                                                                                            gcTCAACTTTATCGTCAATACTCTTCTGGCTCAAAAACCAGAACTGTCTGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTAtt                                                                                                                     <  1:224988/133‑1 (MQ=255)   tCAACTTTATCGTTAATACTATTCTGGCACAAAAACCAGAACTGTCTGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCa                                                                                                                   <  1:443635/133‑1 (MQ=255)     aaCTTTATCGTTAATACTATTCTGGCACAAAAACCAGAACTGTCTGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAgg                                                                                                                 >  1:74554/1‑133 (MQ=255)                 aaTACTATTCTGGCACAAAAACCAGAACTGTCTGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTGCGTAGTGAGCGAGGAGGCTGTTATTGAGGGGATCAAAGATg                                                                                                     >  2:166047/1‑133 (MQ=255)                      tATTCTGGCACAAAAACCAGAACTGTCTGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGct                                                                                                >  2:329661/1‑133 (MQ=255)                        ttCTGGCACAAAAACCAGAACTGTCTGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTAAGGATACCAAAGAAGCAGcttc                                                                                              >  1:284647/1‑133 (MQ=255)                             gCACAAAAACCAGAACTGTCTGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTg                                                                                         >  2:437757/1‑133 (MQ=255)                             gCACAAAAACCAGAACTGTCTGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTg                                                                                         >  1:22612/1‑133 (MQ=255)                             gCACAAAAACCAGAACTGTCTGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTg                                                                                         >  1:506936/1‑133 (MQ=255)                                caAAAACCAGCACTGTCTGCGCAGCTGACTGCGATCGGTCATCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCAc                                                                                      <  2:286542/133‑1 (MQ=255)                                                 tGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCt                                                                     <  1:387195/133‑1 (MQ=255)                                                 tGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATACAGATCCGTTTGCACCACTGCAAAACCCGGCt                                                                     >  2:275936/1‑133 (MQ=255)                                                    gcAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCAc                                                                  <  1:437757/133‑1 (MQ=255)                                                         tGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGAt                                                             <  2:74554/133‑1 (MQ=255)                                                         tGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTCTGCACCGTTGCTCATCCCGGCTCACCTGAt                                                             >  1:213518/1‑133 (MQ=255)                                                                      tCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCgaagaa                                                <  2:22612/133‑1 (MQ=255)                                                                                      ggcggcGAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCAGGCTCACCTGATAGGTCTCGAAGAAGCTCTGAAATCTTTcc                                >  2:279492/1‑133 (MQ=255)                                                                                        cggcgAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCca                              >  1:470774/1‑133 (MQ=255)                                                                                        cggcgAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCca                              >  2:438452/1‑133 (MQ=255)                                                                                          gcgAAAAGTATACCAGCTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCcaca                            >  2:567495/1‑133 (MQ=255)                                                                                                          cTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAAc            >  1:588316/1‑133 (MQ=255)                                                                                                          cTCCGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAAc            >  1:114623/1‑133 (MQ=255)                                                                                                             cGTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGtt         >  1:44795/1‑133 (MQ=255)                                                                                                              gTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTg        <  1:283853/133‑1 (MQ=255)                                                                                                              gTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTg        <  1:275936/133‑1 (MQ=255)                                                                                                              gTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTg        <  2:420600/133‑1 (MQ=255)                                                                                                              gTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTg        <  2:470774/133‑1 (MQ=255)                                                                                                              gTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTg        <  2:506936/133‑1 (MQ=255)                                                                                                              gTAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTg        <  2:541295/133‑1 (MQ=255)                                                                                                               tAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTGc       >  2:211396/1‑133 (MQ=255)                                                                                                               tAGTGATCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATATTTCCCACAGCTGACAGACAAAAACGTTGc       >  1:289714/1‑133 (MQ=255)                                                                                                                    aTCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTGCTGTa   <  2:432849/132‑1 (MQ=255)                                                                                                                    aTCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTGCTGTa   >  1:432849/1‑132 (MQ=255)                                                                                                                    aTCGATGAGCCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTGCTGTAt  >  1:62323/1‑133 (MQ=255)                                                                                                                             |                                                                                                                            GCTCAACTTTATCGTTAATACTATTCTGGCACAAAAACCAGAACTGTCTGCGCAGCTGACTGCTATCGGTCACCGTATCGTACACGGCGGCGAAAAGTATACCAGCTCCGTAGTGATCGATGAGTCTGTTATTCAGGGTATCAAAGATGCAGCTTCTTTTGCACCGCTGCACAACCCGGCTCACCTGATCGGTATCGAAGAAGCTCTGAAATCTTTCCCACAGCTGAAAGACAAAAACGTTGCTGTAT  >  NC_000913/2413667‑2413914

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