Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,936,262 A→G 5.2% G371G (GGT→GGC recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,936,2620AG 5.2% 264.4 / 5.0 96G371G (GGT→GGCrecJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base A (43/48);  new base G (2/3);  total (45/51)
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.90e-01

CCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTG  >  minE/1936194‑1936328
                                                                    |                                                                  
ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATAcc                                                                  >  1:2235271/1‑71 (MQ=255)
ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATAcc                                                                  >  1:4594059/1‑71 (MQ=255)
ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATAcc                                                                  >  1:3489499/1‑71 (MQ=255)
   tCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGa                                                                >  1:2384058/1‑70 (MQ=255)
   tCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGAc                                                               <  1:84228/71‑1 (MQ=255)
     aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACaa                                                             >  1:2802002/1‑71 (MQ=255)
     aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACaa                                                             >  1:1560042/1‑71 (MQ=255)
     aCCTGCTGGCGAAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACaa                                                             >  1:3696156/1‑71 (MQ=255)
      ccTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAAc                                                            <  1:1249402/71‑1 (MQ=255)
        tgctgGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACaa                                                             >  1:777631/1‑68 (MQ=255)
        tgctgGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGc                                                          >  1:72816/1‑71 (MQ=255)
         gctgGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGcc                                                         >  1:3404605/1‑71 (MQ=255)
          ctgGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGcc                                                         <  1:3760568/70‑1 (MQ=255)
             gcgcAAACGCGATAACCGGACGGTGAAAACGCTTTTTGATGCGCGAAGCCAGAATACCGACAACGCCCt                                                       <  1:3030602/69‑1 (MQ=255)
             gcgcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCt                                                       <  1:1348590/69‑1 (MQ=255)
             gcgcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCt                                                       <  1:4677479/69‑1 (MQ=255)
             gcgcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCt                                                       <  1:2258066/69‑1 (MQ=255)
               gcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATg                                                   >  1:469383/1‑71 (MQ=255)
               gcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATg                                                   <  1:1862926/71‑1 (MQ=255)
                cAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGc                                                  <  1:1012487/71‑1 (MQ=255)
                cAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGc                                                  <  1:2173753/71‑1 (MQ=255)
                cAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGc                                                  >  1:1566342/1‑71 (MQ=255)
                 aaaGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGcc                                                 <  1:1006588/71‑1 (MQ=255)
                 aaaGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGcc                                                 <  1:1607817/71‑1 (MQ=255)
                 aaaCGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGc                                                  >  1:3413475/1‑70 (MQ=255)
                  aaCGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGc                                                  <  1:4397949/69‑1 (MQ=255)
                   aCGCGATAACCGGACGGTGAAACCGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCAt                                               >  1:1970958/1‑71 (MQ=255)
                   aCGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCa                                                >  1:3316508/1‑70 (MQ=255)
                    cgcgATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCAtt                                              <  1:1777652/71‑1 (MQ=255)
                     gcgATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATAcc                                                                  <  1:2028664/50‑1 (MQ=255)
                     gcgATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTc                                             >  1:4390214/1‑71 (MQ=255)
                     gcgATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTc                                             <  1:413474/71‑1 (MQ=255)
                        aTAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCggg                                          <  1:4064388/71‑1 (MQ=255)
                         tAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCgggg                                         <  1:2518901/71‑1 (MQ=255)
                           aCCTGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGt                                        <  1:1713670/70‑1 (MQ=255)
                            ccGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGa                                      <  1:2347892/71‑1 (MQ=255)
                              ggACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATa                                    <  1:1823114/71‑1 (MQ=255)
                              ggACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATa                                    <  1:3440101/71‑1 (MQ=255)
                               gACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATAc                                   <  1:3181319/71‑1 (MQ=255)
                                aCGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGa                                                     >  1:2615641/1‑52 (MQ=255)
                                  ggTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACAtt                                <  1:4301011/71‑1 (MQ=255)
                                   gTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCt                                                       <  1:3243164/47‑1 (MQ=255)
                                   gTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCt                                                       <  1:2707206/47‑1 (MQ=255)
                                    tGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGa                                      <  1:3790536/63‑1 (MQ=255)
                                    tGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGAt                                     <  1:729779/64‑1 (MQ=255)
                                    tGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGc                              <  1:1654039/71‑1 (MQ=255)
                                      aaaaCGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCg                                            <  1:4006942/55‑1 (MQ=255)
                                      aaaaCGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACAtt                                >  1:1114265/1‑67 (MQ=255)
                                      aaaaCGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGcc                             <  1:1583279/70‑1 (MQ=255)
                                      aaaaCGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGc                              >  1:2776513/1‑69 (MQ=255)
                                       aaaCGCTCTTTGATGCGCGAAGCCAGAATACCGAc                                                               <  1:3228804/35‑1 (MQ=255)
                                        aaCGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGc                                                  <  1:1157516/47‑1 (MQ=255)
                                          cGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGa                                      >  1:2667700/1‑57 (MQ=255)
                                          cGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGccc                        >  1:1973284/1‑71 (MQ=255)
                                          cGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGccc                        >  1:210671/1‑71 (MQ=255)
                                           gCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACAtt                                <  1:3781513/62‑1 (MQ=255)
                                           gCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccg                       >  1:1335089/1‑71 (MQ=255)
                                            ctctTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgc                      >  1:492313/1‑71 (MQ=255)
                                            ctctTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGCCATTCAGGGTGATACATTGCCAGcccac                      <  1:2228018/71‑3 (MQ=255)
                                             tctTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACa                                  >  1:1211429/1‑58 (MQ=255)
                                             tctTTGATGCGCGAAGCCAGAATACCGACAACGCACTGATGCCATTCGGGGTGATACa                                  >  1:2543836/1‑58 (MQ=255)
                                               tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCAtt                                              >  1:368876/1‑44 (MQ=255)
                                               tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCAt                                               <  1:1035930/43‑1 (MQ=255)
                                               tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGa                                      >  1:4469350/1‑52 (MQ=255)
                                               tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                               <  1:4190407/59‑1 (MQ=255)
                                               tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCgg                   >  1:2601104/1‑71 (MQ=255)
                                               tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCgg                   >  1:721625/1‑71 (MQ=255)
                                               tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCgg                   >  1:3693228/1‑71 (MQ=255)
                                               tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCgg                   >  1:3509719/1‑71 (MQ=255)
                                               tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCAAGCCCGCCgg                   >  1:2943506/1‑71 (MQ=255)
                                               tttGAAGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCgg                   >  1:2784141/1‑71 (MQ=255)
                                                 tGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATa                                    >  1:1554590/1‑52 (MQ=255)
                                                 tGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACa                                  >  1:1945234/1‑54 (MQ=255)
                                                 tGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACa                                  >  1:3041472/1‑54 (MQ=255)
                                                 tGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACa                                  >  1:4047699/1‑54 (MQ=255)
                                                 tGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgccg                    <  1:210214/68‑1 (MQ=255)
                                                 tGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCgg                   <  1:2067016/69‑1 (MQ=255)
                                                 tGATGCGCGAAGCCAGAATACCGAAAACGCCCTGATGCCATTCGGGGTGATACATTGCCCGc                          >  1:2547244/1‑62 (MQ=255)
                                                   aTGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAg               <  1:1054943/71‑1 (MQ=255)
                                                      cgcgAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACa                                  <  1:4402823/49‑1 (MQ=255)
                                                      cgcgAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCg             <  1:169101/70‑1 (MQ=255)
                                                        cgAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGt                                        <  1:4341002/41‑1 (MQ=255)
                                                        cgAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCgtgt          <  1:2067757/71‑1 (MQ=255)
                                                        cgAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCgtgt          >  1:1758279/1‑71 (MQ=255)
                                                         gAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTCCAGCCCGCCGGGTAGCGTGTCa        <  1:2544932/71‑1 (MQ=255)
                                                          aaGCCAGAATACCGACAACGCCCTGATGCCATTCGTGGTGATACa                                  >  1:2706926/1‑45 (MQ=255)
                                                          aaGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCa        >  1:2417378/1‑71 (MQ=255)
                                                           aGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCAc       <  1:1913254/71‑1 (MQ=255)
                                                             ccAGAATACCGACAACGCCCTGATGCCATTCgggg                                         <  1:439656/35‑1 (MQ=255)
                                                             ccAGAATACCGACAACGCCCTGATGCCATTCGGGGt                                        <  1:725565/36‑1 (MQ=255)
                                                             ccAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACAtt                                <  1:3037881/44‑1 (MQ=255)
                                                                gAATGCCAACACCGCCCTGATGCCATTCAGGGTGATACATTGCCAGccca                       >  1:2093469/1‑49 (MQ=255)
                                                                gAATGCCAACAACGCCCTGATGCCATTCAGGGTGATACATTGCCAGCCCACCTGGCAGCGTGTCGCGGCta  >  1:2472693/1‑70 (MQ=255)
                                                                  aTACCGACAACGCCCTGATGCCATTCGGGGTGATa                                    >  1:2403989/1‑35 (MQ=255)
                                                                   tACCGACAACGCCCTGATGCCATTCGGGGTGATACAt                                 <  1:1428526/37‑1 (MQ=255)
                                                                   tACCGACAACGCCCTGATGCCATTCGGGGTGATACAt                                 <  1:2123838/37‑1 (MQ=255)
                                                                    aCCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCAc       <  1:1981876/62‑1 (MQ=255)
                                                                    |                                                                  
CCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTG  >  minE/1936194‑1936328

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: