Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,457,8950CT21.4% 54.0 / 10.4 28G133G (GGC→GGT) paaF2,3‑dehydroadipyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base C (10/12);  new base T (3/3);  total (13/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 = 9.66e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGGCGCGATTGCAGGCCTTTAATAAACCACTTATCGCAGCCGTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGA  >  NC_000913/1457749‑1458038
                                                                                                                                                  |                                                                                                                                               
tGGGCGCGATTGCAGGCCTTTAATAAACCACTTATCGCAGCCGTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGcat                                                                                                                                               >  3:105260/1‑149 (MQ=255)
             aGGCCTTTAATAAACCACTTATCGCAGCCGTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGTATCATGCCTGgcg                                                                                                                                    >  1:135531/1‑147 (MQ=255)
             aGGCCTTTAATAAACCACTTATCGCAGCCGTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGTATCATGCCTGgcg                                                                                                                                    >  1:138040/1‑147 (MQ=255)
             aGGCCTTTAATAAACCACTTATCGCAGCCGTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTCTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGTATCATGCCTGgcg                                                                                                                                    >  7:139646/1‑147 (MQ=255)
                      ataaACCACTTATCGCAGCCGTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACg                                                                                                                         >  5:91206/1‑149 (MQ=255)
                             aCTTATCGCAGCCGTGAATGGTTACGCGCTTTTGGCGGGTTTCGAACTGGGATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGTATCATGCCTGGCGCAGGCGGAACGCAACGtt                                                                                                                  <  8:139646/149‑1 (MQ=255)
                             aCTTATCGCAGCCGTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGTATCATGCCTGGCGCAGGCGGAACGCAACGtt                                                                                                                  <  2:135531/149‑1 (MQ=255)
                             aCTTATCGCAGCCGTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGTATCATGCCTGGCGCAGGCGGAACGCAACGtt                                                                                                                  <  2:138040/149‑1 (MQ=255)
                                 aTCGCAGCCGTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAAt                                                                                                              >  2:58298/1‑149 (MQ=255)
                                         cGTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAgtg                                                                                                      >  6:71636/1‑149 (MQ=255)
                                          gTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAgtgt                                                                                                     <  3:70055/149‑1 (MQ=255)
                                                         ctgggggCGGGTTGCGAACTGGCATTGGTGTGCGATGTGGTGGTTGCAGGAGAGAACGCGCGTTTGGGGGTGCCGGAAAGCACTCTCGGCATCATGCCTGGCGCAGGAGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGc                                                                                      <  4:105260/146‑1 (MQ=255)
                                                          ttGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGc                                                                                      >  5:18429/1‑148 (MQ=255)
                                                            ggggCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGcc                                                                                     <  5:71636/147‑1 (MQ=255)
                                                                   gTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATg                                                                            <  1:58298/149‑1 (MQ=255)
                                                                       cGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTT‑GGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCt                                                                       >  1:67333/1‑149 (MQ=255)
                                                                         aaCTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTT‑GGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGa                                                                     <  2:67333/149‑1 (MQ=255)
                                                                              gCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCg                                                                  <  6:18429/148‑1 (MQ=255)
                                                                                aTTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTg                                                                         >  5:63180/1‑139 (MQ=255)
                                                                                aTTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTg                                                                         <  6:63180/139‑1 (MQ=255)
                                                                                           aTGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCAcc                                                    <  6:91206/149‑1 (MQ=255)
                                                                                                          gagagaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcag                                     >  7:7360/1‑149 (MQ=255)
                                                                                                                  gcgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCgg                             <  8:7360/149‑1 (MQ=255)
                                                                                                                                      aaTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTcc         >  1:119537/1‑149 (MQ=255)
                                                                                                                                       aTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCAcc                                                    <  4:82398/105‑1 (MQ=255)
                                                                                                                                       aTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCAcc                                                    >  3:82398/1‑105 (MQ=255)
                                                                                                                                           ctctCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTc          <  4:36714/143‑1 (MQ=255)
                                                                                                                                           ctctCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTc          >  3:36714/1‑143 (MQ=255)
                                                                                                                                                ggCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGa  <  2:119537/146‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TGGGCGCGATTGCAGGCCTTTAATAAACCACTTATCGCAGCCGTCAATGGTTACGCGCTTGGGGCGGGTTGCGAACTGGCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGA  >  NC_000913/1457749‑1458038

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