Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,6750TG51.6% 33.9 / 42.8 31T218P (ACC→CCC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base T (12/3);  new base G (16/0);  total (28/3)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-05

ACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGT  >  NC_000913/4346528‑4346821
                                                                                                                                                   |                                                                                                                                                  
aCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGGACCGAGGGGACg                                                                                                                                                 >  1:236269/1‑151 (MQ=255)
aCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACg                                                                                                                                                 >  2:151516/1‑151 (MQ=255)
   ccTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGAGGGGACGCaa                                                                                                                                              >  2:493404/1‑150 (MQ=255)
                                 gTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGTTTTTTTGTTTGc                                                                                                                 >  2:354427/1‑150 (MQ=255)
                                      tCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGGTGCAGGACCGCGGGGAAGCCTTTTTTCGCCGAGGAGGGTTTTCAGTTTTTCggggg                                                                                                           >  1:78185/1‑149 (MQ=255)
                                      tCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGCTGTAGGACCGAGGGGACCCATTTTCTCGTCGAGGGAGGTTTTCAGGTTGCCagggg                                                                                                           >  1:578854/1‑149 (MQ=255)
                                      tCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTCGTTTGCCggggg                                                                                                           >  2:592928/1‑149 (MQ=255)
                                      tCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTTAGTTTGCCggggg                                                                                                           >  2:52230/1‑149 (MQ=255)
                                      tCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTCAGTTTGCCGGGAg                                                                                                           >  2:298464/1‑151 (MQ=255)
                                      tCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTAGTTTGCCGGGAg                                                                                                           >  2:426008/1‑151 (MQ=255)
                                        gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGCATTTTTTTGACCAGGAAGGTTTTCAGTTTTTCgggggg                                                                                                          >  1:189740/1‑147 (MQ=255)
                                        gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTCAGTTTGCCGGGGGt                                                                                                          >  2:673857/1‑150 (MQ=255)
                                             aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCATGCTGCAATACCGGAGGTACGCGTTTTCTCGACGAGGAAAATTTTTCGTTTTCCCGGAGGcagcag                                                                                                    >  2:251339/1‑151 (MQ=255)
                                                                            gttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGAGTCAGCAGGGCTTTGGGTTCCTGGGAGAGATACGTTttg                                                                     >  2:261439/1‑151 (MQ=255)
                                                                                       cAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTTCCGGGGGTCCGCAGGGGTTTGGTTTCCTGGTAGAGATAAGGTTTGTTGGCAGAgc                                                           >  2:511500/1‑150 (MQ=255)
                                                                                                 gCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTCAGTTTTCCGGGGGGCAGCAGGGCTTTGGTTTCCTGGGAGAGATACGTTTTGTTGGCAGAGCCGCCGCCttt                                                 >  2:625936/1‑150 (MQ=255)
                                                                                                                 atatGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGGAGAGATACGGTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAtt                                >  2:311331/1‑151 (MQ=255)
                                                                                                                 atatGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGTCAGCAGGGCTTTGGTTTTCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAtt                                >  1:322715/1‑151 (MQ=255)
                                                                                                                                  ggCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                 >  2:219013/1‑149 (MQ=255)
                                                                                                                                   gCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTa              >  1:551429/1‑151 (MQ=255)
                                                                                                                                    cTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTa              >  1:694076/1‑150 (MQ=255)
                                                                                                                                    cTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTa              >  2:535477/1‑150 (MQ=255)
                                                                                                                                    cTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAc             <  2:390703/151‑1 (MQ=255)
                                                                                                                                    cTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAc             <  1:23503/151‑1 (MQ=255)
                                                                                                                                      gCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACcg           >  2:98351/1‑151 (MQ=255)
                                                                                                                                           aCCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTa       >  2:94047/1‑150 (MQ=255)
                                                                                                                                           aCCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTa       >  1:529781/1‑150 (MQ=255)
                                                                                                                                           aCCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTa       >  2:121569/1‑150 (MQ=255)
                                                                                                                                           aCCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTAc      <  1:298463/151‑1 (MQ=255)
                                                                                                                                                gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGt  >  2:529833/1‑150 (MQ=255)
                                                                                                                                                gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGt  >  2:7351/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGT  >  NC_000913/4346528‑4346821

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