Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,073,966 C→A 5.8% H207Q (CAC→CAA ebgA → cryptic beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,073,9660CA 5.8% 325.8 / 3.1 100H207Q (CAC→CAAebgAcryptic beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (54/40);  new base A (3/3);  total (57/43)
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.99e-01

Only 100 of 106 total aligned reads displayed.

ACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGCCACGCTTTC  >  Exported/2073899‑2074034
                                                                   |                                                                    
acGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCt                                                                    >  1:829457/1‑70 (MQ=255)
 cGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCt                                                                    <  1:3948230/69‑1 (MQ=255)
 cGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTaa                                                                  >  1:1198604/1‑71 (MQ=255)
 cGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTaa                                                                  >  1:3439484/1‑71 (MQ=255)
     gAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCa                                                              >  1:5678395/1‑71 (MQ=255)
        gACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCata                                                            >  1:3255080/1‑70 (MQ=255)
          ccAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATAtt                                                          >  1:1936165/1‑70 (MQ=255)
           cAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCata                                                            >  1:2154221/1‑67 (MQ=255)
            aGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAAc                                                       <  1:2294261/71‑1 (MQ=255)
             ggATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAAc                                                       >  1:1389205/1‑70 (MQ=255)
             ggATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAAc                                                       >  1:1723047/1‑70 (MQ=255)
               atatGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCata                                                            >  1:4513104/1‑63 (MQ=255)
                tatGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGAtt                                                   >  1:2886445/1‑71 (MQ=255)
                  tgtgGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACg                                                      <  1:2453607/66‑1 (MQ=255)
                    tggtggTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCa                                                <  1:77616/70‑1 (MQ=255)
                    tggtggTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCa                                                <  1:1981342/70‑1 (MQ=255)
                      gtggtCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGAttt                                                  <  1:3278962/66‑1 (MQ=255)
                       tggtCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACtg                                             >  1:2525863/1‑70 (MQ=255)
                         gtCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAacac                                                                      >  1:5575770/1‑43 (MQ=255)
                         gtCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCt                                                                    <  1:1539173/45‑1 (MQ=255)
                         gtCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGc                                          >  1:395267/1‑71 (MQ=255)
                            aGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTa                                       >  1:292464/1‑71 (MQ=255)
                              cGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTAcc                                     <  1:3681086/71‑1 (MQ=255)
                                gggATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAAc                                                       <  1:168354/51‑1 (MQ=255)
                                 ggATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGAc                                  <  1:3132722/71‑1 (MQ=255)
                                 ggATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGAc                                  >  1:1713250/1‑71 (MQ=255)
                                  gATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGa                                                     >  1:5213767/1‑51 (MQ=255)
                                   aTCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACt                                 <  1:1029236/70‑1 (MQ=255)
                                      ttCCGCGATGTTTATCTGGTCGGAAAACACCTAACGGATATTAACGATTTCa                                                >  1:4331588/1‑52 (MQ=255)
                                      ttCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCa                                                >  1:3461273/1‑52 (MQ=255)
                                       tCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGa                             >  1:4651357/1‑70 (MQ=255)
                                        ccGCGATTTTTATCTGGTCGGAAAACACCTAACGCata                                                            >  1:4828960/1‑38 (MQ=255)
                                        ccGCGATGTTTATCTGGTCGGAAAACACCTAACGCata                                                            >  1:1606797/1‑38 (MQ=255)
                                        ccGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAAc                                                       >  1:1130270/1‑43 (MQ=255)
                                        ccGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGAttt                                                  >  1:4138034/1‑48 (MQ=255)
                                        ccGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGAttt                                                  >  1:653203/1‑48 (MQ=255)
                                        ccGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCa                                                >  1:634237/1‑50 (MQ=255)
                                        ccGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACttt                               <  1:4284029/67‑1 (MQ=255)
                                        ccGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACttt                               <  1:3201223/67‑1 (MQ=255)
                                        ccGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGAc                            >  1:4151343/1‑70 (MQ=255)
                                        ccGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGAc                            <  1:5428645/70‑1 (MQ=255)
                                        ccGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGAc                            <  1:4709083/70‑1 (MQ=255)
                                        ccGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGAc                            >  1:5041309/1‑70 (MQ=255)
                                         cgcgATGTTTATCTGGTCGGAAAACACCTAACGCATATTAAc                                                       >  1:305445/1‑42 (MQ=255)
                                         cgcgATGTTTATCTGGTCGGAAAACACCTAACGCATATTAAc                                                       >  1:5466928/1‑42 (MQ=255)
                                         cgcgATGTTTATCTGGTCGGAAAACACCTAACGCATATTAAc                                                       >  1:2320361/1‑42 (MQ=255)
                                         cgcgATGTTTATCTGGTCGGAAAACACCTAACGCATATTAAc                                                       >  1:815047/1‑42 (MQ=255)
                                            gATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGAc                                  >  1:1264038/1‑60 (MQ=255)
                                            gATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGAc                                  >  1:2589824/1‑60 (MQ=255)
                                            gATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGAc                                  >  1:2421286/1‑60 (MQ=255)
                                            gATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGa                             >  1:1201593/1‑65 (MQ=255)
                                            gATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGa                             >  1:3684648/1‑65 (MQ=255)
                                            gATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGa                          >  1:3926465/1‑68 (MQ=255)
                                              tGTTTATCTTGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCt                     <  1:5625361/71‑1 (MQ=255)
                                              tGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTAcc                                     <  1:2279403/55‑1 (MQ=255)
                                              tGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACttt                               <  1:3397401/61‑1 (MQ=255)
                                              tGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCt                     <  1:4443063/71‑1 (MQ=255)
                                                tttATCTGGTCGGAAAACACCTAACGCATATTAAc                                                       <  1:1722337/35‑1 (MQ=255)
                                                tttATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACtgtg                                           <  1:67031/47‑1 (MQ=255)
                                                tttATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTAt                   <  1:5114093/71‑1 (MQ=255)
                                                tttATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTAt                   <  1:3762761/71‑1 (MQ=255)
                                                  tATCTGGTCGGAAAACACCTAACGCATATTAACGAttt                                                  <  1:2925203/38‑1 (MQ=255)
                                                  tATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTAtt                  >  1:542087/1‑70 (MQ=255)
                                                  tATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTg                 >  1:4225194/1‑71 (MQ=255)
                                                  tATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTg                 >  1:4588310/1‑71 (MQ=255)
                                                    tCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGc                <  1:1888743/70‑1 (MQ=255)
                                                    tCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGc                <  1:1759830/70‑1 (MQ=255)
                                                    tCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGc                <  1:3266544/70‑1 (MQ=255)
                                                    tCTGATCGGAAAACAACTCACGCACATTAACGATTTCACCGTGCGCACCGACTTTGACGAAGCCTATTGc                <  1:1957621/70‑1 (MQ=255)
                                                    tCTGATCGGAAAACAACTCACGCACATTAACGATTTCACCGTGCGCACCGACTTTGACGAAGCCTATTGc                <  1:2095253/70‑1 (MQ=255)
                                                    tCTGATCGGAAAACAACTCACGCACATTAACGATTTCACCGTGCGCACCGACTTTGACGAAGCCTATTGc                <  1:300068/70‑1 (MQ=255)
                                                      tGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTa                                       >  1:1232471/1‑45 (MQ=255)
                                                       ggTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTAttcga               >  1:2641969/1‑65 (MQ=255)
                                                       ggTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGa              >  1:1863789/1‑69 (MQ=255)
                                                       ggTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGa              >  1:2740612/1‑69 (MQ=255)
                                                        gTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAg                        >  1:5222327/1‑58 (MQ=255)
                                                        gTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGcc                      <  1:865727/60‑1 (MQ=255)
                                                        gTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGc                <  1:955188/66‑1 (MQ=255)
                                                        gTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCg               <  1:4859575/67‑1 (MQ=255)
                                                        gTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCg               <  1:1340072/67‑1 (MQ=255)
                                                        gTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCg               <  1:784480/67‑1 (MQ=255)
                                                          cGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGAc                                  <  1:1718101/46‑1 (MQ=255)
                                                          cGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGCCa         >  1:4692536/1‑71 (MQ=255)
                                                           ggAAAACACCTAACGCATATTACCGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGcc          <  1:3561103/69‑1 (MQ=255)
                                                           ggAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGcc                      <  1:4515636/57‑1 (MQ=255)
                                                           ggAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGcc                      <  1:4470045/57‑1 (MQ=255)
                                                           ggAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGcc                      <  1:2875200/57‑1 (MQ=255)
                                                           ggAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGCCAc        >  1:3647121/1‑71 (MQ=255)
                                                             aaaaCACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCt                     >  1:4038652/1‑56 (MQ=255)
                                                             aaaaCACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTa                    >  1:3234138/1‑57 (MQ=255)
                                                             aaaaCACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGCCACGc      >  1:2644696/1‑71 (MQ=255)
                                                             aaaaCAACTCACGCACATTAACGATTTCACCGTGCGCACCGa                                   >  1:5116277/1‑42 (MQ=19)
                                                             aaaaCAACTCACGCACATTAACGATTTCACCGTGCGCACCGa                                   >  1:3234705/1‑42 (MQ=19)
                                                             aaaaCAACTCACGCACATTAACGATTTCACCGTGCGCACCGa                                   >  1:5014876/1‑42 (MQ=19)
                                                              aaaCACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTa                    >  1:3913278/1‑56 (MQ=255)
                                                                acacCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGa                             >  1:4229548/1‑45 (MQ=255)
                                                                acacCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGCCACGCttt   >  1:1790712/1‑71 (MQ=255)
                                                                 cacCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGCCACGCTTTc  <  1:4456771/71‑1 (MQ=255)
                                                                  acCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGa                             >  1:3779181/1‑43 (MQ=255)
                                                                  acCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGa                             >  1:5379298/1‑43 (MQ=255)
                                                                   |                                                                    
ACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGCCACGCTTTC  >  Exported/2073899‑2074034

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