Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,227,1560CG31.6% 68.8 / 12.2 38A59G (GCG→GGG) fecI_10putative RNA polymerase sigma factor FecI
Reads supporting (aligned to +/- strand):  ref base C (6/20);  new base G (12/0);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 6.86e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCACAACTGGCTCACCGGCTGGTTGCGTCGCCGTCTGGGCTGCCCGCACAGCGCCGCCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACTGGAGCGTGCCTACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCG  >  KT2440/5227027‑5227283
                                                                                                                                 |                                                                                                                                
accacAACTGGCTCACCGGCTGGTTGCGTCGCCGTCTGGGCTGCCCGCACAGCGCCGCCGATCTGGCCCAGGCCACCTACCTGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGGTTc                                                                                                                             >  4:209217/1‑135 (MQ=255)
      aCTGGCTCACCGGCTGGTTGCGTCGCCGTCTGGGCTGCCCGCACAGCGCCGCCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCa                                                                                                                       <  8:159032/135‑1 (MQ=255)
        tGGCTCACCGGCTGGTTGCGTCGCCGTCTGGGCTGCCCGCACAGCGCCGCCGATCTGGCCCAGGCCACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGGGCGGTTGATGCcccc                                                                                                                     >  4:84673/1‑135 (MQ=255)
        tGGCTCACCGGCTGGTTGCGTCGCCGTCTGGGCTGCCCGCACAGCGCCGCCGATCTGGCCCAGGACACCTACAAGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGGGTTTGGTGCCAcc                                                                                                                     >  6:75566/1‑135 (MQ=255)
          gcgcACCGGCGGGTTGCGGCGCCGTCTGGGCTGCCCGCACCGCGCCGCCGATCTGGCCCAGGCCACCTACATGCGCCTGCTCCAGGGGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGt                                                                                                                   <  3:187790/132‑1 (MQ=255)
          gCTCACCGGCTGGTTGCGTCGCCGTCTGGGCTGCCCGCACAGCGCCGCCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGGTTCTTGCcccc                                                                                                                     >  8:47592/1‑133 (MQ=255)
            ccACCGGCTGGTGGCGTCGCGGTCTGGGCTGCCCGCACCGCGCCGCGGATCTGGCCCAGTACCCCTACATGCGCCTGCGCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAg                                                                                                                 <  3:294111/134‑1 (MQ=255)
                  ggtGGTTGCGGCGCGGGCTGGGCTGCCCGCACAGCGCCTCCGATCTGGCCCAGGACACCTACCTGCGCCGGATCCAGGCGCGCTAGGCGCCCCAGTTGGTGGAGCCGCGGGCGTTTCTTGCCACCGTAGCCAAAc                                                                                                           <  2:48141/133‑1 (MQ=255)
                   cTGGTTGCGTCGCCGTCTGGGCTGCCCGCACAGCGCCGCCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAAcg                                                                                                          <  7:275239/135‑1 (MQ=255)
                    tGGTTGCGTCGCCGTCTGGGCTGCCCGCACAGCGCCGCCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAAcgc                                                                                                         <  2:252830/135‑1 (MQ=255)
                         gCGTCGCCGTCTGGGCTGCCCGCACAGCGCCGCCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGGTGGTGGAGCCGCGGGGGGTTCTTGCCCCAGGCGCCAAAcgcggttt                                                                                                    >  4:187602/1‑131 (MQ=255)
                                          gCCCGCACCGCGCCGCCGAGCTGGCCCAGGCCACCTACAGGCCCCGGCGCCAGGGGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTccgcc                                                                                  <  5:166885/135‑1 (MQ=255)
                                          gCCCGCACAGCGCCGCCGATCGGGCCCCGTACACCTACATGGGCCGGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTccgcc                                                                                  <  5:75566/135‑1 (MQ=255)
                                          gCCCGCACAGCGCCGCCGAGCGGGCCCTGGACACCTACATGCGCCTGCTGCAGGCGCGCGAGGCGCCCCAGGTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTccgcc                                                                                  <  2:22069/135‑1 (MQ=255)
                                            ccGCACAGCGCCGCCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTccgcc                                                                                  <  7:88153/133‑1 (MQ=255)
                                            ccGCACAGCGCCGCCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGGGGGGTTCTTGGCACCGGTGCCAAACGCGGA‑TTGGGTTACCAATTCcgccgc                                                                                >  5:52224/1‑135 (MQ=255)
                                             cGCACAGCGCCGCCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGGGGGGTTGTTGGCCCCCGAGCCCCACCCGGATTTTTGTAACC‑CATCCgccgcc                                                                               >  2:281754/1‑135 (MQ=255)
                                              gCACAGCGCCGCCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGGGCGGTTCTTGCCCCCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCa                                                                              >  8:88153/1‑135 (MQ=255)
                                                    cgccgcCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACt                                                                        <  3:287808/135‑1 (MQ=255)
                                                    cgccgcCGATCGGCCCCAGGACACCACAAGGCGCCTGCTCCCGCCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACt                                                                        <  6:203239/135‑1 (MQ=255)
                                                               tGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGGGGGTGTGGTGCCCCCGGAGCCCCACCCGTA‑TTTTGGTTCCCCATCCGCCGCCCGGCACTGGCGCGTGccg                                                             >  4:18628/1‑134 (MQ=255)
                                                                   cccgtacacCTACATGCGCGTGCTCCCGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACTGGAGCGTGCCTACCt                                                         <  6:15895/130‑1 (MQ=255)
                                                                   cccggACACCTCCATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACTGGAGCGTGCCTACCt                                                         <  6:25068/132‑1 (MQ=255)
                                                                         cacCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGGGGGGTTGTTGGCCCCGTAGCCCAACCCGTT‑TTGTGTTAACACATCCCCCGCCACGAACAGGCGCGTGCCGACCAGCAAgc                                                   >  5:289548/1‑135 (MQ=255)
                                                                          acCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCCCCGGAGCCAAACGCGTA‑TTGTGTAACCCCATCCGCCGCCAGGGACTGGAGCGTGCCTACCTGCCAgcg                                                  >  5:117375/1‑135 (MQ=255)
                                                                              aCATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGGGGGTGTTGCCCCCGTCGCCAACCGCGTA‑TTGTGTTACCACTTCCGCCGCCCGGAACTGGGGCGGGCCTACCTGCAAGCGCTCg                                              >  1:173455/1‑135 (MQ=255)
                                                                               cATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGGGGTTCTTTCCCCCGGAGGCCAACGCGGT‑TTGTGTTACCACGTCCGCCGCCCGGCACGGGGGCGTGGCGACCGGCCAGCGCTCGc                                             >  3:47974/1‑135 (MQ=255)
                                                                                aTGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACTGGAGCGTGCCTACCTGCAAGCGCTCGc                                             <  3:84673/134‑1 (MQ=255)
                                                                                aTGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACTGGAGCGTGCCTACCTGCAAGCGCTCGc                                             <  2:7025/134‑1 (MQ=255)
                                                                                    gcCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACTGGAGCGTGCCTACCTGCAAGCGCTCGCACAg                                         <  6:140332/134‑1 (MQ=255)
                                                                                    gcCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGGGGGGGTGTTGGCCCCGGTGGCCAACGCGGA‑TTGTGTTTCCACATCCCGCCCCCGGGACTGGGGCGGGGCGCCCCGCCAgcgcgcgcccccg                                        >  6:291728/1‑130 (MQ=255)
                                                                                    gcCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGGGCGGATCTTGCCACCGGAGGCAAACGCGGA‑TTGTGTAACCACGTCCGCCGCCAGGAACAGGCGCGGGCCGACCTGCCACCgcgcgccccgg                                        >  1:213587/1‑130 (MQ=255)
                                                                                         cTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGGGTTTGTTGGCACCGGAGGCCAACCCGGA‑TTGTGTTAACACTTACGCCGCCACGAACGGGATCGTGCCTTCCTGCCAGCGCTCGCACAGGGCCCg                                   >  2:45759/1‑135 (MQ=255)
                                                                                          tCCAGGGGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCACACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACTGGAGCGTGCCTACCTGCAAGCGCTCGCACAGGTccc                                    <  7:122697/133‑1 (MQ=255)
                                                                                                       ggCGCCCCAGTTGGTGGAGCCGCGTGGGTTTGGTGCCCCCGGAGCCCAACGCGTA‑TTGTGTTACCCCTTGCGCCGCCAGGAACTGGTGGGTGCCTACCTGCAAGCGCACGCCCAGGTCCCGGAAGAGGGGGCGcc                     >  4:226415/1‑135 (MQ=255)
                                                                                                       ggCGCCCCAGTTGGTGGAGCCGCGTGGGGTTGGTGGCCCCGGAGCCCAACGCGGA‑TTGTGTTACCCCTTCCGCCGCCCGGCACTGGGGCGTGGCGACCTGCCAGCGCGCGCACAGGGCCCGGGAGAGGTCGGGcc                     >  6:299250/1‑135 (MQ=255)
                                                                                                         cgcCCCAGTTGGTGGAGCCGCGTGCGTTGCTTGCCACCGTAGCCAAACGCGTA‑ATGTGTAACCACTTCCGCCGCCACGAACTGGAGCGTGCCTACCTGCAAGCGCTCGCACAGGACCCGGAAGAGGTCGCGCCGt                   <  7:110851/135‑1 (MQ=255)
                                                                                                                          ccGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACTGGAGCGTGCCTACCTGCAAGCGCTCGCACAGGTCCTGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCg  <  7:47592/135‑1 (MQ=255)
                                                                                                                             cgtgcgtTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACTGGAGCGTGCCTACCTGCAAGCGCTc                                               <  2:51122/87‑1 (MQ=255)
                                                                                                                             cgtgcgtTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACTGGAGCGTGCCTACCTGCAAGCGCTc                                               >  1:51122/1‑87 (MQ=255)
                                                                                                                                 |                                                                                                                                
ACCACAACTGGCTCACCGGCTGGTTGCGTCGCCGTCTGGGCTGCCCGCACAGCGCCGCCGATCTGGCCCAGGACACCTACATGCGCCTGCTCCAGGCGCGCGAGGCGCCCCAGTTGGTGGAGCCGCGTGCGTTTCTTGCCACCGTAGCCAAACGCGTA‑TTGTGTAACCACTTCCGCCGCCAGGAACTGGAGCGTGCCTACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCG  >  KT2440/5227027‑5227283

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