Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,7030CT44.1% 36.5 / 13.8 35L208L (CTG→CTAfumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base C (3/16);  new base T (15/0);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 4.48e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTG  >  NC_000913/4346563‑4346850
                                                                                                                                             |                                                                                                                                                   
tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCgg                                                                                                                                           >  2:543681/1‑151 (MQ=255)
  aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGGGGGCGCAGTACCGAGGGGACGCATTTTTTTGACGAGGGAGGTTTTTAAGTTGGCgggg                                                                                                                                          >  1:448285/1‑151 (MQ=255)
  aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGa                                                                                                                                         <  1:525649/151‑1 (MQ=255)
    cGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGGAGTACCGCGGGGACGCATTTTTTCGACGAGGGAGGTTTTTCGTTTGGCggggg                                                                                                                                        >  1:166737/1‑148 (MQ=255)
       aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCTACGAGGAAGGTTTTTAGTTTGCCggggg                                                                                                                                        >  2:54609/1‑145 (MQ=255)
        aTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGGACCCAGGGGAAGCATTTTTTCGACGAGGGAGGTTTTTAGTTTTCCggggg                                                                                                                                        >  2:539788/1‑144 (MQ=255)
              aGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCGGGGGGAAGCATTTTTTTGGCGAGGGAGGTTTTTAGTTTTCCGGGGGGCAGCAggg                                                                                                                               >  1:473450/1‑149 (MQ=255)
                 gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGCAGGGCTGCAGGGACGAGGGGGCGCCTTTTTTCGAGGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCgggggctt                                                                                                                           >  2:322470/1‑146 (MQ=255)
                  cTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTAACGGGGGGACGGCTTTTTTTGGCGAGGGAGGTTTTTAGTTTTCCGGGGGTcagacggg                                                                                                                               >  1:643120/1‑140 (MQ=255)
                  cTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTATCGAGGGGACGCCATTTTTCGACGAGGGAGGTTTTCAGTTTGCCGGGGGTCAGCAGGGGttt                                                                                                                           >  2:684663/1‑149 (MQ=255)
                    tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCCTTTTTTCGACGAGGGAGGTTTTTAGTTTTCCGGGGGGCAGCCGGGCTTTGGtt                                                                                                                       >  2:91077/1‑151 (MQ=255)
                      cTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGAGGGAGGTTTTTAGTTTTCCGGGAGTCAGCAGGGGTTTTGttt                                                                                                                      >  2:443060/1‑150 (MQ=255)
                        aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt                                                                                                                   <  1:441710/151‑1 (MQ=255)
                        aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt                                                                                                                   <  2:553710/151‑1 (MQ=255)
                         aCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTg                                                                                                                  <  2:88699/151‑1 (MQ=255)
                            tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTgg                                                                                                                 <  1:535738/149‑1 (MQ=255)
                              aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCCGCAGTACCGAGGGGAAGCATTTTTTCCACGAGGAAGGTTTTTCGTTTGCCGGGGGGCAGCAGGGCCTTGGTTTTCTGGTag                                                                                                              >  2:193176/1‑150 (MQ=255)
                              aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTGCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTaga                                                                                                             <  1:366511/151‑1 (MQ=255)
                                ggatttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGTTGCAGCACCGAGGGGGGGGACTTTTTTGATGAGGGAGTTTTTTAAATTGCCGGGGGTGCGCCGGGGTTTTGTTTTCTg                                                                                                                  >  2:661702/4‑144 (MQ=255)
                                  ttttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCAGTAAAGAGGGTACGCATTTTTTTGGCGAGGGAGGTTTTTCGTTTGCCGGGGGGCCGCAGGGCTTTGGTTTCCTGGgagaga                                                                                                           >  1:284947/1‑149 (MQ=255)
                                   tttCATGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCATGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTCGTTCCCTGTTAGAGATa                                                                                                         <  1:88745/150‑1 (MQ=255)
                                       aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGGCTTTTCTCGACGAGGGAGATTTTTAGTTTGCCGGGGGGCAGCAGGGGTTTTGTTTTCTGGgagagaaatggtt                                                                                                    >  2:70845/1‑144 (MQ=255)
                                          ttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGCAGTACCGAGGGGACGCATTTTCTCGCCGAGGGAGATTTTTAGGTTGCCGGGAGTCCGCAGGGGTTTGGTTTCCTTGGAGAGATAAGttt                                                                                                    >  1:411366/1‑148 (MQ=255)
                                                 ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCaga                                                                                          <  2:200344/151‑1 (MQ=255)
                                                   gCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGGTGCAGTACCGAGGGTACGCACTTTCTCGACGAGGAAGGTTTTCAGTTTGCCGGGGGTGAGCAGGGGGTTTGGTTCCTGGGAGAGATACGTTTTTTTGGCagag                                                                                         >  1:112320/1‑150 (MQ=255)
                                                      gACGTACCGCCAAACACAAACGCGATATGGTACGGCGGGCAGGGCGCCGAAACCGGGGGACGCATTTTTTCGAAGAGGGAGTTTTTTTGGTTGTCGGGGGGTAGCAGGGGTTTGGGTTCCTGGTAGGGATACGttttt                                                                                                  >  1:564209/1‑137 (MQ=255)
                                                              gCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGGACCGCGGGTACCCATTTTTTCGACGAGGAAGTTTTTTAGGTTGCCGGGGGGCGGCAGGGCTTTTGTTTCCTGGTAGAGATACGGTTTGTTGGGAGAGCCGCCGCCctt                                                                              >  2:230753/1‑150 (MQ=255)
                                                                      caAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgca                                                                     <  1:573055/151‑1 (MQ=255)
                                                                              atatGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAtt                                                             <  1:459174/151‑1 (MQ=255)
                                                                              atatGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAtt                                                             <  1:304212/151‑1 (MQ=255)
                                                                                      cggcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCa                                                     <  1:641176/151‑1 (MQ=255)
                                                                                                     aGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACcgcg                                      <  1:579006/151‑1 (MQ=255)
                                                                                                      gTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGt                                     <  2:442706/151‑1 (MQ=255)
                                                                                                                                   ggAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGt         <  1:598788/151‑1 (MQ=255)
                                                                                                                                           ttCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGgtg  <  2:473469/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTG  >  NC_000913/4346563‑4346850

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