Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 276,008:1 +AC 100% coding (167/891 nt) cyoE ← protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE276,0081.A94.7% 131.3 / 1.3 38coding (167/891 nt)cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (36/0);  total (38/0)
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.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*minE276,0082.C94.7% 136.4 / 0.6 38coding (167/891 nt)cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (36/0);  total (38/0)
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.

TCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAA‑‑ACACACAACCCGACGCCA  >  minE/275958‑276026
                                                   ||                  
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGcc   >  1:450852/1‑70 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:480081/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:113609/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:371811/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:398531/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:426453/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:444687/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:45819/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:47926/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:345929/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:506895/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:510401/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:533283/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:551668/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:611147/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:6671/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:81428/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:245324/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:133261/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:140756/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:132305/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:160266/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:191295/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:197209/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:218571/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:231214/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:339701/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:261131/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:281176/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:296597/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:323125/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:32627/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:125225/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAA‑‑Acacacacaa          >  1:410028/1‑58 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAA‑‑Acacaca             >  1:378590/1‑58 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGAAGTAGTTGTTAAACACACACAACCCGACGCca  >  1:454279/1‑71 (MQ=255)
tcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTAGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:149779/1‑71 (MQ=255)
tcttcGGCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACACAACCCGACGCca  >  1:343022/1‑71 (MQ=255)
                                                   ||                  
TCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAA‑‑ACACACAACCCGACGCCA  >  minE/275958‑276026

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: