Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,6850CT44.1% 26.7 / 15.4 34E214E (GAG→GAA) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base C (4/15);  new base T (15/0);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 2.24e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGT  >  NC_000913/4346537‑4346815
                                                                                                                                                    |                                                                                                                                   
cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTtctc                                                                                                                                   <  2:180018/151‑1 (MQ=255)
              tCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGttt                                                                                                                     <  2:417497/150‑1 (MQ=255)
                          tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGAAGGTTTTTAGTTTGCCgg                                                                                                        >  2:543681/1‑151 (MQ=255)
                            aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGGGGGCGCAGTACCGAGGGGACGCATTTTTTTGACGAGGGAGGTTTTTAAGTTGGCgggg                                                                                                       >  1:448285/1‑151 (MQ=255)
                            aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGa                                                                                                      <  1:525649/151‑1 (MQ=255)
                              cGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGGAGTACCGCGGGGACGCATTTTTTCGACGA‑GGGAGGTTTTTCGTTTGGCggggg                                                                                                     >  1:166737/1‑148 (MQ=255)
                                 aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCTACGA‑GGAAGGTTTTTAGTTTGCCggggg                                                                                                     >  2:54609/1‑145 (MQ=255)
                                  aTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGGACCCAGGGGAAGCATTTTTTCGACGA‑GGGAGGTTTTTAGTTTTCCggggg                                                                                                     >  2:539788/1‑144 (MQ=255)
                                        aGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCGGGGGGAAGCATTTTTTTGGCGA‑GGGAGGTTTTTAGTTTTCCGGGGGGCAGCAggg                                                                                            >  1:473450/1‑149 (MQ=255)
                                           gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGCAGGGCTGCAGGGACGAGGGGGCGCCTTTTTTCGAGGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCgggggctt                                                                                        >  2:322470/1‑146 (MQ=255)
                                            cTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTAACGGGGGGACGGCTTTTTTTGGCGA‑GGGAGGTTTTTAGTTTTCCGGGGGTcagacggg                                                                                            >  1:643120/1‑140 (MQ=255)
                                            cTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTATCGAGGGGACGCCATTTTTCGACGA‑GGGAGGTTTTCAGTTTGCCGGGGGTCAGCAGGGGttt                                                                                        >  2:684663/1‑149 (MQ=255)
                                              tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCCTTTTTTCGACGA‑GGGAGGTTTTTAGTTTTCCGGGGGGCAGCCGGGCTTTGGtt                                                                                    >  2:91077/1‑151 (MQ=255)
                                                cTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGA‑GGGAGGTTTTTAGTTTTCCGGGAGTCAGCAGGGGTTTTGttt                                                                                   >  2:443060/1‑150 (MQ=255)
                                                  aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt                                                                                <  2:553710/151‑1 (MQ=255)
                                                  aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt                                                                                <  1:441710/151‑1 (MQ=255)
                                                   aCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTg                                                                               <  2:88699/151‑1 (MQ=255)
                                                      tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTgg                                                                              <  1:535738/149‑1 (MQ=255)
                                                        aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCCGCAGTACCGAGGGGAAGCATTTTTTCCACGA‑GGAAGGTTTTTCGTTTGCCGGGGGGCAGCAGGGCCTTGGTTTTCTGGTag                                                                           >  2:193176/1‑150 (MQ=255)
                                                        aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTGCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTaga                                                                          <  1:366511/151‑1 (MQ=255)
                                                          ggatttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGTTGCAGCACCGAGGGGGGGGACTTTTTTGATGA‑GGGAGTTTTTTAAATTGCCGGGGGTGCGCCGGGGTTTTGTTTTCTg                                                                               >  2:661702/4‑144 (MQ=255)
                                                            ttttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCAGTAAAGAGGGTACGCATTTTTTTGGCGA‑GGGAGGTTTTTCGTTTGCCGGGGGGCCGCAGGGCTTTGGTTTCCTGGgagaga                                                                        >  1:284947/1‑149 (MQ=255)
                                                             tttCATGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCATGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTCGTTCCCTGTTAGAGATa                                                                      <  1:88745/150‑1 (MQ=255)
                                                                 aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGGCTTTTCTCGACGA‑GGGAGATTTTTAGTTTGCCGGGGGGCAGCAGGGGTTTTGTTTTCTGGgagagaaatggtt                                                                 >  2:70845/1‑144 (MQ=255)
                                                                    ttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGCAGTACCGAGGGGACGCATTTTCTCGCCGA‑GGGAGATTTTTAGGTTGCCGGGAGTCCGCAGGGGTTTGGTTTCCTTGGAGAGATAAGttt                                                                 >  1:411366/1‑148 (MQ=255)
                                                                           ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCaga                                                       <  2:200344/151‑1 (MQ=255)
                                                                             gCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGGTGCAGTACCGAGGGTACGCACTTTCTCGACGA‑GGAAGGTTTTCAGTTTGCCGGGGGTGAGCAGGGGGTTTGGTTCCTGGGAGAGATACGTTTTTTTGGCagag                                                      >  1:112320/1‑150 (MQ=255)
                                                                                gACGTACCGCCAAACACAAACGCGATATGGTACGGCGGGCAGGGCGCCGAAACCGGGGGACGCATTTTTTCGAAGA‑GGGAGTTTTTTTGGTTGTCGGGGGGTAGCAGGGGTTTGGGTTCCTGGTAGGGATACGttttt                                                               >  1:564209/1‑137 (MQ=255)
                                                                                        gCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGGACCGCGGGTACCCATTTTTTCGACGA‑GGAAGTTTTTTAGGTTGCCGGGGGGCGGCAGGGCTTTTGTTTCCTGGTAGAGATACGGTTTGTTGGGAGAGCCGCCGCCctt                                           >  2:230753/1‑150 (MQ=255)
                                                                                                caAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgca                                  <  1:573055/151‑1 (MQ=255)
                                                                                                        atatGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAtt                          <  1:459174/151‑1 (MQ=255)
                                                                                                        atatGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAtt                          <  1:304212/151‑1 (MQ=255)
                                                                                                                cggcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCa                  <  1:641176/151‑1 (MQ=255)
                                                                                                                               aGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACcgcg   <  1:579006/151‑1 (MQ=255)
                                                                                                                                gTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGt  <  2:442706/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                   
CATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGT  >  NC_000913/4346537‑4346815

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: ‑