Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,326,669 C→T 17.1% L114L (CTG→CTA wcaK ← predicted pyruvyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,326,6690CT17.1% 209.9 / 31.9 82L114L (CTG→CTAwcaKpredicted pyruvyl transferase
Reads supporting (aligned to +/- strand):  ref base C (28/40);  new base T (3/11);  total (31/51)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTG  >  minE/1326603‑1326739
                                                                  |                                                                      
cGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAcagcag                                                                      <  1:3466805/69‑1 (MQ=255)
  gCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGTAGgcg                                                                   >  1:3809160/1‑70 (MQ=255)
       ccATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACg                                                               >  1:1976393/1‑69 (MQ=255)
       ccATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGaa                                                             <  1:2293826/71‑1 (MQ=255)
        cATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAg                                                            <  1:1246386/71‑1 (MQ=255)
         aTAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGt                                                           <  1:773520/71‑1 (MQ=255)
         aTAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGt                                                           <  1:3606789/71‑1 (MQ=255)
          tAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACg                                                               <  1:2286460/66‑1 (MQ=255)
           agagATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGTAGGCGCACGAAGTCg                                                         >  1:1842776/1‑71 (MQ=255)
            gagaTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCgg                                                        >  1:80350/1‑71 (MQ=255)
            gagaTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCg                                                         <  1:3655468/70‑1 (MQ=255)
            gagaTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCg                                                         <  1:2319582/70‑1 (MQ=255)
            gagaTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCg                                                         <  1:2764978/70‑1 (MQ=255)
            gagaTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCg                                                         <  1:1910264/70‑1 (MQ=255)
            gagaTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCg                                                         <  1:1932450/70‑1 (MQ=255)
            gagaTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCg                                                         <  1:2672872/70‑1 (MQ=255)
            gagaTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCg                                                         <  1:2710778/70‑1 (MQ=255)
             agaTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGt                                                       <  1:2544471/71‑1 (MQ=255)
              gaTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGTAGGCGCACGAAGTCGGt                                                       <  1:3399297/70‑1 (MQ=255)
              gaTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTg                                                      <  1:2394554/71‑1 (MQ=255)
              caTCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGTAGGCGCACGAAGTCGGt                                                       <  1:1405090/69‑1 (MQ=255)
                tCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGaa                                                    >  1:647687/1‑71 (MQ=255)
                 cGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCg                                                         >  1:2744268/1‑65 (MQ=255)
                 cGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGa                                                     >  1:202334/1‑69 (MQ=255)
                 cGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAAt                                                   >  1:3875931/1‑71 (MQ=255)
                   aCAAAAAACTATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATcc                                                 >  1:3811973/1‑71 (MQ=255)
                   aCAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTg                                                      <  1:239775/66‑1 (MQ=255)
                   aCAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATcc                                                 >  1:3622577/1‑71 (MQ=255)
                   aCAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATcc                                                 >  1:2763073/1‑71 (MQ=255)
                     aaaaaaCGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCt                                               <  1:2026565/71‑1 (MQ=255)
                      aaaaaCGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCa                                                                 >  1:3070489/1‑52 (MQ=255)
                      aaaaaCGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGa                                                              <  1:976734/55‑1 (MQ=255)
                      aaaaaCGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCt                                               <  1:1038519/70‑1 (MQ=255)
                       aaaaCGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATcc                                                 >  1:2709847/1‑67 (MQ=255)
                       aaaaCGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTgg                                             >  1:2120599/1‑71 (MQ=255)
                        aaaCGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGaa                                                    >  1:1297307/1‑63 (MQ=255)
                        aaaCGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGaa                                                    >  1:2257518/1‑63 (MQ=255)
                         aaCGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTg                                                      >  1:1682417/1‑60 (MQ=255)
                          aCGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGggc                                           >  1:3687420/1‑70 (MQ=255)
                           cGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGggcga                                         <  1:1779545/71‑1 (MQ=255)
                           cGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGggcga                                         <  1:1426051/71‑1 (MQ=255)
                                cgccgACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGggcgatggc                                     >  1:1431615/1‑70 (MQ=255)
                                cgccgACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGggcgatggc                                     >  1:519124/1‑70 (MQ=255)
                                 gccgACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGggcgatggcga                                   <  1:3133473/71‑1 (MQ=255)
                                  ccgACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGggcgatggcga                                   >  1:3750079/1‑70 (MQ=255)
                                  ccgACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGggcgatggcga                                   >  1:1236325/1‑70 (MQ=255)
                                  ccgACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGggcgatggcga                                   <  1:763128/70‑1 (MQ=255)
                                   cgACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGata                                 <  1:205675/71‑1 (MQ=255)
                                    gACCTGGATAATGGCGTCGTAACCTGACAGTAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGatat                                <  1:115765/71‑1 (MQ=255)
                                    gACCTGGATAATGGCGTCGTAACCTGACAGTAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGatat                                <  1:2122260/71‑1 (MQ=255)
                                    gACCTGGATAATGGCGTCGTAACCTGACAGTAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGatat                                <  1:3256011/71‑1 (MQ=255)
                                    gACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACTAAGTCGGTGAATCCCTGGGCGATGGCGatat                                >  1:3804770/1‑71 (MQ=255)
                                    gACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGatat                                >  1:442933/1‑71 (MQ=255)
                                    gACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGatat                                >  1:3781369/1‑71 (MQ=255)
                                     aCCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTg                                                      <  1:1797269/48‑1 (MQ=255)
                                     aCCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATAtt                               <  1:449797/71‑1 (MQ=255)
                                          gATAATGGCGTCGTAACCTGACAGTAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCa                          <  1:3876228/71‑1 (MQ=255)
                                          gATAATGGCGTCGTAACCTGACAGTAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCa                          <  1:3871501/71‑1 (MQ=255)
                                          gATAATGGCGTCGTAACCTGACAGTAGGCGCACGAAGTCGGGGAATCCCTGGGCGATGGCGATATTGCGCa                          <  1:711674/71‑1 (MQ=255)
                                          gATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTgcgc                           <  1:455958/70‑1 (MQ=255)
                                          gATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTgcgc                           <  1:740513/70‑1 (MQ=255)
                                              aTGGCGTCGTAACCTGACAGTAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCt                       <  1:3006064/70‑1 (MQ=255)
                                              aTGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCtt                      <  1:2040315/71‑1 (MQ=255)
                                                 gCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGc                    <  1:2219166/70‑1 (MQ=255)
                                                    tcgtAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAgt                <  1:956911/71‑1 (MQ=255)
                                                    tcgtAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAgt                <  1:3378969/71‑1 (MQ=255)
                                                      gtAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAgtgt              <  1:1407304/71‑1 (MQ=255)
                                                      gtAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAgtgt              <  1:1348372/71‑1 (MQ=255)
                                                       tAACCTGACAGTAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGc                        <  1:2970671/60‑1 (MQ=255)
                                                       tAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGggcgatggcga                                   <  1:3473643/49‑1 (MQ=255)
                                                       tAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGggcgatggcga                                   <  1:2697893/49‑1 (MQ=255)
                                                       tAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTgcgc                           >  1:1684958/1‑57 (MQ=255)
                                                       tAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAgtggc             <  1:2396289/71‑3 (MQ=255)
                                                       tAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTc             <  1:4083824/71‑1 (MQ=255)
                                                       tAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTc             <  1:2992438/71‑1 (MQ=255)
                                                       tAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTc             <  1:2940360/71‑1 (MQ=255)
                                                       tAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTc             <  1:2065268/71‑1 (MQ=255)
                                                          ccTGACAGTAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGt          <  1:1688736/71‑1 (MQ=255)
                                                          ccTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGt          >  1:2873947/1‑71 (MQ=255)
                                                              aCAGTAGGCGCACGAAGTCGATGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGAcg      >  1:1433214/1‑71 (MQ=255)
                                                              aCAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGAcg      >  1:3880485/1‑71 (MQ=255)
                                                              aCAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGAcg      >  1:3303571/1‑71 (MQ=255)
                                                              aCAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGAc       >  1:1693656/1‑70 (MQ=255)
                                                                  cagGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGa        >  1:3083663/1‑65 (MQ=255)
                                                                  cagGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTg  <  1:705557/71‑1 (MQ=255)
                                                                  |                                                                      
CGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTG  >  minE/1326603‑1326739

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