Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 13,094 A→T 42.1% D311V (GAT→GTT)  dnaK → chaperone Hsp70, co‑chaperone with DnaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE13,0940AT42.1% 38.3 / 75.6 64D311V (GAT→GTT) dnaKchaperone Hsp70, co‑chaperone with DnaJ
Reads supporting (aligned to +/- strand):  ref base A (16/21);  new base T (15/12);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATAT  >  minE/13026‑13163
                                                                    |                                                                     
cgACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTc                                                                     <  1:3168501/70‑1 (MQ=255)
aGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTc                                                                     <  1:867047/71‑1 (MQ=255)
 gACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCt                                                                    <  1:2969641/71‑1 (MQ=255)
    gcgACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGt                                                                 <  1:488142/71‑1 (MQ=255)
     cgACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTa                                                                <  1:1962237/71‑1 (MQ=255)
        ccGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAAc                                                             >  1:1668595/1‑71 (MQ=255)
           gTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGt                                                          <  1:2857606/71‑1 (MQ=255)
           gTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGt                                                          >  1:102384/1‑71 (MQ=255)
           gTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGt                                                          >  1:3428857/1‑71 (MQ=255)
             ccGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAAcc                                                            <  1:1539028/67‑1 (MQ=255)
              cGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTg                                                                   >  1:386296/1‑59 (MQ=255)
              cGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTg                                                                   >  1:1377900/1‑59 (MQ=255)
                 aaCACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCAtt                                                    >  1:1472583/1‑71 (MQ=255)
                 aaCACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCAtt                                                    <  1:2848674/71‑1 (MQ=255)
                 aaCACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCAtt                                                    >  1:2634449/1‑71 (MQ=255)
                 aaCACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCAt                                                     >  1:1551778/1‑70 (MQ=255)
                  acacATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTg                                                   <  1:1976498/71‑1 (MQ=255)
                      aTGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCATTGAg                                                 >  1:1318979/1‑69 (MQ=255)
                      aTGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCATTGAGcc                                               >  1:1666927/1‑71 (MQ=255)
                       tGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCg                                              <  1:2955258/71‑1 (MQ=255)
                       tGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCg                                              <  1:2787770/71‑1 (MQ=255)
                       tGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCg                                              <  1:1061704/71‑1 (MQ=255)
                        gAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTg                                                   <  1:2158066/65‑1 (MQ=255)
                         aaCATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCt                                            >  1:3711279/1‑71 (MQ=255)
                          aCATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCATTGAGCCGCt                                            >  1:851106/1‑70 (MQ=255)
                                aaGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAg                                       <  1:2489664/69‑1 (MQ=255)
                                 aGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTg                                    >  1:2057838/1‑71 (MQ=255)
                                     aCTCGTGCGAAACTGGAAAGCCTGTTTGAAGTTCTGGTAAACCGTTCCAtt                                                    <  1:312808/51‑1 (MQ=39)
                                     aCTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACt                                <  1:2809303/71‑1 (MQ=255)
                                      cTCGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTg                               <  1:2178947/71‑1 (MQ=255)
                                        cGTGCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCATTGAGCCGCTGaaa                                        <  1:1886859/60‑1 (MQ=255)
                                        cGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTg                                           <  1:2248659/57‑1 (MQ=255)
                                        cGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTg                                           <  1:531207/57‑1 (MQ=255)
                                         gTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAg                            <  1:1389751/71‑1 (MQ=255)
                                         gTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAg                            <  1:1788155/71‑1 (MQ=255)
                                         gTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAg                            <  1:1808803/71‑1 (MQ=255)
                                          tGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAAc                                                             <  1:2409426/37‑1 (MQ=255)
                                          tGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAAc                                                             <  1:2384484/37‑1 (MQ=255)
                                           gCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTc                                                        >  1:1787033/1‑41 (MQ=255)
                                           gCGAAACTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTc                                                        >  1:3487856/1‑41 (MQ=255)
                                            cGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTg                                    <  1:1954509/60‑1 (MQ=255)
                                            cGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGAc                         >  1:2969437/1‑71 (MQ=255)
                                             gAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTcc                                                       <  1:272293/40‑1 (MQ=255)
                                             gAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTcc                                                       <  1:2380225/40‑1 (MQ=255)
                                             gAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACg                        >  1:545664/1‑71 (MQ=255)
                                              aaaCTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGc                       <  1:2650164/71‑1 (MQ=255)
                                               aaCTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCt                      >  1:3685175/1‑71 (MQ=255)
                                               aaCTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCt                      >  1:63472/1‑71 (MQ=255)
                                                aCTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCAt                                                     >  1:2292928/1‑39 (MQ=255)
                                                aCTGGAAAGCCTGGTTGAAGTTCTGGTAAACCGTTCCAt                                                     >  1:3167161/1‑39 (MQ=255)
                                                     aaaGCCTGGTTGAAGTTCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTg                >  1:1666756/1‑71 (MQ=255)
                                                         ccTGGTTGAAGTTCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCg            >  1:1642706/1‑71 (MQ=255)
                                                            ggTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTa          >  1:2057774/1‑70 (MQ=255)
                                                            ggTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTAt         >  1:3705494/1‑71 (MQ=255)
                                                             gTTGAAGTTCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATc        >  1:3622389/1‑71 (MQ=255)
                                                              gtGAAGTTCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTg                     <  1:3181576/56‑1 (MQ=255)
                                                                gAAGTTCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGa     >  1:737387/1‑71 (MQ=255)
                                                                gAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCa                             >  1:1434344/1‑47 (MQ=255)
                                                                gAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCa                             >  1:284808/1‑47 (MQ=255)
                                                                  aGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGata   >  1:2007711/1‑71 (MQ=255)
                                                                  aGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGat    >  1:1366204/1‑70 (MQ=255)
                                                                   gATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGatat  <  1:3320123/71‑1 (MQ=255)
                                                                   gATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGatat  <  1:1648480/71‑1 (MQ=255)
                                                                   gATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGatat  <  1:2275473/71‑1 (MQ=255)
                                                                    |                                                                     
AGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATAT  >  minE/13026‑13163

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: