Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 15,386 C→T 20.5% 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,3860CT20.5% 104.6 / 22.4 54intergenic (+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 (22/19);  new base T (8/3);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

GTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCTCTTC  >  minE/15318‑15454
                                                                    |                                                                    
gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACt                                                                     >  1:2462978/1‑70 (MQ=255)
  gggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTcc                                                                   >  1:3841102/1‑70 (MQ=255)
  agtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTcc                                                                   >  1:825505‑M1/66‑70 (MQ=255)
    tgtgctatctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCAt                                                                >  1:564670‑M1/64‑71 (MQ=255)
      tgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCATaa                                                              >  1:3359500‑M1/62‑71 (MQ=255)
      tgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCATaa                                                              >  1:3503225‑M1/62‑71 (MQ=255)
      aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATaa                                                              <  1:1472933/71‑1 (MQ=255)
        gCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAgc                                                            >  1:3355855/1‑71 (MQ=255)
          cTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCAc                                                                      <  1:1190796/59‑1 (MQ=255)
          cTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAgcgc                                                          >  1:3076065/1‑71 (MQ=255)
             ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCt                                                         <  1:2698701/69‑1 (MQ=255)
               tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTc                                                                    <  1:1750825/56‑1 (MQ=255)
               tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTa                                                        <  1:3209705/68‑1 (MQ=255)
               tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt                                                     >  1:3299374/1‑71 (MQ=255)
               cgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCATAAGCGCTaa                                                       >  1:2377006‑M1/53‑69 (MQ=255)
                gtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCATAAGCGCTAACtt                                                    >  1:388121‑M1/52‑71 (MQ=255)
                aaaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt                                                     >  1:3013502/1‑70 (MQ=255)
                  caatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCATAAgcgc                                                          <  1:2074114‑M1/14‑1 (MQ=255)
                   aaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAg                                                 >  1:3109025/1‑71 (MQ=255)
                       gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt                                             <  1:2301506/71‑1 (MQ=255)
                           cactggtgctgtaagaccatacataaaaaagccgggtataATTCCCATAAGCGCTAACTTAAGGGTTGTgg                                         >  1:3335450‑M1/41‑71 (MQ=255)
                              ttGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                      <  1:584890/71‑1 (MQ=255)
                              tggtgctgtaagaccatacataaaaaagccggggataATTCCCATAAGCGCTa                                                        >  1:3298036‑M1/38‑53 (MQ=255)
                                gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCAc                                                                      >  1:71128/1‑37 (MQ=255)
                                gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                    >  1:1099710/1‑71 (MQ=255)
                                   gATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                 >  1:1724997/1‑71 (MQ=255)
                                    atatGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATaa                                                              >  1:1312009/1‑41 (MQ=255)
                                      atGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTa                                                   >  1:3551033/1‑50 (MQ=255)
                                       tGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                              <  1:2917456/70‑1 (MQ=255)
                                       tGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                              <  1:3897397/70‑1 (MQ=255)
                                         cGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTaa                                                  >  1:1283318/1‑48 (MQ=255)
                                         cGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                           >  1:3952050/1‑71 (MQ=255)
                                          gAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                   >  1:1519043/1‑62 (MQ=255)
                                          gAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                          <  1:243991/71‑1 (MQ=255)
                                               ccTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                     <  1:1035708/71‑1 (MQ=255)
                                                cTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTgg                                         <  1:736040/50‑1 (MQ=255)
                                                 ataaaaaagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                   <  1:234280‑M1/53‑1 (MQ=255)
                                                  gTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                  >  1:3339663/1‑71 (MQ=255)
                                                   tAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt                                             <  1:340530/43‑1 (MQ=255)
                                                   tAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                       <  1:3302265/49‑1 (MQ=255)
                                                    aaaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgt                >  1:3893103/1‑71 (MQ=255)
                                                     aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAgg                                                >  1:1679547/1‑38 (MQ=255)
                                                     aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtg               >  1:809605/1‑71 (MQ=255)
                                                     aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtg               >  1:2771355/1‑71 (MQ=255)
                                                       gTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                      <  1:1213525/46‑1 (MQ=255)
                                                       aagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgcc             <  1:1899622‑M1/59‑1 (MQ=255)
                                                        tGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgccg            <  1:305575/71‑1 (MQ=255)
                                                         tgCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgccg            <  1:1060814/69‑1 (MQ=25)
                                                         ggCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgccgt           >  1:1943939/1‑71 (MQ=255)
                                                          gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgccgtg          <  1:2857892/71‑1 (MQ=255)
                                                          gCGGGGAGTAATCCCATAAGCGCTAACTAAGGGTTGTGGTATTACGACTGATATGATTTAAcgtgccgat          >  1:3106755/1‑68 (MQ=25)
                                                            ggggATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                     >  1:3929615/1‑42 (MQ=255)
                                                             gggATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                          <  1:1580968/52‑1 (MQ=255)
                                                               gATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGctctc     <  1:3335497/71‑1 (MQ=255)
                                                                 tCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCtctt   <  1:2917546/71‑1 (MQ=255)
                                                                  cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCtcttc  >  1:1601470/1‑71 (MQ=255)
                                                                    cTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                    >  1:2305853/1‑35 (MQ=255)
                                                                    cTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgccg            <  1:2956801/59‑1 (MQ=255)
                                                                    |                                                                    
GTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCTCTTC  >  minE/15318‑15454

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: