Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,147,1330TG56.7% ‑5.3 / 11.2 30intergenic (+187/+30)ybdG_2/bioPMiniconductance mechanosensitive channel YbdG/Biotin transporter
Reads supporting (aligned to +/- strand):  ref base T (1/12);  new base G (4/13);  total (5/25)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTC  >  KT2440/5147002‑5147267
                                                                                                                                   |                                                                                                                                      
cGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCGGCCCt                                                                                                                                   >  8:478105/1‑137 (MQ=255)
          cTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGGGGGCCCGCCGcccccc                                                                                                                         >  8:1792/1‑136 (MQ=17)
                 ggttATCCGCAGCGGAGTGGAAATGACATGGGGGTGGGGGTGGCCGGGCGCCCGGCGGCAGAGCCGATGCGGCGTGGCATTGAAGATATGTGTGAGGCTTGATATTGCGGGGGCGGCCCTGCAGCCCCCGGATCaaa                                                                                                                  <  4:170260/134‑1 (MQ=255)
                 cTGTATCAGCAGCCGCGGGCCAATGACATGCGCGTGGGGTGGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGGCGCATTGAAGATATGAGTGAGGCTGGAGATTGCCGGGGCTGCCCGGCAGCCCCCGGATCaaa                                                                                                                  <  6:76333/137‑1 (MQ=18)
                      acagctgCCGTGTGGAAATGACTGGCGGGGGGGGTGGGCCGGGCGCACGGCGGCAGAGCCGATGCGGCGTCGATTTCAAGATATGAGGGTGGGTTGTGATTGCGGGGGCGGCCCTGCAGCCCCCGGATCAAAGGGCg                                                                                                             <  6:13753/136‑1 (MQ=255)
                       cagcagCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGGGGCGCATCAACCCCCGGATCAAACGGCgg                                                                                                            >  5:91581/1‑136 (MQ=18)
                             ccGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGGGGCCCTGCCGCCCCCCGATAAAAGGGCGAGTCacg                                                                                                      >  1:3466/1‑135 (MQ=18)
                                         gACATGCGCGTGGGGGTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCGGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCa                                                                                          <  6:474262/137‑1 (MQ=255)
                                          aCATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCGGCAGCCCCCCGCTCAAAGGGCGAGGCCGGCCTTGCCGGcccc                                                                                         >  6:91130/1‑135 (MQ=21)
                                          aCATGCGCGTGGGGTGGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGGCGCATTGAAGATATGAGTGCGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAg                                                                                         <  4:432812/137‑1 (MQ=255)
                                             tGCGCGGGGGGTGGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGCGTGAGGCTTGAGATTGCCGGGGCGGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAgcgg                                                                                      <  5:507088/137‑1 (MQ=255)
                                                  ggggggTGGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGAATTGAACATATGAGTGAGGCTTGAGATTCCCGGGGCGGCCCTGCCCCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCt                                                                                 <  1:154021/135‑1 (MQ=255)
                                                     gggTGGGCCGGCGGCACGGCGGCAGACCCGATGCCGCGTCGCTTTTAACATATGAGTGAGGCTTGAGATGGCCGGGGCGCCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCtgtg                                                                              <  2:180193/137‑1 (MQ=255)
                                                      gtttGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGGCGCATTGAAGATATGAGGGCGGCTTGAGATGGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                             <  2:324711/135‑1 (MQ=255)
                                                      ggTTGGCCGGGCGCAGGGCGGCAGAGCGGTGGCGGCGTCGCATTGAAGATATGAGTGAGGGTTGCGATGGCGGGGGGGGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                             <  4:472905/137‑1 (MQ=255)
                                                      ggTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCGGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                             <  4:66139/137‑1 (MQ=255)
                                                      ggTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGGCGCATTGAAGATATGAGTGCGGCTTGAGATTGCCGGGGCGGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                             <  5:139453/137‑1 (MQ=255)
                                                      ggTTGGCCGGGCGCACGGCGGAAGAGCCGATCCCGCGTCGCATTGAATATATGAGTGCGGCTTGAGATTGCCGGGGCTGCCCGGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                             <  4:454992/137‑1 (MQ=38)
                                                      ggTGGGGCGGGGGCAGTGGGGCATCGCGGATGCGGGGTCGCTTTTAATATTTGAGTGCGGCTTGTGATTGCCGGGGCGGCCCTCAAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                             <  4:164131/137‑1 (MQ=255)
                                                      ggTGGGCGGGACGCACGGCGGCAGAGCCTATGCCGCGTCGCATTGAATATTTGAGTGAGGCTTGAGATGGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                             <  1:91791/137‑1 (MQ=255)
                                                                 cgcACGGCGGCAGCGCCGATGACGCGTCTCATTGAAGATATCAGGGAGGCTTGCGATTCCCGGGGCTGCCCTGCCCCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGcccaccca                                                                  <  7:219717/137‑1 (MQ=255)
                                                                          gcAGAGCCGATGCGGCGTCGATTTCAAGATATGAGGGCGGCTGGCGACGGCCGGGGCGGCCCTGCAGCCCCCGGATCAAAGGGCGAGGCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCa                                                         <  3:47363/137‑1 (MQ=255)
                                                                               gCCGTTGCCGCGTCGCATTTAAGATTTGAGTGAGGCTTTGGATTGCCGGGGCGGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCcac                                                    <  4:145510/137‑1 (MQ=255)
                                                                               gCCGATGCCGCGTCTCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCcac                                                    <  4:450888/137‑1 (MQ=255)
                                                                               gCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGACTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCcac                                                    <  8:117640/137‑1 (MQ=255)
                                                                               gCCGATGCCGCGGCGCATTGAAGATTTGAGTGAGTCTTGCGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCcac                                                    <  8:257717/137‑1 (MQ=255)
                                                                                          gTCGCATTGAAGATATGAGTGAGGCTTGAGGTTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAg                                         <  8:123065/137‑1 (MQ=255)
                                                                                                                gtgttGAGATTGCCGGGGCGGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAAc                   <  7:420021/134‑1 (MQ=255)
                                                                                                                       taTTGCGGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAgccgcc            <  5:176356/136‑1 (MQ=255)
                                                                                                                                 gCTGCCCTGCACCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTc  <  3:62150/137‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                      
CGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTC  >  KT2440/5147002‑5147267

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