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

TTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAAC  >  KT2440/5146998‑5147262
                                                                                                                                       |                                                                                                                                  
ttGCCGGAGTTCGGCTTGCGCCGGTAGCAGCAGCCGAGTGGCAATGAGATGCGGGTGGGGTTGGCCGGGCGCACGGCGGCTGAGCCGATGCCGGGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTg                                                                                                                                   >  1:267125/1‑137 (MQ=17)
          tCGGGTTGAGCCTGTATCAGCAGCCGAGTGGCATTGACATGCGCGTGGGGTTGGCCGGGCGCAGGGGGGAAGAGCCTATGCCGCGTCTCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGcc                                                                                                                         <  4:412047/137‑1 (MQ=17)
          ggggCTTGAGCGTGTATCAGCCCCCGTGTGGCATTGACTTGCGGGTGGGGGTGGCCGGGCGCCCGGGGGCAGAGCCGATGCCGCGGCGAATTGAAGATATGAGGGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGcc                                                                                                                         <  2:381130/135‑1 (MQ=17)
                                           ttGACATGCGCGGGGGGGGGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATCCGTGCGGCTTGCGATTGCCGGGGCTGCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGc                                                                                       <  8:343130/136‑1 (MQ=255)
                                                 ggcggGTGGGGTGGGCCGGGCGCACGGCGGCAGAGCCGATGCCGGGTCGCATTGAATATCTGAGGGCGGCTTGAGATTGCCGGGGCGGCCCGGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAgcgg                                                                                 <  6:379137/136‑1 (MQ=255)
                                                 ggcgcGGGGGGGGGGCCGGGCGCCCGGCGGCAGAGCCGCTGCCGCGGCGAATTGAAGATATGCGTGCGGCTTGCGATTGCCGGGGCGGCCCGGCCGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAgcgg                                                                                 <  6:372386/136‑1 (MQ=255)
                                                          ggggggCCGGGCGCCGGGGGGCAGAGCCGATGCCGCGGCGCATTGAAGATAGGCGTCCGGCTGGCGATGGCCGGGGCGGCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                        <  3:146286/133‑1 (MQ=255)
                                                          ggTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                        <  6:353722/137‑1 (MQ=255)
                                                          ggTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCGGCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                        <  4:63871/137‑1 (MQ=255)
                                                          ggTGGGCCGGGCGCACGGCGCCAGAGCCGATGCCGCGGCGCATTGAAGATATGAGTGCGGCTTGAGATGGCCGGGGCGGCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                        <  7:127660/137‑1 (MQ=255)
                                                            ttGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCGGGCCTGCCGCCCC‑CGGCTGAAAGGGCGCGTCAGGCCTTGGCGGCCCGCGGCGCCGGTgggc                                                                      >  5:432638/1‑137 (MQ=255)
                                                            ggggCCGGGCGCACGGGGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgcgc                                                                      <  8:122272/135‑1 (MQ=255)
                                                                 cGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGGTGGCCGGCCGCCCC‑CCGGGCAAACGGGGGGGGAGGGCTTGCCGGCCCGCGGCGCCGGTGGGCGcccc                                                                 >  4:52461/1‑136 (MQ=255)
                                                                 cGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGGGGCCCAGCCCCCCC‑CCGAGTAAAAGGCGCGGACGGCCATGCCGGCCCCCGGCGCCAGTgcccgcccc                                                                 >  6:93925/1‑136 (MQ=255)
                                                                    gcgcACGGCGGCATAGCCGATGCCGCGGCGCATTGAATATATGAGTGAGGCTTGAGATGGCGGGGGCGGCCCTGCCGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGcccaccc                                                              <  7:179585/137‑1 (MQ=255)
                                                                         cggcgtCAGCGCCGTTGCCGCGTCTCATTGATGTTATGAGGGAGGCTGGTGATTGGCGGGGCTGCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAggcg                                                         <  6:91207/137‑1 (MQ=255)
                                                                                   gcctttGCCGCGTCGAATTTAAGATATGAGTGAGGCTTGAGATGGCGGGGGCGGCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCcac                                               <  3:81906/132‑1 (MQ=255)
                                                                                                cGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTTCCCTGCCGCCCC‑CCGATGAAAGGGGGCGTAACGCCTTTCCCGCCACCGGCCCCAGTGCGGGCCCCCCCCGGCGGGGCAGCCCCACCCCCGCCAggg                                  >  8:356474/1‑136 (MQ=255)
                                                                                                      gAAGATATGAGTGCGGCTTGCGATTCCCGGGGCGGCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATcc                            <  8:71771/137‑1 (MQ=255)
                                                                                                        aGATATCAGTCAGGCTAGAGATTGCGGGGGCTCCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATcccc                          <  8:399358/137‑1 (MQ=255)
                                                                                                                                 gggggCGGCCCCGCCCCCCCGCGGATCAAAGGGCGAGGCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAAc  <  2:234348/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
TTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCC‑CGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAAC  >  KT2440/5146998‑5147262

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