Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,542,683 Δ1 bp intergenic (+50/‑432) 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,6830T.93.3% 44.3 / ‑0.9 15intergenic (+50/‑432)fimE/fimAregulator for fimA/type 1 fimbriae major subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (9/5);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

TCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAA‑C‑TGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAA  >  NC_000913/4542441‑4542878
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                     
ccATATGTAAAGGCATGCTTCCGGTTATGCATAAGCCGAGCGTGCAGCAGATACTTGTTCAATTGAGGATTGTCTCGAGCCTCCAATCATTCACCATACTGTGCGTTGTCCCGCCCGTAATGCTCATCGTTATCCAGGACTCTTGGAAAGAAATAATCAAATAATAGAAAAATTCAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATATAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCaa                                                                                                                                                                                               >  2:31567/2‑250 (MQ=255)
      gTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCCTACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTACAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAAAAAATACAAGACAAT‑GGGGCCAAAAC‑TGTc                                                                                                                                                                                       >  2:36957/1‑251 (MQ=255)
           ggCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAAC‑TGTCCatat                                                                                                                                                                                  >  2:426549/1‑251 (MQ=255)
                    gcagTTATGATGTGGCGGAGCGATGTACAGATACTCGTTTAATTCAGGTATACCGCGGGAATCGAAATATTCGCTATACTGTGCTTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATGGATAATAAAGTTAAAAAAACAAATAAATACAAGTCAATTGGGGCCAAA‑T‑TGTCCATATCataaataa                                                                                                                                                                         <  1:135975/248‑1 (MQ=255)
                                                    aCTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAAC‑TGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAtt                                                                                                                                         >  2:307908/1‑251 (MQ=255)
                                                                 aGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAA‑AATGTCCATATCATAAATAAGTTACGTATTTTTTCTCAATCATAAAAATATTAAAAAAcaacaa                                                                                                                             >  2:343119/1‑250 (MQ=255)
                                                                      tATCTCGGGCATCGAAATATTCGCCAAACTGTGCGTTATACCGACAGTAGTGCTGATCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAAC‑TGTCCATATCATAAATAAGTTAAGAATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAGAGc                                                                                                                        >  2:326007/1‑250 (MQ=255)
                                                                         ctcGGGCATCGACATATTCGCCCTACTGTGCGTTATACCGCAAGTAATGCTGCTCGTTTTGCCGGATTATGGGACAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATACTAAAGATAAAAAAACCAATAAATACAAGACAAT‑GGGGCCAAAAC‑TGCCCATATCATAAATAAGTTCCGTTTTTTTTCTCAATCTTAAACATATTAAAAAACGACAAAAAGCATc                                                                                                                     >  2:13146/1‑250 (MQ=255)
                                                                                                         ttATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAAC‑TGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGt                                                                                     >  1:39563/1‑250 (MQ=255)
                                                                                                           ataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAAC‑TGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTaaa                                                                                  <  1:307908/251‑1 (MQ=255)
                                                                                                                            gctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAAC‑TGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCt                                                                 <  1:104122/251‑1 (MQ=255)
                                                                                                                                                                      gAAAAATTAAAAAGAGAAGAGGGTTGATGTAACTTATTGATAAGAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCAAAAAC‑TGTCCAGATCATAAATAAGTTACGTATTTTGTCTCAAGCATAAAAATATTAAAAAACGATAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCGCCGCTATATGTAAAGCTAACGTTTCTGTTGCTTGACGCATCt                       <  2:260884/251‑1 (MQ=255)
                                                                                                                                                                             tAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAAC‑TGTCCATATCATAAATAAGTTACGTCTTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGt                                                      >  1:320627/1‑213 (MQ=255)
                                                                                                                                                                             tAAAAAGAGAAGAGGATTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAAC‑TGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATGTCACATCACCTCCGCtataggt                                                      <  2:320627/213‑4 (MQ=255)
                                                                                                                                                                                           gTTTGATTTAACGTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAAAC‑TTTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTTTAAATTAATCTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATGCTTCTCTCCaaaa  <  1:469575/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                     
TCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAA‑C‑TGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAA  >  NC_000913/4542441‑4542878

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