Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,842,680 C→T K146K (AAG→AAA 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,6800CT84.6% 86.4 / 8.6 39K146K (AAG→AAAynjFpredicted phosphatidyl transferase, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base T (0/33);  total (0/39)
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.63e-01

CACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAG  >  W3110S.gb/1842623‑1842734
                                                         |                                                      
cacaAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTAt                                                    <  1:1106710/62‑1 (MQ=255)
cacaAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTAt                                                    <  1:964615/62‑1 (MQ=255)
cacaAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTAt                                                    <  1:232165/62‑1 (MQ=255)
cacaAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTAt                                                    <  1:672830/62‑1 (MQ=255)
cacaAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTAt                                                    <  1:49466/62‑1 (MQ=255)
cacaAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTAt                                                    <  1:1955581/62‑1 (MQ=255)
                                                  aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa  <  1:550325/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:414930/62‑2 (MQ=255)
                                                  aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa  <  1:509305/62‑2 (MQ=255)
                                                  aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa  <  1:2495285/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:1067332/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:1738041/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)
                                                  aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa  <  1:2254759/62‑2 (MQ=255)
                                                  aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa  <  1:2263935/62‑2 (MQ=255)
                                                  aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa  <  1:2296859/62‑2 (MQ=255)
                                                  aaaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa  <  1:2522104/62‑2 (MQ=255)
                                                   aaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa  <  1:1912078/61‑2 (MQ=255)
                                                   aaaCGATTTATGCGCATAACCGGGATTATCAATCTGATGTTTCGCCGCCAGCGCAGCaaaa  <  1:2751992/61‑2 (MQ=255)
                                                         |                                                      
CACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAG  >  W3110S.gb/1842623‑1842734

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: