Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,936,265 G→A 5.2% V370V (GTC→GTT recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,936,2650GA 5.2% 304.3 / 4.2 95V370V (GTC→GTTrecJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base G (41/49);  new base A (2/3);  total (43/52)
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.28e-01

TCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCA  >  minE/1936197‑1936335
                                                                    |                                                                      
tCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGa                                                                       >  1:2384058/1‑70 (MQ=255)
tCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGAc                                                                      <  1:84228/71‑1 (MQ=255)
  aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACaa                                                                    >  1:1560042/1‑71 (MQ=255)
  aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACaa                                                                    >  1:2802002/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:2258066/69‑1 (MQ=255)
          gcgcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCt                                                              <  1:4677479/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)
                  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:3440101/71‑1 (MQ=255)
                           ggACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATa                                           <  1:1823114/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:2707206/47‑1 (MQ=255)
                                gTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCt                                                              <  1:3243164/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:210671/1‑71 (MQ=255)
                                       cGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGccc                               >  1:1973284/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:3509719/1‑71 (MQ=255)
                                            tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCgg                          >  1:721625/1‑71 (MQ=255)
                                            tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCgg                          >  1:3693228/1‑71 (MQ=255)
                                            tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCgg                          >  1:2601104/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)
                                                                  ccGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGc           >  1:1502947/1‑64 (MQ=255)
                                                                   cGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTcc   <  1:784016/71‑1 (MQ=255)
                                                                    gACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCa  <  1:3414912/71‑1 (MQ=255)
                                                                    |                                                                      
TCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCA  >  minE/1936197‑1936335

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: