Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,178 A→C 15.7% G692G (GGT→GGG pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,1780AC15.7% 223.9 / 46.0 115G692G (GGT→GGGpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (62/35);  new base C (10/8);  total (72/43)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-01

Only 100 of 117 total aligned reads displayed.

CATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATA  >  minE/608109‑608247
                                                                     |                                                                     
caTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                                                                      <  1:3379660/71‑1 (MQ=255)
caTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                                                                      >  1:1334170/1‑71 (MQ=255)
  tttCACGGTTCATCACGTTAACTTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                                                    >  1:4915494/1‑71 (MQ=255)
  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                                                    >  1:2504342/1‑71 (MQ=255)
  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc                                                                     >  1:4281342/1‑70 (MQ=255)
  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAaca                                                                    >  1:3366388/1‑71 (MQ=255)
    tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccat                                                                   <  1:700995/70‑1 (MQ=255)
      aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc                                                                 <  1:2767935/70‑1 (MQ=255)
       cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                                               <  1:4829242/71‑1 (MQ=255)
       cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                                               <  1:5077221/71‑1 (MQ=255)
       cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                                               <  1:3110125/71‑1 (MQ=255)
        ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                                               >  1:5401046/1‑70 (MQ=255)
        ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                                               >  1:3961175/1‑70 (MQ=255)
        ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                                               >  1:5186034/1‑70 (MQ=255)
            cattaCGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCcagcaggc                                                          >  1:2194135/1‑69 (MQ=255)
            catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                                                          >  1:4958019/1‑71 (MQ=255)
            catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                                                          >  1:5157747/1‑71 (MQ=255)
             attacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGcc                                                         <  1:4386120/68‑1 (MQ=9)
              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                                             <  1:2189330/64‑1 (MQ=9)
              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                                             <  1:5129726/64‑1 (MQ=9)
              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg                                                        <  1:3477946/69‑2 (MQ=255)
              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg                                                        <  1:410356/69‑2 (MQ=255)
              tcaCGTTAACGTTCAGGTCCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAg                                                            <  1:244548/67‑1 (MQ=255)
                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                                                                >  1:1947077/1‑61 (MQ=255)
                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                                                                >  1:3737946/1‑61 (MQ=255)
                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                                                       <  1:1231279/70‑1 (MQ=255)
                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                                                       <  1:348123/70‑1 (MQ=255)
                 cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc                                                      >  1:1497237/1‑70 (MQ=255)
                 cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccaga                                                      >  1:3046544/1‑69 (MQ=255)
                   ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagccag                                                   >  1:5154923/1‑71 (MQ=255)
                     gaCGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                                             >  1:5356049/2‑59 (MQ=9)
                     gaCGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                                             >  1:847201/2‑59 (MQ=9)
                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATcagcccagcc                                                     >  1:1906388/1‑67 (MQ=255)
                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAcc                                                         >  1:1661064/1‑63 (MQ=255)
                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAgg                                                  <  1:2450950/70‑1 (MQ=255)
                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                                                 >  1:4374389/1‑71 (MQ=255)
                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                                                 >  1:3673025/1‑71 (MQ=255)
                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                                                 >  1:3487062/1‑71 (MQ=255)
                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                                                 >  1:5228483/1‑71 (MQ=255)
                       cGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGATCAGCCAGGTTg                                               >  1:1063503/1‑71 (MQ=255)
                        gTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTgg                                              >  1:1924003/1‑71 (MQ=255)
                        gTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTgg                                              <  1:4275958/71‑1 (MQ=255)
                           cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                           >  1:4177110/1‑71 (MQ=255)
                             ggTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                         >  1:2654879/1‑71 (MQ=255)
                              gtgtTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACccatccag                                                               <  1:5696337/48‑2 (MQ=255)
                              gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                         >  1:2324261/1‑70 (MQ=255)
                              gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                        <  1:1212158/71‑1 (MQ=255)
                               tGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                                                    >  1:4328242/1‑42 (MQ=255)
                                gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                       >  1:2530585/1‑70 (MQ=255)
                                gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                       >  1:1493020/1‑70 (MQ=255)
                                 cTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGaa                                     <  1:12246/71‑1 (MQ=255)
                                  tGACCGCCTTCGACATACGCTTCGTGGTGGAAATACACATCCAGCAGGCCGACAAGGTTGGTTTTACGTAc                                    >  1:1881386/1‑71 (MQ=255)
                                  tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                                                          >  1:1664053/1‑49 (MQ=255)
                                  tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                         <  1:4780395/66‑1 (MQ=255)
                                   gACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                       <  1:4104449/67‑1 (MQ=255)
                                    accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                       >  1:935576/1‑66 (MQ=255)
                                    accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                       >  1:3098408/1‑66 (MQ=255)
                                    accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                                   <  1:1572301/70‑1 (MQ=255)
                                    accaccTACGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACtt                                  >  1:4563970/1‑71 (MQ=255)
                                        ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGa                                      >  1:1552807/1‑63 (MQ=255)
                                        ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCAAGGTTGGTCTTACGa                                      >  1:3175673/1‑63 (MQ=255)
                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                >  1:5594752/1‑50 (MQ=255)
                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                >  1:5352547/1‑50 (MQ=255)
                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                >  1:2331445/1‑50 (MQ=255)
                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                >  1:2413756/1‑50 (MQ=255)
                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                >  1:4045519/1‑50 (MQ=255)
                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                >  1:4920653/1‑50 (MQ=255)
                                            cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                            <  1:3116944/69‑1 (MQ=255)
                                            cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                            <  1:4803447/69‑1 (MQ=255)
                                            cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt                          >  1:4551955/1‑71 (MQ=255)
                                            cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt                          >  1:330922/1‑71 (MQ=255)
                                              atggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt                        <  1:1041323/71‑1 (MQ=255)
                                              atggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                         >  1:1830649/1‑70 (MQ=255)
                                               tggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt                        <  1:5088166/70‑1 (MQ=255)
                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGt                                             <  1:571495/47‑1 (MQ=255)
                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc                     >  1:756104/1‑71 (MQ=255)
                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc                     >  1:4544237/1‑71 (MQ=255)
                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc                     <  1:3715514/71‑1 (MQ=255)
                                                  atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc                    <  1:260989/71‑1 (MQ=255)
                                                  atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc                    <  1:3692112/71‑1 (MQ=255)
                                                   tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                           >  1:3721636/1‑47 (MQ=255)
                                                   tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                           >  1:4465836/1‑47 (MQ=255)
                                                   tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                           >  1:4405462/1‑47 (MQ=255)
                                                   tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                           >  1:3246150/1‑47 (MQ=255)
                                                   tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc                     >  1:64198/1‑69 (MQ=255)
                                                    gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                           >  1:1995157/1‑62 (MQ=255)
                                                    gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                           >  1:774782/1‑62 (MQ=255)
                                                    gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                           >  1:2977408/1‑62 (MQ=255)
                                                     cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtc                              >  1:251112/1‑58 (MQ=255)
                                                     cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGt                 >  1:5692686/1‑71 (MQ=255)
                                                      ttCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt                        <  1:3583775/63‑1 (MQ=255)
                                                      ttCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt                        <  1:725121/63‑1 (MQ=255)
                                                       tCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgcg                >  1:4196184/1‑68 (MQ=255)
                                                        cGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                       <  1:460381/46‑1 (MQ=255)
                                                        cGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                       <  1:2983188/46‑1 (MQ=255)
                                                            gtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg          >  1:5161950/1‑71 (MQ=255)
                                                             tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg         <  1:2945434/71‑1 (MQ=255)
                                                                  gTAACCATCCATCAGACCAGCCAGGTTGGTCTTTCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGa    <  1:4596744/71‑1 (MQ=255)
                                                                     aCCATCTATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa  <  1:5236490/70‑1 (MQ=255)
                                                                     aCCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGa    >  1:1746875/1‑68 (MQ=255)
                                                                     |                                                                     
CATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATA  >  minE/608109‑608247

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: