Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,216 C→G 7.0% E680Q (GAA→CAA) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,2160CG 7.0% 261.9 / 6.4 84E680Q (GAA→CAA) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (36/42);  new base G (2/4);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAC  >  minE/608146‑608286
                                                                      |                                                                      
ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                                        <  1:599943/71‑1 (MQ=255)
ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                                        >  1:831423/1‑71 (MQ=255)
ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACAATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                                        <  1:1910382/71‑1 (MQ=255)
ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACAATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                                        <  1:1914503/71‑1 (MQ=255)
 caccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcg                                                                       <  1:2452448/71‑1 (MQ=255)
    cTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtcg                                                                    >  1:1584842/1‑71 (MQ=255)
       cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                                                   <  1:3281764/69‑1 (MQ=255)
       cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                                                   <  1:158074/69‑1 (MQ=255)
       cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt                                                                 >  1:1257412/1‑71 (MQ=255)
        gatggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                                                                <  1:2778301/71‑1 (MQ=255)
         atggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                                        <  1:547807/62‑1 (MQ=255)
           ggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                                             >  1:1084660/1‑71 (MQ=255)
           ggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                                             >  1:978465/1‑71 (MQ=255)
           ggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                                             >  1:1056648/1‑71 (MQ=255)
            gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTAGTCTTTAcc                                                            <  1:1333530/71‑1 (MQ=255)
             atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc                                                           <  1:1049229/71‑1 (MQ=255)
             atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc                                                           <  1:3641242/71‑1 (MQ=255)
                cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGt                                                        <  1:2006578/71‑1 (MQ=255)
                 ttCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgcg                                                       >  1:3317690/1‑69 (MQ=255)
                 ttCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgcg                                                       >  1:555346/1‑69 (MQ=255)
                     tggtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCa                                                          >  1:1810574/1‑64 (MQ=255)
                     tggtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCa                                                          >  1:2123364/1‑64 (MQ=255)
                     ctgtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                                                    <  1:1079716/68‑1 (MQ=255)
                        tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                                                    <  1:1732766/67‑1 (MQ=255)
                        tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                                                    <  1:1129711/67‑1 (MQ=255)
                        tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                                                    <  1:2592002/67‑1 (MQ=255)
                        tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                                                <  1:1470584/71‑1 (MQ=255)
                          gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                                        <  1:2564473/45‑1 (MQ=255)
                          gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                                                <  1:81576/69‑1 (MQ=255)
                            aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACGCAGTGCGTTCGGAAc                                             >  1:3145403/1‑70 (MQ=255)
                            aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                                                <  1:3243649/67‑1 (MQ=255)
                            aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                                                <  1:2212706/67‑1 (MQ=255)
                                aCCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                                                  <  1:2868811/45‑1 (MQ=255)
                                  catccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAc                                             <  1:2684017/64‑1 (MQ=255)
                                   atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag                                      >  1:2396999/1‑70 (MQ=255)
                                   atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag                                      >  1:2197010/1‑70 (MQ=255)
                                   atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga                                     <  1:983901/71‑1 (MQ=255)
                                   atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga                                     <  1:1198443/71‑1 (MQ=255)
                                   atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga                                     <  1:2909350/71‑1 (MQ=255)
                                     ccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAg                                   >  1:2451130/1‑71 (MQ=255)
                                      catcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAgg                                  >  1:3194138/1‑71 (MQ=255)
                                         cAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg                                >  1:2079794/1‑70 (MQ=255)
                                          aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt                               <  1:2219934/70‑1 (MQ=255)
                                          aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt                               <  1:541844/70‑1 (MQ=255)
                                          aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt                               <  1:1355717/70‑1 (MQ=255)
                                           gACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAg                             <  1:189502/71‑1 (MQ=255)
                                             ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                                                 >  1:179306/1‑49 (MQ=255)
                                             ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                                                 >  1:2560959/1‑49 (MQ=255)
                                             ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                                                 >  1:2847790/1‑49 (MQ=255)
                                             ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAg                             >  1:470641/1‑69 (MQ=255)
                                              cagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGGACGATAGAGAAGGTGTAg                             >  1:341976/1‑68 (MQ=255)
                                               agccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAg                             >  1:3269415/1‑67 (MQ=255)
                                               agccagGTTGGTCTTACGAACTTCGTCGTCTATACCCAGTGCGTTCGGAACGATAGAGAAGGTGTa                              >  1:2142812/1‑66 (MQ=255)
                                                gccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGgaga                         >  1:2956301/1‑70 (MQ=255)
                                                 ccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag                                      <  1:2029097/56‑1 (MQ=255)
                                                 ccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATa                       <  1:3022282/71‑1 (MQ=255)
                                                  cagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAc                      >  1:2464540/1‑71 (MQ=255)
                                                   agGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGgaga                         <  1:1493229/67‑1 (MQ=255)
                                                   agGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAcc                     <  1:601451/71‑1 (MQ=255)
                                                   agGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAcc                     <  1:1577235/71‑1 (MQ=255)
                                                    ggTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCa                    >  1:1274257/1‑71 (MQ=255)
                                                    ggTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCa                    <  1:2778570/71‑1 (MQ=255)
                                                     gTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCAt                   <  1:2998289/71‑1 (MQ=255)
                                                     gTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCAt                   <  1:632466/71‑1 (MQ=255)
                                                      ttGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGATTGAGAAGGTGTACGAGATCCCATc                  <  1:672893/71‑1 (MQ=255)
                                                      ttGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGATTGAGAAGGTGTACGAGATCCCATc                  <  1:3163479/71‑1 (MQ=255)
                                                      ttGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGATTGAGAAGGTGTACGAGATCCCATc                  <  1:2991555/71‑1 (MQ=255)
                                                      ttGGTCTTACGAACTTGGTCGTCTTTACCCAGTGCTTTCGGAACGATAGAGAACGTGTAGGAGATACCa                    <  1:3379116/69‑1 (MQ=255)
                                                      ttGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATc                  >  1:3214630/1‑71 (MQ=255)
                                                      ttGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATc                  >  1:1002129/1‑71 (MQ=255)
                                                       tGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCt                 >  1:3175971/1‑71 (MQ=255)
                                                        ggTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCt                 <  1:3048909/70‑1 (MQ=255)
                                                         gTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCtt                <  1:3229215/70‑1 (MQ=255)
                                                          tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATaga                                       <  1:225981/46‑1 (MQ=255)
                                                          tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCtt                <  1:2306891/69‑1 (MQ=255)
                                                          tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCt                 <  1:343722/68‑1 (MQ=255)
                                                           cTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATag                                        <  1:989949/44‑1 (MQ=255)
                                                           cTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATag                                        <  1:1673309/44‑1 (MQ=255)
                                                           cTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATag                                        <  1:2567041/44‑1 (MQ=255)
                                                           cTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATag                                        <  1:3274462/44‑1 (MQ=255)
                                                            ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTa                              >  1:3565558/1‑53 (MQ=255)
                                                            ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTa                              >  1:200017/1‑53 (MQ=255)
                                                            ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGga                           >  1:215671/1‑56 (MQ=255)
                                                            ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGc            >  1:2327995/1‑71 (MQ=255)
                                                              aCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg           >  1:3613149/1‑70 (MQ=255)
                                                               cGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt          >  1:1394865/1‑70 (MQ=255)
                                                               cGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt          <  1:3484541/70‑1 (MQ=255)
                                                               cGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGt          >  1:2259443/1‑70 (MQ=255)
                                                               cGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTa         <  1:543070/71‑1 (MQ=255)
                                                                  aCTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAg       >  1:613120/1‑70 (MQ=255)
                                                                  aCTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAg       >  1:719779/1‑70 (MQ=255)
                                                                  aCTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAg       >  1:3429852/1‑70 (MQ=255)
                                                                  aCTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAg       >  1:913473/1‑70 (MQ=255)
                                                                  aCTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGc      <  1:3072800/71‑1 (MQ=255)
                                                                      cgtcgtCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAc  >  1:1681878/1‑71 (MQ=255)
                                                                      |                                                                      
CCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAC  >  minE/608146‑608286

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: ‑