Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 429,759 T→A intergenic (+59/‑70) secF → / → yajD SecYEG protein translocase auxillary subunit/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb429,7590TA81.6% 104.8 / 19.6 49intergenic (+59/‑70)secF/yajDSecYEG protein translocase auxillary subunit/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/2);  new base A (0/40);  total (7/42)
Fisher's exact test for biased strand distribution p-value = 4.19e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

TCCCGTTGATGTTGAAAATCCCGGTCAGAAGATCGGGATTTTTTTTGATGTATGGAGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTC  >  W3110S.gb/429707‑429812
                                                    |                                                     
tCCCGTTGATGTTGAAAATCCCGGTCAGAAGATCGGGATTTTTTTTGATGTATGGAGATTg                                               >  1:2348104/1‑61 (MQ=255)
tCCCGTTGATGTTGAAAATCCCGGTCAGAAGATCGGGATTTTTTTTGATGTATGGAGATTg                                               >  1:527544/1‑61 (MQ=255)
tCCCGTTGATGTTGAAAATCCCGGTCAGAAGATCGGGATTTTTTTTGATGTATGGAGATTg                                               >  1:525627/1‑61 (MQ=255)
tCCCGTTGATGTTGAAAATCCCGGTCAGAAGATCGGGATTTTTTTTGATGTATGGAGATTg                                               >  1:2134999/1‑61 (MQ=255)
tCCCGTTGATGTTGAAAATCCCGGTCAGAAGATCGGGATTTTTTTTGATGTATGGAGATTg                                               >  1:1731472/1‑61 (MQ=255)
tCCCGTTGATGTTGAAAATCCCGGTCAGAAGATCGGGATTTTTTTTGATGTATGGAGATTg                                               >  1:2072952/1‑61 (MQ=255)
tCCCGTTGATGTTGAAAATCCCGGTCAGAAGATCGGGATTTTTTTTGATGTATGGAGATTg                                               >  1:1746963/1‑61 (MQ=255)
                                     aTTTTTTTTGATGTATGGAGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCt         <  1:133079/62‑1 (MQ=255)
                                     aTTTTTTTTGATGTATGGAGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCt         <  1:1635191/62‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:387628/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2321770/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2429130/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2529183/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2607741/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2674720/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2774684/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2781401/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2787180/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2804222/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2244557/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:444440/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:481940/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:523433/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:525994/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:681927/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:914108/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2286030/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:111645/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1188332/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1212526/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1331440/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1514773/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1618053/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1657256/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1736233/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1748739/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:177620/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1794408/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1808801/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1883844/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2006540/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2035986/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2107522/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2143410/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:1094580/61‑1 (MQ=255)
                                             tGATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2319410/61‑1 (MQ=255)
                                              gATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:2911600/60‑1 (MQ=255)
                                              gATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:496643/60‑1 (MQ=255)
                                              gATGTAAGGTGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTc  <  1:824696/60‑1 (MQ=255)
                                                    |                                                     
TCCCGTTGATGTTGAAAATCCCGGTCAGAAGATCGGGATTTTTTTTGATGTATGGAGATTGCGAAAATCCCCGCATCTTGGGAAACTGCGCGTAACCCTACATTTC  >  W3110S.gb/429707‑429812

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