Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,542,692 +A intergenic (+59/‑423) fimE → / → fimA regulator for fimA/type 1 fimbriae major subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,542,6921.A85.0% 48.8 / 5.2 20intergenic (+59/‑423)fimE/fimAregulator for fimA/type 1 fimbriae major subunit
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (14/3);  total (16/4)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

TTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAA‑CTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAAT  >  NC_000913/4542459‑4542926
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          
ttGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATAACataaat                                                                                                                                                                                                                           >  1:230076/1‑251 (MQ=255)
                         ggtgCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCAt                                                                                                                                                                                                  >  1:162268/1‑251 (MQ=255)
                                            ttCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAACAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGAc                                                                                                                                                                               >  1:347510/1‑251 (MQ=255)
                                                          ggtagaataaaTATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAA‑CTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAAc                                                                                                                                                                 >  1:107062/9‑251 (MQ=255)
                                                           ggCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACt                                                                                                                                                                >  1:8556/1‑251 (MQ=255)
                                                            gCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTg                                                                                                                                                               >  1:267057/1‑251 (MQ=255)
                                                                    atatTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGAATAAAAATATTAAAAAAAGACAAAAAGCATCTAACTGTTTGatat                                                                                                                                                       >  2:290073/1‑251 (MQ=255)
                                                                                                     aTGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATCat                                                                                                                                                                                                                               >  2:248076/1‑146 (MQ=255)
                                                                                                     aTGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATCat                                                                                                                                                                                                                               <  1:248076/146‑1 (MQ=255)
                                                                                                                   gCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTatat                                                                                                        >  2:209049/1‑251 (MQ=255)
                                                                                                                                    aaaTACTCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCt                                                                                       >  1:328843/1‑251 (MQ=255)
                                                                                                                                                                      gagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAAAAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCTTAAAAATATTAAAAAACGACAAAAAGAATCTAACTGTTTTATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCACGACGCATCTTCATCAttcttatttca                                                      >  1:333684/1‑244 (MQ=255)
                                                                                                                                                                                  aaCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTc                                          >  1:264845/1‑250 (MQ=255)
                                                                                                                                                                                     ttATTGATAATAAAGTTAAAAAAACAAATAATCACAAGAAAATTGTGCCTAAA‑CTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAACAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATAGAAAATTAATTTCACATCACCTCAGCTATATGTAAAGCTAACGTTTCTGTGGCTAGACCCAGCTTCCTCATTCTTCTCTCAATAAACAACCTAATGc                                      >  2:195186/1‑251 (MQ=255)
                                                                                                                                                                                                  agataacaaaaCAAATAAATACAAAACA‑GTGGGGCCAAAACTATCAATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAGACAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTAATTCTTATCTCCAAAAACCCCCGCATGCAATATAAACATCt                         >  1:316580/8‑251 (MQ=255)
                                                                                                                                                                                                     tAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAACTGTCCATATAATAAATAAGATACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACACCCGCTATATGTAAAGCTACCGTTTCTTTGGATAGACGCCTCTACCTAATTCTTATATCCAAAAACCACCATATGCAATATAAACAtatct                       >  2:148509/1‑246 (MQ=255)
                                                                                                                                                                                                             cAAATAAATACAAGACAATTGGGGCCAAA‑CTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGa              <  1:11635/251‑1 (MQ=255)
                                                                                                                                                                                                                  aaataCAAGACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACa         <  1:383841/251‑1 (MQ=255)
                                                                                                                                                                                                                      aCAAGACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAg     <  1:209049/251‑1 (MQ=255)
                                                                                                                                                                                                                          gACAAT‑GGGGCCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAATTAAAGATAACAATAAAat  >  1:205736/1‑250 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          
TTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAA‑CTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAAT  >  NC_000913/4542459‑4542926

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