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.95.2% 83.8 / ‑2.0 21intergenic (+50/‑432)fimE/fimAregulator for fimA/type 1 fimbriae major subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (7/13);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAAT  >  NC_000913/4542457‑4542909
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
gCTTGCGGTTATGAATTGGCGGAGCGTGGGGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATTTTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCataaa                                                                                                                                                                                                           <  1:11451/251‑1 (MQ=255)
            gAATGGGCGGGGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATAGTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTAt                                                                                                                                                                                               <  1:87831/251‑1 (MQ=255)
                    ggAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                                                                                                                                                       <  1:214234/251‑1 (MQ=255)
                                        gTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCCCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATAATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaa                                                                                                                                                                    >  2:141465/1‑250 (MQ=255)
                                                                    aaaTATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGata                                                                                                                                       <  1:235084/251‑1 (MQ=255)
                                                                              aaTACTGTGCGTTATACCGCCAGTAATGAGTCTCGTTTTGCCGGATTATGGGAAATAAATAATCTCATAAACGAAAAAGTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAATACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAATCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAattat                                                                                                                             <  2:46975/250‑1 (MQ=255)
                                                                                     ggCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTTTTTGATATGTAAATTATTTCTAtt                                                                                                                      <  1:155367/250‑1 (MQ=255)
                                                                                       cGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGt                                                                                                                    <  1:169066/251‑1 (MQ=255)
                                                                                             aCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaatt                                                                                                               >  1:139063/1‑250 (MQ=255)
                                                                                                  cAGNAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGATTTTTAATTTAACTTATTGATAAGAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGTCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAATCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTACTGTAAAGTAATTTc                                                                                                         <  2:418737/251‑1 (MQ=255)
                                                                                                      aaTGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTcaca                                                                                                      <  1:393290/250‑1 (MQ=255)
                                                                                                            gctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCt                                                                                                <  1:4171/250‑1 (MQ=255)
                                                                                                                ggtttGCCGGATAATGGGAAAGTCCTAATATCATAAACGCAAAATTAACAACAGAAGAGATTTGATTTAACTTATTAATAATAAAGTTAAAAAAACAAAGAAATACAAGACA‑GTGGGGCCAAACTGTCCATATCATAAATTAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACATCCGCt                                                                                           >  2:93631/3‑251 (MQ=255)
                                                                                                                     gCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTatat                                                                                       <  1:340534/250‑1 (MQ=255)
                                                                                                                                                                        gagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACACATAAATAAAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATACGTTACGTATTTTTTCTCAAGCATAAAAATATTAACAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCa                                    >  2:208875/1‑250 (MQ=255)
                                                                                                                                                                        gagGTTGGATTTAACTTATGGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATAATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTAACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCa                                    >  2:214171/1‑250 (MQ=255)
                                                                                                                                                                            tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAac                               >  1:253279/1‑251 (MQ=255)
                                                                                                                                                                                          ttGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAGCTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAat                  >  1:433966/1‑250 (MQ=255)
                                                                                                                                                                                             ataataAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATaa              <  1:232231/251‑1 (MQ=255)
                                                                                                                                                                                               aataaAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAAc            <  2:278829/251‑1 (MQ=255)
                                                                                                                                                                                                          aaaaaCAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTataaat  >  2:124004/1‑250 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
GCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAAT  >  NC_000913/4542457‑4542909

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: ‑