Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,548,379 G→A 73.1% intergenic (+34/‑31) fimA → / → fimI major type 1 subunit fimbrin/fimbrial protein involved in type 1 pilus biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,548,3790GA73.1% 39.0 / 15.6 26intergenic (+34/‑31)fimA/fimImajor type 1 subunit fimbrin/fimbrial protein involved in type 1 pilus biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (5/2);  new base A (10/9);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTT  >  W3110S.gb/4548311‑4548444
                                                                    |                                                                 
aTGTGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAgg                                                                   <  1:402740/69‑1 (MQ=255)
aTGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAgg                                                                   <  1:770568/69‑1 (MQ=255)
aTGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAgg                                                                   <  1:1786259/69‑1 (MQ=255)
aTGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAgg                                                                   <  1:1847350/69‑1 (MQ=255)
aTGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAgg                                                                   <  1:463973/69‑1 (MQ=255)
   cGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGc                                                              >  1:1030913/1‑71 (MQ=255)
   cGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGc                                                              >  1:668436/1‑71 (MQ=255)
   cGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGc                                                              >  1:57879/1‑71 (MQ=255)
   cGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGc                                                              >  1:1896620/1‑71 (MQ=255)
                              aaTAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGAGATGCCCACCCTTGTGCGATAAAAATAACGa                                    >  1:768393/1‑70 (MQ=255)
                              aaTAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGAGATGCCCACCCTTGTGCGATAAAAATAACGa                                    >  1:497390/1‑70 (MQ=255)
                              aaTAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGAGATGCCCACCCTTGTGCGATAAAAATAACGAt                                   >  1:1035057/1‑71 (MQ=255)
                              aaTAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGAGATGCCCACCCTTGTGCGATAAAAATAACGAt                                   >  1:832217/1‑71 (MQ=255)
                              aaTAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGAGATGCCCACCCTTGTGCGATAAAAATAACGAt                                   >  1:703773/1‑71 (MQ=255)
                              aaTAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGAGATGCCCACCCTTGTGCGATAAAAATAACGAt                                   >  1:1956581/1‑71 (MQ=255)
                              aaTAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGAGATGCCCACCCTTGTGCGATAAAAATAACGAt                                   >  1:1645389/1‑71 (MQ=255)
                                tAACCTACCCAGGTTCAGGGACGTCATTACGGGCAgg                                                                   <  1:1199964/37‑1 (MQ=255)
                                            gTTCAGGGACGTCATTACGGGCAGAGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGat                     <  1:301781/71‑1 (MQ=255)
                                              tCAGGGACGTCATTACGGGCAGAGATGCCCTCCCTTGTGCGATAAAAAAAACGATGAAAAGGAag                         >  1:78018/1‑65 (MQ=255)
                                              tCAGGGACGTCATTACGGGCAGAGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGatta                   >  1:1724678/1‑71 (MQ=255)
                                              tCAGGGACGTCATTACGGGCAGAGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGatt                    >  1:1965996/1‑70 (MQ=255)
                                                        cATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAg         >  1:604441/1‑71 (MQ=255)
                                                             cGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGcgtc     <  1:939949/70‑1 (MQ=255)
                                                             cGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGcgtc     <  1:1442538/70‑1 (MQ=255)
                                                               ggCAGAGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGtt  <  1:1809478/71‑1 (MQ=255)
                                                               ggCAGAGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGtt  <  1:1891778/71‑1 (MQ=255)
                                                               ggCAGAGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGtt  <  1:735576/71‑1 (MQ=255)
                                                               ggCAGAGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGtt  <  1:557004/71‑1 (MQ=255)
                                                               ggCAGAGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGtt  <  1:506468/71‑1 (MQ=255)
                                                               ggCAGAGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGtt  <  1:503972/71‑1 (MQ=255)
                                                               ggCAGAGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGtt  <  1:221615/71‑1 (MQ=255)
                                                               ggCAGAGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGtt  <  1:37842/71‑1 (MQ=255)
                                                                    |                                                                 
ATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTT  >  W3110S.gb/4548311‑4548444

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: