Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,148,033 2 bp→CA coding (2631‑2632/3318 nt) yegE → predicted diguanylate cyclase, GGDEF domain signalling protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,148,0330TC85.9% 241.7 / 24.3 128G877G (GGT→GGCyegEpredicted diguanylate cyclase, GGDEF domain signalling protein
Reads supporting (aligned to +/- strand):  ref base T (0/18);  new base C (110/0);  total (110/18)
Fisher's exact test for biased strand distribution p-value = 2.64e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01
*W3110S.gb2,148,0340GA85.9% 236.3 / 25.0 128V878I (GTC→ATC) yegEpredicted diguanylate cyclase, GGDEF domain signalling protein
Reads supporting (aligned to +/- strand):  ref base G (0/18);  new base A (110/0);  total (110/18)
Fisher's exact test for biased strand distribution p-value = 2.64e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

Only 100 of 128 total aligned reads displayed.

AATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAA  >  W3110S.gb/2148007‑2148086
                          ||                                                    
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:356867/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:577035/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:575646/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:57381/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:553955/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:539439/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:537954/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:535209/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:500561/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:448685/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:395105/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:580346/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:334616/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:33060/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:298672/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:286521/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:267939/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:264080/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:249189/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:22389/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:201324/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:172710/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:813030/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:964131/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:954498/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:9474/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:939707/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:909646/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:898263/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:894858/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:890597/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:840531/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:831207/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1009667/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:797226/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:760451/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:759376/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:740208/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:722232/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:692154/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:650559/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:632416/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:624983/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:601561/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1173514/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1336394/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1335099/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1313775/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1298676/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1297996/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1265799/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1257500/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1254889/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1230372/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1195932/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1463381/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1167314/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1159872/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1155670/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1094176/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1081671/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1069352/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1068918/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1039156/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:102165/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1004340/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1456802/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1412063/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1399643/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1414/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1450685/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1398044/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1395022/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1392538/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1455737/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1374033/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1373977/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:172703/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:136868/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1355361/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1347945/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggct                                      >  1:1338609/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggc                                       >  1:772921/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACggc                                       >  1:1093483/1‑42 (MQ=25)
aaTCAGTTGATGATGATCGCCCACGGCATCGCTTCGCCACg                                         >  1:440082/1‑41 (MQ=25)
aaTCAGTTGATGATGATCGCCCACCGCATCGCTTCGCCACggc                                       >  1:1262810/1‑42 (MQ=21)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:965243/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:928864/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:1021291/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:1074493/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:549502/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:868439/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:821604/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:305178/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:746748/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:376221/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:1203109/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:620169/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:471057/60‑1 (MQ=255)
                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa  <  1:561672/60‑1 (MQ=255)
                          ||                                                    
AATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAA  >  W3110S.gb/2148007‑2148086

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: