Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,543,958 T→G V77V (GTT→GTG)  fimI → fimbrial protein involved in type 1 pilus biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,543,9580TG98.0% 178.3 / ‑4.1 51V77V (GTT→GTG) fimIfimbrial protein involved in type 1 pilus biosynthesis
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (21/29);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

CGGCGGAAATATGCAATTTCAGGGCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGAT  >  NC_000913/4543811‑4544105
                                                                                                                                                   |                                                                                                                                                   
cggcggAAATATGCAATTTCAGGGCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGAtt                                                                                                                                                  <  1:452389/151‑1 (MQ=255)
 ggcggAAATATGCAATTTCAGGGCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTc                                                                                                                                                 >  1:170834/1‑151 (MQ=255)
   cggAAATATGCAATTTCAGGGCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCAt                                                                                                                                               <  2:508119/151‑1 (MQ=255)
      aaaTATGCAATTTCAGGGCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTa                                                                                                                                            <  2:557762/151‑1 (MQ=255)
       aaTATGCAATTTCAGGGCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTAc                                                                                                                                           <  2:116037/151‑1 (MQ=255)
           tGCAATTTCAGGGCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACggg                                                                                                                                        >  1:529572/1‑150 (MQ=255)
              aaTTTCAGGGCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATg                                                                                                                                    >  2:546286/1‑151 (MQ=255)
                   cAGGGCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGc                                                                                                                                <  1:12038/150‑1 (MQ=255)
                    aGGGCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCAc                                                                                                                              <  1:618672/151‑1 (MQ=255)
                     gggCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACg                                                                                                                             >  2:8959/1‑151 (MQ=255)
                         gTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACggtgg                                                                                                                         >  1:536524/1‑151 (MQ=255)
                            attattGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGgtga                                                                                                                      <  2:396924/151‑1 (MQ=255)
                               attGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGgtgagt                                                                                                                    <  2:459233/150‑1 (MQ=255)
                                 tGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGgtgagtga                                                                                                                  >  2:800093/1‑150 (MQ=255)
                                  gCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAATCGGTTTCATGCGGTTGGGGAAGATATCGCACCGGTGCGTTTTGTGATTCATTTACGGGAATGTAGGACGGTGGTGAGTGaa                                                                                                                 <  1:287916/150‑1 (MQ=255)
                                  gCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAAc                                                                                                                <  1:867561/151‑1 (MQ=255)
                                  gCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAAc                                                                                                                >  1:322417/1‑151 (MQ=255)
                                   cGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACg                                                                                                               >  2:69894/1‑151 (MQ=255)
                                        aCTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACgtgt                                                                                                            >  1:928759/1‑149 (MQ=255)
                                                gATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCg                                                                                                  <  2:314140/151‑1 (MQ=255)
                                                     aaGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCa                                                                                             <  2:810055/151‑1 (MQ=255)
                                                      aGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCAc                                                                                            >  1:180867/1‑151 (MQ=255)
                                                                  aCAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGAt                                                                                <  1:696186/151‑1 (MQ=255)
                                                                         aCGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTaaa                                                                          >  2:450331/1‑150 (MQ=255)
                                                                          cGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTaaaaa                                                                        <  2:575578/151‑1 (MQ=255)
                                                                               aaTATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCgg                                                                   <  2:476808/151‑1 (MQ=255)
                                                                               aaTATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCgg                                                                   >  2:330625/1‑151 (MQ=255)
                                                                                    ggggCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGAtgt                                                               <  2:29339/150‑1 (MQ=255)
                                                                                         aaaTCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTc                                                         <  1:253331/151‑1 (MQ=255)
                                                                                         aaaTCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTc                                                         >  2:416777/1‑151 (MQ=255)
                                                                                         aaaTCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTc                                                         <  1:753536/151‑1 (MQ=255)
                                                                                         aaaTCAGCAGTAACCGGTTTCATGAGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGCGCttt                                                          >  1:645036/1‑150 (MQ=255)
                                                                                             cagcagTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGt                                                      <  2:163343/150‑1 (MQ=255)
                                                                                              agcagTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTgg                                                    <  2:320906/151‑1 (MQ=255)
                                                                                              agcagTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTgg                                                    >  2:30054/1‑151 (MQ=255)
                                                                                               gcagTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTggg                                                   >  2:334423/1‑151 (MQ=255)
                                                                                                      ccGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAgggg                                            >  2:35859/1‑151 (MQ=255)
                                                                                                          tttCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGcc                                          >  2:4052/1‑149 (MQ=255)
                                                                                                             cATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAggg                                      <  1:99415/150‑1 (MQ=255)
                                                                                                                      ggggAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGccacc                             >  2:139941/1‑150 (MQ=255)
                                                                                                                       gggAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCaa                           <  1:546268/151‑1 (MQ=255)
                                                                                                                        ggAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAat                          <  2:529589/151‑1 (MQ=255)
                                                                                                                         gAAGATAGCGCACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAat                          >  1:647392/1‑150 (MQ=255)
                                                                                                                                 cgcACCGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTa                 <  2:925794/151‑1 (MQ=255)
                                                                                                                                     ccGGTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCgt             >  1:610816/1‑151 (MQ=255)
                                                                                                                                       ggTGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCgttg           <  1:915761/151‑1 (MQ=255)
                                                                                                                                         tGCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCgttgtt         >  2:915741/1‑151 (MQ=255)
                                                                                                                                          gCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCgttgtt         <  1:764278/150‑1 (MQ=255)
                                                                                                                                          gCCTTTTGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGttt        <  2:189190/151‑1 (MQ=255)
                                                                                                                                              tttGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTtga      <  2:339115/149‑1 (MQ=255)
                                                                                                                                               ttGTGATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTtgatg    <  2:362445/150‑1 (MQ=255)
                                                                                                                                                tgtgATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTtgatgat  >  2:807257/1‑151 (MQ=255)
                                                                                                                                                tgtgATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTtgatga   <  1:267537/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGGCGGAAATATGCAATTTCAGGGCGTCATTATTGCGGAAACTTGCCGGATTGAAGCCGGTGATAAACAAATGACGGTCAATATGGGGCAAATCAGCAGTAACCGGTTTCATGCGGTTGGGGAAGATAGCGCACCGGTGCCTTTTGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGAT  >  NC_000913/4543811‑4544105

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