Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,549,754 G→A 80.8% intergenic (+43/‑24) 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,7540GA80.8% 53.1 / 8.7 26intergenic (+43/‑24)fimC/fimDchaperone, periplasmic/outer membrane usher protein, type 1 fimbrial synthesis
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (11/10);  total (13/13)
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.29e-01

AAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACA  >  W3110S.gb/4549688‑4549814
                                                                  |                                                            
aaaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGt                                                           >  1:1939447/1‑70 (MQ=255)
        gggCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTCAATAAAAAGAATTGACTGCCAGGTGATTTTaa                                                   <  1:2360584/70‑1 (MQ=255)
                  taaTAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGTCGGGTGATTTTAAGCCGGAGGaat                                        <  1:1878328/70‑1 (MQ=255)
                  gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaat                                        <  1:1845110/71‑1 (MQ=255)
                  gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaat                                        <  1:2737510/71‑1 (MQ=255)
                  gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaat                                        <  1:2698542/71‑1 (MQ=255)
                       aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                   >  1:2352574/1‑71 (MQ=255)
                                 aaTTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTa                                                    <  1:1102103/44‑1 (MQ=255)
                                   tttttCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAg                                                  >  1:392182/1‑44 (MQ=255)
                                     tttCGCCTGAATAAAAAGAATTGACTGTCGGGTGATTTTAAGCCGGAGGaa                                         >  1:2341513/1‑51 (MQ=255)
                                     tttCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAttt                      >  1:986620/1‑70 (MQ=255)
                                     tttCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAttt                      <  1:411987/70‑1 (MQ=255)
                                     tttCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAttt                      >  1:1128879/1‑70 (MQ=255)
                                         gCCTGAATAAAAAGAATTGACTGTCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGAc                 <  1:2565297/71‑1 (MQ=255)
                                              aaTAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAc            >  1:2303952/1‑71 (MQ=255)
                                              aaTAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAc            >  1:2300609/1‑71 (MQ=255)
                                                tAAAAAGAATTGACTGTCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                    <  1:470439/61‑1 (MQ=255)
                                                tAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAc            >  1:865722/1‑69 (MQ=255)
                                                  aaaaGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCatat                              >  1:2355627/1‑49 (MQ=255)
                                                      gAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAg         >  1:2002291/1‑66 (MQ=255)
                                                        aTTGACTGTCGGGTGATTTTAAGCCGGAGGAATAATGTCatat                              >  1:1644060/1‑43 (MQ=255)
                                                        aTTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGATTTTAAGACTTTACCAGCGAAac   >  1:1976249/1‑70 (MQ=255)
                                                        aTTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTa             <  1:82670/60‑1 (MQ=255)
                                                        aTTGACTGCCAGGTGATTTTAAGCCGGAGCAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAaca  <  1:2037733/71‑1 (MQ=255)
                                                         ttGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAttt                      <  1:278195/50‑1 (MQ=255)
                                                         ttGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                    <  1:1357755/52‑1 (MQ=255)
                                                                  |                                                            
AAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACA  >  W3110S.gb/4549688‑4549814

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: