Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,491,067 G→A 37.2% V56V (GTG→GTA yjgP → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,491,0670GA37.2% 36.5 / 41.4 43V56V (GTG→GTAyjgPconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (16/11);  new base A (12/4);  total (28/15)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTCGGCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGACGCTGGGCAAACT  >  W3110S.gb/4491000‑4491135
                                                                   |                                                                    
tCCTCGGCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCg                                                                   >  1:2905316/1‑71 (MQ=255)
tCCTCGGCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCg                                                                   >  1:2897927/1‑71 (MQ=255)
tCCTCGGCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCg                                                                   >  1:3039669/1‑71 (MQ=255)
tCCTCGGCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCg                                                                   >  1:2370796/1‑71 (MQ=255)
tCCTCGGCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCg                                                                   >  1:2229346/1‑71 (MQ=255)
tCCTCGGCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCg                                                                   >  1:1389404/1‑71 (MQ=255)
tCCTCGGCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCg                                                                   >  1:1224922/1‑71 (MQ=255)
     ggCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGaaa                                                               >  1:1941065/1‑70 (MQ=255)
     ggCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAAt                                                              >  1:678203/1‑71 (MQ=255)
     ggCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAAt                                                              >  1:1041597/1‑71 (MQ=255)
     ggCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAAt                                                              >  1:2687168/1‑71 (MQ=255)
     ggCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAAt                                                              >  1:2631459/1‑71 (MQ=255)
     ggCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAAt                                                              >  1:2577734/1‑71 (MQ=255)
     ggCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAAt                                                              >  1:2270899/1‑71 (MQ=255)
     ggCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAAt                                                              >  1:19747/1‑71 (MQ=255)
     ggCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAAt                                                              >  1:1449199/1‑71 (MQ=255)
     ggCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAAt                                                              >  1:1446824/1‑71 (MQ=255)
                      cGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGc                                             <  1:1744623/71‑1 (MQ=255)
                      cGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGc                                             <  1:2988268/71‑1 (MQ=255)
                      cGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGc                                             <  1:680991/71‑1 (MQ=255)
                      cGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGc                                             <  1:318022/71‑1 (MQ=255)
                             ccGGCGAATCTGGTTCTCTCCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAg                                      >  1:1231147/1‑71 (MQ=255)
                             ccGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAg                                      >  1:1188722/1‑71 (MQ=255)
                                              ctcCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTg                                  >  1:2430343/1‑58 (MQ=255)
                                              ctcCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGtt                                >  1:2901802/1‑60 (MQ=255)
                                              ctcCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGtt                                >  1:873206/1‑60 (MQ=255)
                                                  cTTCTCGGGTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCtgat                 <  1:1002060/71‑1 (MQ=255)
                                                  cTTCTCGGGTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCtgat                 <  1:158594/71‑1 (MQ=255)
                                                  cTTCTCGGGTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCtgat                 <  1:1162385/71‑1 (MQ=255)
                                                  cTTCTCGGGTTGGGCGGGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCtgat                 <  1:3326337/71‑1 (MQ=255)
                                                        gggTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAc              <  1:3219572/68‑1 (MQ=255)
                                                        gggTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAc              <  1:953813/68‑1 (MQ=255)
                                                        gggTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAc              <  1:2117694/68‑1 (MQ=255)
                                                        gggTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAc              <  1:1823028/68‑1 (MQ=255)
                                                        gggTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAc              <  1:1044828/68‑1 (MQ=255)
                                                        gggTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAc              <  1:1044818/68‑1 (MQ=255)
                                                         ggTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAc              <  1:485390/67‑1 (MQ=255)
                                                                 gTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGACGCTGGGCAAACt  >  1:2985374/1‑71 (MQ=255)
                                                                 gTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGACGCTGGGCAAACt  >  1:2846853/1‑71 (MQ=255)
                                                                 gTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGACGCTGGGCAAACt  >  1:269608/1‑71 (MQ=255)
                                                                 gTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGACGCTGGGCAAACt  >  1:2686657/1‑71 (MQ=255)
                                                                 gTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGACGCTGGGCAAACt  >  1:1284124/1‑71 (MQ=255)
                                                                 gTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGACGCTGGGCAAACt  >  1:1006927/1‑71 (MQ=255)
                                                                   |                                                                    
TCCTCGGCGCAGCGGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGACGCTGGGCAAACT  >  W3110S.gb/4491000‑4491135

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