Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,362 C→G 24.0% G631A (GGA→GCA) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,3620CG24.0% 86.9 / 26.7 50G631A (GGA→GCA) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (21/17);  new base G (4/8);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAG  >  minE/608296‑608428
                                                                  |                                                                  
gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                                                <  1:2120626/71‑1 (MQ=255)
    cGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc                                                            <  1:5113031/71‑1 (MQ=255)
      gAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgc                                                          <  1:3688099/71‑1 (MQ=255)
        aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc                                                        >  1:980134/1‑71 (MQ=255)
        aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc                                                        >  1:3105067/1‑71 (MQ=255)
              agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac                                                  <  1:1399283/71‑1 (MQ=255)
              agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac                                                  <  1:2490337/71‑1 (MQ=255)
                 ggCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                >  1:2161674/1‑70 (MQ=255)
                  gCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacgacg                                              <  1:3342709/71‑1 (MQ=255)
                      cggCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCg                                          >  1:616796/3‑71 (MQ=255)
                       ggCACCTTTGAGGTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCg                                          <  1:276826/69‑1 (MQ=255)
                        gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                >  1:3379603/1‑63 (MQ=255)
                        gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                >  1:3639095/1‑63 (MQ=255)
                        gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACgacgac                                            >  1:933964/1‑67 (MQ=255)
                         cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACgacgac                                            >  1:5403555/1‑66 (MQ=255)
                         cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                       >  1:5464753/1‑71 (MQ=255)
                         cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                       >  1:5213759/1‑71 (MQ=255)
                           ccTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCttg                                     <  1:2552236/71‑3 (MQ=255)
                              ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg                                                     >  1:2966184/1‑52 (MQ=255)
                              ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                  >  1:2646935/1‑71 (MQ=255)
                              ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                  >  1:4138635/1‑71 (MQ=255)
                               tgcggTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCgt                                  <  1:1831786/67‑1 (MQ=255)
                               tCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                  <  1:2347125/70‑1 (MQ=255)
                                cTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                >  1:2409864/1‑55 (MQ=255)
                                   gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgtgt                                <  1:502969/68‑1 (MQ=255)
                                   gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                             >  1:5106630/1‑71 (MQ=255)
                                   gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                             >  1:3947601/1‑71 (MQ=255)
                                   gTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAc                             <  1:1035538/71‑1 (MQ=255)
                                    tCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCgta                                 >  1:4677054/1‑65 (MQ=255)
                                      aCGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTggg                                    >  1:928333/1‑61 (MQ=255)
                                      aCGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAccgg                          >  1:4037640/1‑69 (MQ=255)
                                      aCGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAccg                           <  1:2079948/70‑2 (MQ=255)
                                        gACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                        <  1:3841163/71‑1 (MQ=255)
                                        gACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                        <  1:1515024/71‑1 (MQ=255)
                                         aCCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGttt                       <  1:4754005/71‑1 (MQ=255)
                                            gTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCt                    <  1:2869956/71‑1 (MQ=255)
                                             tGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                        >  1:5672208/1‑66 (MQ=255)
                                              gCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTa                  >  1:4342200/1‑71 (MQ=255)
                                               cATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCtggc                 <  1:1104124/71‑4 (MQ=255)
                                                aTCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTgg                                     >  1:2202266/1‑50 (MQ=255)
                                                   gggTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATa             <  1:1332168/71‑1 (MQ=255)
                                                      ttAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCg                                          <  1:1314312/39‑1 (MQ=255)
                                                        aGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACaa        >  1:4628740/1‑71 (MQ=255)
                                                          cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg      >  1:5220220/1‑71 (MQ=255)
                                                          cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg      <  1:3273572/71‑1 (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)
                                                                gTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGttt                       <  1:5369907/48‑1 (MQ=255)
                                                                  ccGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTAcc                >  1:321984/1‑53 (MQ=255)
                                                                  |                                                                  
GCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAG  >  minE/608296‑608428

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