Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,842,719 T→C A133A (GCA→GCG ynjF ← predicted phosphatidyl transferase, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,842,7190TC76.7% 75.2 / 19.7 43A133A (GCA→GCGynjFpredicted phosphatidyl transferase, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (10/0);  new base C (0/33);  total (10/33)
Fisher's exact test for biased strand distribution p-value = 5.22e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

AAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAAACAACA  >  W3110S.gb/1842673‑1842774
                                              |                                                       
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:414930/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2263935/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:1067332/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2495285/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2522104/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2651462/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2665161/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2802340/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2882137/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:377002/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2296859/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:509305/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:550325/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:614446/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:665667/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:887119/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:915861/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:962792/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:978248/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:1738041/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:1082057/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:1184289/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:1229842/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:1235477/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:1535216/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:1552446/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2254759/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:1756840/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:1947998/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2135986/62‑2 (MQ=255)
aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2171842/62‑2 (MQ=255)
 aaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:2751992/61‑2 (MQ=255)
 aaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa                                          <  1:1912078/61‑2 (MQ=255)
                                         tttGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAaacaaca  >  1:2721046/1‑61 (MQ=255)
                                         tttGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAaacaaca  >  1:1685388/1‑61 (MQ=255)
                                         tttGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAaacaaca  >  1:2008024/1‑61 (MQ=255)
                                         tttGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAaacaaca  >  1:2520428/1‑61 (MQ=255)
                                         tttGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAaacaaca  >  1:1512823/1‑61 (MQ=255)
                                         tttGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAaacaaca  >  1:1439973/1‑61 (MQ=255)
                                         tttGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAaacaaca  >  1:616534/1‑61 (MQ=255)
                                         tttGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAaacaaca  >  1:124931/1‑61 (MQ=255)
                                         tttGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAaacaaca  >  1:2434966/1‑61 (MQ=255)
                                         tttGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAaacaaca  >  1:2145513/1‑61 (MQ=255)
                                              |                                                       
AAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAAACAACA  >  W3110S.gb/1842673‑1842774

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: