Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,261,3510AG45.5% 11.1 / 11.4 25F51S (TTC→TCC) nicP_10Porin‑like protein NicP
Reads supporting (aligned to +/- strand):  ref base A (4/8);  new base G (10/0);  total (17/8)
Fisher's exact test for biased strand distribution p-value = 1.69e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGCTTCACGCCGAACAGGGCGATACCATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGC  >  KT2440/4261217‑4261478
                                                                                                                                      |                                                                                                                               
ccAGCTTCACTCCGAACAGGGCGATACCATCCAGGCCAAACCCGCTCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGa                                                                                                                                 <  6:154525/135‑1 (MQ=255)
        ccGCCGAACAGGGCGATACCATCCCGGCCAAAACCGATCCGCCCCGGGGTGTCGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATcgcg                                                                                                                          <  5:134244/133‑1 (MQ=255)
        aCGCCGAACAGGGCGATACCATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATcgc                                                                                                                           <  4:42156/133‑1 (MQ=255)
        aCGCCGAACAGGGCGATACCATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATcgc                                                                                                                           <  4:190071/133‑1 (MQ=255)
                      gATACCATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGGCTGGCGGGTGTagtaatagtg                                                                                                           >  5:117368/1‑134 (MQ=38)
                           cATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCa                                                                                                      <  6:140446/135‑1 (MQ=255)
                                aGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTAGCCCTTGGCGGCCCCCGGGTCGCGGCAAACGAGGTTGAACTAGAAGTTGCGCCATGtg                                                                                                 >  6:156500/1‑134 (MQ=38)
                                 ggCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCa                                                                                                <  7:102816/135‑1 (MQ=255)
                                          cGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGAATCGCGGTTGAAGTAGTAGTTTCGCAATGTCAGCCCAccca                                                                                       >  6:79509/1‑134 (MQ=255)
                                          cGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATAGCGGTTTAAGGAGTAGTTTCGCAATGTCAGcccccc                                                                                         >  8:129419/1‑133 (MQ=255)
                                             gcAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTgg                                                                                    <  1:100156/134‑1 (MQ=255)
                                               aGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGCTTAGCGGGTGTAGTTGTAGTTTGGGAATGTCGGGCCCCCCGTGGcg                                                                                  >  3:117303/1‑134 (MQ=255)
                                               aGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGCATCGCGGGTGGTGTAGTAGTTTGGGAGTGTCAGCCCCCCCGTGGcc                                                                                  >  1:155429/1‑134 (MQ=255)
                                               aGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGCAATAGCGGTTGAAGTAGTAGTTGCGCAATTTAAGCCCACCCTTTGCa                                                                                  >  8:15516/1‑135 (MQ=255)
                                               aGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAACTCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCAACCTTGGCa                                                                                  >  8:178816/1‑135 (MQ=255)
                                                              ccAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGAATGGCGGTTGAAGTATTATTTGCGCAATGTGATCCCCCCCGTGGGAGCTTCCAAGAAcccc                                                                   >  6:77326/1‑135 (MQ=255)
                                                              ccAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATAGCGGGTGAAGTAGTAGTTTCGCAATGGAAGCCCACCCTTGGCATCATGCAACAAGccc                                                                   >  4:110490/1‑135 (MQ=255)
                                                                     tGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGCAATCGCGGTTGAAGTAGTTGTTGCGCCAATTCCGCCCACCCATGGGAACTTCCAAGAAGccccccgcgt                                                            >  6:156367/1‑131 (MQ=255)
                                                                      gAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGGAATAGCGGTTGTAGTTGTATTTGGCCGATGTAAGCCCCCCCGTTGGGTCGTCCAAAAACccccccgaagg                                                           >  8:81567/1‑130 (MQ=255)
                                                                                                               ccTTGGCGGCCCCCGGGTCGCGGGAATCGCGGGTGAAGTAGTAGTTGGGCAATGTCCGCCCCCCCATGGGATCTTACAAGAAGCCCCCCGCTTGCGCCCACGCCGGCCAACCACGGCCCCACGCCAGCGACACAg                  >  6:102139/1‑135 (MQ=255)
                                                                                                               ccTTGGCGGCCCCCGGGTCGCGGAAATCGCGGGGGAAGTAGTAGTTTCGCAATGTCAGACCCCCCGTGGCATCGTCCAAGACGCCCCCCGATTTCGCCCCCGCCGGCCCACCCACGCCCAACCCCCGCGaacaac                  >  2:120587/1‑131 (MQ=255)
                                                                                                                 ttGGCGGCCCCCGGGTCGCGGGAATCGCGGGTGAAGTAGTTGTTGCGCAATGTGAGCCCCCCCATGGCCTCTTCCACGAAGCCCCCCTCATGCGCCCACGCCGCCCAACCACGGCCCAACCCCACCGAAACAccc                >  4:141678/1‑135 (MQ=255)
                                                                                                                 ttGGCGGCCCCCGGGTCGCGGGAATCGCGGGGGGAGGTGTAGTTGGGCAATGTCAGCCCCCCCGTGGCCTCTTCCCAGAAGCCCCCCGCCTGCGCCCCGGCCGGCCCACCCCCGCCCAACGCCCCCGCAACAccc                >  4:200011/1‑135 (MQ=255)
                                                                                                                      ggcCCCCGGGTCGCGGGAATCGCGGGTGAAGTAGTAGTTGCGCAATGTCAGCCCCCCCATGGCATCTTACCAGAAGCCCCCCACTTTCGCCCCGGCCCGCCAACCCACGCCCAACGCCCCCGCAACCGCCacgcc           >  2:129143/1‑131 (MQ=255)
                                                                                                                               gTCGGGTAAATCGCGTTTTATGTAGTTGTTGCGAAATGTCCGCCCACCCTTGGCATATTACAATAAGCCCCCCTCTTGCGCCCTGGCGGGCAACCAACGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGc  <  6:117368/135‑1 (MQ=255)
                                                                                                                               gTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGCCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGc  <  2:155429/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                               
CCAGCTTCACGCCGAACAGGGCGATACCATCCAGGCCAAAACCGATCAGCCCCGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGC  >  KT2440/4261217‑4261478

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: