Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,834,9220GT46.2% 1.1 / 17.5 26intergenic (‑44/+208)yqhD/CJ019_02688Alcohol dehydrogenase YqhD/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (6/8);  new base T (12/0);  total (18/8)
Fisher's exact test for biased strand distribution p-value = 2.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGCTTCAGTCGGAATCAGCTCATCCAGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGA  >  ACB122/2834788‑2835046
                                                                                                                                      |                                                                                                                            
gggCTTCAGTCGGAATCAGCTCATCCAGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTg                                                                                                                              >  2:9298/1‑135 (MQ=255)
        gTCGGAATCAGCTCATCCAGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTGGATTTTTTTGTGaaaat                                                                                                                      >  5:176454/1‑133 (MQ=255)
          cGGAATCAGCTCATCCAGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTGTATTTTTTGGTGAAAtttt                                                                                                                    >  5:55105/1‑135 (MQ=255)
              aTCAGCTCATCCAGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTGATTTTTTTTGTGGAATTTTTAtt                                                                                                                >  4:197238/1‑135 (MQ=255)
                              gggCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTTTATTTTTTTTTGTAATTTTTTTTTCAGAAGCTCGTAggg                                                                                                >  1:164509/1‑134 (MQ=255)
                              gggCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGgc                                                                                                >  8:161885/1‑135 (MQ=255)
                                   aTGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTTTATTTTTTTGTGGAATTTTTTTTTCTTGAGCTCGTGGgcgcgcc                                                                                           >  2:133775/1‑134 (MQ=255)
                                   aTGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTGTTTTTTTTTGTGTAATTTTTATTTCAGAAGCTAGTCGgcgcgcc                                                                                           >  3:60745/1‑134 (MQ=255)
                                      gtttggttACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCAtt                                                                                         <  2:54102/134‑1 (MQ=255)
                                           gttACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCa                                                                                    <  7:71750/134‑1 (MQ=255)
                                                       aaTTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGTGCACCAAATAGCTGaa                                                                       <  7:161885/135‑1 (MQ=255)
                                                       aaTTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGa                                                                        <  7:23590/134‑1 (MQ=255)
                                                           cgcgTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTGTATTTTTTGTTTTAATTTTTTTTGGATAAGATCGTCGGCGCGCGGTGGGCACCCAATAGCTgtatctt                                                                   >  5:3461/1‑129 (MQ=255)
                                                            gcgTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTTTATTTTTTTGTGTAATTTTTATTTGAGAAGCTCGTCGGCGCGCATTGCGCACCAAAAAGCAGAAGATTg                                                                  >  5:63546/1‑135 (MQ=255)
                                                                ggggTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTTGTGTAATTTTTTTTTCAGAAGCTCGTCGGCGCGCCTTTCGCACCCAATAGCTGAAGATTTGtt                                                               >  1:38392/1‑134 (MQ=255)
                                                                   ggtgtgGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGaa                                                           <  2:30280/133‑1 (MQ=255)
                                                                          aaGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTTTATTTTTTTGTGTAAATTTTTTTTCAGAAGCTAGTCGGCGCCCATTGAGCACCAAAAAGCAGAATATTTGTTATAAgcttgc                                                     >  7:118881/1‑134 (MQ=255)
                                                                             tCAAAATTGAGCATAGATTTCTCTCCTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCt                                                 <  3:7768/135‑1 (MQ=255)
                                                                                   ttGAGCATAGATTTCTATACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGc                                           <  1:183795/135‑1 (MQ=255)
                                                                                         aTAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTTGTGTAATTTTTTTTTGAGAAGCTCGTCGGCGCGCATTTCGCCCCAAAAAGCTGAATATTTGTTTTAATCTTTGTTATAAAAGCaaa                                        >  5:110853/1‑131 (MQ=255)
                                                                                         aTAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTTCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCATGCTTCTAAAGGAAACCaa                                     >  4:151319/1‑135 (MQ=255)
                                                                                             aTTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTg                                 <  2:146367/135‑1 (MQ=255)
                                                                                                          ttAAAAGTAGACCAGTCGTATATTTTTTTGTGAAATTTTTATTTGAGAAGCGAGTCGGCGCGCATTGCGCCCCAAATAGCTGAAGATTGGTTTGATGCTTGCTTGTAAAGGCAACCAAGTTGGTGTGGCATTCGt                    >  2:112199/1‑135 (MQ=255)
                                                                                                          ttAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTTTTGGAGAAGCTCGTAGGCGCGCATTGCGCACCAAATATCTGAAGATTGGTTTGTAGCTTGCTTGTAAAGGCAACCAACTTGGTGTGGCCTTCGt                    >  4:41218/1‑135 (MQ=255)
                                                                                                                      cAGTCGTATATTTTTTTGTGTAATTTTTATTGCAGAAGCTCGTCGGCGCGCCTTGCGCCCCAAAAACCTGAAGATTGGTTAGAAGCTTTCGTCTAAGGGCAACCAACGTGGGTTGGCATTCGTCTGCAATCCTgg        >  6:97626/1‑135 (MQ=255)
                                                                                                                           gtttatTTTTTGGTGAAATTTTTATTGGAGAAGCTCGTCGGCGCCCATTGCGCACCAAATAGCTGAAGATTTGTTTTTAGCTTGCTTCTAAAGGCCACCCAGTTGGGGTGGGCTTCGGCTGGATTCCGGGTACCg   >  4:204168/4‑135 (MQ=255)
                                                                                                                            ttttttttttTGTGTAATTTTTTTTTGAGAAGCTCGTCGGCGCGCCTTTCGCCCCAAATAGATGTAGATTTGTTTGAATCTTTGTTATAAATGCAACCAAGTTGGTGTGGGATTGGTCTTGCTTCCGGGTCCCgc  >  2:133860/5‑134 (MQ=255)
                                                                                                                                      |                                                                                                                            
GGGCTTCAGTCGGAATCAGCTCATCCAGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGA  >  ACB122/2834788‑2835046

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: