Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,364,3960AT42.1% 55.9 / 15.5 38P530P (CCA→CCT) yhcDputative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base A (18/4);  new base T (0/16);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 2.25e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTG  >  NC_000913/3364250‑3364543
                                                                                                                                                  |                                                                                                                                                   
ggTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTg                                                                                                                                                 <  2:252137/151‑1 (MQ=255)
      aaTGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGc                                                                                                                                               >  2:786740/1‑147 (MQ=255)
         gACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGAGCaa                                                                                                                                        <  2:139683/151‑1 (MQ=255)
               aGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCg                                                                                                                                  <  1:298546/151‑1 (MQ=255)
                   aCCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCcac                                                                                                                              >  1:757470/1‑151 (MQ=255)
                      ttAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACt                                                                                                                           >  1:4594/1‑151 (MQ=255)
                       tAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTg                                                                                                                          >  2:273708/1‑151 (MQ=255)
                                    aaaaaaGCGTAAAACCGCCCCAGGGTTTCCGTAAACTGGAAGCATCACCTTAGTCCTGCAAAAAATAAGAAAGACGATGGGAATCGGTTCTACTTCAACAACAGTTTTCTTTTTGGCAAAAAAAACACCGCCACACTGTATACTCGCCATg                                                                                                             <  2:182299/148‑1 (MQ=255)
                                       aaaGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACg                                                                                                          >  1:441384/1‑151 (MQ=255)
                                           cGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATc                                                                                                       >  1:339671/1‑150 (MQ=255)
                                                aaaCACGCCACGTTTTCCGTAAACTGGCGGCATAACCTTATTGCTTAAAAAAATAATGAAGCAGATGGTGATCTGTTCTACTAAAACATCAGTATTCCTTTTGGAAAAAAAACCACCGCCACACTGTATACTCGCCATGACGATCATaaaa                                                                                                 <  2:91173/151‑1 (MQ=255)
                                                aaaCACGCAGCGGTTTCCGTCAACGGGAAGCATCAGCTTGGTGCTTAAAAAAATAATGAAGACGATGGTGTTCTGTTCTCCGCCAACAACAATATTCCATTTGGAAGACCAAACGCCGCCACACTGTATACTCGCCATGACGATCATaaa                                                                                                  <  1:90346/150‑1 (MQ=255)
                                                    acGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAAccc                                                                                              >  1:490801/1‑150 (MQ=255)
                                                     cGCCACGGTTTCCGTCAACGGGCAGCATAAACTTAGTTCTTAAAAAAAAAATGAAGACGATGGTGATCTGTTCTACGTAAACATCAGTATTCCTTTTGGCAAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAAccc                                                                                              <  2:708278/149‑1 (MQ=255)
                                                       ccACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACt                                                                                           >  1:661899/1‑150 (MQ=255)
                                                       ccACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACACCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACt                                                                                           >  1:729679/1‑150 (MQ=255)
                                                                ccGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTTTTCTACGTCAACATCAGTATTCCTTTTGGAAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACt                                                                                   <  1:6460/149‑1 (MQ=255)
                                                                          gcagcaTCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATggg                                                                        >  1:39271/1‑150 (MQ=255)
                                                                           caacaTCAACTTAGTGCTTAAAAAAATAATGAACCGGAGGGTTATCTTTTCTCCTCCAACTCCTGTATCCCTTTTGGCAGATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGa                                                                       <  1:176199/147‑1 (MQ=255)
                                                                            agcaTCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAg                                                                      <  2:720185/150‑1 (MQ=255)
                                                                                  aatttAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTTTTCTCCGTCAACATCAGTTTTCCTTTTGCCAAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGtct                                                                <  2:189627/147‑1 (MQ=255)
                                                                                                  aaaaaaggaagcggttGGTGATCTGTTTTACGTAAACATCAGTATTTTTTTTGGAAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCctac                                               <  2:448447/136‑1 (MQ=255)
                                                                                                     taatGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTAt                                            >  2:715718/1‑151 (MQ=255)
                                                                                                      aatGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTAtg                                           >  1:786569/1‑151 (MQ=255)
                                                                                                          aaGACGATGGTGATCTGTTCTACGTAAACATCAGTATTCCTTTTGGCAGATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGa                                        <  1:786745/150‑1 (MQ=255)
                                                                                                           aGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAAt                                      >  2:336652/1‑151 (MQ=255)
                                                                                                                 tGGTGATCTTTTCTCCGTAAACATCAGTTTTCATTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGa                                 <  2:349872/150‑1 (MQ=255)
                                                                                                                  ggTGATCTTTTCTACGTCAACATCAGTATCCCTTTTGGAAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCt                                   <  2:637925/147‑1 (MQ=255)
                                                                                                                  ggTGATCTGTTCTACGTCACCATCATTATTCCTTTTGGCAAACCACCCTCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGaa                                <  2:280119/150‑1 (MQ=255)
                                                                                                                         tGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATc                         >  1:346728/1‑150 (MQ=255)
                                                                                                                            tCTACGTCAACATCATTATTCTTTTTGGAAGATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCacg                      <  1:771847/150‑1 (MQ=255)
                                                                                                                             cTACGTCAACACCAGTATCCCTTTTGGAAGATCAACCACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCacg                      <  2:17891/149‑1 (MQ=255)
                                                                                                                              tACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCacgacg                   >  1:406946/1‑151 (MQ=255)
                                                                                                                                cGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCAcgacga                  >  1:147395/1‑150 (MQ=255)
                                                                                                                                      cATCAGTATTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACg               <  1:492042/147‑1 (MQ=255)
                                                                                                                                          aGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACg        >  2:68184/1‑150 (MQ=255)
                                                                                                                                            tATTCCTTTTGCAAGACCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAg      <  2:14310/150‑1 (MQ=255)
                                                                                                                                               tCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTg  >  2:117826/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTG  >  NC_000913/3364250‑3364543

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