Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,222,206 C→T 8.8% N306N (AAC→AAT ebgA → cryptic beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,222,2060CT 8.8% 208.7 / 8.6 68N306N (AAC→AATebgAcryptic beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (34/28);  new base T (4/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATC  >  W3110S.gb/3222138‑3222275
                                                                    |                                                                     
cagcaATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGt                                                                     <  1:953953/71‑1 (MQ=255)
cagcaATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGt                                                                     <  1:2680331/71‑1 (MQ=255)
cagcaATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGt                                                                     <  1:2052336/71‑1 (MQ=255)
 agcaATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGt                                                                     <  1:1494918/70‑1 (MQ=255)
    aaTGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACCCTGAAAGACGCCAACGGCAACGTTCTg                                                                 >  1:4080999/1‑71 (MQ=255)
     aTGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTgg                                                                <  1:2234930/71‑1 (MQ=255)
     aTGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGACAACGTTCTgg                                                                <  1:2876455/71‑1 (MQ=255)
        gTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCttg                                                                >  1:206929/1‑66 (MQ=255)
        gTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTgg                                                                >  1:1311610/1‑68 (MQ=255)
        gTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTgg                                                                >  1:232290/1‑68 (MQ=255)
        gTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAg                                                             >  1:1180702/1‑71 (MQ=255)
              agAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGtt                                                                    >  1:936254/1‑58 (MQ=255)
                aaTCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAgt                                                            <  1:1420708/64‑1 (MQ=255)
                 aTCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCcac                                                    <  1:2637517/71‑1 (MQ=255)
                 aTCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGTGCcac                                                    >  1:1190634/1‑70 (MQ=255)
                  tCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCcac                                                    >  1:2858473/1‑70 (MQ=255)
                   ccccTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGaa                                                              >  1:2593799/1‑59 (MQ=255)
                    cccTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGtt                                                                    >  1:3124564/1‑52 (MQ=255)
                    cccTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGaa                                                              >  1:1758822/1‑58 (MQ=255)
                    cccTTATCTTTACCATCAGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGaa                                                              >  1:3474959/1‑58 (MQ=255)
                       ttATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAAcgcg                                              >  1:1569892/1‑71 (MQ=255)
                        tATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCATCGTTCTGGAATTGGTGCCACAAcgc                                               <  1:144542/69‑1 (MQ=255)
                        tATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGt                                             >  1:3454655/1‑71 (MQ=255)
                        tATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGt                                             >  1:289194/1‑71 (MQ=255)
                         aTCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGtt                                            >  1:1753755/1‑71 (MQ=255)
                         aTCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGaa                                                              >  1:1749313/1‑53 (MQ=255)
                         aTCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAgtgg                                                          >  1:2325323/1‑57 (MQ=255)
                         aTCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGtt                                            >  1:1342195/1‑71 (MQ=255)
                         aTCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGtt                                            >  1:2573641/1‑71 (MQ=255)
                          tCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTg                                           >  1:810656/1‑71 (MQ=255)
                           cTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCcac                                                    <  1:139027/61‑1 (MQ=255)
                            tttACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGc                                         >  1:2721121/1‑71 (MQ=255)
                            tttACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGc                                         >  1:1178421/1‑71 (MQ=255)
                            tttACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGCAAGTGGTGCCACAACGCGTTGGc                                         >  1:652854/1‑71 (MQ=255)
                                  aTCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTTGTGCCACAACGCGTTGGCTTCCGt                                   <  1:744023/71‑1 (MQ=255)
                                  aTCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGt                                   <  1:1810736/71‑1 (MQ=255)
                                  aTCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGt                                   <  1:927378/71‑1 (MQ=255)
                                  aTCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGt                                   <  1:2264431/71‑1 (MQ=255)
                                   tCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGt                                   >  1:2258338/1‑70 (MQ=255)
                                       gTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGa                                 <  1:3305433/68‑1 (MQ=255)
                                       gTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACTCGTTGGCTTCCGTg                                  <  1:3385898/67‑1 (MQ=255)
                                        tCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAgt                                                            <  1:1817340/40‑1 (MQ=255)
                                        tCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAgt                                                            <  1:763325/40‑1 (MQ=255)
                                         cATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATc                             >  1:1274095/1‑70 (MQ=255)
                                          aTGACGCTGAAAGACGCTAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCa                            <  1:2549169/70‑1 (MQ=255)
                                          aTGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCcac                                                    <  1:2763314/46‑1 (MQ=255)
                                          aTGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCcac                                                    <  1:942999/46‑1 (MQ=255)
                                          aTGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCcac                                                    <  1:938134/46‑1 (MQ=255)
                                           tGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCa                            <  1:4060533/69‑1 (MQ=255)
                                            gACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAg                         >  1:4082907/1‑71 (MQ=255)
                                            gACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAg                         >  1:3156714/1‑71 (MQ=255)
                                            gACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAg                         >  1:50206/1‑71 (MQ=255)
                                             aCGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACaa                                                  >  1:431877/1‑45 (MQ=255)
                                             aCGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCaaa                          >  1:3683716/1‑69 (MQ=255)
                                             aCGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCaaa                          >  1:3508892/1‑69 (MQ=255)
                                             aCGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCaaa                          >  1:2089465/1‑69 (MQ=255)
                                                    aaGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGAc                 >  1:3509061/1‑71 (MQ=255)
                                                    aaGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGAc                 >  1:3197849/1‑71 (MQ=255)
                                                          ccAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGa                                 <  1:3269500/49‑1 (MQ=255)
                                                          ccAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCa                            <  1:3283618/54‑1 (MQ=255)
                                                          ccAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTg           <  1:3575399/71‑1 (MQ=255)
                                                          ccAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTg           <  1:1833628/71‑1 (MQ=255)
                                                          ccAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTg           <  1:2385274/71‑1 (MQ=255)
                                                            aaCGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGt          >  1:3061718/1‑70 (MQ=255)
                                                             aCGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTc        <  1:181728/71‑1 (MQ=255)
                                                               ggCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGt              >  1:1211826/1‑63 (MQ=255)
                                                                  aaCGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGtctggtctgg     >  1:229196/1‑69 (MQ=255)
                                                                   aCGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATc  <  1:3181671/71‑1 (MQ=255)
                                                                    |                                                                     
CAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATC  >  W3110S.gb/3222138‑3222275

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: