Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,936,247 G→A 58.7% R376R (CGC→CGT recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,936,2470GA58.7% 25.5 / 53.6 75R376R (CGC→CGTrecJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base G (27/4);  new base A (15/29);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 4.81e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCC  >  minE/1936181‑1936305
                                                                  |                                                          
aCCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATgc                                                           >  1:215038/1‑68 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:1058974/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:925532/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:784069/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:3379977/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:2561897/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:2434772/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:2313532/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:1937056/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:1930865/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:1885089/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:1746966/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:1010619/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:1111433/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:1614508/1‑66 (MQ=255)
      cAGCGTACCGTCGCCCGCAGGCGCAAAGGCGATAACCGGACGGTGAAAACGCTCTTTGATACGCGat                                                      >  1:1521854/1‑66 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:45203/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:2974404/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:3090993/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:1206938/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:2380234/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:2336510/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:110201/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:2077410/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:2064816/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:2017489/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:1080077/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:1378882/1‑65 (MQ=255)
             ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:906538/1‑65 (MQ=255)
                  aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:674059/1‑60 (MQ=255)
                  aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:3076974/1‑60 (MQ=255)
                  aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:979411/1‑60 (MQ=255)
                  aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:362802/1‑60 (MQ=255)
                  aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:3128344/1‑60 (MQ=255)
                  aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:2346232/1‑60 (MQ=255)
                  aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:3322314/1‑60 (MQ=255)
                  aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                                 >  1:3212865/1‑60 (MQ=255)
                            gcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCt                                <  1:1693717/67‑1 (MQ=255)
                                              cGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATAc            <  1:892725/69‑1 (MQ=255)
                                                 tGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACAtt         <  1:33155/69‑1 (MQ=255)
                                                 tGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACAtt         <  1:46611/69‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:3226067/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:3394312/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:2757499/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:466495/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:2810742/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:1854574/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:1827457/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:1809288/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:2101853/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:1263529/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:2220417/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:1771586/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:2405406/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:2579207/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:2734859/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:1710882/65‑1 (MQ=255)
                                                   aaaaCGCTCTTTGATACGCGATGCAAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:574772/65‑1 (MQ=255)
                                                    aaaCGCTCTTTGATACGCGATGCCAGAATGCCAACACTGCCCTGATGCCATTCAGGGTGATACa           <  1:2666724/64‑1 (MQ=255)
                                                    aaaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:3337728/64‑1 (MQ=255)
                                                     aaCGCTCTTTGATACGCGATGCCAGAATGCCAACAATGCCCTGATGCCATTCAGGGTGATACa           <  1:2526623/63‑1 (MQ=255)
                                                     aaCGCTCTTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGCCATTCAGGGTGATACATTGCCa     <  1:266200/69‑1 (MQ=255)
                                                     aaCGCTCTTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGCCATTCAGGGTGATACATTGCCa     <  1:2250229/69‑1 (MQ=255)
                                                     aaCGCTCTTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGCCATTCAGGGTGATACATTGCCa     <  1:3201266/69‑1 (MQ=255)
                                                     aaCGCTCTTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGCCATTCAGGGTGATACATTGCCa     <  1:461993/69‑1 (MQ=255)
                                                     aaCGCTCTTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGCCATTCAGGGTGATACATTGCCa     <  1:2756445/69‑1 (MQ=255)
                                                     aaCGCTCTTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGCCATTCAGGGTGATACATTGCCa     <  1:2058694/69‑1 (MQ=255)
                                                      aCGCTCTTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGCCATTCAGGGTGATACATTGCCa     <  1:3206759/68‑1 (MQ=255)
                                                      aCGCTCTTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGCCATTCAGGGTGATACATTGCCa     <  1:2925541/68‑1 (MQ=255)
                                                       cgcgCTTTGATACGCGATGCCAGAATGCCAACAACGCCCTGATGCCATTCAGGGTGATACATTGCCa     <  1:1248924/67‑1 (MQ=255)
                                                        gCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGcc  >  1:982115/1‑69 (MQ=255)
                                                        gCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGcc  >  1:1986778/1‑69 (MQ=255)
                                                        gCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGcc  >  1:666719/1‑69 (MQ=255)
                                                        gCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGcc  >  1:206632/1‑69 (MQ=255)
                                                        gCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGcc  >  1:2851254/1‑69 (MQ=255)
                                                        gCTATTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGcc  >  1:344011/1‑69 (MQ=255)
                                                                  |                                                          
ACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCC  >  minE/1936181‑1936305

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: