Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,511,656 T→C 100% G308G (GGT→GGC pta → phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,511,6560TC88.0% 68.0 / 2.6 25G308G (GGT→GGCptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (22/0);  total (24/1)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTC  >  minE/1511596‑1511701
                                                            |                                             
ggCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATg                                      >  1:1017238/1‑70 (MQ=255)
ggCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATg                                      >  1:2734380/1‑70 (MQ=255)
                         ggCGTATAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:4512/1‑43 (MQ=38)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:1032893/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:698921/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:669123/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:626162/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:605430/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:407298/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:2878383/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:2814880/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:2759162/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:2650753/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:2576245/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:2431401/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:2218364/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:1799757/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:1762886/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:1634469/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:1583537/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:1509076/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:1462409/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:1366666/1‑43 (MQ=255)
                         ggCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGaaa                                        >  1:1359760/1‑43 (MQ=255)
                                            tgctgctgACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTc  <  1:993881/62‑1 (MQ=255)
                                                            |                                             
GGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTC  >  minE/1511596‑1511701

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