Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,973,851 T→C 44.8% V51A (GTT→GCT)  metK → methionine adenosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,973,8510TC44.8% 0.2 / 20.2 29V51A (GTT→GCT) metKmethionine adenosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (8/8);  new base C (4/9);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

GCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTT  >  minE/1973787‑1973918
                                                                |                                                                   
gcgATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAACCTAAGTAAAAACCGGCATGGTTTTa                                                                 >  1:2430401/1‑69 (MQ=255)
gcgATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGAGAAACCTACGTAAAAACCGGCATGGCTTTa                                                                 >  1:1056150/1‑69 (MQ=255)
                 tCCGAAAGCACGCGTTGCTTGCGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATcac                                                <  1:664756/69‑1 (MQ=255)
                 tCCGAAAGCACGCGTTGCTTGCGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATcac                                                <  1:1931702/69‑1 (MQ=255)
                  ccGAAAGCACGCGTTGCTTGCGAAACCTACGTAAAAACCGGCATGGCTTTAGTTGGCGGCGAAATcacc                                               >  1:1025985/1‑69 (MQ=255)
                  ccGAAAGCACGCGTTGCTTGCGAAACCTACGTAAAAACCGGCATGGCTTTAGTTGGCGGCGAAATcacc                                               >  1:1011749/1‑69 (MQ=255)
                  ccGAAAGCACGCGTTGCTTGCGAAACCTACGTAAAAACCGGCATGGCTTTAGTTGGCGGCGAAATcacc                                               >  1:465170/1‑69 (MQ=255)
                     aaaGCACGCGTTGCTTGCGAAACCTACGTAAAAACTGGCATGGTTTTAGTTGGCGGCGAAATCaccacc                                            >  1:2853604/1‑69 (MQ=255)
                     aaaGCACGCGTTGCTTGCGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCaccacc                                            >  1:55907/1‑69 (MQ=255)
                     aaaGCACGCGTTGCTTGCGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCaccacc                                            >  1:3277008/1‑69 (MQ=255)
                     aaaGCACGCGTTGCTTGCGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCaccacc                                            >  1:2455565/1‑69 (MQ=255)
                             cgTTGCTTGCGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTg                                    <  1:3259929/69‑1 (MQ=255)
                                 gcttgcGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTa                                <  1:3072358/69‑1 (MQ=255)
                                 gcttgcGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTa                                <  1:2777578/69‑1 (MQ=255)
                                 gcttgcGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTa                                >  1:3243125/1‑69 (MQ=255)
                                 gcttgcGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTa                                <  1:1968709/69‑1 (MQ=255)
                                 gcttgcGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTa                                >  1:455645/1‑69 (MQ=255)
                                 gcttgcGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTa                                <  1:516722/69‑1 (MQ=255)
                                                      aCCGGCATGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGt           <  1:812531/69‑1 (MQ=255)
                                                         ggCATGGCTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAac        <  1:2541634/69‑1 (MQ=255)
                                                         ggCATGGCTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAac        <  1:2488214/69‑1 (MQ=255)
                                                         ggCATGGCTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAac        <  1:244055/69‑1 (MQ=255)
                                                         ggCATGGCTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAac        <  1:1961027/69‑1 (MQ=255)
                                                         ggCATGGCTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAac        <  1:1812028/69‑1 (MQ=255)
                                                         ggCATGGCTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAac        <  1:1695379/69‑1 (MQ=255)
                                                         ggCATGGCTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAac        <  1:1185533/69‑1 (MQ=255)
                                                         ggCATGGCTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAac        <  1:743670/69‑1 (MQ=255)
                                                         ggCATGGCTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAac        <  1:1130292/69‑1 (MQ=255)
                                                               gTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGtt  >  1:174399/1‑69 (MQ=255)
                                                                |                                                                   
GCGATCCTCGAACAGGATCCGAAAGCACGCGTTGCTTGCGAAACCTACGTAAAAACCGGCATGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTT  >  minE/1973787‑1973918

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: