Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,780,7060AT57.1% ‑1.5 / 10.2 21intergenic (+96/‑29)ahpC_2/ahpFAlkyl hydroperoxide reductase C/Alkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base T (11/1);  total (15/6)
Fisher's exact test for biased strand distribution p-value = 4.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGC‑GTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGG  >  KT2440/2780587‑2780823
                                                                                                                         |                                                                                                                     
cTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACCCCCGGGCGACCTCGCC‑AGGGC‑GGTTTTTTTTGCGAGTTTTGaaa                                                                                                         >  3:362310/1‑134 (MQ=255)
   taCCTCTTAGGAAATATCTAAGGTCATGCCAATAACGCGGGCGTTAAGCAGCCGCAAAAGCTTAAGTACCTACCGCCCCTTCAAACGCCCGGGCGACCTCGCC‑TGGGC‑GTTTTTTTATCCGAGTTTTTAAAAAgg                                                                                                     <  1:14684/134‑1 (MQ=255)
                     tAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACCCCCGGGGGGCCACGGC‑GGGGC‑GGGTTTGTATGCGCGTTTTGTAAAAGGCATCGCCCGGCTGTTggg                                                                                   >  7:202003/1‑134 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGC‑GGTTTTGTTTCCGAGTTTTGAAAAAGGAATCGCCCGGATGTTGGCCGCCACGCt                                                                          >  6:190102/1‑134 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑GGGGC‑GTGTTTGTTTCCGCGGTTTGGAAAAGGAATCGCCCGGTTGTTGGGCGCCCCGcg                                                                          >  7:104903/1‑133 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑GGGGC‑GGTTTTGTTTCCGAGTTTTGAAAAATGCAAAGCCCGGTTGTTTGGCGCCACgag                                                                          >  6:252768/1‑132 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCCGGCGCCCGCGCC‑TGGGC‑GTTTTTGTTTACCAGGTTTTTAAAAAGAATCGCCCGGATTTTGGACGCCACGcg                                                                          >  8:14803/1‑133 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACCCCCGGGCGACCTCGCC‑TGGGC‑GGTGTTGTAGCCGAGGTTTGAAAAAGGATTCGGCCGGATTGTGGACGCCCCGcct                                                                         >  3:338759/1‑133 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCAGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCCACCTAGCC‑TGGGG‑GTTTTTGTATACGAGTTTTTAAAAATGCATGGCCCCTAAGTTGGGCGCCAccgg                                                                          >  5:242577/1‑131 (MQ=255)
                                ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAAACCCCGGGGGCCCACGCC‑AGGGCGGGTTTTGTTACGGGGTTTTG‑AAAAGGAATCGCCCCGATGTTGGTCGcccccggg                                                                         >  2:253579/1‑128 (MQ=255)
                                      cgcgAGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGC‑GTTTTTGTATCCGAGGTTTGAAAAAGGGATCGCCCGTATGTTGGACGCCACGCTTAAATCg                                                                   >  4:172708/1‑134 (MQ=255)
                                            cggTACGCCGCCGCCAAAGCTGAAGTACCAACCGCCCCGTAAAACGCCCGGGCGACCTCGCC‑TGGGC‑GTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGa                                                            <  1:112699/135‑1 (MQ=255)
                                            cggTAAGCAGCCGCCAATCGTGAAGGACCTACCGCCCGGTCAAACGCCCGGGCGACCGGGGC‑TGGGC‑TTTTTTGTATCCTAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACt                                                              <  6:237761/133‑1 (MQ=255)
                                             ggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGCCCTCGCC‑TGGGC‑GTTTTTGTTTGCGAGTTTTGAAAAAGGAATCGGCCGGTTGGTGGCCGCCCCGCGTAAATCGCAACTGa                                                            >  3:354147/1‑134 (MQ=255)
                                             ggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGCCCTCGCC‑GGGGC‑GGTTTTGTTTGCGCCGTTTGAAAAATGAATAGCCCGGATGTTGGTCGCCCCGCCTAAAACGCAAcagt                                                            >  7:289758/1‑131 (MQ=255)
                                                         ccAAAGCTGATGTACCTCCCGCCCCGGCAAACGCCCGGGCGACCTCGCC‑TGGGC‑GTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGa                                               <  4:240329/135‑1 (MQ=255)
                                                                            cGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGC‑GTTTTTGTATCCGATTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGAt                             <  7:188279/134‑1 (MQ=255)
                                                                               cccGTCAAACGCCCGGGCGCCCTCGCCGGGGGC‑GTTTTT‑TTTAAGAGTTTTGAAAAAGGAATAGCCCGGATGTTGGCCGCCACGCGTAAAAAGGAACAGAAAACCAACATGGAGCGGGTCACCCACCCgcacaa                          >  8:76930/1‑129 (MQ=255)
                                                                                  gTCAAACCCCCGGGCGACCTCGCC‑TGGGC‑GTTTTTGTTTCCGCGTTTTGTAAAAGGAATAGCCCCGATGTTGGGCGCCCCGCGGTAATAGCCACTGAAAACCAACATTGTGCGGGGGACCCCCCCGATCGCGAta                      >  6:64949/1‑134 (MQ=255)
                                                                                                cgcccTCGCC‑GGGGC‑GTTTTTGTTGCCGAGGTTGGCAAAAGGAATCGCCCGGATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcga        <  5:5013/132‑1 (MQ=255)
                                                                                                      cGCC‑TGGGC‑GTTTTTGTTTCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCCCGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGCTCGAGATCGTTGCGTCCCTCGACGACgg  >  8:194556/1‑135 (MQ=255)
                                                                                                                         |                                                                                                                     
CTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCC‑TGGGC‑GTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGG  >  KT2440/2780587‑2780823

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