Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,549,001 G→A C65Y (TGC→TAC)  fimH → minor component of type 1 fimbriae

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,0010GA100.0% 238.4 / NA 73C65Y (TGC→TAC) fimHminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/35);  total (38/35)

GCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCAC  >  NC_000913/4548874‑4549130
                                                                                                                               |                                                                                                                                 
gCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGAtt                                                                                                                        <  2:1660536/139‑1 (MQ=255)
     tAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATaa                                                                                                                             <  1:812255/129‑1 (MQ=255)
     tAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATaa                                                                                                                             >  2:812255/1‑129 (MQ=255)
     tAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCg                                                                                                                   <  2:149272/139‑1 (MQ=255)
     tAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCg                                                                                                                   >  1:149272/1‑139 (MQ=255)
            gCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCa                                                                                                            >  2:1511377/1‑139 (MQ=255)
               aaTGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTa                                                                                                         >  1:820997/1‑139 (MQ=255)
                           ctCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTca                                                                                             >  2:1056545/2‑139 (MQ=255)
                            tCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACa                                                                                                       <  1:793301/128‑1 (MQ=255)
                            tCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACa                                                                                                       >  2:793301/1‑128 (MQ=255)
                              ccTATTGGCGGTGGCCGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGAtt                                                                                                                        >  2:1128892/1‑109 (MQ=255)
                              ccTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGAtt                                                                                                                        <  1:1128892/109‑1 (MQ=255)
                                 aTTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGc                                                                                       >  2:1048929/1‑139 (MQ=255)
                                     gCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGa                                                                                                                 <  1:773402/109‑1 (MQ=255)
                                     gCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGa                                                                                                                 >  2:773402/1‑109 (MQ=255)
                                      cGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATaa                                                                                                                             >  1:938348/1‑96 (MQ=255)
                                      cGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATaa                                                                                                                             <  2:938348/96‑1 (MQ=255)
                                      cGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGa                                                                                                     <  2:1311882/120‑1 (MQ=255)
                                      cGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGa                                                                                                     >  1:1311882/1‑120 (MQ=255)
                                        gTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAgg                                                                                <  2:1372732/139‑1 (MQ=255)
                                          ggCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                              <  2:820997/139‑1 (MQ=255)
                                          ggCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCAACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCt                                                                              <  1:1056545/139‑1 (MQ=255)
                                                               aaCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTca                                                                                             <  2:1110138/103‑1 (MQ=255)
                                                               aaCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTca                                                                                             >  1:1110138/1‑103 (MQ=255)
                                                               aaCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTAt                                                         >  1:1689373/1‑139 (MQ=255)
                                                                aCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACa                                                                                                       >  2:802077/1‑92 (MQ=255)
                                                                aCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACa                                                                                                       <  1:802077/92‑1 (MQ=255)
                                                                aCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                              >  2:437246/1‑117 (MQ=255)
                                                                aCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                              <  1:437246/117‑1 (MQ=255)
                                                                aCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATc                                                        >  2:862297/1‑139 (MQ=255)
                                                                 ccTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCt                                                       >  2:365984/1‑139 (MQ=255)
                                                                   ttGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCAt                                                                                                                               <  1:583158/65‑1 (MQ=255)
                                                                   ttGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCAt                                                                                                                               >  2:583158/1‑65 (MQ=255)
                                                                     gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACa                                                                                                       >  1:102043/1‑87 (MQ=255)
                                                                     gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACa                                                                                                       <  2:102043/87‑1 (MQ=255)
                                                                     gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAAtt                                                   >  2:459540/1‑139 (MQ=255)
                                                                      cgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAAttt                                                  >  2:63740/1‑139 (MQ=255)
                                                                            tcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATg                                                                                                <  1:1481765/87‑1 (MQ=255)
                                                                            tcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATg                                                                                                >  2:1481765/1‑87 (MQ=255)
                                                                            tcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCgg                                            <  1:63740/139‑1 (MQ=255)
                                                                              gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                          <  2:1689373/139‑1 (MQ=255)
                                                                              gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                          <  1:459540/139‑1 (MQ=255)
                                                                              gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                          <  1:1511377/139‑1 (MQ=255)
                                                                              gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                          <  1:117727/139‑1 (MQ=255)
                                                                              gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                          <  2:1345399/139‑1 (MQ=255)
                                                                              gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTACGGCGGCGTGTTATCTAATTTTTCCGGGa                                          <  1:1048929/139‑1 (MQ=255)
                                                                                   tgtgGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTa                                     >  1:1636480/1‑139 (MQ=255)
                                                                                     tgGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGAtt                                                                                                                        >  1:1189856/1‑54 (MQ=255)
                                                                                     tgGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGAtt                                                                                                                        <  2:1189856/54‑1 (MQ=255)
                                                                                     tgGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCg                                                                                                                   >  2:1471857/1‑59 (MQ=255)
                                                                                     tgGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCg                                                                                                                   <  1:1471857/59‑1 (MQ=255)
                                                                                       gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                               >  2:555814/1‑61 (MQ=255)
                                                                                       gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                               <  1:555814/61‑1 (MQ=255)
                                                                                       gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGAcc                                        >  1:678952/1‑132 (MQ=255)
                                                                                       gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGAcc                                        <  2:678952/132‑1 (MQ=255)
                                                                                       gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGAcc                                        >  2:1490509/1‑132 (MQ=255)
                                                                                       gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGAcc                                        <  1:1490509/132‑1 (MQ=255)
                                                                                       gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTaaa                                   <  2:280000/137‑1 (MQ=255)
                                                                                       gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTaaa                                   >  1:280000/1‑137 (MQ=255)
                                                                                        ggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAata                                >  2:1208919/1‑139 (MQ=255)
                                                                                                   gTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                         <  1:1272998/135‑1 (MQ=255)
                                                                                                   gTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                         >  2:1272998/1‑135 (MQ=255)
                                                                                                       tGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCGACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGGAAAATATAGTGGCAGTAGCTa                 >  2:984046/1‑139 (MQ=255)
                                                                                                         gATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATc               >  1:1080353/1‑139 (MQ=255)
                                                                                                             tttCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGgcggcg                                                               <  2:746581/87‑1 (MQ=255)
                                                                                                             tttCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGgcggcg                                                               >  1:746581/1‑87 (MQ=255)
                                                                                                              ttCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTca                                                                                             <  2:1079106/56‑1 (MQ=255)
                                                                                                              ttCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTca                                                                                             >  1:1079106/1‑56 (MQ=255)
                                                                                                              ttCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                              >  2:756141/1‑71 (MQ=255)
                                                                                                              ttCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                              <  1:756141/71‑1 (MQ=255)
                                                                                                                      aaaTCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                          >  1:1236038/1‑99 (MQ=255)
                                                                                                                      aaaTCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                          <  2:1236038/99‑1 (MQ=255)
                                                                                                                      aaaTCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTaccac  >  2:1500451/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                                 
GCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCAC  >  NC_000913/4548874‑4549130

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