Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,655 A→C 7.0% A533A (GCT→GCG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,6550AC 7.0% 218.4 / 8.0 85A533A (GCT→GCG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (41/38);  new base C (2/4);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCA  >  minE/608586‑608722
                                                                     |                                                                   
gATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagcag                                                                    >  1:2594567/1‑71 (MQ=255)
gATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagca                                                                     <  1:1947354/70‑1 (MQ=255)
 aTAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCagcagc                                                                   <  1:4490458/71‑1 (MQ=255)
  tAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagcag                                                                    >  1:2958687/1‑69 (MQ=255)
  tAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCagcagc                                                                   >  1:3295102/1‑70 (MQ=255)
  tAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCagcagc                                                                   >  1:4181604/1‑70 (MQ=255)
     ccAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACg                                                               <  1:827431/71‑1 (MQ=255)
      cAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCaa                                                                 >  1:4241703/1‑68 (MQ=255)
       aGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGa                                                             >  1:4964313/1‑71 (MQ=255)
        gACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGa                                                             >  1:2606141/1‑70 (MQ=255)
          ccGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAg                                                          <  1:3184799/71‑1 (MQ=255)
           cGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagca                                                                     >  1:798317/1‑59 (MQ=255)
           cGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACa                                                           >  1:5132053/1‑69 (MQ=255)
           cGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGa                                                         <  1:4881319/71‑1 (MQ=255)
             tCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCagcagc                                                                   <  1:2299178/59‑1 (MQ=255)
               ttCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAg                                                     >  1:2221473/1‑71 (MQ=255)
               ttCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAg                                                     >  1:5559560/1‑71 (MQ=255)
                 cGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                    <  1:525254/70‑1 (MQ=255)
                 cGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                    <  1:5559538/70‑1 (MQ=255)
                 cGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                    <  1:829012/70‑1 (MQ=255)
                  gTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGc                                                    <  1:5654942/69‑1 (MQ=255)
                   tCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGAt                                                 <  1:4089282/71‑1 (MQ=255)
                     aCGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATaa                                               >  1:4233397/1‑70 (MQ=255)
                     aCGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATAc                                               >  1:1333560/1‑71 (MQ=255)
                       gAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACg                                                               <  1:914416/53‑1 (MQ=255)
                       gAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACca                                             <  1:2480866/71‑1 (MQ=255)
                       gAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACca                                             <  1:463864/71‑1 (MQ=255)
                       gAATCGATTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGAc                                                            <  1:4183711/56‑1 (MQ=255)
                         aTCGGTTTAACTTTCGCATATTTTATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACcaca                                           <  1:3384186/71‑1 (MQ=255)
                         aTCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACcaca                                           <  1:2769197/71‑1 (MQ=255)
                          tCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCacac                                          >  1:5277265/1‑71 (MQ=255)
                             gTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGcc                                       <  1:1521978/71‑1 (MQ=255)
                             gTTAAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGc                                        <  1:3387994/70‑1 (MQ=255)
                               ttAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCAt                                     >  1:5109093/1‑71 (MQ=255)
                                tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                    >  1:3085328/1‑71 (MQ=255)
                                tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                    >  1:5075723/1‑71 (MQ=255)
                                tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                    >  1:115620/1‑71 (MQ=255)
                                 aaCTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                    >  1:2842431/1‑70 (MQ=255)
                                   cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                                      <  1:4223321/50‑1 (MQ=255)
                                   cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                                      <  1:38620/50‑1 (MQ=255)
                                   cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                                      <  1:4423978/50‑1 (MQ=255)
                                   cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                                      <  1:758341/50‑1 (MQ=255)
                                   cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTg                                 <  1:117589/71‑1 (MQ=255)
                                    tttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                                >  1:2713393/1‑71 (MQ=255)
                                     ttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTg                                 >  1:5359604/1‑69 (MQ=255)
                                     ttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCg                               >  1:4107728/1‑71 (MQ=255)
                                      tCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                                <  1:2585078/69‑1 (MQ=255)
                                      tCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                                <  1:2390300/69‑1 (MQ=255)
                                      gggCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCga                              <  1:373464/69‑2 (MQ=255)
                                       cGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATgg                                   >  1:1881742/1‑65 (MQ=255)
                                       cGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATgg                                   >  1:1824049/1‑65 (MQ=255)
                                        gCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGAt                            <  1:937516/71‑1 (MQ=255)
                                         cttatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATa                           <  1:2431491/69‑1 (MQ=255)
                                         cATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATa                           <  1:920176/71‑1 (MQ=255)
                                         cATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGa                                                  >  1:1089330/1‑48 (MQ=255)
                                           tatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACa                                                           <  1:4210037/37‑1 (MQ=255)
                                           tatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACa                                                           <  1:3466546/37‑1 (MQ=255)
                                           tatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGAt                                                 <  1:1127564/47‑1 (MQ=255)
                                           tatTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGcccg                                                     <  1:1613128/43‑3 (MQ=12)
                                            atTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                    <  1:574690/43‑1 (MQ=255)
                                            atTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACg                        >  1:1149481/1‑71 (MQ=255)
                                             tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAg                                                          >  1:3417255/1‑36 (MQ=255)
                                             tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAg                                                          >  1:1767876/1‑36 (MQ=255)
                                             tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGa                                                  <  1:218253/44‑1 (MQ=255)
                                             tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                       >  1:4624953/1‑71 (MQ=255)
                                             tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                       >  1:4120457/1‑71 (MQ=255)
                                             tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                       >  1:3464094/1‑71 (MQ=255)
                                             tttGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGc                                                    >  1:1975401/1‑42 (MQ=255)
                                              ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                       >  1:4482926/1‑70 (MQ=255)
                                              ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                       >  1:3549713/1‑70 (MQ=255)
                                               tGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTc                      <  1:4495706/70‑1 (MQ=255)
                                               tGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACACCATGGTGCGGATAACGTCa                     >  1:1826472/1‑71 (MQ=255)
                                                  ttGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACgg                  <  1:709857/71‑1 (MQ=255)
                                                   tGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGt                 >  1:3881570/1‑71 (MQ=255)
                                                    gCAGACAGGGAGTCAGCAGCGACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                >  1:1409918/1‑71 (MQ=255)
                                                    gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATaa                          >  1:5666001/1‑61 (MQ=255)
                                                    gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATaa                          >  1:2960014/1‑61 (MQ=255)
                                                    gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                >  1:5365467/1‑71 (MQ=255)
                                                    gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                >  1:4431275/1‑71 (MQ=255)
                                                    gCAGACAGGGAGTCAGCAGCAACGGACAGAACAGCGATACCACACGCCATGGTGCGGATaa                          >  1:898897/1‑61 (MQ=255)
                                                     cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                <  1:70206/70‑1 (MQ=255)
                                                     cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCg               >  1:3893114/1‑71 (MQ=255)
                                                       gACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCgcg             <  1:994016/71‑3 (MQ=255)
                                                       gACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTg             <  1:2664318/71‑1 (MQ=255)
                                                        aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTg             >  1:4380746/1‑70 (MQ=255)
                                                        aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGc            <  1:259841/71‑1 (MQ=255)
                                                         cAGGGAGTCAGCAGCACCGGACAGACCAGCGATACcaca                                           >  1:4979632/1‑39 (MQ=255)
                                                         cAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGc            <  1:5063168/70‑1 (MQ=255)
                                                         cAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCa           <  1:2870297/71‑1 (MQ=255)
                                                          aGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAg          >  1:4929491/1‑71 (MQ=255)
                                                            ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCa           >  1:767944/1‑68 (MQ=255)
                                                            ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCa           >  1:4325930/1‑68 (MQ=255)
                                                                tCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt    <  1:5488704/71‑1 (MQ=255)
                                                                  agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAg          >  1:4966893/1‑63 (MQ=255)
                                                                  agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt    >  1:2339797/1‑69 (MQ=255)
                                                                  agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt    >  1:3859894/1‑69 (MQ=255)
                                                                  agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCa  <  1:5218487/71‑1 (MQ=255)
                                                                     |                                                                   
GATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCA  >  minE/608586‑608722

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