Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,390,615 A→G 18.7% 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,6150AG18.7% 175.6 / 36.2 80Q249Q (CAA→CAGnfoendonuclease IV with intrinsic 3'‑5' exonuclease activity
Reads supporting (aligned to +/- strand):  ref base A (30/35);  new base G (6/9);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCG  >  minE/1390545‑1390685
                                                                      |                                                                      
gCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCaa                                                                        <  1:645476/71‑1 (MQ=255)
   gcgTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGga                                                                      >  1:1117587/1‑70 (MQ=255)
   gcgTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAgac                                                                     <  1:1553782/71‑1 (MQ=255)
         aCCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGt                                                               <  1:3757260/71‑1 (MQ=255)
         aCCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGt                                                               <  1:2085752/71‑1 (MQ=255)
            gCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTc                                                            <  1:3238087/71‑1 (MQ=255)
            gCCATCATAGCCTCGGTGAAGGCAATATCGGGCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTc                                                            <  1:2590592/71‑1 (MQ=255)
             ccATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCg                                                           <  1:3588733/71‑1 (MQ=255)
                tcatAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCaa                                                                        <  1:497728/55‑1 (MQ=255)
                  atAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGga                                                                      >  1:3107503/1‑55 (MQ=255)
                  atAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAgacga                                                                   >  1:1060787/1‑58 (MQ=255)
                  atAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACgg                                                       >  1:381405/1‑70 (MQ=255)
                  atAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACgg                                                       >  1:2307481/1‑70 (MQ=255)
                  atAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACgg                                                       <  1:2069635/70‑1 (MQ=255)
                   tAGCCTCGGTGAAGGCAATATCGGTCATGATGCTTTCCGCTGGATCATGCAAgacg                                                                    <  1:3958388/56‑1 (MQ=255)
                    aGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGa                                                          >  1:3254524/1‑65 (MQ=255)
                      ccTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTc                                                            <  1:1955480/61‑1 (MQ=255)
                      ccTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCAt                                                    <  1:3794228/69‑1 (MQ=255)
                       cTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTc                                                            <  1:1684288/60‑1 (MQ=255)
                          ggTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCa                                                     >  1:1175644/1‑64 (MQ=255)
                          ggTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCt                                              >  1:1866889/1‑71 (MQ=255)
                            tGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGAcc                                                                 <  1:1665315/50‑1 (MQ=255)
                             gAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGAt                                           >  1:2521716/1‑71 (MQ=255)
                              aaGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAgacgac                                                                  >  1:1200213/1‑47 (MQ=255)
                              aaGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAgacgac                                                                  >  1:2113851/1‑47 (MQ=255)
                                ggCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCt                                        >  1:3417377/1‑71 (MQ=255)
                                 gCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGt                                                               <  1:2624687/47‑1 (MQ=255)
                                 gCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTc                                       >  1:3767952/1‑71 (MQ=255)
                                   aaTATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACgg                                                       <  1:4139565/53‑1 (MQ=255)
                                   aaTATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACgg                                                       <  1:3738454/53‑1 (MQ=255)
                                   aaTATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCg                                      >  1:105877/1‑70 (MQ=255)
                                   aaTATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                     <  1:1249164/71‑1 (MQ=255)
                                    atatCGGTCATGATGCGTTCCGCTGGATCATGCAAg                                                                       <  1:221349/36‑1 (MQ=255)
                                    atatCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTc                                                            <  1:1834054/47‑1 (MQ=255)
                                    atatCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                     >  1:3902796/1‑70 (MQ=255)
                                     tatCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGaaa                                   <  1:1929967/71‑1 (MQ=255)
                                     tatCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGaaa                                   >  1:3897926/1‑71 (MQ=255)
                                     tatCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGaaa                                   >  1:2245248/1‑71 (MQ=255)
                                     tatCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGaaa                                   >  1:1473654/1‑71 (MQ=255)
                                       tCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGa                                            >  1:2918319/1‑60 (MQ=255)
                                       tCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAAc                                  >  1:3285774/1‑70 (MQ=255)
                                       tCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAAc                                  <  1:2399049/70‑1 (MQ=255)
                                       tCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAAc                                  >  1:3841776/1‑70 (MQ=255)
                                         ggTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCa                                <  1:3296979/70‑1 (MQ=255)
                                           tCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTc                                                  <  1:3813151/50‑1 (MQ=255)
                                            cATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCaa                            >  1:3878099/1‑71 (MQ=255)
                                            cATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCaa                            >  1:2104545/1‑71 (MQ=255)
                                             atgatgCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAAc                           >  1:771757/1‑71 (MQ=255)
                                                atgCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                        <  1:4139049/71‑1 (MQ=255)
                                                atgCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                        >  1:227514/1‑71 (MQ=255)
                                                atgCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                        <  1:948327/71‑1 (MQ=255)
                                                 tgCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCgg                       <  1:3169352/71‑1 (MQ=255)
                                                 tgCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                        >  1:2553965/1‑70 (MQ=255)
                                                   cGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGat                     <  1:1715520/71‑1 (MQ=255)
                                                    gTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCt                                              >  1:294473/1‑45 (MQ=255)
                                                    gTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGata                    <  1:977069/71‑1 (MQ=255)
                                                     ttCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                        <  1:2145915/66‑1 (MQ=255)
                                                     ttCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGata                    >  1:449879/1‑70 (MQ=255)
                                                       ccGCTGGATCATGCAGGACGACCTTTTCGACTGCATTc                                                  <  1:1858001/38‑1 (MQ=37)
                                                       ccGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt                 >  1:3629683/1‑71 (MQ=255)
                                                       ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt                 >  1:182310/1‑71 (MQ=255)
                                                        cGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt                 <  1:3455922/70‑1 (MQ=255)
                                                          cTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAAc                                  >  1:442998/1‑51 (MQ=255)
                                                             gATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCag           <  1:527727/71‑1 (MQ=255)
                                                               tCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaa         <  1:1588061/71‑1 (MQ=255)
                                                               tCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaa         <  1:2247148/71‑1 (MQ=255)
                                                                 aTGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGatat                   >  1:2720381/1‑59 (MQ=255)
                                                                 aTGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGatat                   >  1:2291608/1‑59 (MQ=255)
                                                                 aTGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAAcc                          >  1:1088036/1‑52 (MQ=255)
                                                                 aTGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCt                 <  1:1752999/61‑1 (MQ=255)
                                                                 aTGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaaga       <  1:32295/71‑1 (MQ=255)
                                                                 aTGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaa         <  1:3429077/69‑1 (MQ=255)
                                                                  tGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag      <  1:3177270/71‑1 (MQ=255)
                                                                  tGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag      <  1:686299/71‑1 (MQ=255)
                                                                  tGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATc                  <  1:194004/59‑1 (MQ=255)
                                                                  tGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag      <  1:2360038/71‑1 (MQ=255)
                                                                  tGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag      <  1:4143538/71‑1 (MQ=255)
                                                                  tGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag      <  1:3461783/71‑1 (MQ=255)
                                                                   gCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaaga       <  1:2674417/69‑1 (MQ=255)
                                                                   gCAAGACGACCGTTTCGCCGGCATTCCGCTGATCCTCGAAACCAt                               >  1:3815758/1‑45 (MQ=255)
                                                                    cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCagaa         >  1:2821025/1‑66 (MQ=255)
                                                                    cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAgaga     <  1:1905503/70‑1 (MQ=255)
                                                                      aGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCg  >  1:3088039/1‑71 (MQ=255)
                                                                      aGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCg  >  1:16968/1‑71 (MQ=255)
                                                                      |                                                                      
GCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCG  >  minE/1390545‑1390685

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: