Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,189,8740CA38.1% 65.5 / 47.2 42R346L (CGT→CTT) yehBputative fimbrial usher protein YehB
Reads supporting (aligned to +/- strand):  ref base C (10/16);  new base A (16/0);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 4.09e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGC  >  NC_000913/2189778‑2189968
                                                                                                   |                                                                                               
accGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAc                                                                                                 <  1:18789208/95‑1 (MQ=255)
aTCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAc                                                                                                 <  1:10410706/97‑1 (MQ=255)
aTCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAc                                                                                                 <  1:3608473/97‑1 (MQ=255)
 tCGAGCCCCCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACg                                                                                                <  1:14901303/97‑1 (MQ=255)
 tCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACg                                                                                                <  1:6541882/97‑1 (MQ=255)
 tCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACg                                                                                                <  1:2208886/97‑1 (MQ=255)
 cggAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACg                                                                                                <  1:6152216/95‑1 (MQ=255)
   gacccaccaTACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGAc                                                                                              <  1:11584636/94‑1 (MQ=255)
   gAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGAc                                                                                              <  1:13100350/97‑1 (MQ=255)
   gAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGAc                                                                                              <  1:10456833/97‑1 (MQ=255)
   gAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGAc                                                                                              <  1:17636435/97‑1 (MQ=255)
    aGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGAcc                                                                                             <  1:4551967/97‑1 (MQ=255)
    aGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGAcc                                                                                             <  1:4297780/97‑1 (MQ=255)
    aGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGAcc                                                                                             <  1:10663778/97‑1 (MQ=255)
     gCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGAccc                                                                                            <  1:19684859/97‑1 (MQ=255)
      tcaccaTACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACccg                                                                                           <  1:10428311/96‑1 (MQ=255)
      ccaccaTACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACccg                                                                                           <  1:4338323/97‑1 (MQ=255)
      ccaccaTACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACccg                                                                                           <  1:834849/97‑1 (MQ=255)
       caccaTACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAcccac                                                                                          >  1:1134227/1‑95 (MQ=255)
        accaTACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACccgcc                                                                                         <  1:2930702/97‑1 (MQ=255)
            tACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAAccgccccca                                                                                     >  1:353976/1‑93 (MQ=255)
            tACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAAACGCCGCtc                                                                                     >  1:18281140/1‑96 (MQ=255)
                gCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAAACGCCGCCAAATc                                                                                 >  1:17594000/1‑97 (MQ=255)
                     aaTAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAAACGCCGCTACAACataat                                                                            >  1:12065623/1‑95 (MQ=255)
                       tAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACCAGGCACGCTATCATATTTc                                                                          >  1:14948848/1‑97 (MQ=255)
                          attattAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAACACCCCGCCAAATCATATTTTGac                                                                       >  1:1887782/1‑97 (MQ=255)
                             attaAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACCCCCCGCCAAATCATATTTTGACACg                                                                    >  1:3478736/1‑97 (MQ=255)
                               taAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAACCCCCCGCGGAATCATATTTTTacaccac                                                                  >  1:4731575/1‑94 (MQ=255)
                                   aCCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAGAACCGCCGCTAAATCATATTTCGACACCCCcggc                                                              >  1:19341998/1‑96 (MQ=255)
                                      aTACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAACCACCGCCACAACATAATTTTACACGCCCGGtgtc                                                           >  1:20166202/1‑94 (MQ=255)
                                         cTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAACCCCCCCCCAAATCATATTTTCACACCCCGGGGTggggc                                                        >  1:17097655/1‑93 (MQ=255)
                                         cTGACACCCCGCCTGGAAAAAATCTTTTTGTTTGCCCCCCCCTTCAATATGGCTCCGCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTgcagc                                                        <  1:10035557/97‑1 (MQ=255)
                                          tGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAACCCCCCCCTAAAACATATTTTTACACGCCGCGGTgcagcg                                                       >  1:16129377/1‑96 (MQ=255)
                                            aTAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAACCCCCGCTCAAACATATTTTTACACCCCGCGGTGGAGCatc                                                     >  1:17745253/1‑96 (MQ=255)
                                             cacccccccTGAAAAAAATCTTTTTTTTTCCCCCCCCCTTAAATATGCCTCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCatat                                                   <  1:10274706/94‑1 (MQ=255)
                                              acccccccTGGAAAAAATCCCTTTTTTTTGCCCCCCCCTTCAATATGGCCCCGCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCatat                                                   <  1:10752520/95‑1 (MQ=255)
                                              accccGCCGGAAAAAAATCCTTTTGTTTCCCCCCCCCTTCAATAGGGCTACGCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATAtt                                                  <  1:12502290/95‑1 (MQ=255)
                                               ccccGCCTGGACAAAATCTTTTTGTTTCCTCCCCCCTTCAATAGGGCTACGCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATAttt                                                 <  1:3178985/96‑1 (MQ=255)
                                               ccccGCCTGCAAAAAATCTTTTTGTGTGCTCCCCCCTTAAATAGGGCCCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATAttt                                                 <  1:162716/96‑1 (MQ=255)
                                               ccccGCCTCGAGAAAATCATTTTTTTTTCCCGCCCTTTTATTACGGCAACGCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATAttt                                                 <  1:20515438/96‑1 (MQ=255)
                                                cccGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAACCCCCCCCCCAAACAAAATTTGTCACGCCCGGGCGCAGCACAtttt                                                >  1:18634573/1‑96 (MQ=255)
                                                  cGCCTGAAAAAAATCTTTTTGTTGCCCCCCCCCTACAATAGGGCTCCGCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGc                                              <  1:14156085/97‑1 (MQ=255)
                                                   gCCTGAAAAAAATCCCTTGTTTGCCCCCCCCCTTAAATATGCCACCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCa                                             <  1:9929613/97‑1 (MQ=255)
                                                      tGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAGACACCCCGCTCAAACATATTTCTACACGCCCGGGGGCAGCACATTTTTCAcca                                          >  1:8406281/1‑96 (MQ=255)
                                                      tGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAACCCCCGCCAAATAATATTTTTTCACCCCCGGGTGCCGCATATTTTGCACCg                                          >  1:601857/1‑97 (MQ=255)
                                                       ggACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAACCCCCCCCAAAATCAATTTTTACACGCCCGGGTGCAGCATAATTTGCAcccc                                         >  1:20361017/1‑95 (MQ=255)
                                                       ggACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAACCCCCGCCAAAAAAAATTTTCACACGCCCGGGTGCAGCATATTTTGCACCGc                                         >  1:20961285/1‑97 (MQ=255)
                                                       ggACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAACCCCCCCCCAAACATATTTTTTAACGCCCGGGGGCAGCCTATTTTGGACCGc                                         >  1:13344777/1‑97 (MQ=255)
                                                        caaaaaaTCATTTTGTTTCCCCCCCCCTACAAACTGGCAACCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCt                                        <  1:9552194/94‑1 (MQ=255)
                                                         aaaaaaTCGTTTTGTTGCCCCCCCCCTTAAATAGGGCTCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTg                                       <  1:13527125/95‑1 (MQ=255)
                                                          cAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAACCCCCCCCCAATCATAATTTCACACGCCGGGGTGGCGCCAATTTTGGTCCGCTGc                                      >  1:18154331/1‑97 (MQ=255)
                                                          cAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAAACCCCGCCCAAACATATTTTTTCACCCCGGGGTGCGGCATAATTTGCACCGCTGc                                      >  1:7649895/1‑97 (MQ=255)
                                                             aaTCACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAAACCGCCGCTAAAAACTATTTTGACACGCCCGGGGGCAGCACATTTTGCAACTCTGCATa                                   >  1:13835344/1‑97 (MQ=255)
                                                                tctttttgtttcCCCCCCCCTTAAATAGGGCTCCGCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAg                                 <  1:8399143/93‑1 (MQ=255)
                                                                cACTTTGTTTGCTCGCCCCTTCAATATGGCTAAAACCCCCCGCCAAAAAATATTTTTACAACCCCGGGTGCGGCATAATTTGGCCCGCTGCACAAgg                                >  1:15052251/1‑97 (MQ=255)
                                                                    tgtgttcgCCCGCCCCTTCAATATGCCCCCGCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCAc                             <  1:1677299/95‑1 (MQ=255)
                                                                      tgtttcCCCCCCCCTTCAATATGGCTCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCa                           <  1:8796241/97‑1 (MQ=255)
                                                                      gttttgCCCCCCCCTTAAATAGGCCTCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCa                           <  1:1809445/95‑1 (MQ=255)
                                                                              cccccctTTAAATAGGGCTCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAggt                    <  1:271068/94‑1 (MQ=255)
                                                                                 ccccTTCAATAGGGCCACCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAggtggt                 <  1:7972154/97‑1 (MQ=255)
                                                                                  cccTTCATTAGGGCCCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGtt                <  1:15072047/97‑1 (MQ=255)
                                                                                    tttCAATAGGGCTCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTAc              <  1:18203941/96‑1 (MQ=255)
                                                                                    gttCAATGCGGCTCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTAc              <  1:506004/96‑1 (MQ=255)
                                                                                    cTTGAATATGCCTCCGCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTAc              <  1:6843661/97‑1 (MQ=255)
                                                                                      tCAATATCGCTCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCg            <  1:6604941/97‑1 (MQ=255)
                                                                                      caaaTAGGGCTCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCg            <  1:4789930/95‑1 (MQ=255)
                                                                                       aaaTGTGGCTCCCCCCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGa           <  1:7033373/96‑1 (MQ=255)
                                                                                            tGGCTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGt      >  1:11678065/1‑97 (MQ=255)
                                                                                            tGGCTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGGTACCCAGCCGt      >  1:12017303/1‑97 (MQ=255)
                                                                                             ggCTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTc     >  1:20477716/1‑97 (MQ=255)
                                                                                              gCTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGCGCCGTcc    >  1:1948705/1‑97 (MQ=255)
                                                                                              gCTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTcc    >  1:20173163/1‑97 (MQ=255)
                                                                                               cTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCg   >  1:18914823/1‑97 (MQ=255)
                                                                                                tACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGc  >  1:1308715/1‑97 (MQ=255)
                                                                                                tACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGc  >  1:14984948/1‑97 (MQ=255)
                                                                                                tACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGc  >  1:20180745/1‑97 (MQ=255)
                                                                                                tACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGc  >  1:19531464/1‑97 (MQ=255)
                                                                                                   |                                                                                               
ATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGC  >  NC_000913/2189778‑2189968

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: