Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,547,549 T→C 13.8% 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,5490TC13.8% 166.0 / 17.4 65intergenic (+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 (22/34);  new base C (3/6);  total (25/40)
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

TAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTA  >  W3110S.gb/4547482‑4547612
                                                                   |                                                               
taattTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATtcttct                                                              >  1:3610644/1‑71 (MQ=255)
taattTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATtcttct                                                              >  1:1952214/1‑71 (MQ=255)
taattTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATtcttct                                                              >  1:2137893/1‑71 (MQ=255)
  attTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTtctc                                                             <  1:3305687/70‑1 (MQ=255)
  attTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTtctc                                                             <  1:2719723/70‑1 (MQ=255)
       acaTCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaa                                                       <  1:973367/71‑1 (MQ=255)
       acaTCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaa                                                       <  1:174050/71‑1 (MQ=255)
       acaTCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaa                                                       <  1:1892105/71‑1 (MQ=255)
       acaTCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaa                                                       <  1:1583510/71‑1 (MQ=255)
       acaTCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaa                                                       <  1:1487704/71‑1 (MQ=255)
       acaTCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaa                                                       <  1:2759163/71‑1 (MQ=255)
       acaTCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaa                                                       <  1:1319240/71‑1 (MQ=255)
            aCCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCTACGCATCTTCCTCATTCTTCTCTCCAAAAacca                                                  >  1:1608144/1‑71 (MQ=255)
            aCCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAacca                                                  >  1:794515/1‑71 (MQ=255)
             ccTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCaa                                                        <  1:2108612/64‑1 (MQ=255)
             ccTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCaa                                                        <  1:247824/64‑1 (MQ=255)
             ccTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCaa                                                        <  1:3213609/64‑1 (MQ=255)
             ccTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCaa                                                        <  1:1250305/64‑1 (MQ=255)
             ccTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCCTCTCTCCaa                                                        <  1:1883710/64‑1 (MQ=255)
                     atatGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGc                                          >  1:2524112/1‑70 (MQ=255)
                        tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTtctct                                                            >  1:3657009/1‑49 (MQ=255)
                        tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTtctct                                                            >  1:1695558/1‑49 (MQ=255)
                        tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAata                                      >  1:2874090/1‑71 (MQ=255)
                        tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAata                                      >  1:2173334/1‑71 (MQ=255)
                        tGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTACAAAAACCACCTCATGCAata                                      >  1:3685473/1‑71 (MQ=255)
                          tAAAGCTAACGTTTCTGTGGCTCGACGCAGCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAAtata                                    <  1:594571/71‑1 (MQ=255)
                            aaGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATaaa                                  <  1:1401261/71‑1 (MQ=255)
                                    gTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCtata                          <  1:1870635/71‑1 (MQ=255)
                                     tttCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCaa                                        <  1:1985796/56‑1 (MQ=255)
                                     tttCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCaa                                        <  1:412832/56‑1 (MQ=255)
                                     tttCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCaa                                        <  1:846894/56‑1 (MQ=255)
                                     tttCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCaa                                        <  1:2742275/56‑1 (MQ=255)
                                       tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCta                            >  1:2893623/1‑66 (MQ=255)
                                       tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTataaa                        >  1:2478430/1‑70 (MQ=255)
                                       tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTataaa                        >  1:46294/1‑70 (MQ=255)
                                       tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTataaa                        >  1:2813643/1‑70 (MQ=255)
                                       tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTataaa                        >  1:3307913/1‑70 (MQ=255)
                                       tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTataaa                        >  1:1713570/1‑70 (MQ=255)
                                       tCTGTGGCTAGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTataaa                        >  1:1827306/1‑70 (MQ=255)
                                         tgtgGCTCGACGCATCTTCCTCATTCCTCTCTCCaa                                                        <  1:1989142/36‑1 (MQ=255)
                                             gCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGAt                 <  1:2603701/71‑1 (MQ=255)
                                             gCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGAt                 <  1:675015/71‑1 (MQ=255)
                                             gCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGAt                 <  1:1546142/71‑1 (MQ=255)
                                             gCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGAt                 <  1:2039538/71‑1 (MQ=255)
                                             gCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGAt                 <  1:2439166/71‑1 (MQ=255)
                                             gCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGAt                 <  1:2731547/71‑1 (MQ=255)
                                               tCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATaa               >  1:1051529/1‑71 (MQ=255)
                                                  aCGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACaa            <  1:532418/71‑1 (MQ=255)
                                                  aCGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACaa            <  1:2122957/71‑1 (MQ=255)
                                                       tCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGAt                 <  1:1916746/61‑1 (MQ=255)
                                                       tCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGAt                 <  1:640505/61‑1 (MQ=255)
                                                       tCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGAt                 <  1:609385/61‑1 (MQ=255)
                                                       tCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAg         <  1:1793774/69‑1 (MQ=255)
                                                       tCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAg         <  1:1353838/69‑1 (MQ=255)
                                                       tCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAg         <  1:2536285/69‑1 (MQ=255)
                                                       tCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAg         <  1:3044425/69‑1 (MQ=255)
                                                       tCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAg         <  1:1828601/69‑1 (MQ=255)
                                                         ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCaa                                        >  1:426058/1‑36 (MQ=255)
                                                         ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATaa               >  1:1366218/1‑61 (MQ=255)
                                                         ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATaa               >  1:700170/1‑61 (MQ=255)
                                                         ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATaa               >  1:1046737/1‑61 (MQ=255)
                                                         ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAata     >  1:913178/1‑71 (MQ=255)
                                                         ttCCTCATTCCTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGaa       >  1:157602/1‑69 (MQ=255)
                                                         ttCCTCATTCCTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGaa       >  1:296823/1‑69 (MQ=255)
                                                         ttCCTCATTCCTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGaa       >  1:1139507/1‑69 (MQ=255)
                                                           ccTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGaa       <  1:2513692/67‑1 (MQ=255)
                                                             tCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAg         <  1:618078/63‑1 (MQ=255)
                                                              cATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTa  >  1:1702774/1‑69 (MQ=255)
                                                                   |                                                               
TAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTA  >  W3110S.gb/4547482‑4547612

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: