Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913238,5270AG46.2% 3.0 / 49.7 39I70I (ATT→ATCykfMlethality reduction protein, putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base G (10/8);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATG  >  NC_000913/238385‑238669
                                                                                                                                              |                                                                                                                                              
atCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACa                                                                                                                                          >  1:116968/1‑149 (MQ=255)
   ctttgACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTTTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGa                                                                                                                                       <  7:100706/146‑1 (MQ=255)
                        aaGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGAc                                                                                                                  <  3:62871/149‑1 (MQ=255)
                           gTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACAt                                                                                                                >  8:105753/1‑148 (MQ=255)
                                    cGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGaa                                                                                                       <  2:116968/148‑1 (MQ=255)
                                                   aaaaTATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACAt                                                                                        >  1:148704/1‑148 (MQ=255)
                                                      atatAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATtctc                                                                                    >  4:30243/1‑149 (MQ=255)
                                                       tataAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTcc                                                                                   <  2:148704/149‑1 (MQ=255)
                                                       tataAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTcc                                                                                   <  7:105753/149‑1 (MQ=255)
                                                                  gCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATaaaa                                                                         >  6:3447/1‑148 (MQ=255)
                                                                      tATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAatat                                                                      <  5:3447/147‑1 (MQ=255)
                                                                           ttATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACAt                                                                                                                >  3:10151/1‑100 (MQ=255)
                                                                           ttATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACAt                                                                                                                <  4:10151/100‑1 (MQ=255)
                                                                                          aaTTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTattatt                                                <  3:30243/149‑1 (MQ=255)
                                                                                              gACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAAAATTTCCGTGGAGGATTTTATTATTGat                                             >  8:133098/1‑147 (MQ=255)
                                                                                                     ttCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACa                                                                                         >  6:39564/1‑97 (MQ=255)
                                                                                                     ttCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACa                                                                                         <  5:39564/97‑1 (MQ=255)
                                                                                                     ttCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTattatt                                                <  3:68926/138‑1 (MQ=255)
                                                                                                     ttCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTattatt                                                >  4:68926/1‑138 (MQ=255)
                                                                                                     ttCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATtctc                                                                                    <  2:16042/102‑1 (MQ=255)
                                                                                                     ttCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATtctc                                                                                    >  1:16042/1‑102 (MQ=255)
                                                                                                      tCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGt                                    >  8:71535/1‑149 (MQ=255)
                                                                                                      tCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTACAATCATAAAATATTTCCGAGGAGCATTTTATTATTGAAAATAGAGGt                                    >  5:153625/1‑149 (MQ=255)
                                                                                                       cGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACCGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGtt                                   >  4:18118/1‑149 (MQ=255)
                                                                                                            ttttGATTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACt                              >  8:84699/1‑149 (MQ=255)
                                                                                                                 aTTAAGAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGt                         <  7:84699/149‑1 (MQ=255)
                                                                                                                  ttAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAata                                           >  7:31076/1‑130 (MQ=255)
                                                                                                                  ttAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAata                                           <  8:31076/130‑1 (MQ=255)
                                                                                                                      gAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATagag                                      <  7:118498/131‑1 (MQ=255)
                                                                                                                      gAATGATTTTATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATagag                                      >  8:118498/1‑131 (MQ=255)
                                                                                                                            ttttATTATCGTAAGTAAGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACa                 >  2:110635/1‑146 (MQ=255)
                                                                                                                               tattatCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAaagaag           >  4:58879/1‑149 (MQ=255)
                                                                                                                                   atCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCAt        <  3:58879/148‑1 (MQ=255)
                                                                                                                                    tCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCAtt       >  3:9514/1‑148 (MQ=255)
                                                                                                                                      gtaagtaaAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCTTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGaa    <  7:71535/149‑1 (MQ=255)
                                                                                                                                         agtaaGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACCGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATg  <  3:18118/148‑1 (MQ=255)
                                                                                                                                         agtaaGATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATg  >  3:137799/1‑148 (MQ=255)
                                                                                                                                           taaAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAata                                           <  4:31299/105‑1 (MQ=255)
                                                                                                                                           taaAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAata                                           >  3:31299/1‑105 (MQ=255)
                                                                                                                                              |                                                                                                                                              
ATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATG  >  NC_000913/238385‑238669

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