Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,547,402 Δ1 bp 100% intergenic (+87/‑395) fimE → / → fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,547,4020T.93.3% 51.4 / ‑1.9 15intergenic (+87/‑395)fimE/fimAtyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (14/0);  total (15/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATG  >  W3110S.gb/4547376‑4547460
                          |                                                          
tGTCCATATCATAAATAAGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCt                >  1:1874395/1‑70 (MQ=255)
tGTCCATATCATAAATAAGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCt                >  1:312215/1‑70 (MQ=255)
tGTCCATATCATAAATAAGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCt                >  1:640882/1‑70 (MQ=255)
tGTCCATATCATAAATAAGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCt                >  1:960289/1‑70 (MQ=255)
        tCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGttt        >  1:1443304/1‑71 (MQ=255)
        tCATAAATAAGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGttt        >  1:1170819/1‑70 (MQ=255)
        tCATAAATAAGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGttt        >  1:1723601/1‑70 (MQ=255)
             aataaGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGa                            >  1:876559/1‑45 (MQ=39)
             aataaGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATg  >  1:1150744/1‑71 (MQ=255)
             aataaGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATg  >  1:1549934/1‑71 (MQ=255)
             aataaGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATg  >  1:1582155/1‑71 (MQ=255)
             aataaGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATg  >  1:1686039/1‑71 (MQ=255)
             aataaGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATg  >  1:1789406/1‑71 (MQ=255)
             aataaGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATg  >  1:349282/1‑71 (MQ=255)
             aataaGTTACGTATTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATg  >  1:773423/1‑71 (MQ=255)
                          |                                                          
TGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATG  >  W3110S.gb/4547376‑4547460

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