Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 13,583 T→G 6.6% V474G (GTT→GGT)  dnaK → chaperone Hsp70, co‑chaperone with DnaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE13,5830TG 6.6% 415.1 / 22.0 151V474G (GTT→GGT) dnaKchaperone Hsp70, co‑chaperone with DnaJ
Reads supporting (aligned to +/- strand):  ref base T (69/72);  new base G (8/2);  total (77/74)
Fisher's exact test for biased strand distribution p-value = 9.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

ACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAG  >  minE/13514‑13653
                                                                     |                                                                      
aCAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGt                                                                        <  1:141291/70‑1 (MQ=255)
  aaaTCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTAc                                                                     >  1:1603674/1‑71 (MQ=255)
    aTCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCt                                                                   >  1:2977479/1‑71 (MQ=255)
          gggTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGata                                                             >  1:6261313/1‑71 (MQ=255)
            gTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                           <  1:4324541/71‑1 (MQ=255)
             tCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                           >  1:357250/1‑70 (MQ=255)
                gTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                       <  1:3606120/71‑1 (MQ=255)
                gTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                       <  1:5629640/71‑1 (MQ=255)
                  tCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                           >  1:855763/1‑65 (MQ=255)
                  tCAACCTAGATGGTATCAACCCGGCACCGCGCGGCAAGCCGCAGATCGAAGTTAcca                                                                   >  1:1937802/1‑56 (MQ=255)
                   cAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGata                                                             >  1:447446/1‑62 (MQ=255)
                   cAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                       >  1:2522621/1‑68 (MQ=255)
                   cAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCt                                                     >  1:5902923/1‑70 (MQ=255)
                   cAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCt                                                     >  1:5186079/1‑70 (MQ=255)
                      ccTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                           <  1:2908979/61‑1 (MQ=255)
                      ccTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACg                                                 >  1:6614363/1‑71 (MQ=255)
                           aTGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                           <  1:4802027/56‑1 (MQ=255)
                           aTGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                           <  1:4503383/56‑1 (MQ=255)
                           aTGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                           <  1:5430815/56‑1 (MQ=255)
                           aTGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                           <  1:166964/56‑1 (MQ=255)
                            tGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATcc                                           >  1:6568638/1‑71 (MQ=255)
                             ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGt                                                                        >  1:3801897/1‑41 (MQ=255)
                             ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGa                                                         >  1:2265025/1‑56 (MQ=255)
                             ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGa                                                         >  1:2500743/1‑56 (MQ=255)
                             ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                       <  1:2033667/58‑1 (MQ=255)
                             ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                       <  1:5559418/58‑1 (MQ=255)
                             ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                       <  1:1805074/58‑1 (MQ=255)
                             ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                       <  1:4185801/58‑1 (MQ=255)
                                aTCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGa                                                               >  1:5636874/1‑47 (MQ=255)
                                aTCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGa                                                               >  1:1241767/1‑47 (MQ=255)
                                aTCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGa                                                               >  1:603605/1‑47 (MQ=255)
                                aTCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGa                                                               >  1:527144/1‑47 (MQ=255)
                                aTCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGa                                                               >  1:1685520/1‑47 (MQ=255)
                                 tCAACCCGTCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCAc                                      <  1:170479/71‑1 (MQ=255)
                                 tCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCa                                       <  1:6876650/70‑1 (MQ=255)
                                 tCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCAc                                      <  1:3016907/71‑1 (MQ=255)
                                 tCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCAc                                      >  1:1825752/1‑71 (MQ=255)
                                 tCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCAc                                      <  1:597100/71‑1 (MQ=255)
                                     cccGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGttt                                  <  1:22488/71‑1 (MQ=255)
                                     cccGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGttt                                  <  1:357776/71‑1 (MQ=255)
                                     cccGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGttt                                  >  1:5999024/1‑71 (MQ=255)
                                      ccGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTc                                 <  1:1931801/71‑1 (MQ=255)
                                               cgcgGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGa                            >  1:775265/1‑67 (MQ=255)
                                               cgcgGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGa                            >  1:6926914/1‑67 (MQ=255)
                                                gcgGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGa                        >  1:6727518/1‑70 (MQ=255)
                                                 cgGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATa                      <  1:3690011/71‑1 (MQ=255)
                                                 cgGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATa                      <  1:6794176/71‑1 (MQ=255)
                                                  ggCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGa                                                   >  1:2598092/1‑41 (MQ=255)
                                                  ggCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGa                                                   >  1:4993896/1‑41 (MQ=255)
                                                    cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGt                                               <  1:6820599/43‑1 (MQ=255)
                                                    cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGt                                               <  1:5059197/43‑1 (MQ=255)
                                                    cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGaa                           <  1:3073577/63‑1 (MQ=255)
                                                    cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGaa                           <  1:6701886/63‑1 (MQ=255)
                                                    cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                    >  1:1712572/1‑70 (MQ=255)
                                                    cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                    >  1:6206015/1‑70 (MQ=255)
                                                    cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                    >  1:6089841/1‑70 (MQ=255)
                                                     aTGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                    <  1:2637055/69‑1 (MQ=255)
                                                      tGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCTAAAGATaaaaa                  >  1:2704660/1‑70 (MQ=255)
                                                      tGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaaaa                  >  1:1107382/1‑70 (MQ=255)
                                                      tGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAAc                 <  1:4166689/71‑1 (MQ=255)
                                                      tGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAAc                 <  1:72073/71‑1 (MQ=255)
                                                      tGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAAc                 <  1:467937/71‑1 (MQ=255)
                                                      tGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAAc                 <  1:6595622/71‑1 (MQ=255)
                                                            aGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCg             >  1:956456/1‑69 (MQ=255)
                                                             taTCGAAGGTAACTTCGATCTCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCgg            >  1:2266511/2‑69 (MQ=38)
                                                             taTCGAAGGTAACTTCGATCTCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCgg            >  1:322143/2‑69 (MQ=38)
                                                             taTCGAAGGTAACTTCGATCTCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGt           >  1:1304122/2‑70 (MQ=38)
                                                             taTCGAAGGTAACTTCGATCTCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGt           >  1:3985335/2‑70 (MQ=38)
                                                             taTCGAAGGTAACTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGt           >  1:1383877/2‑70 (MQ=255)
                                                             gATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTa          <  1:6774982/71‑1 (MQ=255)
                                                              aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                         <  1:4438990/55‑1 (MQ=255)
                                                              aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                         <  1:244803/55‑1 (MQ=255)
                                                              aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                         <  1:2607715/55‑1 (MQ=255)
                                                              aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                         <  1:2622870/55‑1 (MQ=255)
                                                              aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                         <  1:6159224/55‑1 (MQ=255)
                                                              aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                    <  1:5552084/60‑1 (MQ=255)
                                                              aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                    <  1:3096015/60‑1 (MQ=255)
                                                              aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                    <  1:5041751/60‑1 (MQ=255)
                                                              aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                    <  1:1901447/60‑1 (MQ=255)
                                                              aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAAc                 <  1:5286095/63‑1 (MQ=255)
                                                               tCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGt           >  1:5000742/1‑68 (MQ=255)
                                                               tCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGt           >  1:4455031/1‑68 (MQ=255)
                                                               tCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGt           >  1:6613648/1‑68 (MQ=255)
                                                               tCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa        <  1:4565934/71‑1 (MQ=255)
                                                               tCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa        <  1:6057743/71‑1 (MQ=255)
                                                                cGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTa          >  1:3307567/1‑68 (MQ=255)
                                                                cGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTa          >  1:326927/1‑68 (MQ=255)
                                                                 gAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa        >  1:4019769/1‑69 (MQ=255)
                                                                 gAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa        >  1:1434893/1‑69 (MQ=255)
                                                                 gAAGGTAACTTCGATCTCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAag       <  1:476033/70‑1 (MQ=38)
                                                                  aaGTTACCTTCGATATCGATGCTGACGGTATCCTGCCCGTTTCCGCGAAAGATAAAAACAGCGGTaaa        >  1:4243467/1‑68 (MQ=255)
                                                                  aaGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa        >  1:2492849/1‑68 (MQ=255)
                                                                  aaGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa        >  1:442240/1‑68 (MQ=255)
                                                                  aaGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa        >  1:901120/1‑68 (MQ=255)
                                                                  aaGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa        >  1:804194/1‑68 (MQ=255)
                                                                  aaGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaa         >  1:2008821/1‑67 (MQ=255)
                                                                   aGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAagag     >  1:2340652/1‑70 (MQ=255)
                                                                   aGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTACCGCGAAAGATAAAAACAGCGGTAAagag     >  1:6350295/1‑70 (MQ=255)
                                                                     ttACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGATAGATAAAAACAGCGGTa          <  1:4535994/63‑1 (MQ=255)
                                                                     ttACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCag  <  1:2702548/71‑1 (MQ=255)
                                                                     |                                                                      
ACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAG  >  minE/13514‑13653

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: