Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,946 A→G 25.2% R436R (CGT→CGC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,9460AG25.2% 250.2 / 88.9 147R436R (CGT→CGC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (50/60);  new base G (21/16);  total (71/76)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

Only 100 of 156 total aligned reads displayed.

CAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGT  >  minE/608876‑609010
                                                                      |                                                                
cAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                                  >  1:2086234/1‑71 (MQ=255)
  aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                                <  1:1813025/71‑1 (MQ=255)
  aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                                <  1:2942887/71‑1 (MQ=255)
  aCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                                                                <  1:62444/71‑1 (MQ=255)
    cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                                                              <  1:3899536/71‑1 (MQ=255)
      gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                                            >  1:2833316/1‑71 (MQ=255)
       cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAc                                                           >  1:664257/1‑71 (MQ=255)
         ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAcc                                                          >  1:3300909/1‑70 (MQ=255)
         ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAcc                                                          >  1:3640368/1‑70 (MQ=255)
            tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                         >  1:4919861/1‑68 (MQ=255)
             cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcgtg                                                              >  1:5664726/1‑60 (MQ=255)
             cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                      <  1:5057062/70‑1 (MQ=255)
              aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                    <  1:2485728/71‑1 (MQ=255)
                tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACt                                                  <  1:3820837/71‑1 (MQ=255)
                 ttctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTg                                                 >  1:5586913/4‑71 (MQ=255)
                  tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                 >  1:3031028/1‑70 (MQ=255)
                  tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                 >  1:2813492/1‑70 (MQ=255)
                   ctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                <  1:4795152/69‑1 (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)
                     cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                              <  1:3273245/70‑1 (MQ=255)
                         aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                         >  1:2402722/1‑71 (MQ=255)
                         aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                         <  1:898099/71‑1 (MQ=255)
                          aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACTCGCACCGAAGAACTGCATTTGtt                                         <  1:2118390/70‑1 (MQ=255)
                          aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                         >  1:3459613/1‑70 (MQ=255)
                           ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                               <  1:1155144/61‑1 (MQ=255)
                            cccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                                      >  1:3756992/2‑71 (MQ=255)
                             ccgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAc                                      >  1:5557805/1‑70 (MQ=255)
                              cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                      <  1:1526179/53‑1 (MQ=255)
                              cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                      <  1:4665827/53‑1 (MQ=255)
                              cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCa                                    <  1:2065457/71‑1 (MQ=255)
                              cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTgct                                         >  1:1081659/1‑64 (MQ=255)
                                ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                       >  1:5516307/1‑66 (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)
                                    tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaa                                                       >  1:3720903/1‑46 (MQ=255)
                                     gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                              <  1:4320014/70‑1 (MQ=255)
                                     aaTTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                              >  1:4177436/2‑54 (MQ=255)
                                      aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                             <  1:261160/70‑1 (MQ=255)
                                       ttGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                      <  1:3908913/44‑1 (MQ=12)
                                        tGCGTACAGCATGGTTTTCCCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTa                                       >  1:3366110/1‑58 (MQ=255)
                                         gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                                         >  1:1466690/1‑38 (MQ=9)
                                         gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaagaa                                                    >  1:4326718/1‑44 (MQ=12)
                                         gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                         >  1:5340523/1‑71 (MQ=255)
                                            tACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                      <  1:4018357/71‑1 (MQ=255)
                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                                   >  1:3498844/1‑40 (MQ=255)
                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                    >  1:5572364/1‑71 (MQ=255)
                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                    >  1:2240311/1‑71 (MQ=255)
                                                 cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                    <  1:5202882/68‑1 (MQ=255)
                                                 cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                    <  1:4089256/68‑1 (MQ=255)
                                                  aTGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                    >  1:3132667/1‑67 (MQ=255)
                                                   tGGTTTTCGGCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                  <  1:4912499/68‑1 (MQ=255)
                                                   tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag               >  1:814231/1‑71 (MQ=255)
                                                   tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACAATCATCGGGCTTACg                  <  1:508464/68‑1 (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:576990/54‑1 (MQ=255)
                                                      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                             <  1:3517010/54‑1 (MQ=255)
                                                      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCtt                     <  1:1943784/62‑1 (MQ=255)
                                                      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc              >  1:4480111/1‑69 (MQ=255)
                                                      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaa            <  1:4812571/71‑1 (MQ=255)
                                                      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaa            <  1:5109958/71‑1 (MQ=255)
                                                       tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                              >  1:5134384/1‑52 (MQ=255)
                                                       tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag           >  1:3055222/1‑71 (MQ=255)
                                                        ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag           >  1:2391923/1‑70 (MQ=255)
                                                        ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag           >  1:2663018/1‑70 (MQ=255)
                                                        ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg           >  1:5332640/1‑68 (MQ=255)
                                                        ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg           >  1:5207214/1‑68 (MQ=255)
                                                        ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg           >  1:5184458/1‑68 (MQ=255)
                                                         tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                      <  1:5587630/58‑1 (MQ=255)
                                                         tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc          <  1:5274704/70‑1 (MQ=255)
                                                         tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc          <  1:202592/70‑1 (MQ=255)
                                                         tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc          <  1:2032018/70‑1 (MQ=255)
                                                           gCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa        <  1:4653418/70‑1 (MQ=255)
                                                           gCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa        <  1:4224293/70‑1 (MQ=255)
                                                             cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaatc                                 <  1:1912273/43‑3 (MQ=19)
                                                              aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc    >  1:4947648/1‑71 (MQ=255)
                                                              aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc    >  1:2480333/1‑71 (MQ=255)
                                                              aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc    >  1:1211034/1‑71 (MQ=255)
                                                               ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                              >  1:1341037/1‑44 (MQ=255)
                                                               ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc          >  1:4642752/1‑64 (MQ=255)
                                                                gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                      <  1:4772410/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:4919366/44‑1 (MQ=255)
                                                                    gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                  >  1:3165380/1‑51 (MQ=255)
                                                                      aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                            >  1:1657053/1‑39 (MQ=255)
                                                                      aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                             <  1:5496042/38‑1 (MQ=255)
                                                                      aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                             <  1:252889/38‑1 (MQ=255)
                                                                      |                                                                
CAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGT  >  minE/608876‑609010

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