Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,224,816 C→T 43.3% 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,8160CT43.3% 10.5 / 34.6 30H147Y (CAT→TAT) ebgCcryptic beta‑D‑galactosidase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (9/8);  new base T (8/5);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGAT  >  W3110S.gb/3224749‑3224884
                                                                   |                                                                    
tCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTat                                                                    <  1:22117/70‑1 (MQ=255)
     ttATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaa                                                              >  1:717292/1‑71 (MQ=255)
     ttATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaa                                                               >  1:643945/1‑70 (MQ=255)
       aTCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaat                                                             <  1:1363507/70‑1 (MQ=255)
        tCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaat                                                             <  1:1025480/69‑1 (MQ=255)
           gCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataa                                                                  >  1:928642/1‑61 (MQ=255)
           gCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataa                                                                  >  1:1216516/1‑61 (MQ=255)
                 aCGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTAcgg                                                  >  1:55054/1‑71 (MQ=255)
                       tCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaataac                                                          <  1:1274419/57‑1 (MQ=255)
                           aaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTa                                                     >  1:1755521/1‑58 (MQ=255)
                           aaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAgg                                        <  1:602569/71‑1 (MQ=255)
                            aaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCaaaa                                          >  1:1635200/1‑68 (MQ=255)
                            aaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAg                                         >  1:838618/1‑69 (MQ=255)
                            aaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGa                                       <  1:1594071/71‑1 (MQ=255)
                             aaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaa                                                              >  1:1485118/1‑47 (MQ=255)
                             aaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaataacaa                                                        >  1:1724762/1‑53 (MQ=255)
                              aGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGt                                     <  1:1026254/71‑1 (MQ=255)
                              aGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGt                                     <  1:740488/71‑1 (MQ=255)
                                  gTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTg                                  <  1:1189683/70‑1 (MQ=255)
                                    tctcAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTg                                  >  1:68211/1‑68 (MQ=255)
                                     ctcAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTAcgg                                                  <  1:1602598/51‑1 (MQ=255)
                                           gTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGc                        >  1:1065863/1‑71 (MQ=255)
                                               ccATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCa                            >  1:330783/1‑63 (MQ=255)
                                                  tCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgc                                 <  1:418750/55‑1 (MQ=255)
                                                   cGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTAc                >  1:216061/1‑71 (MQ=255)
                                                              aTTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCg       >  1:11821/1‑69 (MQ=255)
                                                              aTTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGttt                                   <  1:368158/41‑1 (MQ=255)
                                                                 tCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgcc                                <  1:1009267/41‑1 (MQ=255)
                                                                 tCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGAt  >  1:1286358/1‑71 (MQ=255)
                                                                 tCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGAt  >  1:965388/1‑71 (MQ=255)
                                                                   |                                                                    
TCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGAT  >  W3110S.gb/3224749‑3224884

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: