Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,602,049 +C coding (1700/2043 nt) prlC → oligopeptidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,602,0461.C96.7% 100.0 / ‑5.5 30R566P (CGC→CCC) prlColigopeptidase A
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (29/0);  minor base A (1/0);  total (30/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

TATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCGGATCAGGGG  >  minE/2601990‑2602059
                                                         |             
tATTCTGCTTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:1260699/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCg                >  1:1173952/1‑57 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:809395/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:299747/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:363005/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:507331/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:517329/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:560998/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:586728/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:642616/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:732245/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:1561429/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:83172/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:890885/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:916450/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:968382/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:968588/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:1015180/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:147470/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:1416297/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:1382641/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:1377465/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:137679/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:1375774/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:1354984/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:1187169/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:117647/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:1175822/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:110944/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCCGGATCAgggg  >  1:110213/1‑71 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGACCCGGATCAggg   >  1:18430/1‑70 (MQ=255)
tATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGccccgc         >  1:1477015/1‑60 (MQ=255)
                                                         |             
TATTCTGCGTCAGCTGGAGTTCGGCCTGTTTGATTTCCGCCTTCATGCCGAGTTCCGCCCGGATCAGGGG  >  minE/2601990‑2602059

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: