Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,390,615 A→G 57.0% 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,6150AG57.0% 36.4 / 102.4 86Q249Q (CAA→CAGnfoendonuclease IV with intrinsic 3'‑5' exonuclease activity
Reads supporting (aligned to +/- strand):  ref base A (17/20);  new base G (26/23);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

GCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAG  >  minE/1390545‑1390681
                                                                      |                                                                  
gCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCaa                                                                    <  1:2931334/71‑1 (MQ=255)
   gcgTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAga                                                                  >  1:3066116/1‑70 (MQ=255)
     gTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAgacga                                                               >  1:3275348/1‑71 (MQ=255)
       tGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGAcc                                                             >  1:3768036/1‑71 (MQ=255)
         aCCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGt                                                           <  1:3028768/71‑1 (MQ=255)
           cGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAgacg                                                                >  1:2721176/1‑64 (MQ=255)
            gCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGgacga                                                               >  1:779883/1‑64 (MQ=255)
            gCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTc                                                        <  1:2183682/71‑1 (MQ=255)
             ccATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTc                                                        >  1:1901954/1‑70 (MQ=255)
             ccATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCg                                                       <  1:2366201/71‑1 (MQ=255)
             ccATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCg                                                       <  1:2197719/71‑1 (MQ=255)
             ccATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCg                                                       <  1:1491157/71‑1 (MQ=255)
              catcatAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGa                                                      >  1:2640918/1‑71 (MQ=255)
              catcatAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGa                                                      >  1:3892565/1‑71 (MQ=255)
               atcatATCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGAc                                                     >  1:1220679/1‑71 (MQ=255)
                tcatAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGgacg                                                                <  1:1155640/59‑1 (MQ=255)
                tcatAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAga                                                                  >  1:1804050/1‑57 (MQ=255)
                tcatAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGa                                                      >  1:2094272/1‑69 (MQ=255)
                 catAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACgg                                                   >  1:1407377/1‑71 (MQ=255)
                 catAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACgg                                                   >  1:1402780/1‑71 (MQ=255)
                 catAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGt                                                           <  1:3063296/63‑1 (MQ=255)
                  atAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGgac                                                                 <  1:2014106/56‑1 (MQ=255)
                  atAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACgg                                                   <  1:3532548/70‑1 (MQ=255)
                  atAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACgg                                                   <  1:3938568/70‑1 (MQ=255)
                   tAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCa                                                 <  1:2935888/71‑1 (MQ=255)
                   tAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCa                                                 <  1:2503202/71‑1 (MQ=255)
                   tAGCCTCGGTGAAGGCAATATCGGGCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCa                                                 <  1:546053/71‑1 (MQ=255)
                    aGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACgg                                                   >  1:1188023/1‑68 (MQ=255)
                     gCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCAtt                                               >  1:1620671/1‑71 (MQ=255)
                            tGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGa                                        <  1:85450/71‑1 (MQ=255)
                            tGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGa                                        <  1:3064890/71‑1 (MQ=255)
                             gAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTcc                                             <  1:3624652/65‑1 (MQ=255)
                              aaGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTg                                         <  1:3884937/68‑1 (MQ=255)
                              aaGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATc                                      <  1:130028/71‑1 (MQ=255)
                                 gCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTc                                   <  1:729966/71‑1 (MQ=255)
                                   aaTATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                 <  1:1674208/71‑1 (MQ=255)
                                   aaTATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                 <  1:2757892/71‑1 (MQ=255)
                                   aaTATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                 <  1:2213882/71‑1 (MQ=255)
                                    atatCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGAc                                                     <  1:1793285/50‑1 (MQ=255)
                                     tatCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGaaa                               >  1:3742567/1‑71 (MQ=255)
                                     tatCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGaaa                               >  1:44019/1‑71 (MQ=255)
                                        cGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGAt                                       >  1:363270/1‑60 (MQ=255)
                                        cGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGt                                                           <  1:1491358/40‑1 (MQ=255)
                                         ggTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGatcatc                                   >  1:3751276/1‑63 (MQ=255)
                                         ggTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTc                                   >  1:1109611/1‑63 (MQ=255)
                                         ggTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCa                            >  1:3495097/1‑70 (MQ=255)
                                         ggTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCAt                           >  1:1408344/1‑71 (MQ=255)
                                          gTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCa                            >  1:2847678/1‑69 (MQ=255)
                                           tCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCa                         <  1:1662084/71‑1 (MQ=255)
                                           tCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCa                         <  1:474266/71‑1 (MQ=255)
                                             atgatgCGTTCCGCTGGATCATGCAGGACGACCGTTTCg                                                       <  1:2462096/39‑1 (MQ=255)
                                             atgatgCGTTCCGCTGGATCATGCAGGACGACCGTTTCg                                                       <  1:541017/39‑1 (MQ=255)
                                             atgatgCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCa                            <  1:3504231/66‑1 (MQ=255)
                                             atgatgCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAAc                       >  1:3488214/1‑71 (MQ=255)
                                             atgatgCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAAc                       >  1:83233/1‑71 (MQ=255)
                                              tgatgCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAAcc                      <  1:2562674/71‑1 (MQ=255)
                                               gatgCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACg                                                    >  1:1163798/1‑40 (MQ=255)
                                                atgCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCAt                                                <  1:1173976/43‑1 (MQ=255)
                                                 tgCGTTCCGCTGGATCATGCAAGACGCCCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                    >  1:3125333/1‑70 (MQ=255)
                                                 tgCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                    >  1:1216162/1‑70 (MQ=255)
                                                 tgCGTTCCGCTGGATCATGCAAGACGACAGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                    >  1:2500734/1‑70 (MQ=255)
                                                     tgccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATc                          <  1:95186/58‑1 (MQ=255)
                                                      tCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATc              >  1:1733302/1‑71 (MQ=255)
                                                       ccGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt             >  1:1982138/1‑71 (MQ=255)
                                                       ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAAccc                     >  1:1706866/1‑63 (MQ=255)
                                                       ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAAccc                     >  1:3513282/1‑63 (MQ=255)
                                                        cGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt             <  1:1934487/70‑1 (MQ=255)
                                                        cGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt             <  1:1447086/70‑1 (MQ=255)
                                                        cGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACAATCAACCCGGATATCt             <  1:772766/70‑1 (MQ=255)
                                                        cGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTg            <  1:2465697/71‑1 (MQ=255)
                                                             gATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGc         <  1:3353529/69‑1 (MQ=255)
                                                             gATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                 <  1:1411147/45‑1 (MQ=255)
                                                             gATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                 <  1:813857/45‑1 (MQ=255)
                                                             gATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                 <  1:3946063/45‑1 (MQ=255)
                                                             gATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTggg          >  1:3046799/1‑68 (MQ=255)
                                                              aTCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt             >  1:2346471/1‑64 (MQ=255)
                                                               tCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATc                          >  1:3393154/1‑50 (MQ=255)
                                                               tCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATc                          >  1:2330021/1‑50 (MQ=255)
                                                               tCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaa     <  1:1927253/71‑1 (MQ=255)
                                                                 aTGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaa     >  1:3304392/1‑69 (MQ=255)
                                                                  tGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag  <  1:821385/71‑1 (MQ=255)
                                                                  tGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag  <  1:530208/71‑1 (MQ=255)
                                                                   gCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAAc                       >  1:3311206/1‑49 (MQ=255)
                                                                   gCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAAc                       >  1:1662436/1‑49 (MQ=255)
                                                                   gCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaa     >  1:1993363/1‑67 (MQ=255)
                                                                   gCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaa     >  1:3155964/1‑67 (MQ=255)
                                                                    cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAAccc                     >  1:899921/1‑50 (MQ=255)
                                                                      |                                                                  
GCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAG  >  minE/1390545‑1390681

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: