Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,549,223 T→G V139G (GTG→GGG)  fimH → type 1 fimbriae D‑mannose specific adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,2230TG96.9% 93.8 / ‑6.2 32V139G (GTG→GGG) fimHtype 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (17/14);  minor base C (1/0);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCC  >  NC_000913/4548991‑4549464
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
cAAATCTTTTGCCATAACGATTATCCGGCAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTTACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTTGCTc                                                                                                                                                                                                                                   >  2:55624/1‑251 (MQ=255)
       tttGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTa                                                                                                                                                                                                                               <  1:55624/248‑1 (MQ=255)
        ttGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGCGGCTCGGCTTATGGAGGCGTGATATCAAATTTTTCCGGGACCGTAAAATAAAGTACAAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGCGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGATCAGTGAGGGCGGGGGGG‑CGATTAAA‑GCTGGGCCATTAATTg                                                                                                                                                                                                                           >  1:114629/1‑251 (MQ=255)
                 cGATTATCCGGAAACCATTACAGACTATGTCACTCTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTa                                                                                                                                                                                                                  <  1:12923/251‑1 (MQ=255)
                          ggAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATACATTTCCTACCACCCGCGAAACGCCGCGAGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGCTTAAA‑GCTGCCTCATTAATTGCCGTGCTTATTTTGCGAc                                                                                                                                                                                                         >  1:24770/1‑251 (MQ=255)
                                            tGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGGGTTATCTAATTTTTCCGGGACCGTAACATATAGTGGAGGTAGCTATCCATTTCATACCACCAGCGAAACGCCGCGCGTTGATTATAATTCGAGAACGGATAAGCACTGGCCGGGGGCGCTTTGTTTTACGCCTATGTGCAGTGAGGGCCGGGGGG‑CGATTAAAAGC‑GGCTCCTTAATTGCCGGGGTTATTTTGCGACAGACCCAACA‑CTATaaa                                                                                                                                                                                       >  2:48169/1‑249 (MQ=255)
                                               cacaCTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCg                                                                                                                                                                                   >  1:211766/1‑251 (MQ=255)
                                                            ggCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGtt                                                                                                                                                                      <  1:111484/251‑1 (MQ=255)
                                                                  gCTTATGGCGGCTTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTtgtgtg                                                                                                                                                                <  1:209069/251‑1 (MQ=255)
                                                                     tATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGTG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTgg                                                                                                                                                               <  1:136755/249‑1 (MQ=255)
                                                                                       aaTTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCaata                                                                                                                                               <  2:61858/247‑1 (MQ=255)
                                                                                           tttCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGAtgt                                                                                                                                        <  1:211718/250‑1 (MQ=255)
                                                                                                      gTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTAc                                                                                                                            >  1:94425/1‑251 (MQ=255)
                                                                                                       tAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACt                                                                                                                           >  1:124688/1‑251 (MQ=255)
                                                                                                       tAAAATATAGCGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGACGCACGTTGTTTATAATTCGACACCGGATAACCAGCGGCCGGTGGCGCTTTCTTTGACGCCTGTGCACAGTGCGGGCGGGGGAG‑CGATTAAA‑GCTGGATCATTAATTGCCGTGATTATTTTACGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGCATAATTACTCCAATAATGATGTGGTGGTGCCTACt                                                                                                                           >  1:179596/1‑251 (MQ=255)
                                                                                                                                    ccTACCCCCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGATGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGCCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTTTGCACTTGATg                                                                                               <  2:24770/250‑1 (MQ=255)
                                                                                                                                                gAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTCGAGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGCGGCGG‑CGACTACA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAAAAGCGATGTTTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACCGGCGGCTGGGATGTTTCTGCTCGAGATGTCACCGTTActct                                                                                  >  1:250274/1‑251 (MQ=255)
                                                                                                                                                       cgcgcgTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTTACCGTTACTCTGCCGGAc                                                                           >  2:5927/1‑251 (MQ=255)
                                                                                                                                                                  tataAGTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGtt                                                                <  2:124688/250‑1 (MQ=255)
                                                                                                                                                                         cGAGAACGGATAAGCCGTGGCCGGTTGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGTTCAGTGcc                                                         <  1:161855/250‑1 (MQ=255)
                                                                                                                                                                              aCGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGGTTAAA‑GCTGGCTCATTACTTGCCGTGCTTATTTTGAGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGAACT‑CCCCTGGTTCCGTGCCCATTcc                                                   >  1:175306/1‑251 (MQ=255)
                                                                                                                                                                               cGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGGGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCAAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTTTGCTCGAGATGACACCGTTACTCTGCCGGACCCACCC‑GGTTCAGTGCCAATTccc                                                  >  2:25647/1‑250 (MQ=255)
                                                                                                                                                                                            gCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTAtt                                     <  2:211766/251‑1 (MQ=255)
                                                                                                                                                                                             ccGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGGCGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGGCT‑ACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTAtt                                     >  2:122781/1‑251 (MQ=255)
                                                                                                                                                                                               ggTGGCGCTTTATTAGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTACAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTTTtgtg                                  >  1:214338/1‑251 (MQ=255)
                                                                                                                                                                                                          aTTTGACGCCTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCa                        >  1:29167/1‑250 (MQ=255)
                                                                                                                                                                                                            ttGACGCCTGTGAGCTGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACG‑ACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCaaa                      <  1:107232/250‑1 (MQ=255)
                                                                                                                                                                                                                  ccTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGTTCAGTGCCAATTCCTCTTAc                                             <  2:180570/221‑1 (MQ=255)
                                                                                                                                                                                                                  ccTGTGAGCAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGTTCAGTGCCAATTCCTCTTAc                                             >  1:180570/1‑221 (MQ=255)
                                                                                                                                                                                                                          cAGTGCGGGCGGGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTTCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGGTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTAcc       >  1:139474/1‑251 (MQ=255)
                                                                                                                                                                                                                              gcgggcggGGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGGATTACCtctc   >  1:198721/1‑251 (MQ=255)
                                                                                                                                                                                                                               cgggcgggGGGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTcc  <  1:5927/251‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
CAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGG‑CGATTAAA‑GCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACC‑AACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACT‑ACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCC  >  NC_000913/4548991‑4549464

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