Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,549,061 A→G Y85C (TAT→TGT)  fimH → minor component of type 1 fimbriae

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,0610AG98.3% 202.5 / NA 58Y85C (TAT→TGT) fimHminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (29/28);  minor base T (1/0);  total (30/28)
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

GTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGG  >  NC_000913/4548928‑4549183
                                                                                                                                     |                                                                                                                          
gTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                       <  1:1248668/139‑1 (MQ=255)
              tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAAt                                                                                                         >  2:1260158/1‑139 (MQ=255)
              tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAAt                                                                                                         >  1:489638/1‑139 (MQ=255)
                      tcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                 <  1:218923/139‑1 (MQ=255)
                        gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                               <  1:1137448/139‑1 (MQ=255)
                        gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                               <  2:489638/139‑1 (MQ=255)
                                 gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGATCGTAAAat                                                                                      >  1:610830/1‑139 (MQ=255)
                                 gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAat                                                                                      >  2:505714/1‑139 (MQ=255)
                                    cAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                   <  1:1107579/139‑1 (MQ=255)
                                      aaaCCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGt                                                                                 <  1:542285/139‑1 (MQ=255)
                                                   gATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                       <  2:816813/88‑1 (MQ=255)
                                                   gATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                       >  1:816813/1‑88 (MQ=255)
                                                   gATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATc                                                                    >  2:1062489/1‑139 (MQ=255)
                                                   gATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATc                                                                    >  2:1210072/1‑139 (MQ=255)
                                                                       ttGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAg                                                                            <  1:489353/111‑1 (MQ=255)
                                                                       ttGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAg                                                                            >  2:489353/1‑111 (MQ=255)
                                                                        tGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                 <  2:1253713/121‑1 (MQ=255)
                                                                        tGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                 >  1:798067/1‑121 (MQ=255)
                                                                        tGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                 <  2:798067/121‑1 (MQ=255)
                                                                        tGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                 >  1:1253713/1‑121 (MQ=255)
                                                                              aaCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAcgccgc                                         <  1:856182/139‑1 (MQ=255)
                                                                                 gATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                      <  2:688347/139‑1 (MQ=255)
                                                                                 gATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                      <  1:505714/139‑1 (MQ=255)
                                                                                    tATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTTTggcg                                                                                                                       >  1:325836/1‑55 (MQ=255)
                                                                                    tATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                       <  2:325836/55‑1 (MQ=255)
                                                                                    tATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                   <  1:928208/91‑1 (MQ=255)
                                                                                    tATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                   >  2:928208/1‑91 (MQ=255)
                                                                                    tATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTgg                                                                               <  1:198138/95‑1 (MQ=255)
                                                                                    tATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTgg                                                                               >  2:198138/1‑95 (MQ=255)
                                                                                     aTCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                                                                              <  2:849666/95‑1 (MQ=255)
                                                                                     aTCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                                                                              >  1:849666/1‑95 (MQ=255)
                                                                                         ggAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTtat                              <  1:1260158/139‑1 (MQ=255)
                                                                                                 ttACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                               <  1:638503/66‑1 (MQ=255)
                                                                                                 ttACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                               >  2:638503/1‑66 (MQ=255)
                                                                                                  tACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                 <  2:901897/95‑1 (MQ=255)
                                                                                                  tACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                 >  1:901897/1‑95 (MQ=255)
                                                                                                         tATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                               <  2:834974/58‑1 (MQ=255)
                                                                                                         tATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                               >  1:834974/1‑58 (MQ=255)
                                                                                                         tATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                   <  1:934449/70‑1 (MQ=255)
                                                                                                         tATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                   >  2:934449/1‑70 (MQ=255)
                                                                                                           tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                               >  1:667235/1‑56 (MQ=255)
                                                                                                           tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                               <  2:667235/56‑1 (MQ=255)
                                                                                                           tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                  <  1:645653/85‑1 (MQ=255)
                                                                                                           tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                  >  2:645653/1‑85 (MQ=255)
                                                                                                           tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAg            >  2:1323029/1‑139 (MQ=255)
                                                                                                           tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAg            >  1:4885/1‑139 (MQ=255)
                                                                                                                 acTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTgg      >  2:970993/1‑139 (MQ=255)
                                                                                                                     cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                 <  2:403633/44‑1 (MQ=255)
                                                                                                                     cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                 >  1:403633/1‑44 (MQ=255)
                                                                                                                     cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                                                                              >  1:191662/1‑63 (MQ=255)
                                                                                                                     cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                                                                              <  2:191662/63‑1 (MQ=255)
                                                                                                                     cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCgg  >  2:684633/1‑139 (MQ=255)
                                                                                                                     cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCgg  >  2:36502/1‑139 (MQ=255)
                                                                                                                     cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCgg  >  2:153412/1‑139 (MQ=255)
                                                                                                                           ggCTCGGCTTGTGGTGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                  <  2:1090501/69‑1 (MQ=255)
                                                                                                                           ggCTCGGCTTGTGGTGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                  >  1:1090501/1‑69 (MQ=255)
                                                                                                                               cGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcc                                      <  2:839111/93‑2 (MQ=255)
                                                                                                                               cGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTCCCTACCACCAGCGAAACGCcgcgcg                                      >  1:839111/1‑93 (MQ=255)
                                                                                                                                     |                                                                                                                          
GTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGG  >  NC_000913/4548928‑4549183

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