Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,236,063 G→A 7.3% V403V (GTG→GTA purH → fused IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,236,0630GA 7.3% 246.3 / 7.3 82V403V (GTG→GTApurHfused IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base G (43/33);  new base A (4/2);  total (47/35)
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.53e-01

ACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCGAAGTTT  >  minE/2235994‑2236132
                                                                     |                                                                     
aCGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCgtgg                                                                      >  1:3357476/1‑71 (MQ=255)
 cgcgTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCgtggt                                                                     >  1:50957/1‑71 (MQ=255)
 cgcgTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCgtggt                                                                     <  1:741185/71‑1 (MQ=255)
  gcgTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCgtggtg                                                                    <  1:3752138/71‑1 (MQ=255)
    gTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGAc                                                                  >  1:2821905/1‑71 (MQ=255)
    gTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGAc                                                                  >  1:87156/1‑71 (MQ=255)
    gTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGAc                                                                  >  1:550322/1‑71 (MQ=255)
      gAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGAcc                                                                 >  1:2887376/1‑70 (MQ=255)
      gAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGAcc                                                                 >  1:216843/1‑70 (MQ=255)
      gAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGAcc                                                                 >  1:212180/1‑70 (MQ=255)
       aaCGGCGGTCTGCTGGTTCAGGATCGTGACCTTGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCaa                                                               <  1:3385548/71‑1 (MQ=255)
       aaCGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCaa                                                               <  1:2693917/71‑1 (MQ=255)
       aaCGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCa                                                                <  1:2040977/70‑1 (MQ=255)
       aaCGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCa                                                                <  1:2967828/70‑1 (MQ=255)
         cggcggTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAAc                                                             >  1:1105850/1‑71 (MQ=255)
         cggcggTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTAGTGACCaaa                                                              >  1:522959/1‑70 (MQ=255)
           gcggTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGt                                                           <  1:3872289/71‑1 (MQ=255)
            cggTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCgtgg                                                                      <  1:1269835/59‑1 (MQ=255)
             ggTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCa                                                                <  1:2504313/64‑1 (MQ=255)
             ggTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCa                                                                <  1:1228681/64‑1 (MQ=255)
              gTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAg                                                        >  1:1945837/1‑71 (MQ=255)
               tCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGc                                                       >  1:1193593/1‑71 (MQ=255)
                ctgctgGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCgtggt                                                                     >  1:1858228/1‑56 (MQ=255)
                 tgctgGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGcc                                                      >  1:3549680/1‑70 (MQ=255)
                 tgctgGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGcc                                                      >  1:3127356/1‑70 (MQ=255)
                   ctgGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGa                                                                   <  1:3789495/55‑1 (MQ=255)
                   ctgGTTCAGGATCGTGACCTGGGGATGGTCGGTGCAGAAGAACTGCGCGTGGTGa                                                                   <  1:3448039/55‑1 (MQ=255)
                      gTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCg                                                 <  1:1937503/70‑1 (MQ=255)
                      gTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCg                                                 >  1:897621/1‑70 (MQ=255)
                         cAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCgtggt                                                                     >  1:1528728/1‑47 (MQ=255)
                         cAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGaa                                               >  1:1789760/1‑69 (MQ=255)
                         cAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGaa                                               >  1:3092543/1‑69 (MQ=255)
                         cAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGaa                                               >  1:1127567/1‑69 (MQ=255)
                          aGGATCGTGACCTGTGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAg                                            >  1:1949934/1‑71 (MQ=255)
                          aGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAAGGTCAGCCGAGCGAAc                                              <  1:2723678/69‑1 (MQ=255)
                          aGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAg                                            <  1:315723/71‑1 (MQ=255)
                           ggATCGTGACCTGGGTATGGTCGGTGCGGAAGAACTGCGCGTAGTGACCAAACGTCAGCCGAGCGAACAgg                                           <  1:1232419/71‑1 (MQ=255)
                           ggATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAg                                                        <  1:2146291/58‑1 (MQ=255)
                            gATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAAc                                                             >  1:1522764/1‑52 (MQ=255)
                            gATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAgg                                           >  1:2621728/1‑70 (MQ=255)
                             aTCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCaaa                                                              >  1:3344144/1‑50 (MQ=255)
                             aTCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCaaa                                                              >  1:3691751/1‑50 (MQ=255)
                             aTCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCaaa                                                              >  1:2794586/1‑50 (MQ=255)
                              tCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAAc                                        >  1:991487/1‑71 (MQ=255)
                                 tGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACt                                       <  1:379868/69‑1 (MQ=255)
                                 tGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGc                                     >  1:3408754/1‑71 (MQ=255)
                                  gACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGc                                     >  1:145741/1‑70 (MQ=255)
                                   aCCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGc                                     >  1:1606149/1‑69 (MQ=255)
                                   aCCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCg                                    <  1:218925/70‑1 (MQ=255)
                                   aCCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGt                                   >  1:207896/1‑71 (MQ=255)
                                   aCCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGt                                   <  1:1881565/71‑1 (MQ=255)
                                     cTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCa                                                                >  1:3455739/1‑40 (MQ=255)
                                     cTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTc                                                          >  1:2001788/1‑46 (MQ=255)
                                        gggATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATgc                              >  1:925057/1‑71 (MQ=255)
                                         ggATGGTCGGTGCGGAAGAACTGCGCGTAGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATgcg                             >  1:1587681/1‑71 (MQ=255)
                                          gATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATgcgc                            >  1:1466455/1‑71 (MQ=255)
                                          gATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATgcgc                            >  1:3738092/1‑71 (MQ=255)
                                              gTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGtt                        <  1:3675226/71‑1 (MQ=255)
                                              gTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGt                         <  1:3248666/70‑1 (MQ=255)
                                               tCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGaa                                               <  1:1991946/47‑1 (MQ=255)
                                               tCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTg                          >  1:2431803/1‑68 (MQ=255)
                                               tCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTc                       <  1:1861229/71‑1 (MQ=255)
                                                cGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCg                                    >  1:1653622/1‑57 (MQ=255)
                                                cGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTct                      >  1:1476088/1‑71 (MQ=255)
                                                 ggTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTct                      <  1:1942758/70‑1 (MQ=255)
                                                    gCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGaa                                               >  1:2816024/1‑42 (MQ=255)
                                                    gCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTctgctg                  >  1:1508152/1‑71 (MQ=255)
                                                       gaagaaCTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGa                <  1:2069400/70‑1 (MQ=255)
                                                        aagaaCTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAg              >  1:2969952/1‑71 (MQ=255)
                                                        aagaaCTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAg              >  1:114068/1‑71 (MQ=255)
                                                        aagaaCTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAg              <  1:1734000/71‑1 (MQ=255)
                                                          gaaCTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTTGAAGGt            <  1:2445000/71‑1 (MQ=255)
                                                          gaaCTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGt            >  1:476551/1‑71 (MQ=255)
                                                            aCTGCGCGTGGTGACCAAACGTCAGCCTAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTgg          <  1:3715445/71‑1 (MQ=255)
                                                            aCTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTgg          >  1:708272/1‑71 (MQ=255)
                                                            aCTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTgg          >  1:635947/1‑71 (MQ=255)
                                                            aCTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTgg          <  1:1662669/71‑1 (MQ=255)
                                                            aCTGCGCGTAGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTgg          >  1:87885/1‑71 (MQ=255)
                                                            aCTGCGCGTAGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTgg          >  1:2461984/1‑71 (MQ=255)
                                                             cTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGc         <  1:689216/71‑1 (MQ=255)
                                                             cTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGc         <  1:2347472/71‑1 (MQ=255)
                                                              tGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGc         >  1:1660803/1‑70 (MQ=255)
                                                                cgcgTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCg        <  1:3828809/69‑1 (MQ=255)
                                                                cgcgTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCg        >  1:711674/1‑69 (MQ=255)
                                                                cgcgTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCGa       <  1:830089/70‑1 (MQ=255)
                                                                cgcgTAGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCg        <  1:2542261/69‑1 (MQ=255)
                                                                  cgTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCGAAGt    <  1:2412103/71‑1 (MQ=255)
                                                                  cgTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCGAAGt    <  1:1413554/71‑1 (MQ=255)
                                                                  cgTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCGAAGt    <  1:3531321/71‑1 (MQ=255)
                                                                   gtggtgACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCGAAGtt   <  1:90884/71‑1 (MQ=255)
                                                                    tggtgACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCGAAGttt  <  1:2384842/71‑1 (MQ=255)
                                                                    tggtgACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCGAAGttt  <  1:1799848/71‑1 (MQ=255)
                                                                     |                                                                     
ACGCGTGAACGGCGGTCTGCTGGTTCAGGATCGTGACCTGGGGATGGTCGGTGCGGAAGAACTGCGCGTGGTGACCAAACGTCAGCCGAGCGAACAGGAACTGCGTGATGCGCTGTTCTGCTGGAAGGTGGCGAAGTTT  >  minE/2235994‑2236132

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: