Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,414,037	A→C	T190P (ACC→CCC)	ackA →	acetate kinase A and propionate kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,414,037	0	A	C	97.8%	123.6 / ‑6.6	45	T190P (ACC→CCC)	ackA	acetate kinase A and propionate kinase 2
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (20/24);  total (20/25)
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

TGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATCCGCAACGGTAAATGCGTTGACACCTCTATGGGCCTGA  >  NC_000913/2413907‑2414166                                                                                                                                   |                                                                                                                                  cgCTGTTTTTGACACCGCGTTCAACCATACGATGCCGGAAGAGGCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCgg                                                                                                                       <  2:383477/142‑1 (MQ=255)    tGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTCCCTTACAACCTGTACAAAGAGCACGGCATCCATCGTTACGGCGCGCACGGCACCAGCCCCCTCTATGTACCCCAGGCAGCCGCac                                                                                                                    >  2:179228/1‑142 (MQ=255)         ttGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATg                                                                                                               <  1:635498/143‑1 (MQ=255)           gACATCGCGTTCCACCATACTATGCCGTAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCt                                                                                                             <  2:490419/143‑1 (MQ=255)            acacCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTg                                                                                                            <  2:450476/143‑1 (MQ=255)              acCGCGTTCCACCAGACTATGCCGGAAGAATCTTACATCTACGCCCTGCCTTACAACCTGTACAAAGAGCACAGCATCCGTCGTTACAGCGCGCACGTCACCAGACACTTCTATGTACCCCAGAAAGCGGCAAAAATGCTaaa                                                                                                          >  2:288608/1‑143 (MQ=255)                  tgTTCCATCAGTCTATGTCTTAATAGTTTTTTCTCTACTTCCTGTTTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACaaa                                                                                                      <  1:543782/142‑1 (MQ=255)                  cgTCCCACCAGACTATGCCGGAAGCGTCTTACCTCTACGCCCCGCCTGACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCATGAAGCGGCAAAAATGCTGAACaaa                                                                                                      <  1:690345/143‑1 (MQ=255)                      ccatcaGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCgg                                                                                                  <  1:72846/143‑1 (MQ=255)                      ccaccaTACTATGCCGCAACACTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCgg                                                                                                  <  1:712874/143‑1 (MQ=255)                      ccaccaGACTAGGCCGTAGGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCgg                                                                                                  <  2:199747/143‑1 (MQ=255)                        accaGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGAAATAAGTAATAACGGCGCGCACGGCACCAGCCACTTCTATGTACACCAGGAAGCGGCAAAAATGCTGAACAGACCGGTa                                                                                                >  1:111376/1‑143 (MQ=255)                          caGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCCCCAGCCCCTTCTATGTACCCCAGGACACGGCAAAAATGCTGAAAAAACCGGAaga                                                                                              >  1:660557/1‑143 (MQ=255)                          caGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACCTCTATGTACCCCCGGACGCGGCGACAGTGCTGAACAAACCGGTaaa                                                                                              >  1:289650/1‑141 (MQ=255)                          caGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGAGCGCACGGCACCAGCAACTACTATGTACCCCAGGAAGCGGCAAAAATGCAGAACAAAACGGTaga                                                                                              >  1:367458/1‑143 (MQ=255)                               tgtGCCGTAAGAGTCGTACCTCTGCGGCCTGCCTTACAACCTGTGCAAAGAGCACGGCATCCGTCGTTACGGCGCCCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAGATGCTGAACAAACCGGTAGAAGAAc                                                                                         <  2:555262/141‑1 (MQ=255)                               tATGCCGGAAGAGTCTTACCACTACGCCCTGCCTTACAACCTGTAAAAAGAGCACGGCATCCGTCGTTACGGAGCGCACGGCACCAGCCACTTCTATGTACCCCAGGACGCGACAAAGATGCGGAAAAAACCGCACGAAGAAc                                                                                         >  2:303679/1‑143 (MQ=255)                                 tGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTg                                                                                       <  1:295116/143‑1 (MQ=255)                                  gtcagAAGAGTCTTACCTCTACGCCCTGCCGTACAACCTGGACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGa                                                                                      <  1:432731/139‑1 (MQ=255)                                      gAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGCATTTAGAACAAACTAATAGAAGGGCACCGGCCACATCTATGTACCCCAGGAAGCGGCAAAAAGGCTGAACAAACACGTAGAAAAACTGAAcac                                                                                  >  2:516188/1‑142 (MQ=255)                                           gTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGCCCCCAGCCACTTCTATGTACCCCAGGAAACAACAAAAATACAGAACAAAAAGATAGAAGAACTGAACAtcatca                                                                             >  1:87364/1‑143 (MQ=255)                                             cTTACCTCTACGCCCTGCCTTACTACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACCGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCAcc                                                                           <  2:526579/143‑1 (MQ=255)                                              ttACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCCTTCCGGCGCGCACGCCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGATGAACAAAACGGTAGAAGAACCGAACAACATCAccc                                                                          >  2:83173/1‑142 (MQ=255)                                                       cGCCCTGCCTTACTACCTGTCCAAATAGCACAGTATTCTTTTTTACGTCGCGCACGGCACCAGCCACTTATATGGACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTgg                                                                 <  1:602378/143‑1 (MQ=255)                                                       cGCCCTGCCGTACAACCTGTAGAAGGAGCAAGGCAGCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTgg                                                                 <  1:563182/143‑1 (MQ=255)                                                           cTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCGGAACAAACCGGTAGAAGAACTGAACATCATCACCTGACACCCgggagc                                                             >  2:358815/1‑140 (MQ=255)                                                           cTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTAAGTACCCCAGGAAGCGGCACAAATGCTGAACAACCAGGTAGCAGAGCTGAACACCATCACCAGCCAccacgtcaa                                                             >  2:358819/1‑136 (MQ=255)                                                            tGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACCTCTACGTACCCCAGGAAGCGGCACAAATGCTGAACAAACCGATAGAAGAACTGAACATCATAACCTGCCACCAGAGCAAc                                                            >  1:581883/1‑143 (MQ=255)                                                               cTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAACCATGCTGAACAAACCGGTAGAAGCACTGAACATCATCACCTGCCACCTGGGCAACggc                                                         >  2:34022/1‑142 (MQ=255)                                                                ttACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCAGCAGCCACTTCTATGTACCCCAGGAAGCGGCAAACATGCTGAACAACCCGGTCGAACAACTGCACATCATCCACTGCCACCTGTGCACTggtg                                                        >  1:589538/1‑143 (MQ=255)                                                                  aCAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACggtggt                                                      <  1:156586/143‑1 (MQ=255)                                                                              aaGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCCCCTGGGCAACGGCGGCTCCGTTTCTGCt                                          >  2:454611/1‑143 (MQ=255)                                                                              aaGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGc                                                               <  1:701649/122‑1 (MQ=255)                                                                              aaGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGc                                                               >  2:701649/1‑122 (MQ=255)                                                                                 ggCACGTCATCCGTCGTTACCGCGCGCACGGCCCCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAGCCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATc                                       <  2:60189/142‑1 (MQ=255)                                                                                 agCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGCGCAACGGCCGTTCCGTTTCTggcata                                       >  2:64505/1‑138 (MQ=255)                                                                                  gCACGGCCTCCGTCGTTACGGGGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATcc                                      <  2:111376/143‑1 (MQ=255)                                                                                       gCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTAGGCAACGGAGGATCCGTATCTGCTATCCGCAAc                                 >  1:687496/1‑143 (MQ=255)                                                                                            cgtggtTACGGCGCGTACGGCACCTGCGGCTTCTATGTACCCCAGGAAGCGGCAAAAGTGCTGAACAAGCCGGTAGAAGAACTGAACATCAGCATCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATCCGCAACGGTaa                            <  2:351157/143‑1 (MQ=255)                                                                                             gtagtTACGGCGCGCACGCAACCTGCCACATCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATCCGCAACGGTaaa                           <  2:596332/140‑1 (MQ=255)                                                                                                   aCGGCGCGCACGGCACCAGCCACTTCTATGTACCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGCAGAACTGAACATCCTCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCCGCTATCCGCAACGGTCAATGCGtc                     >  2:124972/1‑142 (MQ=255)                                                                                                                cACCAGCCACTTTTATGTACCCCAGGAAGCGTCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATCCGCAACGGTAAATGCGTTGACACCTCTATgg        <  1:316010/143‑1 (MQ=255)                                                                                                                  ccAGCCACTTCTATGTGCCCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATCCGCAACGGTAAATGCGTTGACACCTCTATGGGc      <  2:660557/143‑1 (MQ=255)                                                                                                                  ccAGCCACGTCTATGTACCCTAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATCCGCAACGGTAAATGCGTTGACACCTCTATGGGc      <  2:367458/143‑1 (MQ=255)                                                                                                                      ccACTTCTATGGTCCCCAGGAAGCGGCAAAAATGCTGAACAAACCTGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATCCGCAACGGTAAATGCGTTGACACCTCTATGGGCCTGa  <  2:366522/143‑1 (MQ=255)                                                                                                                                   |                                                                                                                                  TGCTGTATTTGACACCGCGTTCCACCAGACTATGCCGGAAGAGTCTTACCTCTACGCCCTGCCTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACCGGTAGAAGAACTGAACATCATCACCTGCCACCTGGGCAACGGTGGTTCCGTTTCTGCTATCCGCAACGGTAAATGCGTTGACACCTCTATGGGCCTGA  >  NC_000913/2413907‑2414166

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