Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,824,280 T→C pseudogene (413/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,824,2800TC94.3% 113.7 / ‑1.4 35pseudogene (413/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (15/18);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGA  >  W3110S.gb/3824224‑3824333
                                                        |                                                     
tCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGAt                                                   <  1:1228723/61‑1 (MQ=255)
tCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGAt                                                   <  1:2470638/61‑1 (MQ=255)
tCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGAt                                                   <  1:620248/61‑1 (MQ=255)
tCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGAt                                                   <  1:4475482/61‑1 (MQ=255)
tCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGAt                                                   <  1:4594106/61‑1 (MQ=255)
tCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGAt                                                   <  1:1971511/61‑1 (MQ=255)
     cTGATGGTGGCACGCGTACCGCGTCGATTACGGGGGCCTGCGTGGCGCTGGCAGAt                                                   <  1:4629118/56‑1 (MQ=39)
             ggCACGCGTACCGCGTCGATTACGTGTGCCTGCGGGGCGCTGGCAGATGCGCTACAGAAGCt                                     <  1:1478199/62‑1 (MQ=255)
             ggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                     <  1:3268744/62‑1 (MQ=255)
             ggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                     <  1:3842325/62‑1 (MQ=255)
             ggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                     <  1:185029/62‑1 (MQ=255)
             ggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                     <  1:1603598/62‑1 (MQ=255)
             ggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                     <  1:1436695/62‑1 (MQ=255)
             ggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                     <  1:1049737/62‑1 (MQ=255)
                     tACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                              >  1:2157177/1‑61 (MQ=255)
                              gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                     >  1:4546468/1‑61 (MQ=255)
                              gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                     >  1:4278134/1‑61 (MQ=255)
                              gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                     >  1:3276071/1‑61 (MQ=255)
                              gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                     >  1:880862/1‑61 (MQ=255)
                              gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                     >  1:3117281/1‑61 (MQ=255)
                              gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                     >  1:290203/1‑61 (MQ=255)
                              gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                     >  1:2381022/1‑61 (MQ=255)
                              gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                     >  1:2083746/1‑61 (MQ=255)
                              gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                     >  1:1548476/1‑61 (MQ=255)
                              gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                     >  1:1404871/1‑61 (MQ=255)
                                     ggTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa              >  1:3527451/1‑61 (MQ=255)
                                         ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaac         >  1:2135517/1‑61 (MQ=255)
                                         ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaac         >  1:4377165/1‑61 (MQ=255)
                                             cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCg      >  1:54993/1‑61 (MQ=255)
                                                ggCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:2031304/62‑1 (MQ=255)
                                                ggCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:1827597/62‑1 (MQ=255)
                                                 gcgcTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                   <  1:3198958/44‑1 (MQ=255)
                                                 gcgcTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                   <  1:413250/44‑1 (MQ=255)
                                                 gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:3464151/61‑1 (MQ=255)
                                                 gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:1495501/61‑1 (MQ=255)
                                                        |                                                     
TCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGA  >  W3110S.gb/3824224‑3824333

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: