Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,543,047 T→G intergenic (+414/‑68) 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,543,0470TG94.7% 46.1 / NA 19intergenic (+414/‑68)fimE/fimAregulator for fimA/type 1 fimbriae major subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (7/11);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

CACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTGTTGATC  >  NC_000913/4542806‑4543271
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                
cacaTCACCTCCGCTATATGTAAAGCTAAAGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTATCTCCAAAAACCACCTCATGCAATATAAACACCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGATTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCCAATGTTTCATCTTTTGTGGGAACACTGTGCAGTGTTGGCAGTCAACCTCGTTGacaaaaca                                                                                                                                                                                                                          >  1:13293/1‑250 (MQ=255)
          ccGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAg                                                                                                                                                                                                                <  2:194825/250‑1 (MQ=255)
               atatGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACg                                                                                                                                                                                                            <  1:68983/249‑1 (MQ=255)
                            aCGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAt                                                                                                                                                                                             <  2:230121/251‑1 (MQ=255)
                              gTTTCTGTGGCTCGACGCCTCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGATTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCCAACTCGTTGACAAAACAAAGTGAACAGAACGACTGCCCATGTCGCttt                                                                                                                                                                                           >  2:112224/1‑251 (MQ=255)
                                    gtgGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGGTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAAt                                                                                                                                                                                     <  2:144044/251‑1 (MQ=255)
                                           gACGCAACTTCCTCATTCTTCTCTCCAAACACCACCTCCTGCACTATAACCTTATATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGACCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTCCAGAACGACTGCCCATGTCGATTTAGAAATAGttttt                                                                                                                                                                              >  2:176875/1‑251 (MQ=255)
                                            aCGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAGGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACGCGTTGACAAAACGAAGTGTACAGAACGATTGCGCATGTCGATTTAGATATAGtttttt                                                                                                                                                                             <  2:128854/251‑1 (MQ=255)
                                                      cTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaaggaa                                                                                                                                                                     >  2:263436/1‑249 (MQ=255)
                                                                     aaaaaCCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTGCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaa                                                                                                                                                    <  1:30082/251‑1 (MQ=255)
                                                                        aaCCACCTCATGCAATATAAACAACTATAAATAAAGATAACAATAGAATATTAAGCCAACAAAGAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTAGTAGTCAAAATGGCCCCAAATGTTGCATCTTTTGGGTGAAAACTGTGCATTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAAc                                                                                                                                                 <  1:117682/251‑1 (MQ=255)
                                                                                                                         atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGct                                                                                                >  1:139740/1‑251 (MQ=255)
                                                                                                                                          aCTGAAAAAGTGTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTGTGAAAGGAAAGCAGCATGAAAAGTAAAACTCTGGCAATCGTGGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCacgac                                                                               <  2:12906/251‑1 (MQ=255)
                                                                                                                                              aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGtt                                                                           <  1:20496/251‑1 (MQ=255)
                                                                                                                                                                   ttCCTCTAGTAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCGGTCCCTCAGTTCTACAGCGGCTCTGGCCGATGCCACGACGGTTAATGGTGGGACCGTTCAcc                                                        <  2:214294/249‑2 (MQ=255)
                                                                                                                                                                          aggAGTCAAAATTGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCGGTCCCTCAGTTCTACAGCGGCTCTGTCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAggc                                                <  2:156405/248‑2 (MQ=255)
                                                                                                                                                                                        ccccAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATGTAGAAATAGTTTTTTGAAAGGAAAGCAGCAGGAAAATTAAAACGCTGTCAATCGTTGTTCTGTCGGCTCTGTCCCTCATTTCTGCAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACGTTAAAGGTGAAGTTGTTAAcgcc                                 <  2:265481/251‑1 (MQ=255)
                                                                                                                                                                                                                 aaaaCTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGCCGGTTACTGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTg        >  2:176240/1‑251 (MQ=255)
                                                                                                                                                                                                                       gtgCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTGTTGATc  >  1:93687/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                
CACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTGTTGATC  >  NC_000913/4542806‑4543271

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