Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 183,899 G→A 45.3% L194L (CTG→CTA rnhB → ribonuclease HII, degrades RNA of DNA‑RNA hybrids

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE183,8990GA45.3% 18.4 / 82.8 64L194L (CTG→CTArnhBribonuclease HII, degrades RNA of DNA‑RNA hybrids
Reads supporting (aligned to +/- strand):  ref base G (19/16);  new base A (19/10);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

ACTGGCTGAACACGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTT  >  minE/183836‑183967
                                                               |                                                                    
aCTGGCTGAACACGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTg                                                               >  1:2108922/1‑71 (MQ=255)
aCTGGCTGAACACGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTg                                                               >  1:1971166/1‑71 (MQ=255)
 cTGGCTGAACACGGCGCGCCCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGc                                                              <  1:2367234/71‑1 (MQ=255)
    gCTGAACACGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTc                                                           >  1:4026787/1‑71 (MQ=255)
    gCTGAACACGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTAGGa                                                                   >  1:3905470/1‑63 (MQ=255)
        aaCACGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGa                                                       <  1:2716480/71‑1 (MQ=255)
        aaCACGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGa                                                       <  1:2523056/71‑1 (MQ=255)
        aaCACGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCt                                                         <  1:2697777/69‑1 (MQ=255)
        aaCACGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTAGGACTTGAGTcc                                                          >  1:3570250/1‑68 (MQ=255)
          cacGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGAtt                                                     >  1:1736530/1‑71 (MQ=255)
          cacGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGa                                                       >  1:3137124/1‑69 (MQ=255)
          cacGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGAtt                                                     >  1:410574/1‑71 (MQ=255)
          cacGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGAtt                                                     >  1:1114019/1‑71 (MQ=255)
                gcgACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTtgtg                                               >  1:2557607/1‑71 (MQ=255)
                    ccGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCg                                            <  1:1689227/70‑1 (MQ=255)
                     cGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGc                                                              <  1:1538382/51‑1 (MQ=255)
                      gAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGc                                                              >  1:3315717/1‑50 (MQ=255)
                        acacCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGAtt                                       <  1:1238201/71‑1 (MQ=255)
                          acCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTAGGACTTGCGTc                                                           <  1:173403/49‑1 (MQ=255)
                           ccATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAg                                    >  1:2847782/1‑71 (MQ=255)
                             aTCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTa                                  <  1:1331563/71‑1 (MQ=255)
                               cGGCGCAGCTTTGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCg                                            <  1:274066/59‑1 (MQ=255)
                                ggCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAAc                               >  1:452787/1‑71 (MQ=255)
                                ggCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACATGCGTCCTGATTCTtg                                                 >  1:3848578/1‑53 (MQ=255)
                                ggCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACATGCGTCCTGATTCTTGTGTCGAGATTAAGTAAAc                               >  1:1624471/1‑71 (MQ=255)
                                ggCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACATGCGTCCTGATTCTTGTGTCGAGATTAAGTAAAc                               >  1:439155/1‑71 (MQ=255)
                                ggCGCAGCTTTGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTTCTGTCGAGAt                                        <  1:4076781/62‑1 (MQ=255)
                                   gcAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCg                             >  1:3994332/1‑70 (MQ=255)
                                    cAGCTTTGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGa                           <  1:304314/71‑1 (MQ=255)
                                      gCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTc                                                           >  1:160583/1‑37 (MQ=255)
                                      gCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTc                                             >  1:442989/1‑51 (MQ=255)
                                      gCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAAt                         <  1:232336/71‑1 (MQ=255)
                                       cTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATc                        <  1:4013391/71‑1 (MQ=255)
                                       cTTTGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCt                                                         <  1:2745678/38‑1 (MQ=255)
                                         ttGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCt                       >  1:2597466/1‑70 (MQ=255)
                                         ttGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCt                       >  1:2558302/1‑70 (MQ=255)
                                         ttGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCt                       >  1:3282564/1‑70 (MQ=255)
                                         ttGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCt                       >  1:2996187/1‑70 (MQ=255)
                                          tGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTg                      >  1:763912/1‑70 (MQ=255)
                                          tGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTg                      >  1:2953477/1‑70 (MQ=255)
                                          tGGGCCTGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTg                      <  1:2547031/70‑1 (MQ=255)
                                            ggCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTaaa                                >  1:2893262/1‑58 (MQ=255)
                                            ggCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAg                   <  1:1014229/71‑1 (MQ=255)
                                               cTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCga                                           >  1:1512263/1‑44 (MQ=255)
                                                tGTCAAACGCGCACTGGGACTTGCGTCCTGATTCtt                                                  >  1:2207368/1‑36 (MQ=255)
                                                tGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCgg                            <  1:2748144/58‑1 (MQ=255)
                                                tGTCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTaa                                     <  1:1994631/49‑1 (MQ=255)
                                                  tCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCga                                           <  1:3602035/41‑1 (MQ=255)
                                                  tCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTaaa                                <  1:2155456/52‑1 (MQ=255)
                                                  tCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGTAATCTGAAGATg                <  1:325561/68‑1 (MQ=255)
                                                  tCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCt             >  1:2902220/1‑71 (MQ=255)
                                                  tCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCgaga                                         >  1:3110429/1‑43 (MQ=255)
                                                  tCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTa                                      <  1:3010122/46‑1 (MQ=255)
                                                  tCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTa                                  >  1:2925016/1‑50 (MQ=255)
                                                  tCAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTg                      <  1:2189479/62‑1 (MQ=255)
                                                   cAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTaaa                                >  1:2244278/1‑51 (MQ=255)
                                                   cAAACGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAAc                               >  1:3454213/1‑52 (MQ=255)
                                                    aaaCGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGAt                 >  1:3936647/1‑65 (MQ=255)
                                                    aaaCGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGAt                 >  1:1077376/1‑65 (MQ=255)
                                                    aaaCGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTaaa                                >  1:2528988/1‑50 (MQ=255)
                                                    aaaCGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCt             >  1:1269547/1‑69 (MQ=255)
                                                    aaaCGCGCACTAGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTg            >  1:2653398/1‑70 (MQ=255)
                                                      aCGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAAc         <  1:3514454/71‑1 (MQ=255)
                                                             cTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGttt  <  1:2063254/71‑1 (MQ=255)
                                                               |                                                                    
ACTGGCTGAACACGGCGCGACCGAACACCATCGGCGCAGCTTTGGGCCTGTCAAACGCGCACTGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTT  >  minE/183836‑183967

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: