Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 970,899 A→G P363P (CCA→CCG ycaQ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb970,8990AG86.0% 131.4 / 13.4 50P363P (CCA→CCGycaQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (7/0);  new base G (43/0);  total (50/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.43e-01

GCCAGACAGGCATCCTTGAAGTTATCTCTCTGTGGTTACAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGC  >  W3110S.gb/970844‑970954
                                                       |                                                       
gCCAGACAGGCATCCTTGAAGTTATCTCTCTGTGGTTACAGGAAGGTATTAAACCAACGACAACg                                                >  1:556380/1‑65 (MQ=255)
gCCAGACAGGCATCCTTGAAGTTATCTCTCTGTGGTTACAGGAAGGTATTAAACCAACGACAACg                                                >  1:2432520/1‑65 (MQ=255)
gCCAGACAGGCATCCTTGAAGTTATCTCTCTGTGGTTACAGGAAGGTATTAAACCAACGACAACg                                                >  1:2173747/1‑65 (MQ=255)
 ccAGACAGGCATCCTTGAAGTTATCTCTCTGTGGTTACAGGAAGGTATTAAACcaacgccaacgc                                               >  1:2302936/1‑65 (MQ=255)
 ccAGACAGGCATCCTTGAAGTTATCTCTCTGTGGTTACAGGAAGGTATTAAACCAACGACAACGc                                               >  1:2223757/1‑65 (MQ=255)
 ccAGACAGGCATCCTTGAAGTTATCTCTCTGTGGTTACAGGAAGGTATTAAACCAACGACAACGc                                               >  1:2047966/1‑65 (MQ=255)
 ccAGACAGGCATCCTTGAAGTTATCTCTCTGTGGTTACAGGAAGGTATTAAACCAACGACAACGc                                               >  1:1822563/1‑65 (MQ=255)
 ccAGACAGGCATCCTTGAAGTTATCTCTCTGTGGTTACAGGAAGGTATTAAACCAAAGACAACGc                                               >  1:1262715/1‑65 (MQ=255)
                                              cagtAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:324943/4‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGAt               >  1:1683681/2‑52 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTgg   >  1:1390711/2‑64 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:555571/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:995744/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:2472626/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:249909/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:2501968/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:2506405/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:349685/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:375707/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:478255/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:919327/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:2250034/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:632266/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:665564/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:714430/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:790245/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:860870/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:867836/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:872425/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:882495/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1834482/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1084485/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1263001/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1267563/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1276465/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1530198/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1531726/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1547200/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1752629/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1812582/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1818099/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:2370859/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1924342/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1925913/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1927302/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1931494/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1945419/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1999734/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:1005402/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:2251279/2‑65 (MQ=255)
                                              caTTAAACCGACGACAACGCTGCAAAAAGGGTTACGTCAGGCGATTACTGATTTCGCTAACTGgc  >  1:2264735/2‑65 (MQ=255)
                                                       |                                                       
GCCAGACAGGCATCCTTGAAGTTATCTCTCTGTGGTTACAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGC  >  W3110S.gb/970844‑970954

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