Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 647,788 G→T S336* (TCG→TAG)  citF ← citrate lyase, citrate‑ACP transferase (alpha) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb647,7880GT96.2% 66.4 / ‑4.0 26S336* (TCG→TAG) citFcitrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (25/0);  total (26/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

CGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAG  >  W3110S.gb/647779‑647840
         |                                                    
cGGGCCAGCTATTGCGCAGCATGTCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:323841/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCa   >  1:756926/1‑61 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCa   >  1:195482/1‑61 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCa   >  1:319255/1‑61 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:890204/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:1195693/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:805353/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:772000/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:744789/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:70004/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:598226/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:591432/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:451947/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:442764/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:420505/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:251302/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:1623256/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:1591151/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:1562233/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:1424205/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:1415313/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:1402512/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:1349812/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:1263617/1‑62 (MQ=255)
cGGGCCAGCTATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAg  >  1:1217478/1‑62 (MQ=255)
      aGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGttt           >  1:308212/1‑47 (MQ=255)
         |                                                    
CGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAG  >  W3110S.gb/647779‑647840

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: