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,0370TC93.8% 87.5 / ‑3.1 32V77A (GTC→GCC) fimHtype 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (14/16);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

AATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGC‑GGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACT  >  NC_000913/4548807‑4549279
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
aaTGAAACGAGTTATTACCCTGTTTGCTGTCCTGTGGAGGGGCTGGTCGGTAAATGCTTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGA‑GGCTCgg                                                                                                                                                                                                                                  <  1:66429/251‑1 (MQ=255)
     aaCGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATTTAAACCTTTCGCCCGTCGCGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCTACGCAAATCTTTTTCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCCCGGCTTAt                                                                                                                                                                                                                             >  1:329167/1‑251 (MQ=255)
           ttattaCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGGGGCACCACCCATGTTTATGTAAACCTTGCGCCCCTCGGGAATGTGGGGCAAAACCTGGGCGGGGATCTTTCGACGCAAATTTTTTGCCAAAACGATTTTCCGGGAAACATTAAAGAATATTCCACACTGCAACGAAGGC‑CGGGTTATggcagc                                                                                                                                                                                                                       >  1:118445/1‑248 (MQ=255)
                       ttccTGTCCTGCTGAGGGGCTGGTCGGTAAATGCCTGGTCTTTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATc                                                                                                                                                                                                               <  1:259352/244‑1 (MQ=255)
                       ttGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCAAGTTTTTGTAAACCCTGCGCCCGTCGTGAATGTGGGGCAAAAACTGGGCGTGGATCTTTCGGCGCAAATCTTTTGCCATAACGGTTATCCGGGAACAATTTAAGACTATGCCACACTGCAACGA‑GGCCCGGCTTATGGCGGCGGGTTATCtcat                                                                                                                                                                                                           >  1:140517/1‑248 (MQ=255)
                                    ggtgggtgggTCGGTAAATGCGTGGTCTTCCGCCTGAAAACCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCg                                                                                                                                                                                              <  2:89819/243‑1 (MQ=255)
                                                   aaaTCCCTGTTCATCCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCa                                                                                                                                                                               <  1:117481/251‑1 (MQ=255)
                                                                         aaaaCCGCCAATTGTACCGCTAACCCTATTGGCGGTGTCAGTGCCAATGTTTATGTAAACCTTGCGTCCGTCGTGAATGTGGGGCAAATCCTGGTCGTTGATCTTTCGACGCAAATCTGTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGTCATTAGCTATCCATTTCCTaccac                                                                                                                                                         <  2:187913/251‑1 (MQ=255)
                                                                                     ggTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAAcgcc                                                                                                                                            <  2:329167/251‑1 (MQ=255)
                                                                                              aTCCCTATTGGCGGTGGCAGCGCCAAGTTTTATGTAAACCTTGCGCCCGTCGTGAAGGTGGGGCAAAACCTGGTCGAGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCgttg                                                                                                                                    <  1:502094/250‑1 (MQ=255)
                                                                                                            tGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTACAATATAGTGGCAGTAGCTATCCATTTCCTAACACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCgtgg                                                                                                                      >  1:71766/1‑247 (MQ=255)
                                                                                                             ggCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAAc                                                                                                                    >  1:436167/1‑251 (MQ=255)
                                                                                                                      aaTGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGCTCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTTACCATTTCCTACCACACGCG‑AAACGCCGCGCGTTGTTTTTTATTCTTGAACGGATaaacc                                                                                                           >  2:213268/1‑248 (MQ=255)
                                                                                                                               gTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGAGGGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGGTGCTATCCATTTTCTACACCCCGCG‑GAACCCCGCGCGTTGTTTAAAATTCGCGTACGGATAAACCGGGGccagg                                                                                                  >  2:172221/1‑248 (MQ=255)
                                                                                                                               gTAAACCTTGCGCCAGCCGTGTATGGGGGGCAAAACCTGGTCGTGGATCTTTCGACGAAAATCTTTGGCCATAACGATTATCCGGAAACCATTACAGGCTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGGGGCTGTAGCTATCCATTTCCGCCCACCAGCG‑AAACGCCGCGCGTTGTTTATAATCTGAGAACGGATAAGCCGTGGCCGGt                                                                                                  <  2:463277/251‑1 (MQ=255)
                                                                                                                                 aaatctttcgcCCTTCGTGAAGGTGGGGCAAAACCTGGTCGTGGATCTTTCGCCGCCAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGGTTATAATTCGAGAACGGATAAGCCGTGGCCTGTgg                                                                                                <  2:530129/243‑1 (MQ=255)
                                                                                                                                          gcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGGAGAGGCTGACCATTTCCTACCACCAGCG‑AAACGCCGCGCGGTGTTTTTAATTCGAGAATGGAAAACCCGTGGCCGGTGGGGCtttttt                                                                                       >  2:121140/1‑248 (MQ=255)
                                                                                                                                                   gAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGAAGCTATCCATTTCCTACCACCAGCG‑AAAAGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCtg                                                                              >  1:52380/1‑251 (MQ=255)
                                                                                                                                                            gCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACAGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGACGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTg                                                                     <  1:147682/251‑1 (MQ=255)
                                                                                                                                                              aaaaCCGGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTgcg                                                                   <  2:140517/251‑1 (MQ=255)
                                                                                                                                                                     ggTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGcgggc‑ggg                                                            <  1:172221/250‑1 (MQ=255)
                                                                                                                                                                                 ttCGACGCAAATCTTTTGCCATAACGAGTCTCCGGAACCCATTACCGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAGACAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTGTATTTGACGCCTGTGAGCTGTGCGGGC‑GGGGTGGCGag                                                    <  2:104008/246‑2 (MQ=255)
                                                                                                                                                                                  tCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGC‑GGGGTGGCGATTAAAg                                               <  2:183802/250‑1 (MQ=255)
                                                                                                                                                                                             cTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAACATATAGTGGCAGTAAGTATCCATTTCCTACCACCCGCGCAAAC‑CCGCGCGTTGTTTTTAATTCCAGAACGGATAAGCCGTGGCCCGGGGCGCTTTATTTGACGCCTGTAAGCAGTGCGGGC‑GGGGGGGCGATTAAAGCGGGcgtatta                                    >  2:424024/1‑244 (MQ=255)
                                                                                                                                                                                                 tGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGC‑GGGGTGGCGATTAAAGCTGGCTCATTAATTg                                <  1:175829/250‑1 (MQ=255)
                                                                                                                                                                                                   ccATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAAGTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCGATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGC‑GGGGTGGCGATTAAAGCTGGCTCATTAATTGCCg                             <  2:172075/251‑1 (MQ=255)
                                                                                                                                                                                                    cATAACGATTATCCGGAAACCATTACGGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGAGTTATCTAATTTTTCCGGGAACGTAAAATATAGGGGCAGTAGCTATCCATTTCCTATCACCAGCG‑AAACGCCGAGCGTTGTTTTTAATTCGAAAACGGATAAGCCCTGGACGGGGGCGCTTTATTTGACGCCTGTGAGCCATGCGGGCAGGGTTGGTG‑TTTAAGCTGGCTCATTAAtttcctg                            >  2:22322/1‑246 (MQ=255)
                                                                                                                                                                                                         cgtttATCCGGAAACCATTACAGTCTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGTCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGC‑GGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTa                       <  2:475435/248‑1 (MQ=255)
                                                                                                                                                                                                               gccGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGTGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGC‑GGGGTGGCGAGTAAAGTTGGCTCATTAATTGCCGTGCTTATTTTg                  <  1:483243/249‑1 (MQ=255)
                                                                                                                                                                                                                 cGGAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACCCCGCGCGTTGTTTATAATTCGAGAACGGGTAAGCCGTGGCCGGGGGCGCTTTATTTGACGCCTGTGGGCAGTGCGGGC‑GGGGTGGCGGTTAAAGCTGGCGCATTAAtttgcgtggtttttttgggg               >  2:50186/1‑233 (MQ=255)
                                                                                                                                                                                                                  ggAAACCATTACAGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGCGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGC‑GGGGTGGCGAATAATGCTGGCTCATTAATTGCCGTGCTTGTTTTGCGa               >  1:81818/1‑250 (MQ=255)
                                                                                                                                                                                                                              aGACTATGCCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGC‑GGTGGGGCGGTTAAAGCTGGCTCATTAATTGCCGTGCTTTTTTTGCGACAGACCAACAACt  >  1:8621/1‑251 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
AATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGA‑GGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCG‑AAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGC‑GGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACT  >  NC_000913/4548807‑4549279

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