Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,547,436 A→G T781A (ACT→GCT)  fimD → type I fimbriae usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,547,4360AG100.0% 167.7 / NA 47T781A (ACT→GCT) fimDtype I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/22);  total (25/22)

ACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACATCAGAGAGTAGCCAGAGTAGCGGCATT  >  NC_000913/4547310‑4547570
                                                                                                                              |                                                                                                                                      
aCCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGAt                                                                                                                           <  1:160721/140‑1 (MQ=255)
    aCTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTg                                                                                                                       >  2:532692/1‑140 (MQ=255)
            ggTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCTGCTCGTGGGGCGATCGTGCGAGCaga                                                                                                               >  2:401942/1‑140 (MQ=255)
                       gctgcCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAgcgc                                                                                                    >  2:428605/1‑140 (MQ=255)
                        ctgcCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGc                                                                                                                  >  2:739232/1‑125 (MQ=255)
                        ctgcCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGc                                                                                                                  <  1:739232/125‑1 (MQ=255)
                        ctgcCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAgcgcg                                                                                                   >  2:581772/1‑140 (MQ=255)
                           ccTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGTTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGt                                                                                                <  1:710733/140‑1 (MQ=255)
                           ccTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGt                                                                                                <  2:406373/140‑1 (MQ=255)
                           ccTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGt                                                                                                <  1:482616/140‑1 (MQ=255)
                                     cTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATaaaa                                                                                      >  1:430020/1‑140 (MQ=255)
                                     cTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATaaaa                                                                                      >  1:711984/1‑140 (MQ=255)
                                     cTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCGCGTGGGGCGATCGTGCCAGCAGAGTCTAAAGCGCGCGTTGGGATaaaa                                                                                      >  2:286640/1‑140 (MQ=255)
                                       gAATATCGGGAAAATAGAGTGGGGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACt                                                                                    <  1:288624/140‑1 (MQ=255)
                                             cGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAAATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCAt                                                                              <  2:723536/140‑1 (MQ=255)
                                                    aTAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTg                                                                       >  2:82564/1‑140 (MQ=255)
                                                      agagTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGAc                                                                     <  1:577146/140‑1 (MQ=255)
                                                       gagTGGCGCTGGATACCAATACTCTTGCTGATAACGTAGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGAcc                                                                    <  1:286640/140‑1 (MQ=255)
                                                       gagTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGAcc                                                                    <  2:430020/140‑1 (MQ=255)
                                                         gTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCca                                                                  <  1:64654/140‑1 (MQ=255)
                                                         gTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCca                                                                  <  1:532692/140‑1 (MQ=255)
                                                          tGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCcac                                                                 >  2:451824/1‑140 (MQ=255)
                                                          tGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCcac                                                                 >  1:671406/1‑140 (MQ=255)
                                                                       cAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATGAAACTGCGCATGACGCTGACCAACAATAATAAGCCGc                                                    >  1:211515/1‑140 (MQ=255)
                                                                            cccTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAACCTGATCATGACGCTGATCCACGATAATAAGCCGCTGCCg                                               >  2:279965/1‑140 (MQ=255)
                                                                            cccTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTTATAACGCTGACCCACATTAATAAGACGCTGCCg                                               >  2:324008/1‑140 (MQ=255)
                                                                             ccTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCcaca                                                                <  1:400888/122‑1 (MQ=255)
                                                                             ccTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCcaca                                                                >  2:400888/1‑122 (MQ=255)
                                                                                 gCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGATGAAGACCCACAATAATCAGCCGCTGCCGTTTgg                                          >  2:353975/1‑140 (MQ=255)
                                                                                 gCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTgg                                          >  1:720229/1‑140 (MQ=255)
                                                                                 gCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTgg                                          >  2:11618/1‑140 (MQ=255)
                                                                                       aaCGTCGATTTAGATGACGCGGTTGGTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGAt                                    <  1:279965/140‑1 (MQ=255)
                                                                                           tCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTg                                >  2:479880/1‑140 (MQ=255)
                                                                                           tCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTg                                <  2:711984/140‑1 (MQ=255)
                                                                                               tttAGATAACGCGGTTGCTAACGTTGTTCCTGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACAt                            <  1:401942/140‑1 (MQ=255)
                                                                                               tttAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGCCGCTGACCCACAGTAAAAAGCCGCTACCGTCTGGGGCGATGGTGACAt                            >  2:298812/1‑140 (MQ=255)
                                                                                                ttAGATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACATc                           <  1:185777/140‑1 (MQ=255)
                                                                                                   gATAACGCGGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACATCaga                        <  1:11618/140‑1 (MQ=255)
                                                                                                        cgcgGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACATCAGAGAGTa                   >  1:132363/1‑140 (MQ=255)
                                                                                                         gcgGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGt                                 >  2:660087/1‑125 (MQ=255)
                                                                                                         gcgGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGt                                 <  1:660087/125‑1 (MQ=255)
                                                                                                         gcgGTTGCTAACATTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACATCAGAGAGTAg                  >  1:440083/1‑140 (MQ=255)
                                                                                                          cgGTTGCTAACGTTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACATCAGAGAGTAGc                 <  1:183313/140‑1 (MQ=255)
                                                                                                                cTAACATTGTTCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGGGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACATCAGAGAGTAGCCAGAGt           <  2:440083/140‑1 (MQ=255)
                                                                                                                  aaCGTTGTTCCCGCTTGTGGGGCGATCGTTCGAGCAGGGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCGCAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACATCAGAGAGTAGCCAGAGTAg         <  1:298812/140‑1 (MQ=255)
                                                                                                                         ttCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACATCAGAGAGTAGCCAGAGTAGCGGCAtt  >  2:589353/1‑140 (MQ=255)
                                                                                                                         ttCCCGCTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACATCAGAGAGTAGCCAGAGTAGCGGCAtt  >  2:619184/1‑140 (MQ=255)
                                                                                                                              |                                                                                                                                      
ACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCACAATAATAAGCCGCTGCCGTTTGGGGCGATGGTGACATCAGAGAGTAGCCAGAGTAGCGGCATT  >  NC_000913/4547310‑4547570

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