Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,265 A→C 12.6% G663G (GGT→GGG pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,2650AC12.6% 247.4 / 31.2 111G663G (GGT→GGGpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (51/46);  new base C (10/4);  total (61/50)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

AGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAC  >  minE/608197‑608335
                                                                    |                                                                      
agGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAcc                                                                      <  1:5183340/71‑1 (MQ=255)
 ggTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCa                                                                     >  1:924811/1‑71 (MQ=255)
 ggTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCa                                                                     >  1:4216897/1‑71 (MQ=255)
  gTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCAt                                                                    >  1:266083/1‑71 (MQ=255)
       tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTa                                                               <  1:199760/71‑1 (MQ=255)
         ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAg                                                              >  1:1325124/1‑70 (MQ=255)
         ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAg                                                              >  1:1839091/1‑70 (MQ=255)
           aCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg                                                            >  1:1604319/1‑70 (MQ=255)
           aCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg                                                            >  1:4504611/1‑70 (MQ=255)
             gAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAg                                                              <  1:4489631/66‑1 (MQ=255)
               aCTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAg                                                        >  1:4495514/1‑70 (MQ=255)
                cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                      <  1:3063704/71‑1 (MQ=255)
                cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                      <  1:3001695/71‑1 (MQ=255)
                cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                      <  1:2603901/71‑1 (MQ=255)
                cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                      <  1:3577571/71‑1 (MQ=255)
                cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                      <  1:4424500/71‑1 (MQ=255)
                cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                      <  1:601071/71‑1 (MQ=255)
                   cgtcgtCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                    <  1:747403/70‑1 (MQ=255)
                     tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATc                                                                   <  1:5622563/53‑1 (MQ=255)
                     tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACgg                                                 >  1:1411211/1‑71 (MQ=255)
                     tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTATAGCGTAAGCAAACgg                                                 <  1:1381964/71‑1 (MQ=255)
                       gtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCtt                                                                 >  1:4431174/1‑53 (MQ=255)
                       gtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCtt                                                                 >  1:2525290/1‑53 (MQ=255)
                       gtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCtt                                                                 >  1:2512603/1‑53 (MQ=255)
                       gtcTTTACCCAGTGCGTTGGAACGATAGAGAAGGTGTAGGAGATACCATCtt                                                                 >  1:1218642/1‑52 (MQ=255)
                            tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                    >  1:2353965/1‑61 (MQ=255)
                            tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                    >  1:2424067/1‑61 (MQ=255)
                            tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                            <  1:350380/69‑1 (MQ=255)
                             aCCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGt                                             <  1:1894733/67‑1 (MQ=255)
                             aCCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                          <  1:4161336/70‑1 (MQ=255)
                             aCCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                          <  1:5430160/70‑1 (MQ=255)
                               ccAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCt                                                                  <  1:4271542/44‑1 (MQ=255)
                               ccAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGc                                                <  1:4068465/62‑1 (MQ=255)
                                 aGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCttt                                                                >  1:3229642/1‑44 (MQ=255)
                                 aGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCttt                                                                >  1:1822003/1‑44 (MQ=255)
                                 aGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCttt                                                                >  1:2203233/1‑44 (MQ=255)
                                  gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAcc                                                                      <  1:291027/37‑1 (MQ=255)
                                  gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAcc                                                                      <  1:2119633/37‑1 (MQ=255)
                                  gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCAt                                                                    <  1:2279280/39‑1 (MQ=255)
                                   tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                    >  1:2595823/1‑70 (MQ=255)
                                   tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                    >  1:3162435/1‑70 (MQ=255)
                                     cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                   >  1:726501/1‑69 (MQ=255)
                                       ttCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTa                                                          <  1:1140240/44‑1 (MQ=255)
                                          ggAACGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                             >  1:3755348/1‑70 (MQ=255)
                                          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                            <  1:1390292/71‑1 (MQ=255)
                                          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                            <  1:5653632/71‑1 (MQ=255)
                                           gAACGATTGAGAAGGTGTACGAGATCCCATCTTTGGCg                                                            <  1:2045071/38‑1 (MQ=21)
                                            aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                             >  1:4607131/1‑68 (MQ=255)
                                             aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCa                                               >  1:3865439/1‑49 (MQ=255)
                                             aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAg                                              >  1:256866/1‑50 (MQ=255)
                                             aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                           >  1:513774/1‑67 (MQ=255)
                                             aCGATTGAGAAGGTGTACGAGATCCCATATTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                           >  1:4234740/1‑67 (MQ=255)
                                              cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                              >  1:2378580/1‑49 (MQ=255)
                                              cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                              >  1:4943854/1‑49 (MQ=255)
                                              cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                              >  1:4597779/1‑49 (MQ=255)
                                               gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                 <  1:5255274/61‑1 (MQ=255)
                                               gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                 <  1:4060079/61‑1 (MQ=255)
                                                aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                      <  1:1327759/71‑1 (MQ=255)
                                                 ttgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                      >  1:5678493/3‑70 (MQ=255)
                                                 tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                      <  1:2004201/38‑1 (MQ=255)
                                                 tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                      >  1:1762916/1‑70 (MQ=255)
                                                 tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                      <  1:4135293/70‑1 (MQ=255)
                                                 tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                      >  1:3348759/1‑70 (MQ=255)
                                                  agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                    >  1:4378990/1‑71 (MQ=255)
                                                  agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                    >  1:5573347/1‑71 (MQ=255)
                                                  agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                    >  1:2969125/1‑71 (MQ=255)
                                                   gagaAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                    >  1:2964523/1‑38 (MQ=255)
                                                   gagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGc                                        <  1:4579883/50‑1 (MQ=255)
                                                     gaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                 <  1:920718/71‑1 (MQ=255)
                                                      aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                <  1:3033685/71‑1 (MQ=255)
                                                        ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                     >  1:1318659/1‑48 (MQ=255)
                                                        ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                     >  1:1242477/1‑48 (MQ=255)
                                                        ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                     >  1:3013444/1‑48 (MQ=255)
                                                        ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                     >  1:3370390/1‑48 (MQ=255)
                                                         gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                      <  1:15115/62‑1 (MQ=255)
                                                         gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCAcc             >  1:397857/1‑71 (MQ=255)
                                                          tgtAGTAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                            <  1:4083921/39‑1 (MQ=255)
                                                          tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                             <  1:2288181/54‑1 (MQ=255)
                                                          tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                             <  1:5068583/54‑1 (MQ=255)
                                                           gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt           >  1:267150/1‑71 (MQ=255)
                                                            tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                           >  1:324205/1‑54 (MQ=255)
                                                             aGGAGATACCATCTTAAGCGAAAGCAAACGGCAGTTTAGc                                        <  1:177471/40‑1 (MQ=38)
                                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                            <  1:1009061/35‑1 (MQ=255)
                                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                            <  1:1787219/35‑1 (MQ=255)
                                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                            <  1:3020729/35‑1 (MQ=255)
                                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                            <  1:3027479/35‑1 (MQ=255)
                                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                          >  1:1595621/1‑37 (MQ=255)
                                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                         >  1:5311233/1‑54 (MQ=255)
                                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt        >  1:415944/1‑71 (MQ=255)
                                                              cgagaTCCCATCTTTGGCGTAGGTGAACGGCAGttt                                           >  1:3768728/2‑36 (MQ=12)
                                                              cgagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                             <  1:813430/49‑1 (MQ=255)
                                                               gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg       >  1:1693457/1‑68 (MQ=255)
                                                               gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg       <  1:834402/71‑4 (MQ=255)
                                                               gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                    >  1:4826911/1‑58 (MQ=255)
                                                               gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                    >  1:3596066/1‑58 (MQ=255)
                                                                  aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                               <  1:3959556/44‑1 (MQ=255)
                                                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                      >  1:498555/1‑35 (MQ=255)
                                                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                      >  1:2090346/1‑35 (MQ=255)
                                                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGa                                  >  1:585966/1‑39 (MQ=255)
                                                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc  >  1:5691076/1‑71 (MQ=255)
                                                                    |                                                                      
AGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAC  >  minE/608197‑608335

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: