Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,925,224 T→A K147N (AAA→AAT holC ← DNA polymerase III, chi subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,925,2240TA96.5% 88.8 / ‑3.4 29K147N (AAA→AATholCDNA polymerase III, chi subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (28/0);  total (28/1)
Fisher's exact test for biased strand distribution p-value = 3.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

AGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTATTTCCAGGTTGCCGTATTCA  >  minE/2925181‑2925243
                                           |                   
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:552267/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:997615/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:989132/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:910363/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:90696/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:899779/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:898190/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:839415/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:785506/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:672620/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:659568/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:653122/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:621096/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:61611/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:592268/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:1015870/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:467787/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:403980/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:343537/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:232052/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:173803/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:143196/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:114734/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:1084580/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:1079722/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:107423/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:106235/1‑62 (MQ=255)
aGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTAATTCCAGGTTGCCGTATTc   >  1:102232/1‑62 (MQ=255)
 gCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTATTTCCAGGTTGCCGTATTCa  <  1:401581/62‑1 (MQ=255)
                                           |                   
AGCGGCTGTTCGATATCTTGTGGGTTATATGTCTTTTCCATTATTTCCAGGTTGCCGTATTCA  >  minE/2925181‑2925243

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: