Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,949 C→T 25.2% A435A (GCG→GCA)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,9490CT25.2% 259.2 / 87.3 147A435A (GCG→GCA) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (52/58);  new base T (21/16);  total (73/74)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 152 total aligned reads displayed.

CTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGT  >  minE/608880‑609019
                                                                     |                                                                      
cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                                                                       >  1:2418760/1‑71 (MQ=255)
  gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                                                     >  1:2833316/1‑71 (MQ=255)
   cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAc                                                                    >  1:664257/1‑71 (MQ=255)
     ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAcc                                                                   >  1:3300909/1‑70 (MQ=255)
       ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                                  >  1:55847/1‑69 (MQ=255)
        tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                                  >  1:4919861/1‑68 (MQ=255)
         cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                               <  1:3217052/70‑1 (MQ=255)
         cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                               <  1:3515728/70‑1 (MQ=255)
         cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                               <  1:5057062/70‑1 (MQ=255)
         cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                               <  1:5081114/70‑1 (MQ=255)
          aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                             >  1:366141/1‑71 (MQ=255)
            tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACt                                                           <  1:3820837/71‑1 (MQ=255)
             ttctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTg                                                          >  1:5586913/4‑71 (MQ=255)
              tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                          >  1:3031028/1‑70 (MQ=255)
               ctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                         <  1:4795152/69‑1 (MQ=255)
               ctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                         <  1:5301396/69‑1 (MQ=255)
                tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                       >  1:1163239/1‑71 (MQ=255)
                tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                             >  1:4202879/1‑65 (MQ=255)
                tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                             >  1:3078734/1‑65 (MQ=255)
                tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                                        <  1:2517823/70‑1 (MQ=255)
                 cGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                                      >  1:462022/1‑71 (MQ=255)
                   tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTg                                                    <  1:5268998/71‑1 (MQ=255)
                     aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                  >  1:2669789/1‑71 (MQ=255)
                     aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                  <  1:898099/71‑1 (MQ=255)
                      aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                  >  1:3459613/1‑70 (MQ=255)
                       ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                                        <  1:2393552/61‑1 (MQ=255)
                        cccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                                               >  1:3756992/2‑71 (MQ=255)
                         ccgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAc                                               >  1:5557805/1‑70 (MQ=255)
                          cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                               <  1:4665827/53‑1 (MQ=255)
                          cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                               <  1:1526179/53‑1 (MQ=255)
                          cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCa                                             <  1:2065457/71‑1 (MQ=255)
                          cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTgct                                                  >  1:3426410/1‑64 (MQ=255)
                          cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTgct                                                  >  1:1081659/1‑64 (MQ=255)
                          cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                              >  1:5186119/1‑70 (MQ=255)
                          cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                              >  1:4221099/1‑70 (MQ=255)
                            ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaat                                           <  1:2375314/71‑2 (MQ=255)
                             cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                          <  1:2771617/71‑1 (MQ=255)
                             cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                          <  1:775996/71‑1 (MQ=255)
                               ttGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                        >  1:4718237/1‑71 (MQ=255)
                                taaTTGCTTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                                        <  1:438023/52‑1 (MQ=255)
                                 gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                       <  1:4320014/70‑1 (MQ=255)
                                 gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                       <  1:4810051/70‑1 (MQ=255)
                                 gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                       <  1:2818379/70‑1 (MQ=255)
                                  aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                      <  1:261160/70‑1 (MQ=255)
                                  aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                      <  1:3114830/70‑1 (MQ=255)
                                   ttGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                               <  1:3908913/44‑1 (MQ=12)
                                     gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                                                  >  1:4149446/1‑38 (MQ=9)
                                     gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                  >  1:3865652/1‑71 (MQ=255)
                                      cGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                             >  1:2067710/1‑43 (MQ=255)
                                      cGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                  <  1:493001/70‑1 (MQ=255)
                                          cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                                            >  1:136434/1‑40 (MQ=255)
                                          cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                                            >  1:5458941/1‑40 (MQ=255)
                                           aGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTAc                            <  1:1128710/71‑1 (MQ=255)
                                           aGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTAc                            <  1:196137/71‑1 (MQ=255)
                                             cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                             <  1:4089256/68‑1 (MQ=255)
                                             cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                           >  1:2583007/1‑70 (MQ=255)
                                               tGGTTTTCGGCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                           <  1:4912499/68‑1 (MQ=255)
                                               tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag                        >  1:814231/1‑71 (MQ=255)
                                                cgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                        >  1:1143248/2‑70 (MQ=255)
                                                 gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                     <  1:2207997/56‑1 (MQ=255)
                                                 gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                       <  1:1461504/70‑1 (MQ=255)
                                                 gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                       >  1:1569553/1‑70 (MQ=255)
                                                  ttttCGCCAGGTTTGCACGCGCACCGATGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                             <  1:4923260/63‑1 (MQ=255)
                                                  ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                      <  1:399597/54‑1 (MQ=255)
                                                  ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                      <  1:576990/54‑1 (MQ=255)
                                                  ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                      <  1:3517010/54‑1 (MQ=255)
                                                  ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCtt                              <  1:1943784/62‑1 (MQ=255)
                                                  ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaa                     <  1:5109958/71‑1 (MQ=255)
                                                   tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                    >  1:5591309/1‑71 (MQ=255)
                                                   tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                    >  1:5117703/1‑71 (MQ=255)
                                                    ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                    >  1:32789/1‑68 (MQ=255)
                                                    ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                    >  1:5207214/1‑68 (MQ=255)
                                                    ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                    >  1:5332640/1‑68 (MQ=255)
                                                    ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                    >  1:5184458/1‑68 (MQ=255)
                                                     tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                   <  1:3043363/70‑1 (MQ=255)
                                                     tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                   <  1:202592/70‑1 (MQ=255)
                                                     tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                   <  1:2032018/70‑1 (MQ=255)
                                                     tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                   <  1:5547664/70‑1 (MQ=255)
                                                     tCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                   <  1:4334380/70‑1 (MQ=255)
                                                       gCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                 <  1:4653418/70‑1 (MQ=255)
                                                         cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaatc                                          <  1:1905524/43‑3 (MQ=19)
                                                          aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc             >  1:1211034/1‑71 (MQ=255)
                                                          aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc             >  1:2480333/1‑71 (MQ=255)
                                                           ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                       <  1:4240566/44‑1 (MQ=255)
                                                           ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                   >  1:4642752/1‑64 (MQ=255)
                                                           ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg            <  1:224252/71‑1 (MQ=255)
                                                           ggTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaa                                            >  1:1476382/1‑39 (MQ=16)
                                                            gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                               <  1:4766030/51‑1 (MQ=255)
                                                            gTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt           <  1:2925223/71‑1 (MQ=255)
                                                             tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                           <  1:1089070/38‑1 (MQ=255)
                                                             tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                     <  1:1834487/44‑1 (MQ=255)
                                                                gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                       <  1:5612882/39‑1 (MQ=255)
                                                                gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                       <  1:2935523/39‑1 (MQ=255)
                                                                gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                  <  1:2511788/44‑1 (MQ=255)
                                                                gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                             <  1:3726161/49‑1 (MQ=255)
                                                                gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                           >  1:3165380/1‑51 (MQ=255)
                                                                 cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                   <  1:2138586/58‑1 (MQ=255)
                                                                  aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                     >  1:4914727/1‑39 (MQ=255)
                                                                  aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                     >  1:1657053/1‑39 (MQ=255)
                                                                  aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                      <  1:5496042/38‑1 (MQ=255)
                                                                     |                                                                      
CTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGT  >  minE/608880‑609019

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