Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,667 A→G 7.1% G529G (GGT→GGC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,6670AG 7.1% 224.9 / 6.5 84G529G (GGT→GGC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (42/36);  new base G (3/3);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

CGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCG  >  minE/608597‑608733
                                                                      |                                                                   
cGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGa                                                                     <  1:4881319/71‑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:829012/70‑1 (MQ=255)
      cGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                                <  1:5559538/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)
            gAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACca                                                         <  1:2480866/71‑1 (MQ=255)
            gAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACca                                                         <  1:463864/71‑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:758341/50‑1 (MQ=255)
                        cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                                                  <  1:4223321/50‑1 (MQ=255)
                        cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                                                  <  1:4423978/50‑1 (MQ=255)
                        cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                                                  <  1:38620/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:2390300/69‑1 (MQ=255)
                           tCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                                            <  1:2585078/69‑1 (MQ=255)
                           gggCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCga                                          <  1:373464/69‑2 (MQ=255)
                            cGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATgg                                               >  1:1824049/1‑65 (MQ=255)
                            cGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATgg                                               >  1:1881742/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)
                                tatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGAt                                                             <  1:1127564/47‑1 (MQ=255)
                                tatTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGcccg                                                                 <  1:1613128/43‑3 (MQ=12)
                                 atTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                                <  1:574690/43‑1 (MQ=255)
                                 atTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACg                                   >  1:1149481/1‑71 (MQ=255)
                                  tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGa                                                              <  1:218253/44‑1 (MQ=255)
                                  tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGt                                  >  1:3464094/1‑71 (MQ=255)
                                  tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGt                                  >  1:4624953/1‑71 (MQ=255)
                                  tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGt                                  >  1:4120457/1‑71 (MQ=255)
                                  tttGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGc                                                                >  1:1975401/1‑42 (MQ=255)
                                   ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGt                                  >  1:4482926/1‑70 (MQ=255)
                                   ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGt                                  >  1:3549713/1‑70 (MQ=255)
                                    tGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTc                                 <  1:4495706/70‑1 (MQ=255)
                                    tGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACACCATGGTGCGGATA‑ACGTCa                                >  1:1826472/1‑71 (MQ=255)
                                       ttGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACgg                             <  1:709857/71‑1 (MQ=255)
                                        tGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGt                            >  1:3881570/1‑71 (MQ=255)
                                         gCAGACAGGGAGTCAGCAGCGACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTc                           >  1:1409918/1‑71 (MQ=255)
                                         gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATa‑a                                     >  1:2960014/1‑61 (MQ=255)
                                         gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATa‑a                                     >  1:5666001/1‑61 (MQ=255)
                                         gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTc                           >  1:4431275/1‑71 (MQ=255)
                                         gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTc                           >  1:5365467/1‑71 (MQ=255)
                                         gCAGACAGGGAGTCAGCAGCAACGGACAGAACAGCGATACCACACGCCATGGTGCGGATa‑a                                     >  1:898897/1‑61 (MQ=255)
                                          cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTc                           <  1:70206/70‑1 (MQ=255)
                                          cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCg                          >  1:3893114/1‑71 (MQ=255)
                                            gACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCgcg                        <  1:994016/71‑3 (MQ=255)
                                            gACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTg                        <  1:2664318/71‑1 (MQ=255)
                                             aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTg                        >  1:4380746/1‑70 (MQ=255)
                                             aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGc                       <  1:259841/71‑1 (MQ=255)
                                              cAGGGAGTCAGCAGCACCGGACAGACCAGCGATACcaca                                                       >  1:4979632/1‑39 (MQ=255)
                                              cAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGc                       <  1:5063168/70‑1 (MQ=255)
                                              cAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCa                      <  1:2870297/71‑1 (MQ=255)
                                               aGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAg                     >  1:4929491/1‑71 (MQ=255)
                                                 ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCa                      >  1:4325930/1‑68 (MQ=255)
                                                 ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCa                      >  1:767944/1‑68 (MQ=255)
                                                     tCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt               <  1:5488704/71‑1 (MQ=255)
                                                       agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAg                     >  1:4966893/1‑63 (MQ=255)
                                                       agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt               >  1:2339797/1‑69 (MQ=255)
                                                       agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt               >  1:3859894/1‑69 (MQ=255)
                                                       agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCa             <  1:5218487/71‑1 (MQ=255)
                                                            caACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑Ac                                    <  1:1283199/43‑1 (MQ=255)
                                                             aaCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑Ac                                    <  1:1612146/42‑1 (MQ=255)
                                                              aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagaga         >  1:5597367/1‑68 (MQ=255)
                                                              aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagaga         >  1:606469/1‑68 (MQ=255)
                                                              aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGc      >  1:225641/1‑71 (MQ=255)
                                                               ccgACAGGCCCGCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCt     >  1:2536424/3‑71 (MQ=255)
                                                                ggACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt    <  1:203083/71‑1 (MQ=255)
                                                                  aCAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg  >  1:190803/1‑71 (MQ=255)
                                                                      |                                                                   
CGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCG  >  minE/608597‑608733

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