Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,549,768 C→T 13.0% intergenic (+57/‑10) fimC → / → fimD chaperone, periplasmic/outer membrane usher protein, type 1 fimbrial synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,549,7680CT13.0% 124.6 / 9.5 46intergenic (+57/‑10)fimC/fimDchaperone, periplasmic/outer membrane usher protein, type 1 fimbrial synthesis
Reads supporting (aligned to +/- strand):  ref base C (15/25);  new base T (2/4);  total (17/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTA  >  W3110S.gb/4549703‑4549835
                                                                 |                                                                   
aTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCggagg                                                                <  1:2427834/71‑1 (MQ=255)
   gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGaa                                                              >  1:2858956/1‑70 (MQ=255)
   gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaat                                                             <  1:938272/71‑1 (MQ=255)
   gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaat                                                             <  1:3197661/71‑1 (MQ=255)
        aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                                        <  1:1677941/71‑1 (MQ=255)
           cAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCat                                                     <  1:345578/71‑1 (MQ=255)
           cAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCat                                                     <  1:2183052/71‑1 (MQ=255)
              ggggAATTTTTCGCCTGAATAAAAAGAATTGACTGTCGGGTGATTTTAAGCCGGAGGAATAATGTCATATc                                                  <  1:1411765/71‑1 (MQ=255)
              ggggAATTTTTCCCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGaat                                                             <  1:3395254/60‑1 (MQ=255)
               gggAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCt                                                 >  1:33858/1‑71 (MQ=255)
               gggAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTAAAGCCGGAGGAATAATGTCATATCt                                                 >  1:312495/1‑71 (MQ=255)
                ggAATTTTTCGCCTGAATAAAAAGAATTGACTGTCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTg                                                <  1:3269839/71‑1 (MQ=255)
                 gAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGa                                               <  1:1477343/71‑1 (MQ=255)
                 gAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGa                                               <  1:2482821/71‑1 (MQ=255)
                  aaTTTTTCGCCTGAATAAAAAGAATTGACTGTCGGGTGATTTTAAGCCGGAGGaata                                                            <  1:3243236/57‑1 (MQ=255)
                  aaTTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCata                                                    >  1:3924618/1‑65 (MQ=255)
                  aaTTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCata                                                    >  1:3195902/1‑65 (MQ=255)
                   aTTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGt                                                        <  1:3556455/60‑1 (MQ=255)
                     ttttCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAAt                                             >  1:855033/1‑69 (MQ=255)
                      tttCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGaataa                                                           >  1:2938381/1‑54 (MQ=255)
                      tttCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAttt                                           <  1:1448135/70‑1 (MQ=255)
                      tttCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAttt                                           >  1:3437946/1‑70 (MQ=255)
                      tttCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAttt                                           >  1:1639499/1‑70 (MQ=255)
                            cTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGa                                               <  1:1863206/60‑1 (MQ=255)
                             tGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCggagg                                                                >  1:1161664/1‑42 (MQ=255)
                               aaTAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACttt                                   <  1:3463185/69‑1 (MQ=255)
                               aaTAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACttt                                   <  1:1021252/69‑1 (MQ=255)
                               aaTAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACttt                                   <  1:3763624/69‑1 (MQ=255)
                               aaTAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                                        <  1:374880/48‑1 (MQ=255)
                                 tAAAAAGAATTGACTGCCAGGTGATTTTAAGCCgg                                                                   <  1:816948/35‑1 (MQ=255)
                                     aaGAATTGACTGTCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGa                           <  1:2944041/71‑1 (MQ=255)
                                     aaGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                                         >  1:2799701/1‑57 (MQ=255)
                                     aaGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                                         >  1:1285920/1‑57 (MQ=255)
                                              cTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACaa                  <  1:3974341/71‑1 (MQ=255)
                                                  cGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCt              <  1:3301734/71‑1 (MQ=255)
                                                    ggTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                                         <  1:1691230/42‑1 (MQ=255)
                                                    ggTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                                         <  1:1652082/42‑1 (MQ=255)
                                                    ggTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                                         <  1:2733111/42‑1 (MQ=255)
                                                       gATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGaaa                         >  1:807715/1‑55 (MQ=255)
                                                       gATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCa          >  1:774105/1‑70 (MQ=255)
                                                       gATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCa          >  1:2717737/1‑70 (MQ=255)
                                                        aTTTTAAGCCGTAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCata        <  1:1866625/71‑1 (MQ=255)
                                                        aTTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCata        >  1:2293077/1‑71 (MQ=255)
                                                        aTTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCata        >  1:1917931/1‑71 (MQ=255)
                                                         ttttAAGCCGGAGGAATAATGTCATATCTGAAttt                                           <  1:2105464/35‑1 (MQ=255)
                                                              aGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTa  <  1:2988524/71‑1 (MQ=255)
                                                                 |                                                                   
ATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTA  >  W3110S.gb/4549703‑4549835

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