Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,091 T→C 16.3% E721E (GAA→GAG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,0910TC16.3% 135.7 / 17.0 55E721E (GAA→GAG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (26/20);  new base C (2/7);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 7.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGAT  >  minE/608021‑608160
                                                                      |                                                                     
tAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGttt                                                                       >  1:941510/1‑71 (MQ=255)
     aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                                                                  >  1:1828542/1‑71 (MQ=255)
          gatGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTACCGGGtttt                                                                      >  1:1958951/3‑62 (MQ=255)
             gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTc                                                            >  1:2861031/1‑69 (MQ=255)
             gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                                                           <  1:2574956/70‑1 (MQ=255)
              aaaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCCTCGag                                                         <  1:589931/71‑1 (MQ=255)
               aaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                                                           >  1:1936877/1‑68 (MQ=255)
                aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                                         >  1:566372/1‑69 (MQ=255)
                aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                                         >  1:529629/1‑69 (MQ=255)
                  gTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGagca                                                       >  1:3640/1‑69 (MQ=255)
                    aCTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAt                                                   >  1:3378394/1‑71 (MQ=255)
                     cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                                                  <  1:1265563/69‑1 (MQ=12)
                     cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                                                  <  1:3645971/69‑1 (MQ=12)
                       gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                 >  1:2362056/1‑70 (MQ=255)
                        cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa                                               >  1:2963865/1‑71 (MQ=255)
                            gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                           <  1:3744660/71‑1 (MQ=255)
                              cAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga                                            <  1:2729242/68‑2 (MQ=9)
                               agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCg                                             <  1:263141/66‑1 (MQ=9)
                               aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                        <  1:385533/71‑1 (MQ=255)
                               aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                        >  1:396873/1‑71 (MQ=255)
                                taTACACGGATGGTCAGCTGCGGATATTTTTCCGGGttt                                                                       >  1:454934/2‑39 (MQ=255)
                                gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                       <  1:803100/68‑1 (MQ=9)
                                gATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                 <  1:303918/61‑1 (MQ=255)
                                 aTACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                       >  1:3380607/1‑70 (MQ=255)
                                   acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                    >  1:1149639/1‑71 (MQ=255)
                                    cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                    <  1:251882/70‑1 (MQ=255)
                                     acGGATGGTCAGCTGCGGAAATTTTTCCGGGTTTTCCATCGCGTCGagcagc                                                     >  1:25087/1‑52 (MQ=255)
                                       ggATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttg                                <  1:1812077/71‑2 (MQ=255)
                                       ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                 >  1:2612834/1‑70 (MQ=255)
                                            gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                 <  1:731714/49‑1 (MQ=255)
                                             tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                              <  1:160098/51‑1 (MQ=255)
                                             tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                          <  1:928150/71‑1 (MQ=255)
                                              cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                           >  1:254837/1‑69 (MQ=255)
                                               aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                        >  1:1010417/1‑71 (MQ=255)
                                               aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                        >  1:3155460/1‑71 (MQ=255)
                                               aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                        <  1:3241491/71‑1 (MQ=255)
                                                gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                      >  1:2916952/1‑56 (MQ=255)
                                                gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                         <  1:1912623/69‑1 (MQ=255)
                                                 cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                      <  1:2538351/71‑1 (MQ=255)
                                                 cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                      >  1:40621/1‑71 (MQ=255)
                                                 cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                      <  1:3583132/71‑1 (MQ=255)
                                                   taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                    >  1:1934025/3‑70 (MQ=255)
                                                   gCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGc                    >  1:3690568/1‑71 (MQ=255)
                                                     ggATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                              >  1:2967272/1‑43 (MQ=255)
                                                      gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                 <  1:1929460/71‑1 (MQ=255)
                                                      gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                 >  1:616967/1‑71 (MQ=12)
                                                         atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGt                                  <  1:1328372/51‑1 (MQ=255)
                                                         atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacca              >  1:564393/1‑71 (MQ=255)
                                                          tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac             <  1:2149654/71‑1 (MQ=255)
                                                           ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc            >  1:534188/1‑71 (MQ=255)
                                                           ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc            >  1:3251582/1‑71 (MQ=255)
                                                           ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc            >  1:855116/1‑71 (MQ=255)
                                                            tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc            <  1:3278536/70‑1 (MQ=255)
                                                                cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga       <  1:1278827/71‑1 (MQ=255)
                                                                cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga       <  1:868636/71‑1 (MQ=255)
                                                                 gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                       <  1:558239/38‑1 (MQ=255)
                                                                 gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgat      <  1:1627687/71‑1 (MQ=255)
                                                                   gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg     >  1:1553947/1‑70 (MQ=255)
                                                                    ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                        <  1:237535/50‑1 (MQ=255)
                                                                    ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                        <  1:1953432/50‑1 (MQ=255)
                                                                    ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga   >  1:1659397/1‑71 (MQ=255)
                                                                      ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat  >  1:1977977/1‑70 (MQ=255)
                                                                      |                                                                     
TAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGAT  >  minE/608021‑608160

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