Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 520,325 T→C 50.6% V115V (GTT→GTC uvrB → excinulease of nucleotide excision repair, DNA damage recognition component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE520,3250TC50.6% 2.5 / 126.6 89V115V (GTT→GTCuvrBexcinulease of nucleotide excision repair, DNA damage recognition component
Reads supporting (aligned to +/- strand):  ref base T (27/17);  new base C (20/25);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGACTACTATCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGTGATGTGGTTGTG  >  minE/520257‑520394
                                                                    |                                                                     
tacGACTACTATCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTaa                                                                     >  1:3456401/1‑71 (MQ=255)
tacGACTACTATCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCaa                                                                     >  1:3255151/1‑71 (MQ=255)
 acGACTACTATCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTa                                                                      >  1:868403/1‑69 (MQ=255)
      tactaTCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAAc                                                                    >  1:1122390/1‑66 (MQ=255)
       actaTCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACa                                                               >  1:2418585/1‑70 (MQ=255)
        ctaTCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACat                                                              >  1:1282879/1‑70 (MQ=255)
        ctaTCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACata                                                             <  1:3216365/71‑1 (MQ=255)
          aTCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACata                                                             >  1:1476435/1‑69 (MQ=255)
           tCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTa                                                                      >  1:683768/1‑59 (MQ=255)
           tCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATAtt                                                           >  1:462256/1‑70 (MQ=255)
           tCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATAtt                                                           >  1:1527341/1‑70 (MQ=255)
           tCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTg                                                          <  1:1457773/71‑1 (MQ=255)
             aGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATagag                                                        >  1:3618655/1‑71 (MQ=255)
             aGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAg                                                        >  1:3460274/1‑71 (MQ=255)
               cgggAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCa                                                      <  1:2333944/69‑1 (MQ=255)
                cGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACatat                                                            <  1:2257339/64‑1 (MQ=255)
                cGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAg                                                     <  1:2597791/71‑1 (MQ=255)
                  gAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGa                                                    >  1:2376170/1‑70 (MQ=255)
                  gAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGa                                                    <  1:2859909/70‑1 (MQ=255)
                   aaGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGa                                                         >  1:3497590/1‑64 (MQ=255)
                   aaGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGc                                                       >  1:2813666/1‑66 (MQ=255)
                    aGCCTATGTACCGAGTTCCGACACTTTCATTAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGc                                                 <  1:2621926/70‑1 (MQ=255)
                     gCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCg                                                <  1:2735278/71‑1 (MQ=255)
                      ccTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCGt                                               <  1:4440091/71‑1 (MQ=255)
                      ccTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCGt                                               <  1:122000/71‑1 (MQ=255)
                        tATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCa                                                      <  1:2148670/62‑1 (MQ=255)
                         aTGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTg                                            >  1:2031824/1‑71 (MQ=255)
                          tGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGt                                           <  1:3834743/71‑1 (MQ=255)
                          tGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGt                                           >  1:1608026/1‑71 (MQ=255)
                           gTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTc                                          >  1:69232/1‑71 (MQ=255)
                            tACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTcc                                         >  1:3619106/1‑71 (MQ=255)
                            tACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTcc                                         >  1:2372667/1‑71 (MQ=255)
                              ccGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACat                                                              <  1:4619543/48‑1 (MQ=255)
                              ccGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACat                                                              <  1:969612/48‑1 (MQ=255)
                              ccGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACa                                                               <  1:3508780/47‑1 (MQ=255)
                              ccGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACa                                                               <  1:2996136/47‑1 (MQ=255)
                              ccGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGa                                                    >  1:657533/1‑58 (MQ=255)
                                gAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGcca                                     >  1:3480106/1‑71 (MQ=255)
                                 aGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGccac                                    >  1:574630/1‑71 (MQ=255)
                                 aGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGccac                                    >  1:236429/1‑71 (MQ=255)
                                 aGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGa                                                    >  1:918049/1‑55 (MQ=255)
                                 aGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTAGCAGATGCGTTTGTCCGccac                                    >  1:4643209/1‑70 (MQ=255)
                                  gTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGccac                                    >  1:2151622/1‑70 (MQ=255)
                                   ttCCGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACat                                                              <  1:2709511/43‑1 (MQ=255)
                                     ccGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGa                                                    >  1:859296/1‑51 (MQ=255)
                                     ccGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGAt                                                   <  1:3898652/52‑1 (MQ=255)
                                     ccGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCaa                                 <  1:2735664/70‑1 (MQ=255)
                                     ccGACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGa                                                    >  1:1612785/1‑51 (MQ=255)
                                      cGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCg                                                <  1:817439/54‑1 (MQ=255)
                                       gACACTTTCATTGAGAAAGATGCCTCTGTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGc                              <  1:2155513/71‑1 (MQ=255)
                                       gACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGc                              <  1:1385800/71‑1 (MQ=255)
                                       gACACTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCa                                                      <  1:2040536/47‑1 (MQ=255)
                                        acacTTTCATTGAGAAAGATGCCTCGGTTAACGAAc                                                                >  1:2294214/1‑36 (MQ=255)
                                        acacTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGtt                                              >  1:3372001/1‑54 (MQ=255)
                                        acacTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAg                               >  1:2906310/1‑69 (MQ=255)
                                        acacTTTCATTGAGAAAGATGCCTCGGTCAACGAACatat                                                            >  1:2743867/1‑40 (MQ=255)
                                          acTTTCATTGAGAAAGATGCCTCGGTTAACGAACata                                                             >  1:586615/1‑37 (MQ=255)
                                          acTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGc                                                 <  1:4079841/49‑1 (MQ=255)
                                          acTTTCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGc                              <  1:1711754/68‑1 (MQ=255)
                                            tttCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGc                         >  1:3157624/1‑71 (MQ=255)
                                             ttCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGccac                                    >  1:4481501/1‑59 (MQ=255)
                                             ttCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCt                        >  1:1092070/1‑71 (MQ=255)
                                              tCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTg                                            >  1:4627428/1‑50 (MQ=255)
                                              tCATTGAGAAAGATGCCTCGGTTAAAGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGc                         >  1:4420312/1‑69 (MQ=255)
                                              tCATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGc                         <  1:3201597/69‑1 (MQ=255)
                                               cATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTg                       <  1:1390866/70‑1 (MQ=255)
                                               cATTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTg                       <  1:2438030/70‑1 (MQ=255)
                                                aTTGAGAAAGATGCCTCGGTCAACGAACATATTGAGc                                                       >  1:170037/1‑37 (MQ=255)
                                                aTTGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTgg                      <  1:3330365/70‑1 (MQ=255)
                                                 ttGAGAAAGATGCCTCGGTCAACGAACATATTGAGCa                                                      <  1:2330563/37‑1 (MQ=255)
                                                 ttGAGAAAGATGCCTCGGTCAACGAACATATTGAGCa                                                      <  1:2975209/37‑1 (MQ=255)
                                                 ttGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTc                                          >  1:3837382/1‑49 (MQ=255)
                                                 ttGAGAAAGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGc                              <  1:1052846/61‑1 (MQ=255)
                                                     gaAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAgcggc                <  1:3319544/71‑1 (MQ=255)
                                                       aaGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGt              <  1:1653067/71‑1 (MQ=255)
                                                       aaGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGt              <  1:855430/71‑1 (MQ=255)
                                                       aaGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGt              <  1:2743123/71‑1 (MQ=255)
                                                        aGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGt              >  1:3388513/1‑70 (MQ=255)
                                                        aGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGTg             <  1:2411406/71‑1 (MQ=255)
                                                        aGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGt              >  1:1186025/1‑70 (MQ=255)
                                                        aGATGCCTCGGTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGTg             <  1:4418793/71‑1 (MQ=255)
                                                          aTGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGTGAt           >  1:3113287/1‑71 (MQ=255)
                                                           tGCCTCGGTCAACGAACATATTGAGCAGATGCGttt                                             >  1:2075655/1‑36 (MQ=255)
                                                               tCGGTCAACGAACATATTGAGCAGATGCGTTTGTCAGc                                       <  1:3372596/38‑1 (MQ=38)
                                                                 ggTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGTGAtgtggttg    <  1:379797/71‑1 (MQ=255)
                                                                 ggTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGTGAtgtggt      <  1:3282156/69‑1 (MQ=255)
                                                                 ggTCAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGTGAtgtggttg    >  1:4375988/1‑71 (MQ=255)
                                                                  gTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGTGAtgtggtt     >  1:1273239/1‑69 (MQ=255)
                                                                  gTTAACGAACATATTGAGCAGATGCGTTTGACCGCCACCAAAGCGATGCTGGAGCGGCGTGAtgtggttgt   <  1:3573091/71‑1 (MQ=255)
                                                                    tAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGTGAtgtggttgtg  <  1:2984061/70‑1 (MQ=255)
                                                                    |                                                                     
TACGACTACTATCAGCCGGAAGCCTATGTACCGAGTTCCGACACTTTCATTGAGAAAGATGCCTCGGTTAACGAACATATTGAGCAGATGCGTTTGTCCGCCACCAAAGCGATGCTGGAGCGGCGTGATGTGGTTGTG  >  minE/520257‑520394

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: