Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,301,534 C→T 100% L179L (CTG→CTA metL ← fused aspartokinase II and homoserine dehydrogenase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,301,5340CT93.1% 24.1 / ‑4.8 14L179L (CTG→CTAmetLfused aspartokinase II and homoserine dehydrogenase II
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (8/5);  total (9/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTTTCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTGCCCGGATGTTGCACCAGCAGCTGTTGCAGCAACGGGTAAGAAAGCCCTTCATCAACCTGCGGTTG  >  minE/2301466‑2301581
                                                                    |                                               
gTTTCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTGCCCGGATGTTGCACcagc                                                 >  1:2764215/1‑69 (MQ=255)
   tCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTGCCCGGATGTTGCACCagcagc                                              >  1:890958/1‑69 (MQ=255)
   tCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTGCCCGGATGTTGCACCagcagc                                              >  1:849182/1‑69 (MQ=255)
   tCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTGCCCGGATGTTGCACCagcagc                                              >  1:1701493/1‑69 (MQ=255)
   tCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTGCCCGGATGTTGCACCagcagc                                              >  1:1749173/1‑69 (MQ=255)
   tCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTGCCCGGATGTTGCACCagcagc                                              >  1:66095/1‑69 (MQ=255)
   tCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTGCCCGGATGTTGCACCagcagc                                              >  1:1866413/1‑69 (MQ=255)
   tCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTGCCCGGATGTTGCACCagcagc                                              >  1:1926906/1‑69 (MQ=255)
   tCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTGCCCGGATGTTGCACCagcagc                                              >  1:2623304/1‑69 (MQ=255)
   tCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTCCCGGATGTTGCACCAGCAGCt                                             >  1:183073/1‑69 (MQ=255)
                      tGATAAACCCGGTCACCACCAGACGTTTGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTac                            >  1:1147370/1‑67 (MQ=255)
                      tGATAAACCCGGTCACCACCAGACGTTTGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTac                            >  1:773682/1‑67 (MQ=255)
                      tGATAAACCCGGTCACCACCAGACGTTTGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTac                            >  1:348318/1‑67 (MQ=255)
                      tGATAAACCCGGTCACCACCAGACGTTTGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTac                            >  1:2986602/1‑67 (MQ=255)
                      tGATAAACCCGGTCACCACCAGACGTTTGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTac                            >  1:2323829/1‑67 (MQ=255)
                      tGATAAACCCGGTCACCACCAGACGTTTGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTac                            >  1:1912541/1‑67 (MQ=255)
                      tGATAAACCCGGTCACCACCAGACGTTTGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTac                            >  1:1665312/1‑67 (MQ=255)
                      tGATAAACCCGGTCACCACCAGACGTTTGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTac                            >  1:1089442/1‑67 (MQ=255)
                                               tttGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTACGAAAGCCCTTCATCAACCTGCGGTtg  <  1:2574661/69‑1 (MQ=255)
                                               tttGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTACGAAAGCCCTTCATCAACCTGCGGTtg  <  1:1935744/69‑1 (MQ=255)
                                               tttGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTACGAAAGCCCTTCATCAACCTGCGGTtg  <  1:1838756/69‑1 (MQ=255)
                                               tttGCCCGGATGTTGTACCAGTAGCTGTTGCAGCAACGGGTACGAAAGCCCTTCATCAACCTGCGGTtg  <  1:1240148/69‑1 (MQ=255)
                                               tttGCCCGGATGTTGTACCAGTAGCTGTTGCAACAACGGGTACGAAAGCCCTTCATCAACCTGCGGTtg  <  1:324370/69‑1 (MQ=255)
                                                                    |                                               
GTTTCACCGGCGTTGTTGCGGCTGATAAATCCGGTCACCACCAGACGTTTGCCCGGATGTTGCACCAGCAGCTGTTGCAGCAACGGGTAAGAAAGCCCTTCATCAACCTGCGGTTG  >  minE/2301466‑2301581

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: