Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,528,968 A→G 83.0% L242P (CTA→CCA)  yjhH ← predicted lyase/synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,528,9680AG83.0% 127.3 / 18.3 53L242P (CTA→CCA) yjhHpredicted lyase/synthase
Reads supporting (aligned to +/- strand):  ref base A (7/2);  new base G (29/15);  total (36/17)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

CTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCG  >  W3110S.gb/4528901‑4529030
                                                                   |                                                              
cTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGtagtag                                                               >  1:3536115/1‑69 (MQ=255)
  gCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtt                                                             >  1:210057/1‑69 (MQ=255)
   cATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGttt                                                            >  1:3590740/1‑69 (MQ=255)
   cATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGttt                                                            >  1:1947613/1‑69 (MQ=255)
   cATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGttt                                                            >  1:1799862/1‑69 (MQ=255)
     tGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGttttt                                                          <  1:3221005/69‑1 (MQ=255)
     tGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGttttt                                                          <  1:493132/69‑1 (MQ=255)
     tGCAGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGttttt                                                          <  1:78875/69‑1 (MQ=255)
      gCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGtttttt                                                         >  1:423190/1‑69 (MQ=255)
      gCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGtttttt                                                         >  1:2316022/1‑69 (MQ=255)
      gCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGtttttt                                                         >  1:540036/1‑69 (MQ=255)
         gTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTAt                                                       <  1:3332810/68‑1 (MQ=255)
         gTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTAt                                                       <  1:1899795/68‑1 (MQ=255)
         gTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTAt                                                       <  1:1886277/68‑1 (MQ=255)
         gTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTAt                                                       <  1:1258075/68‑1 (MQ=255)
            tttGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTAtt                                                      <  1:1160448/66‑1 (MQ=255)
            tttGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAg                                                   >  1:3065516/1‑69 (MQ=255)
            tttGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAg                                                   >  1:525758/1‑69 (MQ=255)
            tttGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAg                                                   >  1:2842123/1‑69 (MQ=255)
              tGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGc                                                  <  1:1887050/68‑1 (MQ=255)
              tGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCg                                                 >  1:114784/1‑69 (MQ=255)
                         aCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTc                                      >  1:934376/1‑69 (MQ=255)
                           aaaCGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAg                                                   <  1:1484315/54‑1 (MQ=255)
                           aaaCGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAg                                                   >  1:1430549/1‑54 (MQ=255)
                           aaaCGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAg                                                   <  1:3452041/54‑1 (MQ=255)
                           aaaCGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAg                                                   <  1:139259/54‑1 (MQ=255)
                           aaaCGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAg                                                   <  1:3020836/54‑1 (MQ=255)
                           aaaCGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAg                                                   >  1:2361382/1‑54 (MQ=255)
                           aaaCGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAg                                                   <  1:3539744/54‑1 (MQ=255)
                              cGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCg                                              >  1:2465680/1‑56 (MQ=255)
                              cGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCgg                                        >  1:1882288/1‑62 (MQ=255)
                              cGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCgg                                        >  1:2352158/1‑62 (MQ=255)
                              cGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCgg                                        >  1:149961/1‑62 (MQ=255)
                               ggTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTc                                      >  1:1101444/1‑63 (MQ=255)
                                    ttCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCg                                              >  1:2371519/1‑50 (MQ=255)
                                    ttCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAgag                                         >  1:2003383/1‑53 (MQ=255)
                                    ttCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCg                                         >  1:2880392/1‑55 (MQ=255)
                                     tCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGtttttt                                                         <  1:956852/38‑1 (MQ=255)
                                     tCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCa                                  >  1:1960864/1‑61 (MQ=255)
                                     tCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAg                                 >  1:1311028/1‑62 (MQ=255)
                                     tCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAg                                 >  1:1965265/1‑62 (MQ=255)
                                     tCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGcc                           >  1:874764/1‑68 (MQ=255)
                                      cGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCtt                         >  1:929552/1‑69 (MQ=255)
                                          ggCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCAc                      >  1:2977682/1‑68 (MQ=255)
                                           gCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCAc                      >  1:1051914/1‑67 (MQ=255)
                                           gCGTAAATAGAGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCAc                      >  1:1559732/1‑67 (MQ=255)
                                                aaTAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAgg                >  1:3582480/1‑68 (MQ=255)
                                                aaTAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGgc               >  1:2799412/1‑69 (MQ=255)
                                                     cGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCtt                         >  1:381771/1‑54 (MQ=255)
                                                      gggCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa         <  1:3267656/69‑1 (MQ=255)
                                                      gggCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa         <  1:3013345/69‑1 (MQ=255)
                                                         cAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGAt      >  1:1543420/1‑69 (MQ=255)
                                                             tGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGcc   >  1:3026880/1‑68 (MQ=255)
                                                             tGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCg  >  1:956152/1‑69 (MQ=255)
                                                              gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCtt                         >  1:3563317/1‑45 (MQ=255)
                                                              gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa         <  1:485860/61‑1 (MQ=255)
                                                              gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa         <  1:2534485/61‑1 (MQ=255)
                                                                   |                                                              
CTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCG  >  W3110S.gb/4528901‑4529030

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