Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 15,386 C→T 19.0% intergenic (+88/+59) dnaJ → / ← mokC chaperone Hsp40, co‑chaperone with DnaK/regulatory protein for HokC, overlaps CDS of hokC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE15,3860CT19.0% 114.8 / 21.3 53intergenic (+88/+59)dnaJ/mokCchaperone Hsp40, co‑chaperone with DnaK/regulatory protein for HokC, overlaps CDS of hokC
Reads supporting (aligned to +/- strand):  ref base C (23/18);  new base T (5/5);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCTCTTCTC  >  minE/15318‑15456
                                                                    |                                                                      
gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTc                                                                      <  1:1515669/71‑1 (MQ=255)
    gCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAt                                                                  <  1:386441/71‑1 (MQ=255)
       gctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCATAAg                                                               <  1:1957464‑M1/11‑1 (MQ=255)
          cTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAgcgc                                                            >  1:3060669/1‑71 (MQ=255)
             ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAgcgat                                                           <  1:1826134/69‑3 (MQ=255)
             ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTa                                                          <  1:961201/70‑1 (MQ=255)
               tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt                                                       >  1:3309587/1‑71 (MQ=255)
               tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt                                                       >  1:3860129/1‑71 (MQ=255)
                 aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACtt                                                      <  1:787139/70‑1 (MQ=255)
                 aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTa                                                     <  1:2997661/71‑1 (MQ=255)
                 aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTa                                                     <  1:131474/71‑1 (MQ=255)
                 aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTa                                                     <  1:1724117/71‑1 (MQ=255)
                   aaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAg                                                   >  1:4110892/1‑71 (MQ=255)
                     tAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAggg                                                 >  1:1855796/1‑71 (MQ=255)
                     tAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAggg                                                 >  1:3488915/1‑71 (MQ=255)
                       gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACt                                                                       >  1:1553591/1‑47 (MQ=255)
                       gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt                                               >  1:2225032/1‑71 (MQ=255)
                       gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGt                                                <  1:2633244/70‑1 (MQ=255)
                              ttGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                        <  1:2807949/71‑1 (MQ=255)
                              ttGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                        <  1:214050/71‑1 (MQ=255)
                              tggtgctgtaagaccatacataaaaaagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                        <  1:1449256‑M1/34‑1 (MQ=255)
                               tGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCAc                                                                        >  1:2485265/1‑38 (MQ=255)
                                gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAAc                                                        >  1:4055419/1‑53 (MQ=255)
                                gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTgg                                           >  1:1591709/1‑66 (MQ=255)
                                 tgctgtaagaccatacataaaaaagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                     <  1:1164800‑M1/37‑1 (MQ=255)
                                 aaGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTcc                                                                     >  1:660592/1‑39 (MQ=255)
                                  aGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                    >  1:505865/1‑71 (MQ=255)
                                     tatGCGAGCACCTGTAAAGTGGCGGGGATCAATCCCATAAGCGCTAACTTAAGGGTtgt                                             >  1:2223481/1‑59 (MQ=255)
                                      taagaccatacataaaaaagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                 >  1:1726692‑M1/30‑70 (MQ=255)
                                      taagaccatacataaaaaagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                 >  1:903766‑M1/30‑70 (MQ=255)
                                       tGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                        <  1:3290071/62‑1 (MQ=255)
                                        agaccatacataaaaaagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                  <  1:961605‑M1/40‑1 (MQ=255)
                                          gAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                            <  1:1743659/71‑1 (MQ=255)
                                            gCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                        <  1:49565/57‑1 (MQ=255)
                                            catacataaaaaagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTgg                                           >  1:2678584‑M1/24‑54 (MQ=255)
                                             cACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt                                               >  1:1168423/1‑49 (MQ=255)
                                              tacataaaaaagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                        <  1:2336033‑M1/50‑1 (MQ=255)
                                               ccTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtg                                              <  1:123433/48‑1 (MQ=255)
                                                cTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                      <  1:3123788/71‑1 (MQ=255)
                                                 tGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTgg                                           <  1:181238/49‑1 (MQ=255)
                                                 tGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTgg                                           <  1:50531/49‑1 (MQ=255)
                                                  taaaaaagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                               >  1:3364167‑M1/18‑60 (MQ=255)
                                                    aaaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgt                  >  1:2453078/1‑71 (MQ=255)
                                                    aaaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcg                   >  1:2677895/1‑70 (MQ=255)
                                                     aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                   >  1:4111992/1‑53 (MQ=255)
                                                       gTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgcc               >  1:943650/1‑71 (MQ=255)
                                                          gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtg                                              >  1:2475878/1‑37 (MQ=40)
                                                          gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtg                                              >  1:1528640/1‑37 (MQ=40)
                                                           cggggagtAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaacgtgccgat            >  1:1108840‑M1/9‑68 (MQ=255)
                                                               gataATTCCCATAAGCGCTAACTTAAGGGTTGTgg                                           >  1:735675‑M1/5‑35 (MQ=255)
                                                                  cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGctc         >  1:14398/1‑66 (MQ=255)
                                                                  cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGctc         >  1:483374/1‑66 (MQ=255)
                                                                   ACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCtcttc    >  1:2821423‑M1/1‑70 (MQ=255)
                                                                   ACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCtcttc    >  1:3025373‑M1/1‑70 (MQ=255)
                                                                    cTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCTCTtctc  <  1:908813/71‑1 (MQ=255)
                                                                    |                                                                      
GTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCTCTTCTC  >  minE/15318‑15456

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: