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.

CTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGG  >  KT2440/2780587‑2780823
                                                                                                                         |                                                                                                                     
cTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACCCCCGGGCGACCTCGCCAGGGCGGTTTTTTTTGCGAGTTTTGaaa                                                                                                         >  3:362310/1‑134 (MQ=255)
   taCCTCTTAGGAAATATCTAAGGTCATGCCAATAACGCGGGCGTTAAGCAGCCGCAAAAGCTTAAGTACCTACCGCCCCTTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTTTATCCGAGTTTTTAAAAAgg                                                                                                     <  1:14684/134‑1 (MQ=255)
                     tAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACCCCCGGGGGGCCACGGCGGGGCGGGTTTGTATGCGCGTTTTGTAAAAGGCATCGCCCGGCTGTTggg                                                                                   >  7:202003/1‑134 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGGTTTTGTTTCCGAGTTTTGAAAAAGGAATCGCCCGGATGTTGGCCGCCACGCt                                                                          >  6:190102/1‑134 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCGGGGCGTGTTTGTTTCCGCGGTTTGGAAAAGGAATCGCCCGGTTGTTGGGCGCCCCGcg                                                                          >  7:104903/1‑133 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCGGGGCGGTTTTGTTTCCGAGTTTTGAAAAATGCAAAGCCCGGTTGTTTGGCGCCACgag                                                                          >  6:252768/1‑132 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCCGGCGCCCGCGCCTGGGCGTTTTTGTTTACCAGGTTTTTAAAAAGAATCGCCCGGATTTTGGACGCCACGcg                                                                          >  8:14803/1‑133 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACCCCCGGGCGACCTCGCCTGGGCGGTGTTGTAGCCGAGGTTTGAAAAAGGATTCGGCCGGATTGTGGACGCCCCGcct                                                                         >  3:338759/1‑133 (MQ=255)
                               cccATAACGCGGGCGGTAAGCAGCAGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCCACCTAGCCTGGGGGTTTTTGTATACGAGTTTTTAAAAATGCATGGCCCCTAAGTTGGGCGCCAccgg                                                                          >  5:242577/1‑131 (MQ=255)
                                ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAAACCCCGGGGGCCCACGCCAGGGCGGGTTTTGTTACGGGGTTTTGAAAAGGAATCGCCCCGATGTTGGTCGcccccggg                                                                         >  2:253579/1‑128 (MQ=255)
                                      cgcgAGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGGTTTGAAAAAGGGATCGCCCGTATGTTGGACGCCACGCTTAAATCg                                                                   >  4:172708/1‑134 (MQ=255)
                                            cggTACGCCGCCGCCAAAGCTGAAGTACCAACCGCCCCGTAAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGa                                                            <  1:112699/135‑1 (MQ=255)
                                            cggTAAGCAGCCGCCAATCGTGAAGGACCTACCGCCCGGTCAAACGCCCGGGCGACCGGGGCTGGGCTTTTTTGTATCCTAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACt                                                              <  6:237761/133‑1 (MQ=255)
                                             ggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGCCCTCGCCTGGGCGTTTTTGTTTGCGAGTTTTGAAAAAGGAATCGGCCGGTTGGTGGCCGCCCCGCGTAAATCGCAACTGa                                                            >  3:354147/1‑134 (MQ=255)
                                             ggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGCCCTCGCCGGGGCGGTTTTGTTTGCGCCGTTTGAAAAATGAATAGCCCGGATGTTGGTCGCCCCGCCTAAAACGCAAcagt                                                            >  7:289758/1‑131 (MQ=255)
                                                         ccAAAGCTGATGTACCTCCCGCCCCGGCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGa                                               <  4:240329/135‑1 (MQ=255)
                                                                            cGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGATTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGAt                             <  7:188279/134‑1 (MQ=255)
                                                                               cccGTCAAACGCCCGGGCGCCCTCGCCGGGGGCGTTTTTTTTAAGAGTTTTGAAAAAGGAATAGCCCGGATGTTGGCCGCCACGCGTAAAAAGGAACAGAAAACCAACATGGAGCGGGTCACCCACCCgcacaa                          >  8:76930/1‑129 (MQ=255)
                                                                                  gTCAAACCCCCGGGCGACCTCGCCTGGGCGTTTTTGTTTCCGCGTTTTGTAAAAGGAATAGCCCCGATGTTGGGCGCCCCGCGGTAATAGCCACTGAAAACCAACATTGTGCGGGGGACCCCCCCGATCGCGAta                      >  6:64949/1‑134 (MQ=255)
                                                                                                cgcccTCGCCGGGGCGTTTTTGTTGCCGAGGTTGGCAAAAGGAATCGCCCGGATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcga        <  5:5013/132‑1 (MQ=255)
                                                                                                      cGCCTGGGCGTTTTTGTTTCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCCCGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGCTCGAGATCGTTGCGTCCCTCGACGACgg  >  8:194556/1‑135 (MQ=255)
                                                                                                                         |                                                                                                                     
CTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGG  >  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: