Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,222,944 4 bp→AGCA coding (2671‑2674/3093 nt) b3076 → evolved beta‑D‑galactosidase alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,222,9440CA99.2% 394.7 / ‑5.6 134Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (69/64);  total (70/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01
*NC_0009133,222,9450AG98.5% 349.5 / NA 134Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (70/62);  minor base T (0/2);  total (70/64)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,222,9460AC100.0% 362.0 / NA 134Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (70/64);  total (70/64)
*NC_0009133,222,9470GA99.2% 388.0 / NA 134A892T (GCC→ACC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (69/62);  minor base C (1/0);  total (70/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

GTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATG  >  NC_000913/3222882‑3223007
                                                              ||||                                                            
gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCCCCGTCGATGc                                                    >  2:1619947/1‑76 (MQ=255)
   gTGGACGGGGCAAAAACTCCCCCGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                 <  1:3465927/76‑1 (MQ=255)
   gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                 >  1:4436035/1‑76 (MQ=255)
   gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                 >  2:4417261/1‑76 (MQ=255)
   gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                 >  1:1507253/1‑76 (MQ=255)
   gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                 >  1:3113975/1‑76 (MQ=255)
   gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                 >  2:2773185/1‑76 (MQ=255)
   gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                 >  1:294914/1‑76 (MQ=255)
   gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGACGCCAt                                                 >  2:1198558/1‑76 (MQ=255)
     ggACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                               >  2:3183492/1‑76 (MQ=255)
      gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                              >  2:3403470/1‑76 (MQ=255)
      gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGACGCCATGtt                                              >  1:1327898/1‑76 (MQ=255)
       ccggggCAAAACCTACCCGGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTc                                             <  1:1428207/73‑1 (MQ=255)
       aCCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTc                                             >  1:4650992/1‑76 (MQ=255)
       aCCGGCCAAAAACACCCCCAACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTc                                             <  1:2816711/76‑1 (MQ=255)
        ccGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCg                                            >  1:4623622/1‑76 (MQ=255)
        ccGGGCAAAAACCACCCCGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCg                                            <  2:2595893/76‑1 (MQ=255)
         ggggCGAAAGCTCGGCCGCCAGCCAGCAGGCTAACATCATCGATAGCTGGCGCAGCACCGTCGATGCCATGTTCga                                           <  2:1921943/75‑1 (MQ=255)
         ggggCGAAAACTCCGCCGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                           <  2:2382181/75‑1 (MQ=255)
         cGGGCGAAAACTACGCCGACAGCCAGCAGGCTCACATCATCGATCTCTGGCGCCGCACCGTCGATGCCAGGTTCga                                           >  2:3080431/1‑76 (MQ=255)
         cGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                           >  1:2555878/1‑76 (MQ=255)
          gggCGAAAACTCCGCCGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                          <  1:4327821/76‑1 (MQ=255)
          gggCGAAAACTACGCGGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                          <  2:1879499/76‑1 (MQ=255)
          gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                          >  1:4196236/1‑76 (MQ=255)
          gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                          >  2:836880/1‑76 (MQ=255)
          gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                          >  1:3496793/1‑76 (MQ=255)
          gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                          >  1:3509808/1‑76 (MQ=255)
            gaaaaacCTCCGCCGCCAGCCAGCAGGCTATCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                        <  1:3686370/73‑1 (MQ=255)
            gCGAAAACTCCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                        <  1:3672421/76‑1 (MQ=255)
             cGAAACCTCCGCGGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAAc                                       <  2:3000520/76‑1 (MQ=255)
             cGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAAc                                       >  1:2174170/1‑76 (MQ=255)
              gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACt                                      >  2:3493314/1‑76 (MQ=255)
              gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACt                                      >  2:2123443/1‑76 (MQ=255)
              gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACt                                      >  1:4340809/1‑76 (MQ=255)
              gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACt                                      >  1:3062219/1‑76 (MQ=255)
              aaaacCTACCCCGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACt                                      <  1:4700293/75‑1 (MQ=255)
              aaaaaCTCCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACt                                      <  1:3860673/75‑1 (MQ=255)
               aaaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                     >  1:4496256/1‑76 (MQ=255)
                aaccTACCCCGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                    <  2:1204681/73‑1 (MQ=255)
                   cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                 >  2:2908066/1‑76 (MQ=255)
                    tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                >  1:417454/1‑76 (MQ=255)
                    tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACCACCCCt                                >  2:4191302/1‑76 (MQ=255)
                       gCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCCTGTTCGAGAACTATCCCTTcc                             >  2:4739118/1‑76 (MQ=255)
                       gCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTcc                             >  2:3187992/1‑76 (MQ=255)
                        ccGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTccc                            <  1:2563899/76‑1 (MQ=255)
                        ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTccc                            >  1:1068697/1‑76 (MQ=255)
                        ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTccc                            >  1:2028301/1‑76 (MQ=255)
                         cgccagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCg                           <  2:2502409/73‑1 (MQ=255)
                          gACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGATCTATCCCTTCCCGc                          <  2:4602544/76‑1 (MQ=255)
                          cccagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                          <  2:432629/74‑1 (MQ=255)
                          cccagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                          <  1:4145277/74‑1 (MQ=255)
                          cccagccagcAGGCTAACATCATCGATATCGGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                          <  2:4767806/74‑1 (MQ=255)
                           ccagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCa                         <  2:4009030/75‑1 (MQ=255)
                           aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCa                         <  1:3969629/76‑1 (MQ=255)
                            cagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAg                        >  2:4008786/1‑76 (MQ=255)
                            cagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAg                        <  1:2503358/76‑1 (MQ=255)
                             agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGa                       >  1:2815435/1‑76 (MQ=255)
                              gccagcGGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaa                      <  2:1260600/76‑1 (MQ=255)
                              gccagcAGGCTAACATCATCGATATCCGGCGCAGCACCGTCGATGACATGCTCGAGAACTAGTCCTTCCCGCAGaa                      >  1:1944898/1‑76 (MQ=255)
                                cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                    >  2:1823580/1‑76 (MQ=255)
                                cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                    >  2:1123510/1‑76 (MQ=255)
                                cagcagGCTAACATCATCGATATCGGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                    <  2:3945940/76‑1 (MQ=255)
                                cagcagGCTAACAGCATCGATATATGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                    <  2:3804180/76‑1 (MQ=255)
                                 agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                   >  1:2977199/1‑76 (MQ=255)
                                  gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                  <  2:4661081/76‑1 (MQ=255)
                                   cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                 >  2:1384153/1‑76 (MQ=255)
                                   cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                 >  1:2948285/1‑76 (MQ=255)
                                   cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                 <  1:1210195/76‑1 (MQ=255)
                                   cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                 >  2:2895008/1‑76 (MQ=255)
                                    agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                >  1:3906972/1‑76 (MQ=255)
                                    agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                >  1:855544/1‑76 (MQ=255)
                                    agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                <  2:1485294/76‑1 (MQ=255)
                                    agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                <  2:2353513/76‑1 (MQ=255)
                                     ggCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt               >  2:3430023/1‑76 (MQ=255)
                                      gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa              >  2:3202992/1‑76 (MQ=255)
                                      gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa              >  1:2779227/1‑76 (MQ=255)
                                      gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa              <  2:2962540/76‑1 (MQ=255)
                                       cgaaCATCAGCGATATCTGGCGCATCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa             <  2:3759333/74‑1 (MQ=255)
                                       cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTac             >  2:518836/1‑75 (MQ=255)
                                        tgaCATCAGCGCTATCTGGCGCAGCACCGTCGCTGCCATGTTCGAGAACGATCCCGTCCCGCAGAACAACGGTAAc            <  2:4168573/74‑1 (MQ=255)
                                        tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc            <  1:1602022/76‑1 (MQ=255)
                                        tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc            <  1:3250796/76‑1 (MQ=255)
                                        tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc            <  2:2589553/76‑1 (MQ=255)
                                        tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc            >  1:2397622/1‑76 (MQ=255)
                                        tAACATCATCGATATCTGGCCCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGCACAACGGTAAc            >  2:2446220/1‑76 (MQ=255)
                                         aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc           >  1:1856328/1‑76 (MQ=255)
                                            atcatcGCTATCTGGCGCATCGCCGTCGATGCCATGTTCGAGAACTACCCCTCCCCGCATAACACCGGTAACCGTc        <  1:2516153/76‑1 (MQ=255)
                                             tcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa       <  1:1144659/76‑1 (MQ=255)
                                             tcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa       <  1:4557349/76‑1 (MQ=255)
                                             tcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa       <  1:277723/76‑1 (MQ=255)
                                             tcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa       <  1:2066380/76‑1 (MQ=255)
                                              catcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAg      <  1:4121065/76‑1 (MQ=255)
                                               atcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGc     >  1:1949872/1‑76 (MQ=255)
                                                tcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCa    <  1:1250703/76‑1 (MQ=255)
                                                 cGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCag   >  2:881220/1‑75 (MQ=255)
                                                 cGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAt   >  2:3881276/1‑76 (MQ=255)
                                                 cGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAt   <  1:3703091/76‑1 (MQ=255)
                                                 cGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAt   >  1:134140/1‑76 (MQ=255)
                                                 cGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGCCCGCac   >  2:223696/1‑75 (MQ=255)
                                                  gATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGACCAACGGTAACCGTCAGCATg  <  2:2153737/76‑1 (MQ=255)
                                                              ||||                                                            
GTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATG  >  NC_000913/3222882‑3223007

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