Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 334,384 T→G 57.5% L334L (CTT→CTG 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,3840TG57.5% 21.8 / 41.1 47L334L (CTT→CTGushAbifunctional UDP‑sugar hydrolase and 5'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base T (13/7);  new base G (15/12);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

ACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAA  >  minE/334317‑334447
                                                                   |                                                               
aCCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCtt                                                                 >  1:2248114/1‑68 (MQ=255)
aCCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCtt                                                                 >  1:2219880/1‑68 (MQ=255)
aCCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCtt                                                                 >  1:3160778/1‑68 (MQ=255)
 ccAGCTGATTCCGTTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTa                                                               >  1:820666/1‑69 (MQ=255)
   aGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTa                                                               <  1:3298609/67‑1 (MQ=255)
   aGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTa                                                               <  1:2984266/67‑1 (MQ=255)
   aGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTa                                                               <  1:2734297/67‑1 (MQ=255)
   aGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTa                                                               <  1:1685310/67‑1 (MQ=255)
    gCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTacac                                                            <  1:1247524/69‑1 (MQ=255)
           ccGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGaa                                                     >  1:1398112/1‑69 (MQ=255)
           ccGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGaa                                                     >  1:382013/1‑69 (MQ=255)
                            gaaAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCa                                    >  1:3043854/1‑69 (MQ=255)
                            gaaAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCa                                    >  1:616795/1‑69 (MQ=255)
                                gTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTa                                                               <  1:1767168/38‑1 (MQ=255)
                                     cTGTGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGAtctc                           <  1:227217/69‑1 (MQ=255)
                                       gggAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGc                         >  1:1426018/1‑69 (MQ=255)
                                       gggAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGc                         >  1:675789/1‑69 (MQ=255)
                                       gggAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGc                         >  1:2636910/1‑69 (MQ=255)
                                       gggAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGc                         >  1:3332461/1‑69 (MQ=255)
                                            gACGGGAAAAGCGAGCTCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                     >  1:3201201/1‑68 (MQ=255)
                                            gACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                     >  1:3183437/1‑68 (MQ=255)
                                            gACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                     >  1:2915451/1‑68 (MQ=255)
                                            gACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                     >  1:379376/1‑68 (MQ=255)
                                               gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCTAATGATCTCGCTGTTATCa                 <  1:118821/69‑1 (MQ=255)
                                               gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                 <  1:2193929/69‑1 (MQ=255)
                                               gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                 <  1:2494275/69‑1 (MQ=255)
                                               gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                 <  1:2468296/69‑1 (MQ=255)
                                               gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                 <  1:2260275/69‑1 (MQ=255)
                                               gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                 <  1:2190042/69‑1 (MQ=255)
                                               gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                 <  1:570269/69‑1 (MQ=255)
                                               gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                 <  1:1814413/69‑1 (MQ=255)
                                               gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                 <  1:1125340/69‑1 (MQ=255)
                                               gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                 <  1:1096979/69‑1 (MQ=255)
                                               gggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                 <  1:1027583/69‑1 (MQ=255)
                                                ggAAAAGCGAGCGCGTGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCa                 <  1:1090171/68‑1 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTTTTATCACCGTTCCAGAACaaa  >  1:1339463/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:3151806/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:2811550/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:2747424/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:2635631/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:2582019/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:230629/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:2300405/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:2189402/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:2105515/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:611639/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:159825/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:648150/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:1409401/1‑68 (MQ=255)
                                                               tGCTGTACACCCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa  >  1:956164/1‑68 (MQ=255)
                                                                   |                                                               
ACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAA  >  minE/334317‑334447

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: