Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,527,855 G→A 100% S280S (TCG→TCA trmC → fused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,527,8550GA76.2% 53.5 / 12.9 46S280S (TCG→TCAtrmCfused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base G (1/10);  new base A (34/1);  total (35/11)
Fisher's exact test for biased strand distribution p-value = 2.89e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACT  >  minE/1527792‑1527917
                                                               |                                                              
tatgggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTgg                                                            <  1:2390052/65‑1 (MQ=255)
   gggCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCt                                                        <  1:138802/69‑1 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCt                                                        >  1:615609/1‑68 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:2606371/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:928641/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:450724/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:357443/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:3370903/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:3366634/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:3012588/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:2620802/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:254795/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:2503185/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:2457787/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:2169908/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:1253366/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:1585926/1‑69 (MQ=255)
    ggCAGCAACAAACAGGAAGCGGCAATTATTGGAGGCGGTATTGCCAGCGCGTTGTTGTCACTGGCGCTa                                                       >  1:1868333/1‑69 (MQ=255)
     gcagcagcaAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTAt                                                      <  1:57909/69‑1 (MQ=255)
             aaaCGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGgcgc                                              <  1:3078851/69‑1 (MQ=255)
             aaaCGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGgcgc                                              <  1:547191/69‑1 (MQ=255)
             aaaCGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGgcgc                                              <  1:546325/69‑1 (MQ=255)
             aaaCGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGgcgc                                              <  1:2029145/69‑1 (MQ=255)
                                              gCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAg             <  1:2085384/69‑1 (MQ=255)
                                              gCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAg             <  1:902956/69‑1 (MQ=255)
                                              gCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAg             <  1:2788132/69‑1 (MQ=255)
                                              gCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAg             <  1:2797754/69‑1 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAg                                     >  1:1304317/1‑42 (MQ=37)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTGTTGCGCgg                  >  1:3136566/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTAt                         >  1:1907035/1‑54 (MQ=38)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:1215840/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:1300823/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:2469039/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:1343247/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:1490643/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:3279625/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:3144857/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:3098534/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:2124634/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:307018/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:1594879/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:3009267/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:1708951/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGCGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:2602510/1‑61 (MQ=255)
                                                 aGCGCGTTGTTGTCACTGGAGCTATTGCGTCGCGGCTGGCAGGTAACGCTTTATTGCGCgg                  >  1:2439845/1‑61 (MQ=38)
                                                         gttgtCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACt  >  1:3278879/1‑69 (MQ=255)
                                                               |                                                              
CACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACT  >  minE/1527792‑1527917

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: