Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,546,601 C→T I502I (ATC→ATT fimD → type I fimbriae usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,546,6010CT97.3% 118.2 / ‑3.7 37coding (1506/2637 nt)fimDtype I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/15);  minor base . (1/0);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAAC  >  NC_000913/4546381‑4546847
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
aaagcATGGGGGCACTGGGCGCTCTGTCTGCGGATATGGCGCAGGCTGATTCCACTCTTCCCGATGACAGTCAGCATGGCGGACAATCGGTGCGTTTTCTCTATGACAAATCGCGCAATGAATCAGGCACGAATATGCAGTTAGTGGGTTGCCGTTTTCCGACCAGCGGATATGTTAATTGCGCTGATACAACATACAGTCGAATGAATGGCTACAACACTGAAACACAGGACGGAGTTAGTCAGGTTAAg                                                                                                                                                                                                                          <  1:360127/247‑1 (MQ=255)
 aaaCATGGGGGCACTGGGTGCTCTGTCTGTGGATATGACGCAGGCTAAGTCCTCACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGc                                                                                                                                                                                                                         <  2:319718/251‑1 (MQ=255)
     aTGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTGATGCAGGTTAAGCCGaa                                                                                                                                                                                                                     <  1:121351/251‑1 (MQ=255)
        ggggCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGCCGGACAATCGGCGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                                                                                                                                                                                                  >  1:159926/1‑251 (MQ=255)
                    gcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                                                                                                                                                                                       >  1:349331/1‑250 (MQ=255)
                                gATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTta                                                                                                                                                                                          <  1:72458/251‑1 (MQ=255)
                                gATATGACGCAGGCTAATTCCACACTGCCCGATGACAGTCAGCATGACGGACAATCGGTGCGGTTTCTCTATAGCAAATCGCTCAATTAATCAGGTACGAAGAGTCAGTTAGTGGGTGACCGTTATTCGTCCAGCGGATATTTTAATTTCGTTGAGACAACATACAGTCGGATGTATTTCTACAACATTGAAACACAGGACTGAGTTATTGAGTTTATGTCGAAATGCACCGACTATTACAACCTCGCTta                                                                                                                                                                                          <  2:412070/251‑1 (MQ=255)
                                           ggCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGGCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcgcg                                                                                                                                                                               <  2:186089/251‑1 (MQ=255)
                                                                      tCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGACTCAGGCACGACTATTCAGTTAGTGGGTTACCGTTATTCGAACAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGc                                                                                                                                                    >  1:110140/1‑251 (MQ=255)
                                                                        aGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCa                                                                                                                                                   >  1:137161/1‑250 (MQ=255)
                                                                                        ggTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGACACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACCGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGCCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATcc                                                                                                                                   >  2:223134/1‑249 (MQ=255)
                                                                                                ttCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACATCTCCCCGTTACTCAGCAACCCTGGCGCCCATCAACACTGtg                                                                                                                           <  2:89369/250‑2 (MQ=255)
                                                                                                    ctATAACAAATCGCTCAATGAATCAGGCACGAATCTTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCa                                                                                                                                                                                                                >  1:113151/1‑161 (MQ=255)
                                                                                                    ctATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCa                                                                                                                                                                                                                <  2:113151/161‑1 (MQ=255)
                                                                                                    atataACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGGTGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCTGGAAATTTCAACGCACCGTTACTCAGCAACTCTGGCGCACATCATCACGGTAtg                                                                                                                         <  2:382146/247‑2 (MQ=255)
                                                                                                           aaaTCGCTCAATGCAAAAGGCACGAATAATCAGTTAGTGGGTTAACGCTATCAGAACAGCGGATCTTTTAATTTCGCTGATACAACATCCAGTCGAATGAATGGCTACAACATTGAAACACCGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATCACACACGCGGGAAATTACAACTCACCGTTACTCCGCAACTCGGGCGCACATCAACACTGTATTTGAGTgga                                                                                                                 >  2:241398/1‑248 (MQ=255)
                                                                                                            aaTCGCTCAATGACTAAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGCCGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACATCGCTTATAACAAACGAGGAAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGGGTGGTAGc                                                                                                              >  2:181714/1‑251 (MQ=255)
                                                                                                             aTCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAACCACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                                                                             >  2:373568/1‑251 (MQ=255)
                                                                                                              tCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAAATCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                                                                             >  1:120483/1‑250 (MQ=255)
                                                                                                                gCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATc                                                                                                          <  1:10962/251‑1 (MQ=255)
                                                                                                                                          aGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTc                                                                                >  1:30894/1‑251 (MQ=255)
                                                                                                                                               gTGGGTTACCGTTATTCGACCAGCTGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGa                                                                           <  1:62402/251‑1 (MQ=255)
                                                                                                                                                 gggTTACCGATATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGTTGAGc                                                                         >  1:320587/1‑251 (MQ=255)
                                                                                                                                                   gTTACCGTTATTCGACCGGCGGATATTTTAATTTCGCTGATACAACATACAGCCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATACAACCTCGCTTATAACAAACGAGGGCAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGCGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGAAAt                                                                      >  2:191559/1‑251 (MQ=255)
                                                                                                                                                                gACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAt                                                          >  1:6042/1‑251 (MQ=255)
                                                                                                                                                                       ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATAc                                                   >  1:92836/1‑251 (MQ=255)
                                                                                                                                                                           atTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGc                                                <  1:10841/250‑1 (MQ=255)
                                                                                                                                                                           atTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCg                                               <  1:291771/251‑1 (MQ=255)
                                                                                                                                                                                aaTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCg                                           >  1:337527/1‑250 (MQ=255)
                                                                                                                                                                                       cTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCAGAACTTCACAGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATc                                   >  1:165793/1‑251 (MQ=255)
                                                                                                                                                                                         gATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCTAAGATATCa                                  <  1:316993/250‑1 (MQ=255)
                                                                                                                                                                                                     aGTCGAATTAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGGTAAGCCGCAATTCACCGACAATACCAACCTCGCTTATAACAAACGCGAGAAATTACACATCACCGTTACTCAGCAACTCGGGAGAACATCAACACTTAATTTGAGTGGTAGCCATAAGACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCACGGCTGGATTAAATACTGCGATCTATGATATCAACTGGACGCTc                       >  2:5176/1‑249 (MQ=255)
                                                                                                                                                                                                      gTCGAATGAATGGCAACAACGTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTAGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCCTCAAACTTATTGGGGAACGAGTACTGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGc                    >  2:149936/1‑251 (MQ=255)
                                                                                                                                                                                                        cGAATGAATGCCTGCAACCATGAAACACAGGACGGAGATATTCAGGTTAAGCCGAAATTCACCGACTATTACACCATTGCTTATAACAAACGCGGGCAATTAAAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAACCATCAAACTTATTGGGGAACGAGTAATGTCGCTGAGCAATTCCAGGCTGGATTAAATACTGCGTTAGAAGATATCAACTGAAAGCTCAGCt                   >  2:279259/1‑250 (MQ=255)
                                                                                                                                                                                                             gaatgGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCAAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGcc             >  2:68200/1‑251 (MQ=255)
                                                                                                                                                                                                                     acaacaTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaa      >  1:180227/1‑250 (MQ=255)
                                                                                                                                                                                                                        acaTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAAc  <  2:349331/251‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
AAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAAC  >  NC_000913/4546381‑4546847

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