Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,666,1730AC52.9% 5.8 / 13.7 19D54A (GAC→GCC) xerDTyrosine recombinase XerD
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base C (9/0);  total (14/5)
Fisher's exact test for biased strand distribution p-value = 9.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGGAAAAAGGCCTGTCCGACAACACCCGCGTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTCCACGGCGCGTTTCCTCTCCGGCCTGCGCGGCTTCTTCCGTTACCTGCTGCGGGAAAAGCTGGTGGC  >  KT2440/1666065‑1666293
                                                                                                            |                                                                                                                        
tGGAAAAAGGCCTGTCCGACAACACCCGCGTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTTTCGCTGCCGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGgc                                                                                                <  7:100820/135‑1 (MQ=255)
     aaaGGCCTGTCCGACAACACCCGCGTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGCGCCGCGCTCACGTGGGCTGGCGcccg                                                                                           >  1:98476/1‑133 (MQ=255)
     aaaGGCCTGTCCGACAACACCCGCGTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCAGATCCTCGATCACTTGGCCTGGCGCCt                                                                                            >  4:96791/1‑134 (MQ=255)
                 gACAACACCCGCGTGTCGTTTCGCAGCGACCTGGCCCTGATCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTAc                                                                               <  6:23492/135‑1 (MQ=255)
                            cgTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTc                                                                    <  6:89040/135‑1 (MQ=255)
                                              ccTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGGCCAGATCCTCGATCACTTTGCCTGGCGCCTGGATCCGGGCTACAAGCCACGTTCCACGGCGCGGTTCCtctc                                                  >  3:39738/1‑135 (MQ=255)
                                                tGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGGCCTGCTGCTCGATCACTTGGCCTGGCGCCGGGATCAGGGCTACAAGCCACGTTCCACGGCGCGGTTGGTCTCCg                                                >  5:60116/1‑135 (MQ=255)
                                                tGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGCTCCTCGATGACTTGGCCTGGCGCCTGGATCAGGGCGACAAGCCACGTTCCACGGCGCGTTTCCTCTACg                                                >  7:143804/1‑135 (MQ=255)
                                                     cTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGATCCGCGCTCACTTTGCCTGGCGCCGGGGTCAGGGCGACAAGCCCCGTTCCACGGCGCGTTTACTCTCCGGccgg                                           >  1:125806/1‑133 (MQ=255)
                                                           aaTGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCAGATACTCGCTCCAATTGCCTGGCGCCTGGCTCAGGGGGACAACCCCAGTTCCACGGCGCGGTTTCACAGCGGCCAGCCCGGc                                     >  6:147604/1‑135 (MQ=255)
                                                                cTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGCTCCTCGCTCACTTGGCCTGGCGCCTGGATCAGGGCTACCACCCCCCTTCCACGGCGCGGTTGCTCTCCGGCCTGCGCGGcttcgt                                >  2:3720/1‑133 (MQ=255)
                                                                   gTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTCCACGGCGCGTTTCCTCTCCGGCCTGCGCGGCTTCTTCCg                             >  8:30568/1‑135 (MQ=255)
                                                                     tGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCCCGTTCCACGGCGCGTTTCCTCTCCGGCCTGCGCGGCTTCTTCCGtt                           >  4:159668/1‑135 (MQ=255)
                                                                     tGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTCCACGGCGCGGTTCCGCTCCGGCCTGCGCGGCGTCTTCCGtt                           >  5:206290/1‑135 (MQ=255)
                                                                     tGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATcc                                                                                                                  >  4:177011/1‑48 (MQ=255)
                                                                     tGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATcc                                                                                                                  <  3:177011/48‑1 (MQ=255)
                                                                                    tATCGCTGCCGGACGCCGGTCGCGCCCTGCTGCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCCCGGTGCACGGCGCGGTTCCTCGCCGGCCTGCGCCGCTTATTGCGTTTCCTGCTGCgggc              >  8:10017/1‑132 (MQ=255)
                                                                                            ccGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTCCACGGCGCGTTTCCTCTCCGGCCTGCGCGGCTTCTTCCGTTACCTGCTGCGGGAAAAGCtggtg    >  7:19640/1‑135 (MQ=255)
                                                                                              ggACGCGGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCATGGATGAGGGCTACAAGCGAGGGTACACGGCGAGTTTCGTCTCAGGCCTGCGCGGCTTCTTCCGTTACCTGCTGCGTGAAATGGTGGTGGc  <  7:156210/135‑1 (MQ=255)
                                                                                                            |                                                                                                                        
TGGAAAAAGGCCTGTCCGACAACACCCGCGTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTCCACGGCGCGTTTCCTCTCCGGCCTGCGCGGCTTCTTCCGTTACCTGCTGCGGGAAAAGCTGGTGGC  >  KT2440/1666065‑1666293

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