Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,540,984 (A)7→6 coding (28/603 nt) fimB → tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,540,9780A.98.0% 194.8 / ‑2.3 50coding (22/603 nt)fimBtyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (24/25);  total (24/26)
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.97e-01

TAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAA  >  NC_000913/4540849‑4541107
                                                                                                                                 |                                                                                                                                 
tagttaTTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACaaaaaaggaa                                                                                                                          <  2:338371/139‑5 (MQ=255)
    tattaAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTc                                                                                                                     <  1:876014/139‑1 (MQ=255)
            aaaCGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGAcc                                                                                                               <  2:770629/137‑1 (MQ=255)
            aaaCGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGAcc                                                                                                               >  1:770629/1‑137 (MQ=255)
                            aaaCATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                             <  1:1208156/139‑1 (MQ=255)
                                 tAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCt                                                                                        <  2:944536/139‑1 (MQ=255)
                                  aaTCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAAt                                                                                               >  2:302508/1‑131 (MQ=255)
                                  aaTCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAAt                                                                                               <  1:302508/131‑1 (MQ=255)
                                    tCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCt                                                                                        >  2:544858/1‑136 (MQ=255)
                                    tCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCt                                                                                        <  1:544858/136‑1 (MQ=255)
                                     cAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTaaa                                                                                    <  1:506603/139‑1 (MQ=255)
                                        gATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACaaaaaaggaa                                                                                                                          <  2:1475536/99‑5 (MQ=255)
                                        gATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACaaaaaaggaa                                                                                                                          >  1:1475536/1‑95 (MQ=255)
                                         aTAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAGAATAAGGCTGATAACaaaaaagg                                                                                                                            <  1:245646/95‑3 (MQ=255)
                                         aTAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAGAATAAGGCTGATAACaaaaaagg                                                                                                                            >  2:245646/1‑93 (MQ=255)
                                                tGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACaaaaaaggaa                                                                                                                          <  2:1329883/91‑5 (MQ=255)
                                                tGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACaaaaaaggaa                                                                                                                          >  1:1329883/1‑87 (MQ=255)
                                                  ttGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATa                                                                                                           <  2:924372/103‑1 (MQ=255)
                                                  ttGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATa                                                                                                           >  1:924372/1‑103 (MQ=255)
                                                  ttGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCg                                                                       >  1:243837/1‑139 (MQ=255)
                                                  ttGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCg                                                                       >  1:232862/1‑139 (MQ=255)
                                                   tGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCgg                                                                      <  2:130019/139‑1 (MQ=255)
                                                   tGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCgg                                                                      <  1:743521/139‑1 (MQ=255)
                                                            gCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgc                                                             <  2:243837/139‑1 (MQ=255)
                                                              tAACCCAGCACAGCTAGTGCGCGTCTGTAGTTATAAGGGAAAAAAGATGAAGAATAAGGCTGGTAACAAAAAAAGGAACTTCCTTACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgca                                                            <  2:1286729/139‑1 (MQ=255)
                                                              tAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgcag                                                           <  2:258810/139‑1 (MQ=255)
                                                                    aGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAg                                                                                                          <  2:487392/86‑1 (MQ=255)
                                                                    aGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAg                                                                                                          >  1:487392/1‑86 (MQ=255)
                                                                     tcacaGCTAGTGCGCGTCTGTAATTGTAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTg                                                                                                                  <  1:922677/76‑1 (MQ=255)
                                                                     gCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTg                                                                                                                  >  2:922677/1‑77 (MQ=255)
                                                                     gCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCa                                                                                                             <  1:1150433/82‑1 (MQ=255)
                                                                     gCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCa                                                                                                             >  2:1150433/1‑82 (MQ=255)
                                                                     gCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgcagca                                                         >  1:785227/1‑134 (MQ=255)
                                                                     gCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgcagca                                                         <  2:785227/134‑1 (MQ=255)
                                                                           cTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGt                                                                                                         >  2:579010/1‑80 (MQ=255)
                                                                           cTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGt                                                                                                         <  1:579010/80‑1 (MQ=255)
                                                                            tAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGt                                             >  1:1032647/1‑139 (MQ=255)
                                                                                   cgTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGGATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGActt                                      >  2:914579/1‑139 (MQ=255)
                                                                                   cgTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGActt                                      >  1:1315849/1‑139 (MQ=255)
                                                                                     tCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGActttg                                    >  1:643358/1‑139 (MQ=255)
                                                                                         tAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACg                                                       <  1:380084/116‑1 (MQ=255)
                                                                                         tAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACg                                                       >  2:380084/1‑116 (MQ=255)
                                                                                                aaGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgcagca                                                         >  1:948990/1‑107 (MQ=255)
                                                                                                aaGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATgcagca                                                         <  2:948990/107‑1 (MQ=255)
                                                                                                 aGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGc                                                                    >  1:95141/1‑95 (MQ=255)
                                                                                                 aGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGc                                                                    <  2:95141/95‑1 (MQ=255)
                                                                                                   ggAAAAACGTTGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAgcagca                                                                              <  2:13704/83‑1 (MQ=255)
                                                                                                   ggAAAAACGTTGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAgcagca                                                                              >  1:13704/1‑83 (MQ=255)
                                                                                                   ggAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCa                      <  2:232862/139‑1 (MQ=255)
                                                                                                   ggAAAAACGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCa                      >  1:407048/1‑139 (MQ=255)
                                                                                                        aaCGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGaaa                                                                                                     >  2:1401540/1‑55 (MQ=255)
                                                                                                        aaCGATGAAGAATAAGGCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGaaa                                                                                                     <  1:1401540/55‑1 (MQ=255)
                                                                                                                       ggCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                             <  2:1077871/48‑1 (MQ=39)
                                                                                                                       ggCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCa                                                                                             >  1:1077871/1‑48 (MQ=255)
                                                                                                                       ggCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAat                                                   >  2:635198/1‑90 (MQ=255)
                                                                                                                       ggCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAat                                                   <  1:635198/90‑1 (MQ=255)
                                                                                                                       ggCTGATAACAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGaaa  >  2:518024/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                 
TAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAA  >  NC_000913/4540849‑4541107

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