Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,196 C→A 11.5% L686L (CTG→CTT pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,1960CA11.5% 300.3 / 15.6 95L686L (CTG→CTTpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (38/46);  new base A (3/8);  total (41/54)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 102 total aligned reads displayed.

GTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGAT  >  minE/608127‑608264
                                                                     |                                                                    
gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca                                                                     <  1:10364/71‑1 (MQ=255)
gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca                                                                     <  1:1102502/71‑1 (MQ=255)
gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca                                                                     <  1:866714/71‑1 (MQ=255)
gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca                                                                     <  1:1246005/71‑1 (MQ=255)
    aCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                                 <  1:3703610/71‑1 (MQ=255)
    aCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                                 >  1:1799271/1‑71 (MQ=255)
     cGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTg                                                                <  1:2821012/71‑1 (MQ=255)
     cGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTg                                                                <  1:3223304/71‑1 (MQ=255)
     cGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTg                                                                <  1:668188/71‑1 (MQ=255)
     cGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTg                                                                <  1:2482210/71‑1 (MQ=255)
      gTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                                 >  1:3239501/1‑69 (MQ=255)
        tCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTc                                                             >  1:1994647/1‑71 (MQ=255)
         gaGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGttt                                                            >  1:2508660/2‑69 (MQ=255)
         cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTc                                                             <  1:3159160/70‑1 (MQ=255)
         cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTc                                                             >  1:1792913/1‑70 (MQ=255)
         cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                            >  1:401254/1‑71 (MQ=255)
         cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                            >  1:2402537/1‑70 (MQ=255)
          aGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtttt                                                           <  1:942347/71‑1 (MQ=255)
          aGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtttt                                                           <  1:2231230/71‑1 (MQ=255)
          aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                            <  1:2804426/70‑1 (MQ=255)
          aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                            <  1:2700611/70‑1 (MQ=255)
          aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                            <  1:2689224/70‑1 (MQ=255)
          aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                            <  1:180736/70‑1 (MQ=255)
          aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                            <  1:2362760/70‑1 (MQ=255)
          aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                            <  1:2343644/70‑1 (MQ=255)
          aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAACCTGGCTGGTCt                                                            <  1:3119832/70‑1 (MQ=255)
           ggTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTa                                                          <  1:3409193/71‑1 (MQ=255)
           ggTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTa                                                          <  1:3824925/71‑1 (MQ=255)
           ggTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACAATCCATCAGACCAGCCAGGTTGGTCTTa                                                          <  1:3880655/71‑1 (MQ=255)
            gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                         <  1:1614315/71‑1 (MQ=255)
            gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                         <  1:1993034/71‑1 (MQ=255)
            gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                         <  1:1251432/71‑1 (MQ=255)
            gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                         >  1:923689/1‑71 (MQ=255)
              gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                                        >  1:2724014/1‑70 (MQ=255)
                  accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAAc                                                     <  1:3631156/69‑1 (MQ=255)
                  accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                                                    <  1:2984463/70‑1 (MQ=255)
                  accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                                                    >  1:1538295/1‑70 (MQ=255)
                   ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                         >  1:1479166/1‑64 (MQ=255)
                   ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                         >  1:2649502/1‑64 (MQ=255)
                   ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                  >  1:2907297/1‑71 (MQ=255)
                    caccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcg                                                 >  1:2952013/1‑71 (MQ=255)
                      ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtc                                               <  1:2284422/71‑1 (MQ=255)
                      ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtc                                               <  1:924198/71‑1 (MQ=255)
                         tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                            >  1:3731260/1‑71 (MQ=255)
                          cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAGGAACTTcgtcgt                                             <  1:2293151/69‑1 (MQ=255)
                          cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                             <  1:2823118/69‑1 (MQ=255)
                          cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                             <  1:474139/69‑1 (MQ=255)
                          cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                            <  1:1797071/70‑1 (MQ=255)
                          cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                            <  1:3724881/70‑1 (MQ=255)
                          cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt                                           >  1:3120132/1‑71 (MQ=255)
                          cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt                                           >  1:3396910/1‑71 (MQ=255)
                          cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt                                           >  1:1495472/1‑71 (MQ=255)
                           gatggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                                          >  1:3933899/1‑71 (MQ=255)
                             tggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                            <  1:2676588/51‑1 (MQ=255)
                               gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                         >  1:2492549/1‑52 (MQ=255)
                                atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc                                     <  1:1021903/71‑1 (MQ=255)
                                   cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                          <  1:334105/47‑1 (MQ=255)
                                   cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                          <  1:151775/47‑1 (MQ=255)
                                   cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                             <  1:463558/60‑1 (MQ=255)
                                   cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGt                                  <  1:3467922/71‑1 (MQ=255)
                                   cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGt                                  <  1:2497312/71‑1 (MQ=255)
                                   cTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgc                                  <  1:504485/71‑2 (MQ=255)
                                     tCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACtgg                                                  >  1:3358305/1‑51 (MQ=255)
                                     tCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACtgg                                                  >  1:2273165/1‑51 (MQ=255)
                                      cGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGaa                                                      >  1:2059806/1‑48 (MQ=255)
                                       gtggtgGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                              >  1:1672782/1‑71 (MQ=255)
                                        ttgtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                                                    <  1:824844/46‑1 (MQ=255)
                                          gtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTc                            <  1:546230/70‑1 (MQ=255)
                                           tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                          <  1:2630236/71‑1 (MQ=255)
                                            ggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                                          <  1:2843820/54‑1 (MQ=255)
                                            ggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                          >  1:677342/1‑70 (MQ=255)
                                             gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc                                     >  1:637500/1‑58 (MQ=255)
                                              aaGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAc                       >  1:571959/1‑71 (MQ=255)
                                               aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAc                       >  1:2924643/1‑70 (MQ=255)
                                                 tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                  >  1:3751140/1‑41 (MQ=255)
                                                 tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                  >  1:3342509/1‑41 (MQ=255)
                                                 tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                             >  1:622982/1‑46 (MQ=255)
                                                 tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                       >  1:1036054/1‑52 (MQ=255)
                                                  aaCCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                            >  1:856909/1‑46 (MQ=255)
                                                    ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCg                               <  1:1975745/57‑1 (MQ=255)
                                                    ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCg                               <  1:2832041/57‑1 (MQ=255)
                                                    ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa                   >  1:1611263/1‑69 (MQ=255)
                                                    ccatccatCAGACCACCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCg                               <  1:1263972/57‑1 (MQ=255)
                                                     catccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag                <  1:1141675/71‑1 (MQ=255)
                                                       tccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag                >  1:203586/1‑69 (MQ=255)
                                                         catcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa                   >  1:3905764/1‑64 (MQ=255)
                                                          atcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTg                                 <  1:126335/49‑1 (MQ=255)
                                                           tcAGACCAGCCAGGTTGGTCTTACGAACTtcgtcg                                              <  1:2829526/35‑1 (MQ=255)
                                                           tcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAGCGATAGAGAAGgtg          >  1:1607805/1‑71 (MQ=255)
                                                           tcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa                   <  1:1007353/62‑1 (MQ=255)
                                                           tcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg          <  1:1191281/71‑1 (MQ=255)
                                                            cAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTc                            >  1:2140048/1‑52 (MQ=255)
                                                            cAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg          >  1:3031071/1‑70 (MQ=255)
                                                             aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt         <  1:3074812/70‑1 (MQ=255)
                                                             aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg          >  1:151327/1‑69 (MQ=255)
                                                               aCCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAgg      <  1:435437/71‑1 (MQ=255)
                                                               aCCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAgg      <  1:2461678/71‑1 (MQ=255)
                                                                ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGtt                             >  1:1256042/1‑47 (MQ=255)
                                                                   gccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGgaga   >  1:1586559/1‑70 (MQ=255)
                                                                   gccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGAt  <  1:1203182/71‑1 (MQ=255)
                                                                     |                                                                    
GTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGAT  >  minE/608127‑608264

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: