Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,938,2390AC52.2% 10.0 / 20.2 23K537N (AAA→AACpqiBIntermembrane transport protein PqiB
Reads supporting (aligned to +/- strand):  ref base A (6/5);  new base C (0/12);  total (6/17)
Fisher's exact test for biased strand distribution p-value = 4.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-03
Rejected as consensus: Frequency below/above cutoff threshold.

ATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTACCACACGTTGGTAC  >  KT2440/2938105‑2938371
                                                                                                                                      |                                                                                                                                    
aTTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGGTTCGTGGCCGGCCGcaa                                                                                                                                      >  1:112601/1‑135 (MQ=255)
aTTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGaaa                                                                                                                                      >  3:74855/1‑135 (MQ=255)
             cgccgcAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTcc                                                                                                                          <  8:299083/134‑1 (MQ=255)
                 gcAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTGGTGGCCGGCCGCAAGCGGCCCGGGGCGCAGa                                                                                                                     >  5:332982/1‑135 (MQ=255)
                  ccgcgcGAGCTATGGGGTGAAACGCGGGCCGAGCTGGAGCGCTTGTCGCGTGCGCTACGTGACCTTGCCCATTCCCGGGCCCCCCACCCCGAGGCGCTGATGCGGGCCCGGCCGAACGCGGCCGGGTCCGCAGAc                                                                                                                    <  2:112601/133‑1 (MQ=255)
                                      aaCGCGGGACGAGCTGGAGCGCAGGTGGCGTTCGATACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGGGGCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACa                                                                                                 <  8:240097/134‑1 (MQ=255)
                                               gttgCTGGAGCTAATGTCGGGTGCGCTACGGTACCTGGCCTATTACCTGGGCCCCCCCCCCGTGTCGTTGATGCGGGGCCGGCCCAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCaa                                                                                       <  6:332982/132‑1 (MQ=255)
                                               cGAGCTGGCGCGCAGGGCGCGTGCGCTCCGTGACCTTGCCGATACCCGGGCCCGCCACCCCGAGGCGCTGATTCGTGCCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCaa                                                                                       <  7:284450/135‑1 (MQ=255)
                                                    gggAGCGCTTGTCGCGTTCGCTACGTGACGTTGCCGATTACCGGGGCCGCCACCCGGAGTCGGTGATTCGGGCCCGGCCTAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACt                                                                                  <  3:185805/134‑1 (MQ=255)
                                                                   gtgcGCTCCGTCACCTGGCCGATTCCCGGGGCCGCCACCCCGAGCCGCTGATGCGGGGCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCAtgtgt                                                                   <  4:155291/132‑1 (MQ=255)
                                                                     tCGCTACGGCACCTTGCCGATAACCGGGGCCCCCCCCCCGAGTCGCTCATGCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTgg                                                                 <  8:201443/135‑1 (MQ=255)
                                                                            gttcccTTGCCGATAACAGGGCCCCCCACCCCGCTGCGCTGATGGGGGGCCGGCCTAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGc                                                            <  8:52597/129‑1 (MQ=255)
                                                                            gggACCTTGCCGATTCCCTGGGCCGCCACCCCGCGGCGGTGTTGCGGGGCCGGCCTAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGc                                                            <  8:46519/131‑1 (MQ=255)
                                                                               aCCTGGCCGATTCCCGGGGCCGCCCCCCGGAGTCGCTGATGGGGGGCGGCCCTAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCctgct                                                       <  3:236601/135‑1 (MQ=255)
                                                                                      cGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGCCCTTGCACCCTGCGGAGACAGCCGATGCAACCACTGCGCCATGCATGTGTGGCCCGCCGGCTGGTGTgg                                                >  4:35882/1‑135 (MQ=255)
                                                                                      cGATTACCTGGCCGCCACCCCGAGTCGCTGATTAGTGGCCCGCCGAAAGCGGCCGGGTCCGCACACCTTTAAACCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATTTGTGGCCAGCCGGCTGGTGTGGc                                               >  7:242840/1‑135 (MQ=255)
                                                                                              gggcccgccCCCCCGCGGCGCTGATGGGGGGCCGCCCCAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCggggg                                        <  5:179573/131‑1 (MQ=255)
                                                                                                  ccgccACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCt                                      >  7:279486/1‑133 (MQ=255)
                                                                                                   cgccACCCCGAGTCGCGGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCCGCCCGTGTACCCGGTGGAGACACCCGATGGCACGCCTGCGGCATGGATGTGTGGCCCGCCGGCTGGTGTGGGGCggggggggca                                   >  6:103622/1‑130 (MQ=255)
                                                                                                         cccGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGCTGCAACGACTGCGCAATGCATGTGTGGCCCGCCTGCTGGTGTGGCTCGGGGGCTGCCGCAGTa                             >  3:344311/1‑135 (MQ=255)
                                                                                                           ggAGTCGCTGTTTCGGGGCGGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACg                           <  4:74855/134‑1 (MQ=255)
                                                                                                              gTCGCGGTTTCGTGCCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCTACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCg                        <  7:233563/135‑1 (MQ=255)
                                                                                                                      tgcggggccggccGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTa                <  7:330465/130‑1 (MQ=255)
                                                                                                                        cggggcccgccAAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTAcc              <  7:140874/132‑1 (MQ=255)
                                                                                                                         ggggcgggccGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTAcc              <  4:19915/132‑1 (MQ=255)
                                                                                                                                    aaaGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCggggg                                        <  4:31351/97‑1 (MQ=255)
                                                                                                                                    aaaGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCggggg                                        >  3:31351/1‑97 (MQ=255)
                                                                                                                                    aaaGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTCCGCCGCACAACTACCACACGTTGGTAc  >  8:67537/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                    
ATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTACCACACGTTGGTAC  >  KT2440/2938105‑2938371

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