Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 334,390 T→C 100% T336T (ACT→ACC ushA → bifunctional UDP‑sugar hydrolase and 5'‑nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE334,3900TC76.4% 74.9 / 17.2 59T336T (ACT→ACCushAbifunctional UDP‑sugar hydrolase and 5'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base T (13/1);  new base C (33/12);  total (46/13)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCG  >  minE/334328‑334456
                                                              |                                                                  
ccGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGaa                                                              >  1:382013/1‑69 (MQ=255)
ccGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGaa                                                              >  1:1398112/1‑69 (MQ=255)
                 gaaAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCa                                             >  1:3043854/1‑69 (MQ=255)
                 gaaAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCa                                             >  1:616795/1‑69 (MQ=255)
                          cTGTGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGAtctc                                    <  1:227217/69‑1 (MQ=255)
                            gggAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGc                                  >  1:2636910/1‑69 (MQ=255)
                            gggAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGc                                  >  1:675789/1‑69 (MQ=255)
                            gggAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGc                                  >  1:3332461/1‑69 (MQ=255)
                            gggAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGc                                  >  1:1426018/1‑69 (MQ=255)
                                 gACGGGAAAAGCGAGCTCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                              >  1:3201201/1‑68 (MQ=255)
                                 gACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                              >  1:3183437/1‑68 (MQ=255)
                                 gACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                              >  1:379376/1‑68 (MQ=255)
                                 gACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                              >  1:2915451/1‑68 (MQ=255)
                                    gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCTAATGATCTCGCTGTTATCa                          <  1:118821/69‑1 (MQ=255)
                                    gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                          <  1:2468296/69‑1 (MQ=255)
                                    gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                          <  1:2494275/69‑1 (MQ=255)
                                    gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                          <  1:2260275/69‑1 (MQ=255)
                                    gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                          <  1:2190042/69‑1 (MQ=255)
                                    gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                          <  1:2193929/69‑1 (MQ=255)
                                    gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                          <  1:1027583/69‑1 (MQ=255)
                                    gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                          <  1:1096979/69‑1 (MQ=255)
                                    gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                          <  1:1125340/69‑1 (MQ=255)
                                    gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                          <  1:1814413/69‑1 (MQ=255)
                                    gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                          <  1:570269/69‑1 (MQ=255)
                                     ggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                          <  1:1090171/68‑1 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTTTTATCACCGTTCCAGAACaaa           >  1:1339463/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:611639/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:2811550/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:2747424/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:2189402/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:956164/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:2635631/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:2582019/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:2105515/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:1409401/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:159825/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:648150/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:230629/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:2300405/1‑68 (MQ=255)
                                                    tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa           >  1:3151806/1‑68 (MQ=255)
                                                           cacaCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:326007/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGt                               >  1:2526150/1‑41 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:626155/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:981354/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:527672/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:2629150/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:2617834/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:2399065/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:2202302/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:2124727/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:1922773/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:1817375/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:1791661/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:1787769/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:1640846/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:1307672/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:117097/1‑63 (MQ=255)
                                                           cacCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:108861/1‑63 (MQ=255)
                                                           cacCCCAGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAgg         >  1:441821/1‑63 (MQ=255)
                                                             cTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAgcg  >  1:2668989/1‑68 (MQ=255)
                                                              |                                                                  
CCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCG  >  minE/334328‑334456

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