Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,123,2380CA52.4% 7.3 / 31.0 43intergenic (‑240/+64)gcvH_1/CJ019_01062Glycine cleavage system H protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (5/15);  new base A (22/0);  total (28/15)
Fisher's exact test for biased strand distribution p-value = 1.57e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGA  >  KT2440/1123115‑1123360
                                                                                                                           |                                                                                                                          
ttAAATTACCGGAGAAGCCTTATTAATTGTGGGTAATGATATCGTTACAGGTCCCTCTGAGAGCAAAAGGGTTTCAGAATCAAACCCTTGTGCCCCGACGGCGCGGCGGCTGGGTATCGTTTTCGCTACGACGTAGc                                                                                                               <  8:355518/136‑1 (MQ=255)
gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGCCGTAGc                                                                                                               >  3:299937/1‑137 (MQ=255)
   aaTTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCACTACCCCGTAGCGaa                                                                                                            >  3:84121/1‑137 (MQ=255)
        ccGGACAAGCCGTATTTATTGGGGGGAATTATATCGTTACCGGGGCCTCTGAGAGGAAACGGGTTTGAGCATCAACCCCTTGTGCACCGACGGCGCGGCGGCTGTGTAGCGTTTGCGCTACGACGTAGCTAAAACGa                                                                                                       <  8:274179/136‑1 (MQ=255)
             taaGCCTTATTTTTTGTGTGTAATGATATCGTTCCTGGGGCCTCTGAGAGTAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGCGGCGGGTGGGTAGCGTTTTCGCTACGACGTCGCGAAAACGATACAt                                                                                                  <  6:208031/136‑1 (MQ=255)
                           tgtgtgGAATGATATCTTTCCTGGGGCCTCTGAGAGGAAAAGGGTTGGAGCATCAACCCTTTGTGCCCCGACGGGTAGGCGCCGGGTTAGCTTTTTCGCTACCACGTAGCGAAAACGATACATGTGGTGCCTGGACt                                                                                    <  6:114022/137‑1 (MQ=255)
                           tgtgtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGAAACATGTGGTGCCTGGACt                                                                                    >  6:386155/1‑137 (MQ=255)
                             tgtgGAATCATATCGTTACCGGTCCCTCTGACAGTAAAAGGGTTTGACCATCAAACCCTTGTGCACCGACGGCGAGGCGGCTGTGTAGCGTTTCCGCAACCACGTCGCTAAAACGATACATGTGGTGCCTGGACTgg                                                                                  <  4:81877/137‑1 (MQ=255)
                              gtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGACGCCGTAGCGAAAACGATACATGTGGTGCCTGGTCTggg                                                                                 >  1:406602/1‑136 (MQ=255)
                               tgtaaTGATATCGTTCCAGGGGCCTCTGCGAGGAAAAGGGTTTCAGCATCAACCCCTTGTGCACCGACGGCGTGGAGGCGGTGTATCGTTTTCGCTACGACGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCg                                                                                <  4:299937/134‑1 (MQ=255)
                               tgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTAAGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCg                                                                                >  6:172469/1‑137 (MQ=255)
                                ggAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGGGCAAACGATACATGTGGTGGCTGGTCTGGCgg                                                                               >  3:112431/1‑136 (MQ=255)
                                  aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAgg                                                                             >  5:175731/1‑137 (MQ=255)
                                  aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGGAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAgg                                                                             >  7:477802/1‑137 (MQ=255)
                                  aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAgg                                                                             >  8:100267/1‑137 (MQ=255)
                                  aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGACGACGTAGCGAAAACGATACATGTGGTGCCGGGCCTGGCGAgg                                                                             >  3:42697/1‑137 (MQ=255)
                                  aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTACACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTAAGACGTAGCGTAAACGAAAAATGTGGTGGAGGGTCTGGCGAgg                                                                             >  5:326734/1‑137 (MQ=255)
                                    ttaaatCTTTACTGGTGCCTATCATAGAAAACGGGTTGGAGCATAAAACCCTTGTCAACGCACGGCGCGGGGGGGGTTTATCGTTTTCGCTACCCCGTCGATAAAACGATACATGTGGTGCCTGGACTGGCGAGGTg                                                                           <  2:98310/133‑1 (MQ=255)
                                     gATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGATACAAGTGGTGCCTGGACTGGCGAGGTgg                                                                          >  7:415672/1‑136 (MQ=255)
                                       tatCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGATACATGTGGTGCCAGGACTGGCGAGGTGCCa                                                                        >  7:396916/1‑137 (MQ=255)
                                         gcGTTAGTGGTGCCGCTGCGAGTAAAAGGGTTTGTGCATAAAACCTTGGTGCCCCGAGGGCGAGTCGGCTTTGAATCTTTTTCGCTACGCCGTTGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAAt                                                                      <  6:184831/136‑1 (MQ=255)
                                                    gCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGCGGCGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTc                                                           <  8:392465/137‑1 (MQ=255)
                                                     ccTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGGAGCGAAAAAGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCGTAGGtaa                                                          >  1:31672/1‑135 (MQ=255)
                                                       tctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGACGAACTATCGAAAAACAAACAATTTGTGCCGGGCCTGGCGACGTGCCAAATGCCTTAGGtaaaa                                                        >  8:77029/1‑133 (MQ=255)
                                                        ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGGGACCGCGGTGCGAAAACGCTACATGTGGTTGCTGGACTTGCGAGGTGCCAATTTGCTTCGGTAAtcg                                                       >  7:52794/1‑136 (MQ=255)
                                                             tgtaaatGGTTTTGAGCATCAAACCCTTGTGCACCGAGGGCGGGGAGGGTGTGTATGGTTTGCCCTACGACTTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACtt                                                  <  7:100267/134‑1 (MQ=255)
                                                                aaaaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACCTACCGAAAAACAATCATGTGGTGCCTGGCCTGGCGACGTGCCAATTGCCTTCGGTGATCACACTTTct                                               >  8:339196/1‑136 (MQ=255)
                                                                    gggTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGGTACGGCGTATCGAAAAAGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCGCACTTTaagcag                                           >  6:196115/1‑137 (MQ=255)
                                                                      gTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGCCGTAGCGAAAACGATAAATGTGGGGGCTGGACTGGCGAGGTGGCAATGGGCGTCGGTGATGTCACTTTCAGCAgat                                         >  3:388362/1‑137 (MQ=255)
                                                                          gAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGATACATTTTGTGCCTGGACTGGCGAGGTGCCCATGGCCGTCGGTCATCTCACTTTGAGCAGATTATc                                     >  5:264914/1‑137 (MQ=255)
                                                                             cATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGACGAAGTAGCGAAAACGATACAATTGGTGCCTGGACTGGCGTGGTGCCAATGGCCTTCGGTCATATCACTTTCAGCAGATGATCCac                                  >  7:391987/1‑136 (MQ=255)
                                                                                  aaCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTAGGTCATCTCACTTTCATCAGATGATCCAGCGaaa                             >  6:229438/1‑137 (MQ=255)
                                                                                  aaCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACCACGTAGCGAAAACGATACATTTGGTGCCTGGACTGGCGAGGTGACAATTGCCTTAGGTCATCTCACATTCACCAGATGATCCACCGaaa                             >  8:288176/1‑137 (MQ=255)
                                                                                   aCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGTGGTTCCAATGGCCTTCGGTCCTCTCACTTTCAGCAGATGATCCAGCGaaaa                            >  8:489614/1‑137 (MQ=255)
                                                                                   aCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGATAAATGTGGTGCCGGGACTGGCGAGGTTCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGaaaa                            >  7:469171/1‑137 (MQ=255)
                                                                                            cACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGATAAATGTGGTTCCTGGACTGGCGCGGTGCCAATGGCCTTGCGTAATGTCACTTTAAGCAGATGATCCAGCGAAAAAGCGCCGgc                   >  8:373734/1‑137 (MQ=255)
                                                                                              ccGACGGCGCGGCGGCTGTGTAGCTTTTTCGCACGGACGTCGTCAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGgcgc                 <  2:406602/137‑1 (MQ=255)
                                                                                               cGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGAAACATGGGGTGCCGGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCCGCGAAAACGCGCCGGCGcc                >  8:124362/1‑137 (MQ=255)
                                                                                                 aCGGCGCGGAGGCTTTGATAGGTTATCTCTAGCATGAAGCGAAAACGATAAATGTGGTGCCTGTACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTc              <  6:409012/137‑1 (MQ=255)
                                                                                                     cgagctggCTTTGTATCTTTTCCGATCCCAGTTCGCTAAACCGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGc          <  8:78625/131‑1 (MQ=255)
                                                                                                       aggaggCTGTGTATCGTTTTTGCTTACCCGTTTGGAAAACCCAACATGTGGTGGCTGGGCAGGGGAGGTTGCAAAGGCCTTCGGGCATATCACATTGAGCAGATTATACAGCGCAAACCCCCCCGCCCCTCGGGCaa        >  8:173988/1‑137 (MQ=255)
                                                                                                            gCTGTGTATCGTTTGCGATACCACTAAGCTAAAACGATAAATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAgg   <  6:175731/137‑1 (MQ=255)
                                                                                                             gggtttATCGTTTGCGCTACCACGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGa  <  7:408266/132‑1 (MQ=255)
                                                                                                                           |                                                                                                                          
GTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGA  >  KT2440/1123115‑1123360

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: ‑