Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,094,8310TG69.6% ‑3.5 / 10.7 26K178N (AAA→AACCJ019_01995hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base G (16/0);  total (22/4)
Fisher's exact test for biased strand distribution p-value = 3.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATGGGCAGCGCGGGTCA  >  KT2440/2094707‑2094961
                                                                                                                            |                                                                                                                                  
cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCGGGGCGTTTTGGAggg                                                                                                                          >  1:249473/1‑135 (MQ=255)
cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCGGGGCGTTTTGCAggg                                                                                                                          >  4:354325/1‑135 (MQ=255)
cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCCGGCGGGGCGGGGCGGTTTGCgggg                                                                                                                          >  5:203233/1‑135 (MQ=255)
     gCCCGGGCCAACAGTTTCTCAAGATACGTCATCTTCTCTTCAAACACCTTCGGCTCCGGCAGCTTGCGCAAATACTGCACGCCGCGCTCGAACGCCCTGCCCGCCAGGCCGGGCCTGCCCTTTTGCAagggcggt                                                                                                                     >  6:291884/1‑130 (MQ=255)
     gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGCGCGCCAGGCCGGGCTGGCCGTTTTTCAGGGCCAGt                                                                                                                     >  2:257340/1‑135 (MQ=255)
             cAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCAGGCCCGGCTGGCCGTTTTTCACGGCCAGGTGGCCGAg                                                                                                             >  6:167110/1‑135 (MQ=255)
                        cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCGGGTGGCCGAGGGTGTTgaag                                                                                                   >  2:6548/1‑134 (MQ=255)
                        cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCGGGCCGTTTTTGAGGGCCGGGTGGCCGAGGTTGTTgaag                                                                                                   >  5:142130/1‑134 (MQ=255)
                                                acacCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCTGGCCGGGCTGGCCGGTTTGCAGGGCCGGGTGGGCGTCGGTGGTGAAGAACTCTATATGCAGCAGCCCCAgca                                                                          >  7:252092/1‑133 (MQ=255)
                                                     ttCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCGGGCCGGGTTGCAGGGCCGGGTGGGCGAGGGTGTTGTAGCACTCTATATTCAGCCGCACCACGAGATgg                                                                     >  3:362116/1‑135 (MQ=255)
                                                        ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCGGGCCGTTTTGCAGGGCCTGTTGGGCGAGGTTGTTGAAGAACTCTATATTCAGCAGCCCCAGGATATGGc                                                                    >  2:226949/1‑133 (MQ=255)
                                                        ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCACGATATGGCgg                                                                  >  8:167118/1‑135 (MQ=255)
                                                          cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGGAGGGCCGGGTGGCCGCGGGTGTTGAAGCACTCTATATTCAGCAGCCCCCGGCTATGGCGGag                                                                >  8:12440/1‑134 (MQ=255)
                                                          cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGCGCGCCTGGCCGGGCTGGCGGTTTTGCACGGGCGGTTGGGCGCTGGTGTTGTAGACCTCTATAATCAGCAGCACCAGGAGATGGCggct                                                                >  3:59691/1‑133 (MQ=255)
                                                                 tgctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGGCGTTTTGCAGGGGCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATTGCGGATTTCCTTg                                                         >  4:89352/1‑135 (MQ=255)
                                                                 tgctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCGGGCCGTTTTGCCGGGCCTGTTGGCCGCGGGTGTTGAAGAACGCTATATTCCGCCGCCCCAGGTTATGGGGGCtttgcgtt                                                         >  4:321427/1‑130 (MQ=255)
                                                                         cAAATACTCCACGCCGCGCTCGACGCCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATg                                                 <  4:56388/135‑1 (MQ=255)
                                                                         aaattACGGCACGCCGCGCGCGATGGCCGGGCCGGCCTGGCCGGGCTGTCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGTCGTATTTCCTTGACCCaa                                                   <  7:289103/129‑1 (MQ=255)
                                                                                    cgccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTTCCGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGCTTTCCTTGACCCAATGCTTGGCCACg                                       >  8:5430/1‑134 (MQ=255)
                                                                                     gccgcGCTCGAAGGCCCCGCGCGCCTGGCCGGGCTGGCCATTTTGCCGGGGCTGGTGGGCGAGGTGTTTGAAGAACTCGATATGCAGCAGCGCCCGGATATTGGGGATTTACTTGACCCCATGCGTGGCCTCgcc                                     >  6:258317/1‑133 (MQ=255)
                                                                                     gccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCATTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTACTTGACCCAATGCATGGCCTCGAc                                     >  6:335032/1‑135 (MQ=255)
                                                                                         cgcTCGACGCCCCGGCCCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCg                                 <  6:142130/135‑1 (MQ=255)
                                                                                          gcTCGAAGGCCCGGCGCGCCAGGCCGGGGTGGGCGTTTTTCCGGGGCGGTTGGGCGCGGGTGTTGAAGAACTCGATATTCAGCCGCACCAGGATATGGCGGGTTTCCTTGGCCCAATGCTTGGCCGCGACAGCgg                                >  4:2358/1‑135 (MQ=255)
                                                                                          gcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGCCTGGTGGCCGAGGGTGTTGAAGCACTCTATATGCCGCCGCACCCGGCTATGGCGGATTTCCGTGACCCAATGCGTGGCCTCGACACCgg                                >  6:157637/1‑135 (MQ=255)
                                                                                             ctacggCCAGCCGCGGCGGCCCGGGCTGTCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCa                             <  4:325851/131‑1 (MQ=255)
                                                                                                                        ggCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAACAAGCCCTCATGGGCAGCGCGGGTCa  >  2:362987/1‑135 (MQ=255)
                                                                                                                            |                                                                                                                                  
CATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATGGGCAGCGCGGGTCA  >  KT2440/2094707‑2094961

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: