Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,549,037 T→C V77A (GTC→GCC)  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,0370TC96.9% 91.2 / ‑5.5 32V77A (GTC→GCC) fimHtype 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (15/16);  minor base A (1/0);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

TGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACT  >  NC_000913/4548804‑4549279
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
tGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAAACTTGCGCCCGTCGTGGATGTGGGGCAAAACCTTGTCGTGGATCTTTTGACGCAAAACTTTTTCCATAACGATTTTTCGGAAAACCTTACAGACTATGCCACACTGCAACGAGGCt                                                                                                                                                                                                                                   >  1:219508/1‑251 (MQ=255)
  tAATGAAACGATTAATTCCCCTTTTTGCGGTACTGCTGATGGGCGGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCGATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTc                                                                                                                                                                                                                                  <  2:150595/250‑1 (MQ=255)
  tAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCCTCGTGGATGTGGGGCAAAAACTGGGCGTGGGTCTTTCGACGCAAATCTTTTGCCATAACGATTATTCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCg                                                                                                                                                                                                                                 >  1:189938/1‑251 (MQ=255)
      aaaaCGAGTTATTCCCCTGTTTGCTGTCCTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATCCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCtt                                                                                                                                                                                                                             <  2:92607/250‑1 (MQ=255)
            aGTTATTACACTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCAAGTTTATGTAAACCTTTCGCCCGGCGGGAATGTGGGGCAAAACCTGGTCGGGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATTCGGAAAACATTAAAGAATATTCCACACttgaaacaggcgcgggctatagcg                                                                                                                                                                                                                       >  1:65237/1‑228 (MQ=255)
                   aCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTTTGTAAACCTTGCGCCCGGCGGGAATGTGGGGCAAAACCTGGTCGTGGGTCTTTCGACCCCAATCTTTTTCCATAACGATTATTCGGAAACCATTACAGACTATGCCACACTGCAACGAAGCTCGGCTTTTGgcggcggggt                                                                                                                                                                                                                >  1:182897/1‑247 (MQ=255)
                                           ggCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCGCCCATTTTTATTTAAAACTCTCGCCCCCCGTGAATGTGGGGCAAAACCTGGGCGTGGATCTTTCGCCGCAAAACTTTTTCCCTAACGCTTATTCGGGAAACCTTAAAGACCATGCCCCACTGCAACCAGGCTCCGCTTATTGCGGGGTGTTATTTCATTTTTTCGGGGCCggta                                                                                                                                                                                         >  2:116931/1‑247 (MQ=255)
                                           ggCTGGTCGGAAAATGCCTGGTCATTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTaaa                                                                                                                                                                                        <  1:129430/251‑1 (MQ=255)
                                               ggTCGGAAAATGCCGGGTAATTCGCCTTAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAGCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAata                                                                                                                                                                                     <  1:184810/250‑1 (MQ=255)
                                                gTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGACCCAAAGATTTTTTAAAACCTGCGCCCGTCGTGGATGTGGGGGAAAACCTGGGTGGGGGTTTTTCGACGCAAATCTTTTTCCATAACGATTATCCGGGAACCATTAACGAATATGCCACACCGCAACGAAGCTCGGCTTATTGCGGGGGGTTATCTAAtttttttcggagcctgaaaaat                                                                                                                                                                                    >  2:185571/1‑233 (MQ=255)
                                                          gcctggtctttcgctacaaacaCCGCCAATGGTACAGCTAGCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAGCCGGAAACCATTACAGACTATGCCACACTGCAACTATGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAg                                                                                                                                                                         <  2:97301/230‑1 (MQ=255)
                                                               gTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTTTGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAAACTGGGCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAAACATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATTGCGGCGTGTTTTCTAATTTTTTCGGGACCGTAAAATATAGTGGCAGGAGCTAtt                                                                                                                                                                    >  1:186385/1‑250 (MQ=255)
                                                                            aaaaCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTaccac                                                                                                                                                       <  1:167231/251‑1 (MQ=255)
                                                                               aCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTATGGAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAAACTGGTCGGGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAAACATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTTCTGGACCGCAAAATATAGTGGCAGTAGCTCTTCATTTTCTtccaccc                                                                                                                                                     >  1:205022/1‑249 (MQ=255)
                                                                                  gCCAATGGTACCGCTATCCCTATTGGCGGTGGGAACGCCAACGTTTTTTTAAACCTTCCGCCCGTCGTGAATGTGGGGCAAAAACTGGGCGTGGCTCTTTTGCCGCGAAACTTTTTGCCTAAAGATTATCCGGGAACCATTTCAAACCATTACAAACAGCAACGAGGCTCGGCTTATGGCGGCGTGTTTTTTAATTTTTTCGGGGCCGTGAAACATATTGGCAGGAGcgttcccttttttccccccccacc                                                                                                                                                 >  2:103635/1‑228 (MQ=255)
                                                                                                       aTTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTtata                                                                                                                             <  1:70707/250‑1 (MQ=255)
                                                                                                          ggCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTTGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAAtt                                                                                                                          <  1:195933/250‑1 (MQ=255)
                                                                                                                            gTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCGACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGt                                                                                                        <  1:277627/250‑1 (MQ=255)
                                                                                                                                aTGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGcc                                                                                                    <  1:79411/250‑1 (MQ=255)
                                                                                                                                    aaaCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCACATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAACATATAGTGGCAGTAGCTATCCATTTCCTACCACCACCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTgg                                                                                               >  2:121160/1‑251 (MQ=255)
                                                                                                                                    aaaCCGTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACTGATAAGCCGTGGCCGGTgg                                                                                               <  2:65245/251‑1 (MQ=255)
                                                                                                                                                                        ggTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCgggg                                                           <  1:227898/251‑1 (MQ=255)
                                                                                                                                                                               gATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGAAGCCTGTGAGCAGTGCGGGCGGGGTGGCGAt                                                    >  1:22809/1‑251 (MQ=255)
                                                                                                                                                                                          gCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCCGCGAAAAGCCGCGCGTTGTTTATAATTCGAGAACGGAGCAGCCGAGGACGGTGGCGCTTTATTTGACGCCTGTGAGCCGTGCGGGCGGGGGGGCGATTAAAGCTGGcc                                         >  2:92036/1‑250 (MQ=255)
                                                                                                                                                                                                         aTAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTg                           <  1:53329/251‑1 (MQ=255)
                                                                                                                                                                                                               aTTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTCTTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTAtt                     >  1:40957/1‑251 (MQ=255)
                                                                                                                                                                                                                ttATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCAGGGCTTAttt                    >  2:23875/1‑251 (MQ=255)
                                                                                                                                                                                                                             gaTTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGAGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGGAGGGTGGCGATTTAAGCTGGCTCATTCATTGCCCTGCTTATTTTGCGACAGAcc         >  2:150949/2‑250 (MQ=255)
                                                                                                                                                                                                                                tACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCtg                                                                             <  2:209711/179‑1 (MQ=255)
                                                                                                                                                                                                                                tACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGGTGTTTATAATTCGAGAACGGAAAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCtg                                                                             >  1:209711/1‑179 (MQ=255)
                                                                                                                                                                                                                                 aCAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTGTAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACcaacaa    <  2:233120/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACt  <  1:197725/251‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
TGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACT  >  NC_000913/4548804‑4549279

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