Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 108,941 C→T 42.8% Y395Y (TAC→TAT aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE108,9410CT42.8% 10.8 / 27.9 42Y395Y (TAC→TATaceEpyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base C (19/5);  new base T (13/5);  total (32/10)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGT  >  minE/108878‑109007
                                                               |                                                                  
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:1155347/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:3411239/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:3374810/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:3176143/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:2924120/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:282914/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:2667870/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:26087/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:2469679/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:1916418/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:1770419/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATACCATTAAAGGTTATGGCAt                                                               >  1:1638299/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACGGTAATCCTTGCTCATAACATTAAAGGTTATGGCa                                                                >  1:1989052/1‑68 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCAt                                                               >  1:2020340/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCAt                                                               >  1:3382085/1‑69 (MQ=255)
cAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCAt                                                               >  1:2873249/1‑69 (MQ=255)
              aaaCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGc                                                   >  1:3323826/1‑67 (MQ=255)
              aaaCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGc                                                   >  1:317928/1‑67 (MQ=255)
              aaaCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGc                                                   >  1:1171331/1‑67 (MQ=255)
              aaaCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGc                                                   >  1:2888403/1‑67 (MQ=255)
              aaaCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGc                                                   >  1:1420819/1‑67 (MQ=255)
              aaaCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGc                                                   >  1:2505768/1‑67 (MQ=255)
              aaaCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGc                                                   >  1:1564038/1‑67 (MQ=255)
              aaaCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGc                                                   >  1:870585/1‑67 (MQ=255)
              aaaCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGc                                                   >  1:2189789/1‑67 (MQ=255)
              aaaCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGc                                                   >  1:1584721/1‑67 (MQ=255)
                    aaGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTa                                           <  1:2655107/69‑1 (MQ=255)
                                          tGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGtt                       >  1:199327/1‑67 (MQ=255)
                                          tGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTa                      >  1:2833301/1‑68 (MQ=255)
                                             tCATACCATTAAAGGTTATGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAg                    <  1:1104752/67‑1 (MQ=255)
                                             tCATACCATTAAAGGTTATGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAg                    <  1:1364345/67‑1 (MQ=255)
                                             tCATACCATTAAAGGTTATGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAg                    <  1:3372873/67‑1 (MQ=255)
                                              cATACCATTAAAGGTTATGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAg                    <  1:1770563/66‑1 (MQ=255)
                                              cATACCATTAAAGGTTATGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAg                    <  1:3040354/66‑1 (MQ=255)
                                                  ccATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATg              >  1:3246018/1‑68 (MQ=255)
                                                  ccATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATg              >  1:802053/1‑68 (MQ=255)
                                                  ccATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATg              >  1:887236/1‑68 (MQ=255)
                                                           gTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGaa            >  1:1187126/1‑61 (MQ=255)
                                                             tACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGgt  <  1:1311134/69‑1 (MQ=255)
                                                             tACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGgt  <  1:1422207/69‑1 (MQ=255)
                                                             tACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGgt  <  1:505323/69‑1 (MQ=255)
                                                             tACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGgt  <  1:2384167/69‑1 (MQ=255)
                                                               |                                                                  
CAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGT  >  minE/108878‑109007

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