Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,421 A→C 14.6% V611V (GTT→GTG pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,4210AC14.6% 175.3 / 26.1 82V611V (GTT→GTGpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (35/35);  new base C (7/5);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTT  >  minE/608352‑608491
                                                                     |                                                                      
aGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACaa                                                                       >  1:4628740/1‑71 (MQ=255)
  cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                                                                     <  1:3273572/71‑1 (MQ=255)
  cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                                                                     >  1:5220220/1‑71 (MQ=255)
   aCCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                                    <  1:1479695/71‑1 (MQ=255)
    cccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                                    >  1:2916497/1‑70 (MQ=255)
      cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                                 <  1:1196942/71‑1 (MQ=255)
      cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                                 <  1:4386550/71‑1 (MQ=255)
           cGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGt                                                             >  1:3633108/1‑70 (MQ=255)
           cGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg                                                            <  1:899285/71‑1 (MQ=255)
             aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                                 <  1:2862864/64‑1 (MQ=255)
             aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                                 <  1:5010659/64‑1 (MQ=255)
             aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGAt                                                          >  1:3351370/1‑71 (MQ=255)
               cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATgg                                                        >  1:3361145/1‑71 (MQ=255)
                 gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                                       >  1:1669343/1‑70 (MQ=255)
                 gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                                       >  1:2058523/1‑70 (MQ=255)
                 gcgcgcCAGCACGACGACCGTCTGGGGTGTTAACCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                                       >  1:2981590/1‑70 (MQ=255)
                   gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                                    <  1:210484/55‑1 (MQ=255)
                   gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                                    <  1:4115035/55‑1 (MQ=255)
                    cgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg                                                    >  1:2490241/1‑70 (MQ=255)
                     gcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                                  >  1:781163/1‑71 (MQ=255)
                         gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg                                                            <  1:392881/57‑1 (MQ=255)
                         gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                              >  1:1605881/1‑71 (MQ=255)
                         gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                              >  1:1599268/1‑71 (MQ=255)
                         gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                              >  1:2419006/1‑71 (MQ=255)
                          cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCa                                                     >  1:2961006/1‑63 (MQ=255)
                          cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGAc                                             <  1:5056468/71‑1 (MQ=255)
                            cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg                                           <  1:985849/71‑1 (MQ=255)
                            cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg                                           <  1:1131484/71‑1 (MQ=255)
                             gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTATaa                                                               >  1:964193/1‑50 (MQ=255)
                             gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACt                                            >  1:779657/1‑69 (MQ=255)
                              acgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                                 <  1:5523749/47‑1 (MQ=255)
                               cgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGt                                        >  1:1411661/1‑71 (MQ=255)
                                gacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTc                                       <  1:2463970/71‑1 (MQ=255)
                                 acCGGCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCg                                      <  1:998088/71‑1 (MQ=255)
                                  ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTa                                                                  >  1:2208333/1‑42 (MQ=255)
                                  ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCgg                                     <  1:1327464/71‑1 (MQ=255)
                                  ccGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                                     >  1:3422674/1‑71 (MQ=255)
                                   cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                                    <  1:3959433/71‑1 (MQ=255)
                                   cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                                    <  1:858281/71‑1 (MQ=255)
                                   cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                                    <  1:3416414/71‑1 (MQ=255)
                                    gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                                >  1:1857874/1‑58 (MQ=255)
                                    gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                                >  1:5010318/1‑58 (MQ=255)
                                       tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                                  <  1:3858501/53‑1 (MQ=255)
                                       tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                                  <  1:3540171/53‑1 (MQ=255)
                                        ggggTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                                 <  1:1368028/37‑1 (MQ=255)
                                              ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGt                                                             >  1:5545369/1‑35 (MQ=12)
                                              ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGa                                                           <  1:18036/37‑1 (MQ=12)
                                              ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGta                                       >  1:2900815/1‑56 (MQ=255)
                                              ttACCCGTTTTCTTACCATACACAACGTTAGAAGt                                                             <  1:3801679/35‑1 (MQ=255)
                                              ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                                  <  1:1265460/46‑1 (MQ=255)
                                              ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACg                         >  1:4274323/1‑71 (MQ=255)
                                              ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGCAt                                  <  1:4533783/62‑1 (MQ=255)
                                               tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                                     <  1:4226166/58‑1 (MQ=255)
                                               tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGt                             <  1:3955736/66‑1 (MQ=255)
                                                aCCGGTTTTCTGGCCGTACACCACGTTAGAAGTGAta                                                         >  1:4473590/1‑36 (MQ=12)
                                                aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                                           >  1:3691654/1‑35 (MQ=255)
                                                aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                                           >  1:1242288/1‑35 (MQ=255)
                                                aCCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                        >  1:115970/1‑70 (MQ=255)
                                                    gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                        >  1:3634530/1‑66 (MQ=255)
                                                    gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                        >  1:1062633/1‑66 (MQ=255)
                                                    gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                   >  1:4710746/1‑71 (MQ=255)
                                                    gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                   >  1:2879597/1‑71 (MQ=255)
                                                    gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                   >  1:2541054/1‑71 (MQ=255)
                                                    gTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGgt                   >  1:2102100/1‑71 (MQ=255)
                                                     ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                                <  1:2378405/41‑1 (MQ=255)
                                                     ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                                <  1:5667625/41‑1 (MQ=255)
                                                     ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                                <  1:38625/41‑1 (MQ=255)
                                                     ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCAc                          <  1:4317281/63‑1 (MQ=255)
                                                      tttCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACt                                            >  1:4245046/1‑44 (MQ=12)
                                                       ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgtgt                 <  1:4831545/70‑1 (MQ=255)
                                                       ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGtgtg                >  1:1910693/1‑71 (MQ=255)
                                                       ttCTTACCATACACAACGTTAGAAGAGATGGTCAGAACAGACTGAGTCGGGATAg                                >  1:3934166/1‑55 (MQ=255)
                                                        tCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                              <  1:2128709/40‑1 (MQ=255)
                                                            aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCa              <  1:3568046/68‑1 (MQ=255)
                                                            aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt            >  1:1920407/1‑70 (MQ=255)
                                                            aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt            <  1:4449503/70‑1 (MQ=255)
                                                             ccATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATATCGTCACGGTAGGTGTGCAGttt          <  1:82175/71‑1 (MQ=255)
                                                             ccATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAg             <  1:5049511/68‑1 (MQ=255)
                                                               aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGt                       <  1:3255532/56‑1 (MQ=255)
                                                               aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCt        >  1:1650878/1‑71 (MQ=255)
                                                                 acacAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGa                                   >  1:5409775/1‑42 (MQ=255)
                                                                 acacAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGa                                   >  1:711807/1‑42 (MQ=255)
                                                                 acacAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGa                                   >  1:2274774/1‑42 (MQ=255)
                                                                  caccacGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAgg                    >  1:471017/1‑56 (MQ=255)
                                                                  cacaACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCa                           >  1:5395537/1‑49 (MQ=255)
                                                                  cacaACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGaa     >  1:2005261/1‑71 (MQ=255)
                                                                     aaCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt  >  1:1961911/1‑71 (MQ=255)
                                                                     aaCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt  >  1:932886/1‑71 (MQ=255)
                                                                     aaCGTTAGAAGTGATGGTCAGAACAGACTCAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt  <  1:1590691/71‑1 (MQ=255)
                                                                     |                                                                      
AGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTT  >  minE/608352‑608491

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: