Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,259 G→C 9.4% S665S (TCC→TCG pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,2590GC 9.4% 266.7 / 14.2 95S665S (TCC→TCGpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (46/40);  new base C (7/2);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

ACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCT  >  minE/608190‑608329
                                                                     |                                                                      
aCCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAgg                                                                       <  1:2844479/71‑1 (MQ=255)
       agGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGATTGAGAAGGTGTACGAGAtcc                                                                 >  1:2780402/1‑68 (MQ=255)
       agGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAcc                                                                <  1:5183340/71‑1 (MQ=255)
        ggTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAcc                                                                >  1:3288311/1‑70 (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:2734750/1‑70 (MQ=255)
                ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAg                                                        >  1:1839091/1‑70 (MQ=255)
                 tACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg                                                      >  1:4288506/1‑71 (MQ=255)
                  aCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg                                                      >  1:1604319/1‑70 (MQ=255)
                  aCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg                                                      >  1:4504611/1‑70 (MQ=255)
                      aCTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAg                                                  >  1:4495514/1‑70 (MQ=255)
                       cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                <  1:3001695/71‑1 (MQ=255)
                       cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                <  1:3577571/71‑1 (MQ=255)
                       cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                <  1:601071/71‑1 (MQ=255)
                       cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                <  1:4424500/71‑1 (MQ=255)
                       cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                <  1:2603901/71‑1 (MQ=255)
                          cgtcgtCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                              <  1:747403/70‑1 (MQ=255)
                           gtcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACg                                            >  1:4515711/1‑71 (MQ=255)
                            tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATc                                                             <  1:5622563/53‑1 (MQ=255)
                            tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACgg                                           >  1:1411211/1‑71 (MQ=255)
                            tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACgg                                           >  1:384793/1‑71 (MQ=255)
                            tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTATAGCGTAAGCAAACgg                                           <  1:1381964/71‑1 (MQ=255)
                              gtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCtt                                                           >  1:2512603/1‑53 (MQ=255)
                              gtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCtt                                                           >  1:2525290/1‑53 (MQ=255)
                              gtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCtt                                                           >  1:4431174/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)
                                         gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                              >  1:5103641/1‑71 (MQ=255)
                                          tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                             <  1:666583/71‑1 (MQ=255)
                                          tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                              >  1:3162435/1‑70 (MQ=255)
                                          tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                              >  1:2595823/1‑70 (MQ=255)
                                            cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                             >  1:726501/1‑69 (MQ=255)
                                            cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                           <  1:184186/71‑1 (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)
                                                   aaCGATAGAGAATGTGTAGGAGATACCATCTTTAg                                                        <  1:1369106/35‑1 (MQ=255)
                                                   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:4943854/1‑49 (MQ=255)
                                                     cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                        >  1:4597779/1‑49 (MQ=255)
                                                     cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                        >  1:2378580/1‑49 (MQ=255)
                                                      gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                           <  1:4060079/61‑1 (MQ=255)
                                                      gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                           <  1:5255274/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:3348759/1‑70 (MQ=255)
                                                        tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                <  1:4135293/70‑1 (MQ=255)
                                                        tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                >  1:1762916/1‑70 (MQ=255)
                                                         agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa              >  1:4378990/1‑71 (MQ=255)
                                                         agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa              >  1:2969125/1‑71 (MQ=255)
                                                         agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa              >  1:5573347/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)
                                                             aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc          <  1:3026919/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:5603934/35‑1 (MQ=255)
                                                                     ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                      <  1:3020729/35‑1 (MQ=255)
                                                                     ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                      <  1:1009061/35‑1 (MQ=255)
                                                                     ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                      <  1:1787219/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)
                                                                     |                                                                      
ACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCT  >  minE/608190‑608329

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: