Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,534,382 C→T P88S (CCT→TCT)  fimC → type 1 fimbria chaperone FimC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,534,3820CT96.7% 220.4 / ‑2.0 61P88S (CCT→TCT) fimCtype 1 fimbria chaperone FimC
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (30/29);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCA  >  NZ_CP009273/4534246‑4534516
                                                                                                                                        |                                                                                                                                      
gACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGtct                                                                                                                                      >  2:897206/1‑139 (MQ=255)
 aCTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTATCGTGACGtctcc                                                                                                                                    >  2:359232/1‑139 (MQ=255)
     gcgTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGtctcct                                                                                                                                   <  1:1118053/137‑1 (MQ=255)
               tCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGAt                                                                                                                        >  1:798172/1‑138 (MQ=255)
                 cGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGaa                                                                                                                     <  2:883636/139‑1 (MQ=255)
                     aGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAggg                                                                                                                  >  1:1226755/1‑138 (MQ=255)
                      gggCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGaa                                                                                                                <  2:1169240/139‑1 (MQ=255)
                       ggCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGaa                                                                                                                >  1:555715/1‑138 (MQ=255)
                       ggCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGaa                                                                                                                >  1:1075905/1‑138 (MQ=255)
                       ggCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGaa                                                                                                                >  1:888743/1‑138 (MQ=255)
                             aaCAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAaga                                                                                                         <  1:637268/139‑1 (MQ=255)
                               cAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAgaga                                                                                                       <  2:452386/139‑1 (MQ=255)
                               cAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAgaga                                                                                                       <  2:953280/139‑1 (MQ=255)
                                aaGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGaa                                                                                                      <  2:762037/139‑1 (MQ=255)
                                        cTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGaa                                                                                                      >  2:104479/1‑131 (MQ=255)
                                        cTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGaa                                                                                                      <  1:104479/131‑1 (MQ=255)
                                                     ataatGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGc                                                                                 <  2:1208992/139‑1 (MQ=255)
                                                           atgaAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACaaa                                                                           >  2:1210994/1‑139 (MQ=255)
                                                                aaTAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAAc                                                                       >  1:83570/1‑138 (MQ=255)
                                                                          aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAGGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCca                                                              >  2:762090/1‑137 (MQ=255)
                                                                          aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCca                                                              >  2:1023975/1‑137 (MQ=255)
                                                                           tttAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcaca                                                            <  1:1210994/138‑1 (MQ=255)
                                                                           tttAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcac                                                             <  2:1093562/137‑1 (MQ=255)
                                                                               aTTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGAcc                                                       >  1:827908/1‑139 (MQ=255)
                                                                                   aaTCATGGGTGGAAAATGCCGATGGTGTAAGGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCggg                                                    <  1:762090/138‑1 (MQ=255)
                                                                                   aaTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCggg                                                    <  2:1238202/138‑1 (MQ=255)
                                                                                   aaTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCggg                                                    <  2:187525/138‑1 (MQ=255)
                                                                                   aaTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCggg                                                    <  2:902201/138‑1 (MQ=255)
                                                                                     tCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGc                                                                                 >  2:487757/1‑107 (MQ=255)
                                                                                     tCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGc                                                                                 <  1:487757/107‑1 (MQ=255)
                                                                                     tCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGa                                                   >  2:560707/1‑137 (MQ=255)
                                                                                       aTGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGt                                               <  1:952967/139‑1 (MQ=255)
                                                                                              gAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCt                                        >  2:300023/1‑139 (MQ=255)
                                                                                                   tGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                                    >  2:1071384/1‑138 (MQ=255)
                                                                                                   tGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATg                                   <  1:1071384/139‑1 (MQ=255)
                                                                                                         tGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATAc                                                                                                   <  2:1270498/69‑1 (MQ=255)
                                                                                                         tGGTGTAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATAc                                                                                                   >  1:1270498/1‑69 (MQ=255)
                                                                                                            tgtAAAGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTaaa                          <  1:897206/139‑1 (MQ=255)
                                                                                                                 aGGATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGAt                     <  2:440521/139‑1 (MQ=255)
                                                                                                                   gATGGTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTc                   <  1:291897/139‑1 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGCGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGt                >  2:767533/1‑139 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATAc                                                                                                   <  2:578167/56‑1 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATAc                                                                                                   >  1:578167/1‑56 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGc                                                                                 <  2:884306/74‑1 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGc                                                                                 >  1:884306/1‑74 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAAt                                                                   >  2:1116732/1‑88 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAAt                                                                   <  1:1116732/88‑1 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAAt                                                                   >  1:544974/1‑88 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAAt                                                                   <  2:544974/88‑1 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcac                                                             >  1:178862/1‑94 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcac                                                             <  2:178862/94‑1 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCggg                                                    >  2:1023725/1‑103 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCggg                                                    >  1:1011374/1‑103 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCggg                                                    <  2:1011374/103‑1 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCggg                                                    <  1:1023725/103‑1 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGAtt                    >  2:393969/1‑135 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGAtt                    <  1:393969/135‑1 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGt                >  1:672455/1‑139 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGt                >  2:84507/1‑139 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGt                >  1:499382/1‑139 (MQ=255)
                                                                                                                      ggTCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGt                >  1:119043/1‑139 (MQ=255)
                                                                                                                        tCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTc                   >  1:880276/1‑134 (MQ=255)
                                                                                                                        tCGTTTTATCGTGACGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTc                   <  2:880276/134‑1 (MQ=255)
                                                                                                                                     aCGTCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCa  >  1:44058/1‑138 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCA  >  NZ_CP009273/4534246‑4534516

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