Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,624,185 G→A A107A (GCG→GCA)  metR → DNA‑binding transcriptional activator, homocysteine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,624,1850GA98.8% 150.4 / ‑5.6 84A107A (GCG→GCA) metRDNA‑binding transcriptional activator, homocysteine‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (83/0);  minor base T (1/0);  total (84/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCGTCTGCGCATTGCCATTGAGTGCCATAGCTGTATTCAGTGGCTGACACCCGCGTTAGAAA  >  W3110S.gb/3624131‑3624192
                                                      |       
gcgcCTGCGCATTGCCATAGAGTGCCTTAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:2372/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCGGCATTAGaaa  >  1:61673/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCTTTAGaaa  >  1:46575/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:59689/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:68710/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:68032/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:67500/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:67172/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:66275/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:65197/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:62294/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:6035/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:71553/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:57201/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:54396/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:54028/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:53492/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:52175/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:51334/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:51185/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:5007/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:49977/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:81512/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:9180/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:9078/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:87623/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:85992/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:85822/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:85354/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:84414/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:83406/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:81918/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:49119/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:81209/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:77047/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:7636/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:75399/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:74781/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:74516/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:74443/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:73497/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:73206/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:22577/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:11947/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:28643/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:28527/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:28514/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:27712/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:27465/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:2493/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:24414/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:23986/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:28720/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:22388/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:20897/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:20660/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:19855/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:19074/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:14372/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:13236/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:12731/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:12428/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:3998/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:47526/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:47290/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:46792/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:45638/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:45398/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:44293/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:43832/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:42890/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:41474/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:31611/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:39489/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:37796/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:37171/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:37017/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:35697/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:34003/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:33107/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTATTAAGTGGCTGACGCCCGCATTAGaaa  >  1:1574/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCATAGCTGTACTCAGTGGCTGACGCCCGCATTAGaaa  >  1:39338/1‑62 (MQ=255)
gcgcCTGCGCATTGCCATAGAGTGCCAAAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:72741/1‑62 (MQ=255)
gcgcCTGCGCATTGCCAGAGAGTGCCATAGCTGTATTCAGTGGCTGACGCCCGCATTAGaaa  >  1:65090/1‑62 (MQ=255)
                                                      |       
GCGTCTGCGCATTGCCATTGAGTGCCATAGCTGTATTCAGTGGCTGACACCCGCGTTAGAAA  >  W3110S.gb/3624131‑3624192

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: ‑