Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,549,768 C→T 75.9% 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,7680CT75.9% 49.4 / 14.2 29intergenic (+57/‑10)fimC/fimDchaperone, periplasmic/outer membrane usher protein, type 1 fimbrial synthesis
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base T (9/13);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGC  >  W3110S.gb/4549706‑4549822
                                                              |                                                      
gAATAACGCAGGGGTAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGaat                                                <  1:1876607/71‑1 (MQ=255)
gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTggag                                                    >  1:1595558/1‑67 (MQ=255)
gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTggag                                                    >  1:2360394/1‑67 (MQ=255)
gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGaat                                                <  1:749866/71‑1 (MQ=255)
gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGaat                                                <  1:741116/71‑1 (MQ=255)
gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGaat                                                <  1:3720312/71‑1 (MQ=255)
gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGaat                                                <  1:2097263/71‑1 (MQ=255)
gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGaa                                                 >  1:2741962/1‑70 (MQ=255)
gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGCAGGaat                                                <  1:479106/71‑1 (MQ=255)
gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaat                                                <  1:381594/71‑1 (MQ=255)
     aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGt                                           <  1:574596/71‑1 (MQ=255)
     aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGt                                           >  1:1369279/1‑71 (MQ=255)
                       taatGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGAc                         <  1:1033287/68‑1 (MQ=255)
                            aaTAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTaa                            <  1:2290287/63‑1 (MQ=255)
                            aaTAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAc                    >  1:2646627/1‑71 (MQ=255)
                              tAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACTTTAc                    >  1:170180/1‑69 (MQ=255)
                                 aaaGAATTGACTGCCGGGTGATTTTAAGCTGGAGGaat                                                <  1:3126072/38‑1 (MQ=255)
                                  aaGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGa              <  1:2914257/71‑1 (MQ=255)
                                         gACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACTTTAc                    >  1:1715062/1‑58 (MQ=255)
                                            tGCTGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAcc                   <  1:243097/56‑1 (MQ=255)
                                            tGCTGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAcc                   <  1:1476941/56‑1 (MQ=255)
                                            tGCCGGGTGATTTTAAGCTGGAGGAATAATGTCat                                        >  1:1046327/1‑35 (MQ=255)
                                            tGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAAt    >  1:3524948/1‑71 (MQ=255)
                                              ccGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGc  <  1:376774/71‑1 (MQ=255)
                                              ccGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGc  <  1:1334378/71‑1 (MQ=255)
                                                      ttttAAGCTGGAGGAATAATGTCATATCTGAATTTaa                            <  1:2789767/37‑1 (MQ=255)
                                                      ttttAAGCTGGAGGAATAATGTCATATCTGAATTTaa                            <  1:1334513/37‑1 (MQ=255)
                                                        ttAAGCTGGAGGAATAATGTCATATCTGAATTTaa                            >  1:3517024/1‑35 (MQ=255)
                                                        ttAAGCTGGAGGAATAATGTCATATCTGAATTTaa                            >  1:3521405/1‑35 (MQ=255)
                                                              |                                                      
GAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGC  >  W3110S.gb/4549706‑4549822

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