Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 4,149,789 C→T R67R (CGC→CGT pabA → aminodeoxychorismate synthase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,149,7890CT75.0% 58.3 / 19.1 36R67R (CGC→CGTpabAaminodeoxychorismate synthase, subunit II
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (27/0);  total (31/5)
Fisher's exact test for biased strand distribution p-value = 3.34e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTGTACGCCAGATGAAGCCGGGATCTCCCTTGACGTTATTCGCCACTATGCCGGGCGCTTGCCGATTCTTGGCGTCTGCCTCGGTCATCAGGCAATGG  >  W3110S.gb/4149745‑4149844
                                            |                                                       
cccTGTACGCCAGATGAAGCCGGGATCTCCCTTGACGTTATTCGCCACTATGCCGGGCGCt                                         <  1:118974/61‑1 (MQ=255)
cccTGTACGCCAGATGAAGCCGGGATCTCCCTTGACGTTATTCGCCACTATGCCGGGCGCt                                         <  1:2373792/61‑1 (MQ=255)
cccTGTACGCCAGATGAAGCCGGGATCTCCCTTGACGTTATTCGCCACTATGCCGGGCGCt                                         <  1:2351829/61‑1 (MQ=255)
cccTGTACGCCAGATGAAGCCGGGATCTCCCTTGACGTTATTCGCCACTATGCCGGGCGCt                                         <  1:1936957/61‑1 (MQ=255)
cccTGTACGCCAGATGAAGCCGGGAGCTCCCTTGACGTTATTCGCCACTATGCCGGGCGCt                                         <  1:378607/61‑1 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:2957537/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:99656/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:97524/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:857388/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:816134/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:799422/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:589476/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:532309/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:474309/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:383718/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:3089252/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:3042071/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:3014242/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:3006425/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:1300453/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:291619/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:1011940/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:124184/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:2305925/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:2253846/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:2117857/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:2082107/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:2017329/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:1989992/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:1254735/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:1797704/1‑45 (MQ=255)
                          ctcCCTTGACGTTATTCGTCACTATGCCGGACGCTTGCCGATTCt                               >  1:1358862/1‑45 (MQ=255)
                                      tATTCGCCACTATGCCGGGCGCTTGCCGATTCTTGGCGTCTGCCTCGGTCATCAGGc       >  1:174207/1‑57 (MQ=255)
                                      tATTCGCCACTATGCCGGGCGCTTGCCGATTCTTGGCGTCTGCCTCGGTCATCAGGc       >  1:1404835/1‑57 (MQ=255)
                                      tATTCGCCACTATGCCGGGCGCTTGCCGATTCTTGGCGTCTGCCTCGGTCATCAGGCAATgg  >  1:932434/1‑62 (MQ=255)
                                      tATTCGCCACTATGCCGGGCGCTTGCCGATTCTTGGCGTCTGCCTCGGTCATCAAGc       >  1:2006456/1‑57 (MQ=255)
                                            |                                                       
CCCTGTACGCCAGATGAAGCCGGGATCTCCCTTGACGTTATTCGCCACTATGCCGGGCGCTTGCCGATTCTTGGCGTCTGCCTCGGTCATCAGGCAATGG  >  W3110S.gb/4149745‑4149844

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: