Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 565,810 A→G 100% S9P (TCG→CCG)  grxA ← glutaredoxin 1, redox coenzyme for ribonucleotide reductase (RNR1a)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE565,8100AG95.1% 126.0 / ‑1.6 41S9P (TCG→CCG) grxAglutaredoxin 1, redox coenzyme for ribonucleotide reductase (RNR1a)
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (0/39);  total (0/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAA  >  minE/565767‑565865
                                           |                                                       
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:648155/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:2281562/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:2303888/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:254035/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:358157/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:429322/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:442886/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:459431/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:498195/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:583089/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:2199452/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:821367/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:86375/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:902265/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:902450/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:925906/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:934628/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:935936/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:956450/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:987511/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1785740/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1223355/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1323517/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1377593/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1388809/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1531975/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1569270/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1572515/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1668783/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1748898/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:118219/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1791887/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1893842/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:1983547/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:2036645/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:2115451/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:2131606/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:2137886/52‑1 (MQ=255)
ttCTCGGCCAGATCTTTGGCACGCACACAGTAAGGGCAACCCGGACGACCaa                                                 <  1:2139629/52‑1 (MQ=255)
                                aGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTaa  <  1:709305/67‑1 (MQ=255)
                                aGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTaa  <  1:1976671/67‑1 (MQ=255)
                                           |                                                       
TTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAA  >  minE/565767‑565865

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: