Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,570,788 T→C 19.6% T231T (ACT→ACC yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,570,7880TC19.6% 204.6 / 43.6 102T231T (ACT→ACCyibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/41);  new base C (16/4);  total (57/45)
Fisher's exact test for biased strand distribution p-value = 2.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 103 total aligned reads displayed.

CAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCC  >  minE/2570720‑2570855
                                                                    |                                                                   
cAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACt                                                                     >  1:2620716/1‑69 (MQ=255)
  aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgt                                                                 <  1:1579372/71‑2 (MQ=255)
  aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                                                 >  1:766332/1‑71 (MQ=255)
  aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                                                 >  1:5504789/1‑71 (MQ=255)
      aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCATGTa                                                               <  1:3341062/69‑1 (MQ=255)
      aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACt                                                             >  1:129133/1‑71 (MQ=255)
       cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTa                                                               <  1:1242373/68‑1 (MQ=255)
          aaGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATc                                                         >  1:5325897/1‑71 (MQ=255)
          aaGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGAtt                                                         >  1:2660324/1‑70 (MQ=255)
          aaGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGAtt                                                         >  1:5628654/1‑70 (MQ=255)
             cAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCg                                                       <  1:930966/70‑1 (MQ=255)
              aaaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCgt                                                      <  1:212453/70‑2 (MQ=255)
               aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                                                    <  1:3208415/71‑1 (MQ=255)
               aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                                                    >  1:1253642/1‑71 (MQ=255)
                aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCa                                                   <  1:1698509/71‑1 (MQ=255)
                aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTa                                                               >  1:3456040/1‑59 (MQ=255)
                aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                                                   >  1:2880584/1‑70 (MQ=255)
                aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                                                   >  1:2733471/1‑70 (MQ=255)
                 atatAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAg                                                  >  1:1906756/1‑71 (MQ=255)
                  tataACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAgg                                                 >  1:3649775/1‑71 (MQ=255)
                    taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                               >  1:4119278/1‑71 (MQ=255)
                    taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                               <  1:3749721/71‑1 (MQ=255)
                    taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                               >  1:1397815/1‑71 (MQ=255)
                    taACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTa                                               >  1:5628176/1‑71 (MQ=255)
                      aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                                             <  1:2085988/71‑1 (MQ=255)
                       ccTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCgtc                                                     >  1:3467639/1‑60 (MQ=255)
                        cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                                                 <  1:2222215/49‑1 (MQ=255)
                        cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                                                 <  1:2274309/49‑1 (MQ=255)
                        cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                                                 <  1:449845/49‑1 (MQ=255)
                        cTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTa                                           >  1:1271195/1‑71 (MQ=255)
                         ttGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGa                                                           >  1:562770/1‑54 (MQ=255)
                          tGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACg                                         >  1:2174325/1‑71 (MQ=255)
                           gAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                                                           >  1:2387129/1‑52 (MQ=255)
                            aGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGc                                        <  1:2295133/70‑1 (MQ=255)
                            aGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGc                                        <  1:1028155/70‑1 (MQ=255)
                               aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt                                    <  1:3218372/71‑1 (MQ=255)
                               aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt                                    <  1:3225794/71‑1 (MQ=255)
                                gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGc                                                      >  1:359309/1‑52 (MQ=255)
                                gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                                             <  1:1560887/61‑1 (MQ=255)
                                gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                   <  1:2670348/71‑1 (MQ=255)
                                gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                   >  1:3098634/1‑71 (MQ=255)
                                  cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                                     >  1:3899776/1‑67 (MQ=255)
                                  cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                                     >  1:5506930/1‑67 (MQ=255)
                                  cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                                     >  1:917187/1‑67 (MQ=255)
                                     tCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCtt                               <  1:124460/70‑1 (MQ=255)
                                       aTTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc                            >  1:2726988/1‑71 (MQ=255)
                                        ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                           <  1:96965/71‑1 (MQ=255)
                                        ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                           <  1:3303062/71‑1 (MQ=255)
                                        ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                           <  1:3667828/71‑1 (MQ=255)
                                        ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                           <  1:3065483/71‑1 (MQ=255)
                                        ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                           <  1:4652072/71‑1 (MQ=255)
                                        ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                           <  1:5155473/71‑1 (MQ=255)
                                        ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                           <  1:2016074/71‑1 (MQ=255)
                                        ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                           <  1:1277911/71‑1 (MQ=255)
                                        ttCGCGCGCCGAGCAAAGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                           <  1:1593534/71‑1 (MQ=255)
                                         tCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc                          >  1:4405786/1‑71 (MQ=255)
                                          cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                           >  1:4819656/1‑69 (MQ=255)
                                             gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                              >  1:3886713/1‑47 (MQ=255)
                                             gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                              >  1:974992/1‑47 (MQ=255)
                                             gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                              >  1:3157193/1‑47 (MQ=255)
                                             gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                              >  1:2029958/1‑47 (MQ=255)
                                             gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                              >  1:2616414/1‑47 (MQ=255)
                                              cgcCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGt                     >  1:4854762/1‑71 (MQ=255)
                                              cgcCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGt                     >  1:2187208/1‑71 (MQ=255)
                                              cgcCGAGCCATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTatg                                         >  1:5221741/1‑49 (MQ=255)
                                               gcagAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                    >  1:5095792/4‑71 (MQ=255)
                                               gcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                    >  1:5027913/1‑71 (MQ=255)
                                                  gAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg                 <  1:4523482/71‑1 (MQ=255)
                                                  gAGCAATGGCTACGTTACTCAGGTACTGAACCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTcc                   >  1:2462052/1‑69 (MQ=255)
                                                   aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                   >  1:13360/1‑52 (MQ=255)
                                                   aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                   >  1:3905850/1‑52 (MQ=255)
                                                   aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                   >  1:1658442/1‑52 (MQ=255)
                                                    gCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg                 <  1:3771478/69‑1 (MQ=255)
                                                    gCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg                 <  1:2791965/69‑1 (MQ=255)
                                                       aTGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATg            >  1:544441/1‑71 (MQ=255)
                                                       aTGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATg            <  1:4950043/71‑1 (MQ=255)
                                                       aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                    <  1:5461260/63‑1 (MQ=255)
                                                       aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                    <  1:2249279/63‑1 (MQ=255)
                                                        tGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATgg           >  1:622311/1‑71 (MQ=255)
                                                        tGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGa              <  1:198856/68‑1 (MQ=255)
                                                         ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg                  >  1:4733665/1‑63 (MQ=255)
                                                         ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg                  >  1:951436/1‑63 (MQ=255)
                                                         ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATgg           >  1:2486290/1‑70 (MQ=255)
                                                         ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt          >  1:5129572/1‑71 (MQ=255)
                                                         ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt          >  1:4296935/1‑71 (MQ=255)
                                                          gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc         >  1:3392954/1‑71 (MQ=255)
                                                          gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc         <  1:2328269/71‑1 (MQ=255)
                                                           cTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCt        <  1:4172809/71‑1 (MQ=255)
                                                              cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt       <  1:1319461/69‑1 (MQ=255)
                                                              cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc      <  1:4657486/70‑1 (MQ=255)
                                                              cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCa     >  1:2203160/1‑71 (MQ=255)
                                                               gTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc                            <  1:4329799/47‑1 (MQ=255)
                                                                ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc                          <  1:3557122/48‑1 (MQ=255)
                                                                ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                    <  1:2937265/54‑1 (MQ=255)
                                                                ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATc   >  1:508147/1‑71 (MQ=255)
                                                                  aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTg                        >  1:2121406/1‑48 (MQ=255)
                                                                  aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTg                        >  1:2601556/1‑48 (MQ=255)
                                                                  aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATcc  >  1:415778/1‑70 (MQ=255)
                                                                   cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                   <  1:879172/36‑1 (MQ=255)
                                                                   cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                   <  1:3699059/36‑1 (MQ=255)
                                                                    |                                                                   
CAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCC  >  minE/2570720‑2570855

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: