Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,900,258 G→A 29.6% A928V (GCT→GTT)  recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,900,2580GA29.6% 101.4 / 51.6 71A928V (GCT→GTT) recBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (26/24);  new base A (8/13);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAGCT  >  minE/1900188‑1900328
                                                                      |                                                                      
cGGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcag                                                                        >  1:361916/1‑71 (MQ=255)
 ggTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcag                                                                        >  1:1288992/1‑70 (MQ=255)
   tGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCAt                                                                     >  1:1652640/1‑71 (MQ=255)
     aCGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATc                                                                    <  1:469138/70‑1 (MQ=255)
        cACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACa                                                                >  1:1345197/1‑71 (MQ=255)
         aCCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACAt                                                               >  1:2544009/1‑71 (MQ=255)
         aCCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACAt                                                               >  1:2922483/1‑71 (MQ=255)
           cgcgcgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATc                                                              >  1:3127706/1‑70 (MQ=255)
             cgcgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAg                                                           >  1:2956039/1‑71 (MQ=255)
               cgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGcc                                                         <  1:3276163/71‑1 (MQ=255)
               cgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGcc                                                         <  1:2318/71‑1 (MQ=255)
               cgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGcc                                                         <  1:372064/71‑1 (MQ=255)
               cgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGcc                                                         <  1:11998/71‑1 (MQ=255)
                ggAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCg                                                        <  1:2235723/71‑1 (MQ=255)
                 gAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGa                                                       <  1:1169303/71‑1 (MQ=255)
                  aaaCTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAg                                                      >  1:2488333/1‑71 (MQ=255)
                  aaaCTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAg                                                      >  1:3314224/1‑71 (MQ=255)
                    aCTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAgg                                                     >  1:976367/1‑70 (MQ=255)
                    aCTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAgg                                                     >  1:2366136/1‑70 (MQ=255)
                    aCTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAgg                                                     >  1:971639/1‑70 (MQ=255)
                    aCTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAgg                                                     >  1:1462274/1‑70 (MQ=255)
                    aCTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTAGCAACGCCTGCAGCATCGACATCCAGCCGAgg                                                     >  1:3119208/1‑70 (MQ=255)
                     cTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCa                                                   <  1:1685675/71‑1 (MQ=255)
                       gATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCAt                                                  >  1:174729/1‑70 (MQ=255)
                        aTGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATc                                                 >  1:3265652/1‑70 (MQ=255)
                        aTGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATc                                                 >  1:2171646/1‑70 (MQ=255)
                        aTGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGa                                                       >  1:1021553/1‑64 (MQ=255)
                         tgtgGTGTTAACGTCGTTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGa                                                       <  1:2159110/63‑1 (MQ=255)
                         tgtgGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGa                                                       <  1:594362/63‑1 (MQ=255)
                           tggtgTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCaaa                                              <  1:2691484/70‑1 (MQ=255)
                           tggtgTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCaaa                                              <  1:2265630/70‑1 (MQ=255)
                            ggtgTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATc                                            >  1:3908236/1‑71 (MQ=255)
                             gtgtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATcc                                           <  1:220629/71‑1 (MQ=255)
                              tgtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCt                                          <  1:2146882/71‑1 (MQ=255)
                               gtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGa                                                       >  1:203510/1‑57 (MQ=255)
                                ttAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcagca                                                                      <  1:2817580/41‑1 (MQ=255)
                                ttAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTg                                         <  1:3303861/70‑1 (MQ=255)
                                 tAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATc                                            <  1:3097454/66‑1 (MQ=255)
                                  aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTggg                                       >  1:2392536/1‑70 (MQ=255)
                                  aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTggg                                       <  1:782903/70‑1 (MQ=255)
                                     gTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCt                                          <  1:2418723/64‑1 (MQ=255)
                                          ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCg                               <  1:3121309/70‑1 (MQ=255)
                                          ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCg                               <  1:3101828/70‑1 (MQ=255)
                                          ttcttcAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCGAATCCTGGGCGATACCGt                              <  1:3556382/71‑1 (MQ=255)
                                          ttcttcAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt                              <  1:3279373/71‑1 (MQ=255)
                                          ttcttcAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt                              <  1:1715766/71‑1 (MQ=255)
                                           tcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt                              >  1:1637013/1‑70 (MQ=255)
                                           tcttcAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGc                                      <  1:3498359/62‑1 (MQ=255)
                                              tcAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGAc                           >  1:1603009/1‑70 (MQ=255)
                                               cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGa                                                       <  1:2467533/41‑1 (MQ=255)
                                               cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGAcc                          <  1:1350217/70‑1 (MQ=255)
                                               cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGAc                           <  1:3843017/69‑1 (MQ=255)
                                               cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCa                         >  1:1409/1‑71 (MQ=255)
                                               cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCa                         >  1:2632049/1‑71 (MQ=255)
                                               cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCa                         >  1:2059225/1‑71 (MQ=255)
                                               cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCa                         >  1:1608748/1‑71 (MQ=255)
                                                  cgacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGc                      <  1:2105720/71‑1 (MQ=255)
                                                  cgacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGc                      <  1:3850063/71‑1 (MQ=255)
                                                  cgacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGc                      >  1:903526/1‑71 (MQ=255)
                                                    acgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTg                    <  1:3023548/71‑1 (MQ=255)
                                                     cgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGTCGATACCGTGACCACGCTg                    <  1:496617/70‑1 (MQ=255)
                                                     cgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGt                   <  1:1421276/71‑1 (MQ=255)
                                                       cTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCt                                          <  1:2025463/46‑1 (MQ=255)
                                                       cTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCt                                          <  1:2623486/46‑1 (MQ=255)
                                                       cTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTg                 >  1:1185949/1‑71 (MQ=255)
                                                       cTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCaaa                                              <  1:772811/42‑1 (MQ=255)
                                                          gCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGtt                  >  1:42441/1‑67 (MQ=255)
                                                          gCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCa               <  1:1297523/70‑1 (MQ=255)
                                                            aaCGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAAc            >  1:1984538/1‑71 (MQ=255)
                                                             aCGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAAc            >  1:2516469/1‑70 (MQ=255)
                                                              cGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCaatcg           >  1:678450/1‑67 (MQ=255)
                                                               gCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAAc            >  1:1652488/1‑68 (MQ=255)
                                                                 cTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCagag       >  1:2957942/1‑71 (MQ=255)
                                                                 cTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGa                            <  1:699393/50‑1 (MQ=255)
                                                                   gcagcaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTg                 >  1:1404971/1‑59 (MQ=255)
                                                                   gcaacaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGAACACGCTGTTGCAAACCagag       >  1:1168939/1‑69 (MQ=255)
                                                                    cagcaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAAcc           >  1:3456724/1‑64 (MQ=255)
                                                                    cagcaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAg    <  1:2115534/71‑1 (MQ=255)
                                                                     agcaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAGc   >  1:1166110/1‑71 (MQ=255)
                                                                      gcaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAGCt  <  1:1081733/71‑1 (MQ=255)
                                                                      |                                                                      
CGGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAGCT  >  minE/1900188‑1900328

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: