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

TGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACT  >  NC_000913/4548804‑4549279
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
tGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAAACTTGCGCCCGTCGTGGATGTGGGGCAAAACCTTGTCGTGGATCTTTTGACGC‑AAAACTTTTTCCATAACGATTTT‑TCGGAAAACCTTACAGACTATGCCACACTGCAACGAGGCt                                                                                                                                                                                                                                   >  1:219508/1‑251 (MQ=255)
  tAATGAAACGATTAATTCCCCTTTTTGCGGTACTGCTGATGGGCGGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAATCTTTTGCGATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTc                                                                                                                                                                                                                                  <  2:150595/250‑1 (MQ=255)
  tAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCCTCGTGGATGTGGGGCAAAAACTGGGCGTGGGTCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑TCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCg                                                                                                                                                                                                                                 >  1:189938/1‑251 (MQ=255)
      aaaaCGAGTTATTCCCCTGTTTGCTGTCCTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATCCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCtt                                                                                                                                                                                                                             <  2:92607/250‑1 (MQ=255)
            aGTTATTACACTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCAAGTTTATGTAAACCTTTCGCCCGGCGGGAATGTGGGGCAAAACCTGGTCGGGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑TCGGAAAACATTAAAGAATATTCCACACttgaaacaggcgcgggctatagcg                                                                                                                                                                                                                       >  1:65237/1‑228 (MQ=255)
                   aCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTTTGTAAACCTTGCGCCCGGCGGGAATGTGGGGCAAAACCTGGTCGTGGGTCTTTCGACCC‑CAATCTTTTTCCATAACGATTAT‑TCGGAAACCATTACAGACTATGCCACACTGCAACGAAGCTCGGCTTTTGgcggcggggt                                                                                                                                                                                                                >  1:182897/1‑247 (MQ=255)
                                           ggCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCGCCCATTTTTATTTAAAACTCTCGCCCCCCGTGAATGTGGGGCAAAACCTGGGCGTGGATCTTTCGCCGC‑AAAACTTTTTCCCTAACGCTTATTCGGGAAACC‑TTAAAGACCATGCCCCACTGCAACCAGGCTCCGCTTATTGCGGGGTGTTATTTCATTTTTTCGGGGCCggta                                                                                                                                                                                         >  2:116931/1‑247 (MQ=255)
                                           ggCTGGTCGGAAAATGCCTGGTCATTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTaaa                                                                                                                                                                                        <  1:129430/251‑1 (MQ=255)
                                               ggTCGGAAAATGCCGGGTAATTCGCCTTAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAAGCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAata                                                                                                                                                                                     <  1:184810/250‑1 (MQ=255)
                                                gTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGACCCAAAGATTTTTTAAAACCTGCGCCCGTCGTGGATGTGGGGGAAAACCTGGGTGGGGGTTTTTCGACGC‑AAATCTTTTTCCATAACGATTAT‑CCGGGAACCATTAACGAATATGCCACACCGCAACGAAGCTCGGCTTATTGCGGGGGGTTATCTAAtttttttcggagcctgaaaaat                                                                                                                                                                                    >  2:185571/1‑233 (MQ=255)
                                                          gcctggtctttcgctacaaacaCCGCCAATGGTACAGCTAGCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAG‑CCGGAAACCATTACAGACTATGCCACACTGCAACTATGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAg                                                                                                                                                                         <  2:97301/230‑1 (MQ=255)
                                                               gTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTTTGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAAACTGGGCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAAACATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATTGCGGCGTGTTTTCTAATTTTTTCGGGACCGTAAAATATAGTGGCAGGAGCTAtt                                                                                                                                                                    >  1:186385/1‑250 (MQ=255)
                                                                            aaaaCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTaccac                                                                                                                                                       <  1:167231/251‑1 (MQ=255)
                                                                               aCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTATGGAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAAACTGGTCGGGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAAACATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTTCTGGACCGCAAAATATAGTGGCAGTAGCTCTTCATTTTCTtccaccc                                                                                                                                                     >  1:205022/1‑249 (MQ=255)
                                                                                  gCCAATGGTACCGCTATCCCTATTGGCGGTGGGAACGCCAACGTTTTTTTAAACCTTCCGCCCGTCGTGAATGTGGGGCAAAAACTGGGCGTGGCTCTTTTGCCGCGAAACTTTTTGCC‑TAAAGATTAT‑CCGGGAACCATTTCAAACCATTACAAACAGCAACGAGGCTCGGCTTATGGCGGCGTGTTTTTTAATTTTTTCGGGGCCGTGAAACATATTGGCAGGAGcgttcccttttttccccccccacc                                                                                                                                                 >  2:103635/1‑228 (MQ=255)
                                                                                                       aTTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTtata                                                                                                                             <  1:70707/250‑1 (MQ=255)
                                                                                                          ggCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTTGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAAtt                                                                                                                          <  1:195933/250‑1 (MQ=255)
                                                                                                                            gTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCGACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGt                                                                                                        <  1:277627/250‑1 (MQ=255)
                                                                                                                                aTGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGcc                                                                                                    <  1:79411/250‑1 (MQ=255)
                                                                                                                                    aaaCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑ACATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAACATATAGTGGCAGTAGCTATCCATTTCCTACCACCACCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTgg                                                                                               >  2:121160/1‑251 (MQ=255)
                                                                                                                                    aaaCCGTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACTGATAAGCCGTGGCCGGTgg                                                                                               <  2:65245/251‑1 (MQ=255)
                                                                                                                                                                        ggTCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCgggg                                                           <  1:227898/251‑1 (MQ=255)
                                                                                                                                                                               gATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGAAGCCTGTGAGCAGTGCGGGCGGGGTGGCGAt                                                    >  1:22809/1‑251 (MQ=255)
                                                                                                                                                                                          gC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCCGCGAAAAGCCGCGCGTTGTTTATAATTCGAGAACGGAGCAGCCGAGGACGGTGGCGCTTTATTTGACGCCTGTGAGCCGTGCGGGCGGGGGGGCGATTAAAGCTGGcc                                         >  2:92036/1‑250 (MQ=255)
                                                                                                                                                                                                         aTAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTg                           <  1:53329/251‑1 (MQ=255)
                                                                                                                                                                                                               aTTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTCTTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTAtt                     >  1:40957/1‑251 (MQ=255)
                                                                                                                                                                                                                ttAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCAGGGCTTAttt                    >  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)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
TGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGC‑AAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACT  >  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: ‑