Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 15,001 C→T 27.2% G278G (GGC→GGT dnaJ → chaperone Hsp40, co‑chaperone with DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE15,0010CT27.2% 93.9 / 32.3 84G278G (GGC→GGTdnaJchaperone Hsp40, co‑chaperone with DnaK
Reads supporting (aligned to +/- strand):  ref base C (37/24);  new base T (8/15);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 4.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTA  >  minE/14938‑15068
                                                               |                                                                   
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:1217574/69‑1 (MQ=255)
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:2972140/69‑1 (MQ=255)
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:2265323/69‑1 (MQ=255)
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:392671/69‑1 (MQ=255)
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:491239/69‑1 (MQ=255)
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:501025/69‑1 (MQ=255)
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:1688161/69‑1 (MQ=255)
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:171623/69‑1 (MQ=255)
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:1745573/69‑1 (MQ=255)
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:1748598/69‑1 (MQ=255)
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:254453/69‑1 (MQ=255)
gCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:1858185/69‑1 (MQ=255)
 cGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:2496195/68‑1 (MQ=255)
  gTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCg                                                              <  1:893372/69‑1 (MQ=255)
  gTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCg                                                              <  1:2747855/69‑1 (MQ=255)
  gTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCg                                                              <  1:2793712/69‑1 (MQ=255)
  gTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCg                                                              <  1:1787178/69‑1 (MQ=255)
  gTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCg                                                              <  1:756985/69‑1 (MQ=255)
  gTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCg                                                              <  1:753976/69‑1 (MQ=255)
  gTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCg                                                              <  1:1142446/69‑1 (MQ=255)
  gTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCg                                                              <  1:1347844/69‑1 (MQ=255)
  gTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCg                                                              <  1:562683/69‑1 (MQ=255)
  gTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCg                                                              <  1:718259/69‑1 (MQ=255)
         caacaaCCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTAcc                                                       >  1:1925068/1‑69 (MQ=255)
         caacaaCCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTAcc                                                       >  1:760087/1‑69 (MQ=255)
         caacaaCCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTAcc                                                       >  1:571407/1‑69 (MQ=255)
         caacaaCCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTAcc                                                       >  1:3298105/1‑69 (MQ=255)
              acCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:1771347/55‑1 (MQ=255)
                cTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCtt                                                <  1:1715210/69‑1 (MQ=255)
                cTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCtt                                                <  1:422625/69‑1 (MQ=255)
                cTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCtt                                                <  1:2842312/69‑1 (MQ=255)
                 tGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTg                                               <  1:1899344/69‑1 (MQ=255)
                 tGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTg                                               <  1:3143612/69‑1 (MQ=255)
                  gTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAATCGAAGTACCGACCCTTGa                                              >  1:1063428/1‑69 (MQ=255)
                  gTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAATCGAAGTACCGACCCTTGa                                              >  1:2251059/1‑69 (MQ=255)
                  gTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAATCGAAGTACCGACCCTTGa                                              >  1:2681873/1‑69 (MQ=255)
                  gTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAATCGAAGTACCGACCCTTGa                                              >  1:2852426/1‑69 (MQ=255)
                  gTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAATCGAAGTACCGACCCTTGa                                              >  1:2905160/1‑69 (MQ=255)
                  gTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAATCGAAGTACCGACCCTTGa                                              >  1:1445345/1‑69 (MQ=255)
                  gTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAATCGAAGTACCGACCCTTGa                                              >  1:1277523/1‑69 (MQ=255)
                  gTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAATCGAAGTACCGACCCTTGa                                              >  1:1680870/1‑69 (MQ=255)
                    aTTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGTGAAAt                                                                <  1:502928/49‑1 (MQ=255)
                     ttGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTa                                                         >  1:1587985/1‑55 (MQ=255)
                     ttGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGa                                              >  1:878537/1‑66 (MQ=255)
                     ttGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATgg                                           >  1:2564636/1‑69 (MQ=255)
                     ttGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATgg                                           >  1:171322/1‑69 (MQ=255)
                     ttGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATgg                                           >  1:2922479/1‑69 (MQ=255)
                     ttGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATgg                                           >  1:1120894/1‑69 (MQ=255)
                     ttGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATgg                                           >  1:1645337/1‑69 (MQ=255)
                     ttGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATgg                                           >  1:1641139/1‑69 (MQ=255)
                     ttGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATgg                                           >  1:2877147/1‑69 (MQ=255)
                               ccGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaaa                                 >  1:939125/1‑69 (MQ=255)
                               ccGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaaa                                 >  1:785211/1‑69 (MQ=255)
                               ccGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaa                                  >  1:3350918/1‑68 (MQ=255)
                               ccGATCAACTTCGCTATGGCGGCGCTGGGTGGCAAAATCGAAGTACCGACCCTTGATGGTCGCGTCaaa                                 >  1:242676/1‑69 (MQ=255)
                                    cAACTTCTCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGaa                            <  1:1500608/69‑1 (MQ=255)
                                    cAACTTCTCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGaa                            <  1:1496606/69‑1 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaa                                  >  1:1124751/1‑63 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaa                                  >  1:164714/1‑63 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaa                                  >  1:271812/1‑63 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaa                                  >  1:741959/1‑63 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaa                                  >  1:118443/1‑63 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaa                                  >  1:1575838/1‑63 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaa                                  >  1:322980/1‑63 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaa                                  >  1:3375018/1‑63 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCaa                                  >  1:489505/1‑63 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGaa                            <  1:868381/69‑1 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGaa                            <  1:3317125/69‑1 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGaa                            <  1:493484/69‑1 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGaa                            <  1:3036955/69‑1 (MQ=255)
                                    cAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAACTACCGACCCTTGATGGTCGCGTCAAACTGaa                            >  1:908744/1‑69 (MQ=255)
                                        ttCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTg                        <  1:2980189/69‑1 (MQ=255)
                                        ttCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTg                        <  1:1119893/69‑1 (MQ=255)
                                          cGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGcc                      >  1:1174226/1‑69 (MQ=255)
                                          cGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGcc                      >  1:1857833/1‑69 (MQ=255)
                                          cGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGcc                      >  1:2097580/1‑69 (MQ=255)
                                                 gcggcgCTGGGTGGCGAAATCGAAGTACCGCCCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGaa               >  1:2235690/1‑69 (MQ=255)
                                                 gcggcgCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGaa               >  1:1372759/1‑69 (MQ=255)
                                                 gcggcgCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGaa               >  1:1615589/1‑69 (MQ=255)
                                                 gcggcgCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGaa               >  1:280608/1‑69 (MQ=255)
                                                 gcggcgCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGaa               >  1:1954639/1‑69 (MQ=255)
                                                            gggCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTg                        <  1:1345646/48‑1 (MQ=255)
                                                              gCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTa  >  1:2896317/1‑69 (MQ=255)
                                                              gCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTa  >  1:105456/1‑69 (MQ=255)
                                                               |                                                                   
GCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTA  >  minE/14938‑15068

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: