Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,609,458	T→C	R324R (CGT→CGC)	hyfG →	hydrogenase 4, subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,609,458	0	T	C	94.3%	114.8 / ‑0.6	35	R324R (CGT→CGC)	hyfG	hydrogenase 4, subunit
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (12/21);  total (13/22)
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.83e-01

TGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGA  >  NC_000913/2609316‑2609601                                                                                                                                               |                                                                                                                                                tgctgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCt                                                                                                                                          >  1:158821/1‑150 (MQ=255)    tgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACgg                                                                                                                                      >  1:567535/1‑151 (MQ=255)               cGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGa                                                                                                                           <  1:177928/151‑1 (MQ=255)                    acacAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATtg                                                                                                                        <  1:198310/149‑1 (MQ=255)                     cacaGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATtggt                                                                                                                      <  2:103481/150‑1 (MQ=255)                      acaGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATtggtgg                                                                                                                    >  1:266533/1‑151 (MQ=255)                        aGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGgtggtg                                                                                                                  <  2:439173/151‑1 (MQ=255)                                 cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCcgcg                                                                                                         >  2:190791/1‑151 (MQ=255)                                    aaCCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGata                                                                                                      <  1:273835/151‑1 (MQ=255)                                    aaCCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTAATTGGTGGTGTTCGCCGCGata                                                                                                      <  2:196666/151‑1 (MQ=255)                                          ggCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATtctc                                                                                                 <  2:403775/150‑1 (MQ=255)                                          ggCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATtctc                                                                                                 >  1:360901/1‑150 (MQ=255)                                           gCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCaa                                                                                               <  2:267564/151‑1 (MQ=255)                                              tctcCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGTGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCaa                                                                                               >  1:531399/1‑148 (MQ=255)                                                 ccTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCa                                                                                         <  1:541345/151‑1 (MQ=255)                                                      cATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTc                                                                                    <  1:226289/151‑1 (MQ=255)                                                       aTTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCt                                                                                   >  1:231554/1‑151 (MQ=255)                                                       aTTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCt                                                                                   <  1:192080/151‑1 (MQ=255)                                                         ttCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGc                                                                                 >  2:311925/1‑151 (MQ=255)                                                            gTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCaa                                                                               >  2:490513/1‑150 (MQ=255)                                                              tGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAAc                                                                             <  2:228369/150‑1 (MQ=255)                                                               ggTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACg                                                                            <  2:561343/150‑1 (MQ=255)                                                                gTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCt                                                                          <  1:467796/151‑1 (MQ=255)                                                                           ggCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTgc                                                               <  1:445628/151‑1 (MQ=255)                                                                            gCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTgcg                                                              >  2:624051/1‑151 (MQ=255)                                                                               ttATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGcgcg                                                            <  1:190796/150‑1 (MQ=255)                                                                                            ccGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGAc                                              <  1:213327/151‑1 (MQ=255)                                                                                             cgcgTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACg                                             >  1:352455/1‑151 (MQ=255)                                                                                                      gAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGc                                    >  1:173573/1‑151 (MQ=255)                                                                                                                    gATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGAtgctg                      <  1:210253/151‑1 (MQ=255)                                                                                                                              ttGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGc            >  1:442145/1‑151 (MQ=255)                                                                                                                                gCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGcc           <  2:209004/150‑1 (MQ=255)                                                                                                                                   gATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGa         <  1:311932/149‑1 (MQ=255)                                                                                                                                      cGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATg    <  2:442137/151‑1 (MQ=255)                                                                                                                                        ggTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGa  <  1:443613/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                TGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGA  >  NC_000913/2609316‑2609601

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