Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,547,549 T→C 17.4% intergenic (+234/‑248) 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,5490TC17.4% 121.1 / 19.5 53intergenic (+234/‑248)fimE/fimAtyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin
Reads supporting (aligned to +/- strand):  ref base T (30/13);  new base C (4/5);  total (35/18)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

TTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCC  >  W3110S.gb/4547486‑4547616
                                                               |                                                                   
ttCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTtctc                                                                 >  1:3591291/1‑68 (MQ=255)
ttCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTtctc                                                                 >  1:743764/1‑68 (MQ=255)
 tCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTtctct                                                                >  1:2516112/1‑68 (MQ=255)
 tCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTtctct                                                                >  1:605449/1‑68 (MQ=255)
 tCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCAtctct                                                                >  1:1898831/1‑68 (MQ=255)
      tCACCTCCGCTATATGTAAAGCTAACGTTTGTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaaa                                                          <  1:1399229/69‑1 (MQ=255)
       cACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaaa                                                          >  1:3211764/1‑68 (MQ=255)
       cACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaaa                                                          >  1:2897202/1‑68 (MQ=255)
        aCCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaa                                                           >  1:610554/1‑66 (MQ=255)
            ccGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAaccacc                                                    <  1:3502457/69‑1 (MQ=255)
                    tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCa                                             <  1:1710378/68‑1 (MQ=255)
                    tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCa                                             <  1:1594357/68‑1 (MQ=255)
                    tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCa                                             <  1:36955/68‑1 (MQ=255)
                    tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCAAAAac                                                        <  1:1605248/57‑1 (MQ=255)
                    tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCAAAAac                                                        <  1:1891034/57‑1 (MQ=255)
                    tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCAAAAac                                                        <  1:3423674/57‑1 (MQ=255)
                    tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCAAAAac                                                        <  1:2892447/57‑1 (MQ=255)
                      tAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAat                                           >  1:177905/1‑68 (MQ=255)
                      tAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAat                                           >  1:891321/1‑68 (MQ=255)
                              aCGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAAtata                                        <  1:1213624/63‑1 (MQ=255)
                              aCGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAAtata                                        <  1:2847316/63‑1 (MQ=255)
                              aCGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAAtata                                        <  1:238198/63‑1 (MQ=255)
                              aCGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACAt                                   >  1:2010435/1‑68 (MQ=255)
                                gTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCaaa                                                           >  1:2388821/1‑42 (MQ=255)
                                gTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCaaa                                                           >  1:3098749/1‑42 (MQ=255)
                                gTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCaaa                                                           >  1:1800509/1‑42 (MQ=255)
                                 tttCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCta                                >  1:1502928/1‑68 (MQ=255)
                                   tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATaaa                                      <  1:3699622/60‑1 (MQ=255)
                                   tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATaaa                                      <  1:2715955/60‑1 (MQ=255)
                                   tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATaaa                                      <  1:3544779/60‑1 (MQ=255)
                                           tCGACGCATCTTCCTCATTCCTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGAt                     <  1:99817/69‑1 (MQ=255)
                                               cGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATaa                   >  1:1407492/1‑67 (MQ=255)
                                                gCATCTTCCTCATTCTTCTCTCCAAAAACCCCCTCATGCAATATAAACATCTATAAATAAAGATAAc                  >  1:39502/1‑67 (MQ=255)
                                                gCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATa                    >  1:2127181/1‑65 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAAtata                                        >  1:1396980/1‑40 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAAtata                                        >  1:619060/1‑40 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACa                 >  1:2793187/1‑63 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACa                 >  1:707419/1‑63 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACa                 >  1:3339319/1‑63 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACa                 >  1:2517429/1‑63 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAg             <  1:3576322/67‑1 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAg             <  1:58046/67‑1 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGaa           >  1:1258326/1‑69 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGaa           >  1:915204/1‑69 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGaa           >  1:1473066/1‑69 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGaa           >  1:760962/1‑69 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGaa           >  1:1871960/1‑69 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGaa           >  1:605177/1‑69 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGaa           >  1:2614193/1‑69 (MQ=255)
                                                     ttCCTCATTCTTCTCTCCAAAAACCACCTCATGAAAtata                                        >  1:148633/1‑40 (MQ=255)
                                                           aTTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAatat        >  1:2726037/1‑66 (MQ=255)
                                                           aTTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTaa     >  1:2680657/1‑69 (MQ=255)
                                                           aTTCCTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTaa     >  1:2194417/1‑69 (MQ=255)
                                                             tcttctCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTaa     >  1:3172055/1‑67 (MQ=255)
                                                              cttctCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGcc  >  1:563643/1‑69 (MQ=255)
                                                              cttctCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGcc  >  1:415676/1‑69 (MQ=255)
                                                              cttctCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGcc  >  1:2460678/1‑69 (MQ=255)
                                                              cttctCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGcc  >  1:2284175/1‑69 (MQ=255)
                                                              cttctCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGcc  >  1:1543108/1‑69 (MQ=255)
                                                              cttctCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGcc  >  1:710033/1‑69 (MQ=255)
                                                              cttctCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGcc  >  1:2942712/1‑69 (MQ=255)
                                                               |                                                                   
TTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCC  >  W3110S.gb/4547486‑4547616

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