Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,151,9670AG23.3% 47.9 / 13.5 30Q231R (CAG→CGG) eddphosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (11/12);  new base G (4/3);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCAT  >  NC_002947/1151834‑1152102
                                                                                                                                     |                                                                                                                                       
aacaaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGTCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGaa                                                                                                                            <  1:113679/147‑1 (MQ=255)
aacaaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGCCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACTGGTTGGTGTTg                                                                                                                               <  1:130923/144‑1 (MQ=255)
     gCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAATGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATg                                                                                                                      >  6:26503/1‑148 (MQ=255)
        gAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATgggg                                                                                                                   >  4:96801/1‑148 (MQ=255)
             ccGATGTGCGCCAGCGCTACGCCGAAGGCAAGCCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACCGCCCCGGCACCTGCACCTACTACGGCTCCGCCAACCCCAACCAGTGGGTGAGGGAAGGGGTGGTGCTGCa                                                                                                              >  5:100526/1‑148 (MQ=255)
                       ccAGCGCTATGCCGAAGGCAATGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGgc                                                                                                    <  5:26503/148‑1 (MQ=255)
                       ccAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGGCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGGTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGgc                                                                                                    <  2:111220/148‑1 (MQ=255)
                                 gCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttc                                                                                           <  8:2973/147‑1 (MQ=255)
                                  ccGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttcg                                                                                          >  5:91854/1‑147 (MQ=255)
                                          aaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCCCCTGCACCTTATACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGGTCTTCAACCCGt                                                                                 >  7:118363/1‑148 (MQ=255)
                                                aGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACAcgc                                                                           <  1:17955/148‑1 (MQ=255)
                                                             tGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGc                                                              >  7:12877/1‑148 (MQ=255)
                                                                             gAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCt                                              >  2:110356/1‑148 (MQ=255)
                                                                                   cTACCACAGCCCCGACACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGcc                                                                                                        >  1:114134/1‑84 (MQ=255)
                                                                                    tACCACAGCCCCGACACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCg                                                                                                       <  2:114134/84‑1 (MQ=255)
                                                                                       cacaGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCgcagca                                      <  1:71171/146‑1 (MQ=255)
                                                                                       cacaGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAgg                                    <  3:114287/148‑1 (MQ=255)
                                                                                       cacaGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAgg                                    >  1:131173/1‑148 (MQ=255)
                                                                                       cacaGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAgg                                    >  1:113846/1‑148 (MQ=255)
                                                                                                     ctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttcg                                                                                          >  6:126492/1‑80 (MQ=255)
                                                                                                      tgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgttcgt                                                                                         <  5:126492/80‑1 (MQ=255)
                                                                                                        cacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCGGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAg                   <  2:66764/148‑1 (MQ=255)
                                                                                                         acctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAg                   <  7:89289/147‑1 (MQ=255)
                                                                                                                       ttccATCACCAACTGGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTc    >  5:137619/3‑148 (MQ=255)
                                                                                                                        gCCATCACCAACCAGTTGGTGCTGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCa   <  3:88685/148‑1 (MQ=255)
                                                                                                                         gaatcacCAACTGGTTGGTGTTGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTc    <  3:19190/142‑1 (MQ=255)
                                                                                                                         ccATCACCAACTGGTTGGTGTTGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCAt  >  1:36707/1‑148 (MQ=255)
                                                                                                                         ccATCACCAACTGGTTGGTGTTGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCAt  >  1:34576/1‑148 (MQ=255)
                                                                                                                           atcacCAACTGGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGAcgcg                             >  2:55141/3‑119 (MQ=255)
                                                                                                                            tcacCAACTGGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACgcgc                            <  1:55141/118‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
AACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCAT  >  NC_002947/1151834‑1152102

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