Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 531,446 C→T 28.0% intergenic (‑1/‑229) ybbS ← / → allA transcriptional activator of the allD operon/ureidoglycolate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb531,4460CT28.0% 36.9 / 15.4 25intergenic (‑1/‑229)ybbS/allAtranscriptional activator of the allD operon/ureidoglycolate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (12/6);  new base T (3/4);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTT  >  W3110S.gb/531382‑531501
                                                                |                                                       
cTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAgg                                                   <  1:1900243/71‑1 (MQ=255)
     gAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATc                                               <  1:1148877/70‑1 (MQ=255)
                       gCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGt                            <  1:1716360/71‑1 (MQ=255)
                       gCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGt                            >  1:1680903/1‑71 (MQ=255)
                       gCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGt                            >  1:1664926/1‑71 (MQ=255)
                       gCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGt                            <  1:1378802/71‑1 (MQ=255)
                                   aaaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTc                >  1:1847505/1‑71 (MQ=255)
                                   aaaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTc                >  1:963330/1‑71 (MQ=255)
                                   aaaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTc                >  1:869015/1‑71 (MQ=255)
                                   aaaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTc                >  1:674102/1‑71 (MQ=255)
                                   aaaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTc                >  1:357941/1‑71 (MQ=255)
                                   aaaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTc                >  1:1832696/1‑71 (MQ=255)
                                   aaaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTc                >  1:1787729/1‑71 (MQ=255)
                                   aaaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTc                >  1:1076488/1‑71 (MQ=255)
                                   aaaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTc                >  1:1414531/1‑71 (MQ=255)
                                   aaaGTCCGCAAGGTTCCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTc                >  1:1299067/1‑71 (MQ=255)
                                      gTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa             <  1:1707594/71‑1 (MQ=255)
                                      gTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa             <  1:1347154/71‑1 (MQ=255)
                                      gTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa             <  1:628181/71‑1 (MQ=255)
                                       tCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa             <  1:67045/70‑1 (MQ=255)
                                               gTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGAt    >  1:1620658/1‑71 (MQ=255)
                                               gTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGAt    >  1:1533611/1‑71 (MQ=255)
                                               gTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCATCTGAt    >  1:603490/1‑71 (MQ=255)
                                                  tCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGAttt  <  1:1973136/70‑1 (MQ=255)
                                                           aaCATCTTCAGGTATCCCCTTTTAAATCCGCAAGt                            <  1:1280857/35‑1 (MQ=255)
                                                                |                                                       
CTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTT  >  W3110S.gb/531382‑531501

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: