Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,224,816 C→T 31.2% H147Y (CAT→TAT)  ebgC → cryptic beta‑D‑galactosidase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,224,8160CT31.2% 75.2 / 37.8 52H147Y (CAT→TAT) ebgCcryptic beta‑D‑galactosidase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (24/11);  new base T (8/8);  total (33/19)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGG  >  W3110S.gb/3224749‑3224880
                                                                   |                                                                
tCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTat                                                                <  1:503805/70‑1 (MQ=255)
     ttATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaa                                                          >  1:3864767/1‑71 (MQ=255)
     ttATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaa                                                          >  1:4370253/1‑71 (MQ=255)
     ttATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaa                                                          >  1:356019/1‑71 (MQ=255)
            cAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATaacaac                                                   <  1:1291798/71‑1 (MQ=255)
              aTAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTa                                                 >  1:2630586/1‑71 (MQ=255)
                aaCGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTAcg                                               >  1:4689746/1‑71 (MQ=255)
                 aCGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaataa                                                       >  1:4609022/1‑62 (MQ=255)
                 aCGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTAcgg                                              >  1:3200796/1‑71 (MQ=255)
                 aCGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTAcgg                                              >  1:46428/1‑71 (MQ=255)
                      gTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCa                                         <  1:1891800/71‑1 (MQ=255)
                         aaaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaataa                                                       >  1:4894565/1‑54 (MQ=255)
                          aaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAg                                     >  1:4214105/1‑71 (MQ=255)
                             aaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTa                                                 >  1:774016/1‑56 (MQ=255)
                               gTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGtt                                <  1:5130091/71‑1 (MQ=255)
                                 ggTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGttt                               <  1:2623580/70‑1 (MQ=255)
                                    tctcAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAg                                  <  1:2270638/64‑1 (MQ=255)
                                    tctcAAAGTCACCATCGAAGATGGTTATTTCCataacaaataacaa                                                    >  1:1642020/1‑46 (MQ=255)
                                    tctcAAAGTCACCATCGAAGATGGTTATTTCCataacaaataac                                                      >  1:2635458/1‑44 (MQ=255)
                                    tctcAAAGTCACCATCGAAGATGGTTATTTCCataacaaataac                                                      >  1:3861205/1‑44 (MQ=255)
                                     ctcAAAGTCACCATCGAAGATGGTTATTTCTataacaaataac                                                      <  1:3032025/43‑1 (MQ=255)
                                     ctcAAAGTCACCATCGAAGATGGTTATTTCTataacaaataac                                                      <  1:1538712/43‑1 (MQ=255)
                                     ctcAAAGTCACCATCGAAGATGGTTATTTCCataacaaataac                                                      <  1:4076800/43‑1 (MQ=255)
                                     ctcAAAGTCACCATCGAAGATGGTTATTTCCataacaaataac                                                      <  1:1750649/43‑1 (MQ=255)
                                       cAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCa                        >  1:4661731/1‑71 (MQ=255)
                                       cAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCa                        >  1:1888306/1‑71 (MQ=255)
                                        aaaGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgccgcc                         <  1:3014014/69‑1 (MQ=255)
                                            tCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGc                    >  1:549385/1‑70 (MQ=255)
                                                 aTCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTAccc               >  1:1271553/1‑70 (MQ=255)
                                                  tCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCaa                                        <  1:1899189/44‑1 (MQ=255)
                                                  tCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTa             >  1:1174104/1‑71 (MQ=255)
                                                    gAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTg                              >  1:941365/1‑52 (MQ=255)
                                                    gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgcc                            >  1:4495287/1‑54 (MQ=255)
                                                    gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgc                             >  1:2958825/1‑53 (MQ=255)
                                                    gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTa                  >  1:3808948/1‑64 (MQ=255)
                                                    gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTAc            >  1:706306/1‑70 (MQ=255)
                                                    gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACt           >  1:3178330/1‑71 (MQ=255)
                                                    gAAGATGGTTATTTCAATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgccgc                          >  1:3789666/1‑56 (MQ=255)
                                                         tGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGttt                               <  1:1205483/46‑1 (MQ=255)
                                                           gTTATGTCCATAACAAATAACAACTACGGCGGCaaaa                                      >  1:2290676/1‑37 (MQ=255)
                                                            ttATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTg                              <  1:899925/44‑1 (MQ=255)
                                                            ttATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCg   <  1:3247604/71‑1 (MQ=255)
                                                             tATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCAc                       <  1:1115049/50‑1 (MQ=255)
                                                              aTTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgcc                            >  1:3223459/1‑44 (MQ=255)
                                                              aTTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGc                    >  1:2276814/1‑52 (MQ=255)
                                                              aTTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTc    >  1:4893398/1‑68 (MQ=255)
                                                              aTTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCg   >  1:596489/1‑69 (MQ=255)
                                                              aTTTCCATAACAAATAACAACTACGGCGGCAAAAGGAg                                  <  1:4299916/38‑1 (MQ=255)
                                                               tttCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCgg  <  1:3121777/69‑1 (MQ=255)
                                                               tttCCATAACAAATAACAACTACGGCGGCAAAAGGa                                   <  1:4572152/36‑1 (MQ=255)
                                                                ttCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCg   >  1:391645/1‑67 (MQ=255)
                                                                ttCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCg   >  1:3962795/1‑67 (MQ=255)
                                                                  ccATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTAcc                >  1:4479290/1‑52 (MQ=255)
                                                                   |                                                                
TCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGG  >  W3110S.gb/3224749‑3224880

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: