Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,671,9950GC31.2% 18.6 / 12.3 16A763G (GCC→GGC) acnBaconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base G (6/5);  new base C (2/3);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGCGTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCACCCTTGACCTTCTCCAGCAGCTTGC  >  NC_002947/2671871‑2672112
                                                                                                                            |                                                                                                                     
ggTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGGCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGCGTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAg                                                                                                            <  2:409279/136‑1 (MQ=255)
                cGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGCGTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCctt                                                                                            <  1:26415/136‑1 (MQ=255)
                                gggTAGAGGTCGAAACCACCGTTGCACCGTCCGCCACACGCGCCTGGTTACCCATACACAGGGAACAGCCAGGGATCTCGATACGCGCACCACTCTTACCGAAGACGCTGTAGTAGCCTTCTTCGGTCAACTgtgc                                                                            <  2:434611/136‑4 (MQ=255)
                                               ccACGGTCGAACCGGTCTGCACGCGTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCttcttc                                                                                        >  1:196175/1‑109 (MQ=255)
                                               ccACGGTCGAACCGGTCTGCACGCGTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCttcttc                                                                                        <  2:196175/109‑1 (MQ=255)
                                                     ccGAACCGGTCTGCACCCGTGCCTGGATGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCgg                                                       <  2:432596/135‑1 (MQ=255)
                                                                  cacaCGCGCCTGGTTACCCATACACAGGGAACAGCCAGGGATCTCGATACGCGCACCACTCTTACCGAAGACGCTGTAGTAGCCTTCTTCGGTCAACTGTGCGGCGTCCATACGGGTTGGCgg                                                       <  1:28964/123‑1 (MQ=255)
                                                                  cacaCGCGCCTGGTTACCCATACACAGGGAACAGCCAGGGATCTCGATACGCGCACCACTCTTACCGAAGACGCTGTAGTAGCCTTCTTCGGTCAACTGTGCGGCGTCCATACGGGTTGGCgg                                                       <  1:31032/123‑1 (MQ=255)
                                                                  cacaCGCGCCTGGTTACCCATACACAGGGAACAGCCAGGGATCTCGATACGCGCACCACTCTTACCGAAGACGCTGTAGTAGCCTTCTTCGGTCAACTGTGCGGCGTCCATACGGGTTGGCgg                                                       >  2:28964/1‑123 (MQ=255)
                                                                  cacaCGCGCCTGGTTACCCATACACAGGGAACAGCCAGGGATCTCGATACGCGCACCACTCTTACCGAAGACGCTGTAGTAGCCTTCTTCGGTCAACTGTGCGGCGTCCATACGGGTTGGCgg                                                       >  2:31032/1‑123 (MQ=255)
                                                                                    cATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCAccc                        >  1:433951/1‑136 (MQ=255)
                                                                                                gcagcCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCACCCTTGACCTTCTcc            >  1:172564/1‑136 (MQ=255)
                                                                                                gcagcCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCACCCTTGACCTTCTcc            >  1:280272/1‑136 (MQ=255)
                                                                                                    ccTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCACCCTTTACCTTCTCcagcc        >  1:176572/1‑135 (MQ=255)
                                                                                                        gCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCACCCTTAACCTTCTCCAACAACtt    >  1:173627/1‑136 (MQ=255)
                                                                                                          aTTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCACCCTTGACCTTCTCCAGCAGCTTGc  <  1:246686/136‑1 (MQ=255)
                                                                                                                            |                                                                                                                     
GGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGCGTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCACCCTTGACCTTCTCCAGCAGCTTGC  >  NC_002947/2671871‑2672112

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