Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,965,8010CA32.7% 97.0 / 38.2 52P428Q (CCG→CAG) deoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (35/0);  new base A (17/0);  total (52/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 = 9.32e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAA  >  minE/2965757‑2965841
                                            |                                        
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACt                          >  1:100612/1‑61 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACt                          >  1:453989/1‑61 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:220577/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:633596/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:660975/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:474045/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:433155/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:679299/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:690979/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:766774/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:769878/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:197067/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:105307/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:119628/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:1171723/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:932720/1‑62 (MQ=255)
ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg                         >  1:1188750/1‑62 (MQ=255)
                      aTTAACCTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:173085/1‑62 (MQ=255)
                      aTTAAACTTGCGGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:184340/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGc        >  1:481250/1‑57 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:587165/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:974280/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:638374/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:954581/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:67088/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:680792/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:93204/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:730385/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:877450/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:987919/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:799552/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:859015/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:245759/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:1007887/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:1016290/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:1055221/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:1193190/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:1193969/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:136474/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:195693/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:23055/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:565061/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:320220/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:333295/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:340736/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:343164/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:373353/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:490562/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:498802/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:53880/1‑62 (MQ=255)
                      aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa   >  1:543447/1‑62 (MQ=255)
                      aTTAAACTTCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa  >  1:200290/1‑62 (MQ=255)
                                            |                                        
GGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAA  >  minE/2965757‑2965841

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: