Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,428,905 G→A 100% S226S (AGC→AGT wzzE ← Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,428,9050GA92.6% 154.8 / 3.0 54S226S (AGC→AGTwzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (50/0);  total (52/2)
Fisher's exact test for biased strand distribution p-value = 4.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTACATCTGTCGCACTGCGCGAAATATTATGCTGCTCAGCAATTTTCAGCGCCTGCTCAATGCTGTTCATCCGGCGGTCGTAGATGGCTTTCGC  >  minE/2428844‑2428937
                                                             |                                
gTACATCTGTCGCACTGCGCGAAATATTATGCTGCTCAGCAATTTTCAGCGCCTGCTCAATGCTGtt                             <  1:1524838/67‑1 (MQ=255)
gTACATCTGTCGCACTGCGCGAAATATTATGCTGCTCAGCAATTTTCAGCGCCTGCTCAATGCTGtt                             <  1:857742/67‑1 (MQ=255)
     tCTGTCGCACTGCGCGAAATATTATGCTGCTCAGCAATTTTCAGCGCCTGCTCAATGc                                 >  1:1707355/1‑58 (MQ=255)
     tCTGTCGCACTGCGCGAAATATTATGCTGCTCAGCAATTTTCAGCGCCTGCTCAATGc                                 >  1:458237/1‑58 (MQ=255)
                                             ttaGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:836798/3‑49 (MQ=255)
                                             tCAGCGCCTGTTCGATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:830144/1‑49 (MQ=21)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCtt      >  1:963267/1‑45 (MQ=25)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:573385/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:593731/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:586803/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:575786/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:779229/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:460233/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:439271/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:438602/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:368013/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:350977/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:596457/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:637927/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:663229/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:293582/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:898080/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:903605/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:9151/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:921796/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:926141/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:94656/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:997736/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:175129/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1221560/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1228734/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1282614/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1295227/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:129761/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1371505/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1436996/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1513468/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1531006/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1577352/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1610128/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:174076/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:348786/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1780707/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1918999/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1938195/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1959761/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:2000290/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:2092351/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:2129897/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:2213465/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:2311159/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:241721/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:281694/1‑49 (MQ=37)
                                             tCAGCGCCTGTTCAATACTGTTCATCCGGCGATCGTAGATAGCTTTCGc  >  1:1168454/1‑49 (MQ=37)
                                                             |                                
GTACATCTGTCGCACTGCGCGAAATATTATGCTGCTCAGCAATTTTCAGCGCCTGCTCAATGCTGTTCATCCGGCGGTCGTAGATGGCTTTCGC  >  minE/2428844‑2428937

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