Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,390,615 A→G 7.6% Q249Q (CAA→CAG nfo → endonuclease IV with intrinsic 3'‑5' exonuclease activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,390,6150AG 7.6% 386.1 / 17.7 132Q249Q (CAA→CAGnfoendonuclease IV with intrinsic 3'‑5' exonuclease activity
Reads supporting (aligned to +/- strand):  ref base A (60/62);  new base G (2/8);  total (62/70)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 138 total aligned reads displayed.

GCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATC  >  minE/1390545‑1390684
                                                                      |                                                                     
gCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCaa                                                                       <  1:1097718/71‑1 (MQ=255)
 ccGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAg                                                                      <  1:2729056/71‑1 (MQ=255)
 ccGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAg                                                                      <  1:2090965/71‑1 (MQ=255)
 ccGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAg                                                                      <  1:3685389/71‑1 (MQ=255)
     gTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAgacga                                                                  <  1:2706011/71‑1 (MQ=255)
     gTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAgac                                                                    >  1:1360636/1‑69 (MQ=255)
        gACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCg                                                               >  1:3123553/1‑71 (MQ=255)
        gACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCg                                                               >  1:2959623/1‑71 (MQ=255)
         aCCGCCATCATAGCCTCGGTGAAGGCAATATCGTTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGt                                                              >  1:2928315/1‑71 (MQ=255)
         aCCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGt                                                              >  1:1582249/1‑71 (MQ=255)
         aCCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGt                                                              <  1:3113875/71‑1 (MQ=255)
           cGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGttt                                                            >  1:3448315/1‑71 (MQ=255)
           cGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGttt                                                            >  1:519994/1‑71 (MQ=255)
           cGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGtt                                                             >  1:3162511/1‑70 (MQ=255)
            gCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTc                                                           <  1:1253868/71‑1 (MQ=255)
            gCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTc                                                           <  1:672887/71‑1 (MQ=255)
             ccATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTc                                                           >  1:1486726/1‑70 (MQ=255)
             ccATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCg                                                          <  1:1176563/71‑1 (MQ=255)
             ccATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCg                                                          <  1:1640257/71‑1 (MQ=255)
             ccATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCg                                                          <  1:520104/71‑1 (MQ=255)
              catcatAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGa                                                         <  1:1814213/71‑1 (MQ=255)
               atcatAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGAc                                                        <  1:3187645/71‑1 (MQ=255)
               atcatAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGAc                                                        <  1:2572562/71‑1 (MQ=255)
                 catAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACgg                                                      >  1:2677675/1‑71 (MQ=255)
                   tAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCa                                                    <  1:2768596/71‑1 (MQ=255)
                     gcgtCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCAt                                                   <  1:403662/67‑1 (MQ=255)
                     gCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCAt                                                   <  1:2474625/70‑1 (MQ=255)
                     gCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCAt                                                   <  1:2700673/70‑1 (MQ=255)
                     gCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCAt                                                   <  1:3076411/70‑1 (MQ=255)
                     gCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCAt                                                   <  1:2016601/70‑1 (MQ=255)
                      ccTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGa                                                         <  1:1876822/63‑1 (MQ=255)
                      ccTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGa                                                         <  1:1217558/63‑1 (MQ=255)
                      ccTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTc                                                 >  1:2551911/1‑71 (MQ=255)
                       cTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCg                                                          <  1:2161016/61‑1 (MQ=255)
                        tCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGAc                                                        >  1:2113032/1‑62 (MQ=255)
                        tCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACg                                                       >  1:1602330/1‑63 (MQ=255)
                        tCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCg                                               >  1:389904/1‑71 (MQ=255)
                        tCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCg                                               >  1:297626/1‑71 (MQ=255)
                        tCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCg                                               >  1:605427/1‑71 (MQ=255)
                        tCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCg                                               >  1:856022/1‑71 (MQ=255)
                          ggTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCAtt                                                  >  1:2971645/1‑66 (MQ=255)
                          ggTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCg                                               <  1:801347/69‑1 (MQ=255)
                          ggTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCt                                             <  1:1493583/71‑1 (MQ=255)
                           gTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCt                                             <  1:1115071/70‑1 (MQ=255)
                           gTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTg                                            <  1:3409718/71‑1 (MQ=255)
                           gTGAAGGAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCg                                                          <  1:2284072/56‑1 (MQ=255)
                            tGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGa                                           <  1:1215368/71‑1 (MQ=255)
                             gAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCg                                                          <  1:3547037/55‑1 (MQ=255)
                               aGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATcc                                        <  1:2660079/71‑1 (MQ=255)
                               aGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATca                                        <  1:293695/71‑2 (MQ=255)
                                ggCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCt                                       <  1:142821/71‑1 (MQ=255)
                                ggCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGtt                                                             >  1:1613684/1‑49 (MQ=255)
                                ggCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCa                                                    >  1:1629714/1‑58 (MQ=255)
                                 gCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTc                                      <  1:703883/71‑1 (MQ=255)
                                  cAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTc                                      >  1:2016046/1‑70 (MQ=255)
                                   aaTATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                    <  1:1017209/71‑1 (MQ=255)
                                     tatCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCAt                                                   <  1:906043/54‑1 (MQ=255)
                                     tatCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGt                                                              <  1:3703413/43‑1 (MQ=255)
                                      atCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAAc                                 >  1:2419125/1‑71 (MQ=255)
                                      atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTc                                      >  1:2534242/1‑66 (MQ=255)
                                       tCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAAcc                                <  1:2673935/71‑1 (MQ=255)
                                        cGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAAcc                                <  1:1336291/70‑1 (MQ=255)
                                         ggTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCAt                              <  1:1770976/71‑1 (MQ=255)
                                          gTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGa                                                         >  1:1161223/1‑43 (MQ=255)
                                          gTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAAc                                 >  1:1374109/1‑67 (MQ=255)
                                           tCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGa                                           >  1:1746374/1‑56 (MQ=255)
                                           tCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCa                            <  1:2525183/71‑1 (MQ=255)
                                            cATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCg                                               >  1:3210654/1‑51 (MQ=255)
                                                atgCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                       <  1:2934682/71‑1 (MQ=255)
                                                 tgCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCgg                      <  1:397683/71‑1 (MQ=255)
                                                 tgCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCgg                      <  1:937381/71‑1 (MQ=255)
                                                   cGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGa                     >  1:1410493/1‑70 (MQ=255)
                                                     ttCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCgg                      <  1:767221/67‑1 (MQ=255)
                                                     ttCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGatat                  <  1:715274/71‑1 (MQ=255)
                                                       ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGc                                              >  1:85562/1‑41 (MQ=255)
                                                       ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGc                                              >  1:35051/1‑41 (MQ=255)
                                                       ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGc                                              >  1:775101/1‑41 (MQ=255)
                                                       ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGc                                              >  1:1523469/1‑41 (MQ=255)
                                                       ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGc                                              >  1:1307764/1‑41 (MQ=255)
                                                       ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGa                     >  1:2424908/1‑66 (MQ=255)
                                                       ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATc                 <  1:2777701/70‑1 (MQ=255)
                                                        cGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt                <  1:3243179/70‑1 (MQ=255)
                                                        cGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt                <  1:772573/70‑1 (MQ=255)
                                                        cGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt                <  1:1733489/70‑1 (MQ=255)
                                                        cGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt                <  1:839473/70‑1 (MQ=255)
                                                        cGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTg               >  1:1444903/1‑71 (MQ=255)
                                                        cGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTg               >  1:258848/1‑71 (MQ=255)
                                                         gCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCa                            <  1:2498410/57‑1 (MQ=255)
                                                            ggATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTggg             >  1:602820/1‑69 (MQ=255)
                                                             gATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTggg             >  1:1676911/1‑68 (MQ=255)
                                                               tCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaa        <  1:3368714/71‑1 (MQ=255)
                                                               tCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCag          >  1:1860340/1‑69 (MQ=255)
                                                               tCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCag          >  1:3043975/1‑69 (MQ=255)
                                                               tCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTCGGCag          >  1:800313/1‑69 (MQ=255)
                                                                   gCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATc                             >  1:1383369/1‑46 (MQ=255)
                                                                   gCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaa        >  1:1047062/1‑67 (MQ=255)
                                                                    cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAgaga    <  1:52554/70‑1 (MQ=255)
                                                                    cAAGACGACCGTTTCGACGGCATTC‑‑CTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATc  <  1:1048384/70‑1 (MQ=255)
                                                                      aGACGACCGTTTCGACGGCATTCCGCTGATCCTCGaa                                   <  1:165957/37‑1 (MQ=255)
                                                                      aGACGACCGTTTCGACGGCATTCCGCTGATCCTCGaa                                   <  1:1111864/37‑1 (MQ=255)
                                                                      |                                                                     
GCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATC  >  minE/1390545‑1390684

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: