Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,687,465 C→T 13.6% G272S (GGC→AGC)  pck ← phosphoenolpyruvate carboxykinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,687,4650CT13.6% 216.2 / 22.4 81G272S (GGC→AGC) pckphosphoenolpyruvate carboxykinase
Reads supporting (aligned to +/- strand):  ref base C (35/35);  new base T (4/7);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

CTTTCGACAGCTTGATAGTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCGGTGCCGGAA  >  minE/2687396‑2687532
                                                                     |                                                                   
cttTCGACAGCTTGATAGTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCg                                                                    >  1:2691364/1‑71 (MQ=255)
 tttCGACAGCTTGATAGTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCgt                                                                   >  1:1723192/1‑71 (MQ=255)
   tCGACAGCTTGATAGTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCgt                                                                   <  1:3028060/69‑1 (MQ=255)
   tCGACAGCTTGATAGTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCgt                                                                   <  1:2780853/69‑1 (MQ=255)
    cGACAGCTTGATAGTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCgtgt                                                                 <  1:2919532/70‑1 (MQ=255)
    cGACAGCTTGATAGTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCgtgt                                                                 <  1:1640155/70‑1 (MQ=255)
             gATAGTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGc                                                       >  1:21455/1‑71 (MQ=255)
               tAGTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGcc                                                      >  1:1679309/1‑70 (MQ=255)
                aGTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCTGTGTTCGTc                                                            <  1:2014841/63‑1 (MQ=255)
                 gTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCTGTGTTCGTCATCGCCaa                                                    <  1:3270685/70‑1 (MQ=255)
                 gTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCaa                                                    >  1:2881381/1‑70 (MQ=255)
                 gTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCaa                                                    <  1:2442238/70‑1 (MQ=255)
                 gTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAAt                                                   >  1:2591020/1‑71 (MQ=255)
                 gTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAAt                                                   >  1:2415001/1‑71 (MQ=255)
                   ttttGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAAt                                                   >  1:2109869/1‑69 (MQ=255)
                   ttttGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAAt                                                   >  1:1149818/1‑69 (MQ=255)
                    tttGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCTGTGTTCGTCATCGCCAATCAg                                                >  1:2347892/1‑71 (MQ=255)
                    tttGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATc                                                  >  1:3014116/1‑69 (MQ=255)
                    tttGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCa                                                 <  1:3044503/70‑1 (MQ=255)
                     ttGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAg                                                <  1:2850975/70‑1 (MQ=255)
                     ttGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAg                                                <  1:1918582/70‑1 (MQ=255)
                      tGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAgg                                               <  1:884911/70‑1 (MQ=255)
                      tGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAgg                                               <  1:587323/70‑1 (MQ=255)
                       gCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCTGTGTTCGTCATCGCCAATCAGGc                                              >  1:1384732/1‑70 (MQ=255)
                       gCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCTGTGTTCGTCATCGCCAATCAGGc                                              >  1:1638629/1‑70 (MQ=255)
                          tAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCTGTGTTCGTCATCGCCAATCAGGCGACg                                          <  1:2709781/71‑1 (MQ=255)
                             cagcCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGttt                                       <  1:1239351/71‑1 (MQ=255)
                              agcCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCTGTGTTCGTCATCGCCAATCAGGCGACGTTTc                                      <  1:186694/71‑1 (MQ=255)
                              agcCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCTGTGTTCGTCATCGCCAATCAGGCGACGTTTc                                      <  1:2914411/71‑1 (MQ=255)
                              agcCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCTGTGTTCGTCATCGCCAATCAGGCGACGTTTc                                      <  1:1370226/71‑1 (MQ=255)
                              agcCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCTGTGTTCGTCATCGCCAATCAGGCGACGTTTc                                      <  1:1425809/71‑1 (MQ=255)
                              agcCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTc                                      <  1:2459443/71‑1 (MQ=255)
                              agcCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTc                                      >  1:3506998/1‑71 (MQ=255)
                                  gccTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGt                                  <  1:2379382/71‑1 (MQ=255)
                                   ccTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGt                                  >  1:3590692/1‑70 (MQ=255)
                                      tCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAg                                                <  1:2437317/53‑1 (MQ=255)
                                       cgtaGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCgg                               <  1:3425881/66‑1 (MQ=255)
                                        gAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCAt                                                          >  1:3320562/1‑41 (MQ=255)
                                        gAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTc                                      >  1:831352/1‑61 (MQ=255)
                                        gAAGTTAAACACGCCGTCAACGTCCCAGCCGTGTTCg                                                              >  1:2196263/1‑37 (MQ=255)
                                         aaGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGa                           <  1:1752270/71‑1 (MQ=255)
                                          aGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGa                           >  1:370413/1‑70 (MQ=255)
                                          aGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGa                           >  1:3478096/1‑70 (MQ=255)
                                            ttAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCg                                <  1:2992053/63‑1 (MQ=255)
                                             tAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGa                           >  1:1909311/1‑67 (MQ=255)
                                             tAAACACGCCGTCATCGTCCCAGCCGTCTTCGTCATCGACAATCAGGCGACGTTTCg                                     >  1:381782/1‑57 (MQ=255)
                                              aaaCACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGt                                         >  1:2652128/1‑52 (MQ=255)
                                              aaaCACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGt                                         >  1:111201/1‑52 (MQ=255)
                                              aaaCACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGt                                         >  1:299563/1‑52 (MQ=255)
                                              aaaCACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCgg                                    >  1:1699535/1‑57 (MQ=255)
                                              aaaCACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCgg                                    >  1:3649007/1‑57 (MQ=255)
                                              aaaCACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAgg                       <  1:3367215/70‑1 (MQ=255)
                                               aaCACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAggg                      <  1:2011409/70‑1 (MQ=255)
                                                 cacGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGggtgg                   >  1:3028556/1‑71 (MQ=255)
                                                 cacGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGggtgg                   <  1:3656885/71‑1 (MQ=255)
                                                 cacGCCGTCATCGTACCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGggtgg                   <  1:560868/71‑1 (MQ=255)
                                                  acGCCGTCATCGTCCCATCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGaa                          <  1:1673405/63‑1 (MQ=255)
                                                  acGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGaa                          <  1:1024009/63‑1 (MQ=255)
                                                  acGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGaa                          <  1:2335303/63‑1 (MQ=255)
                                                  acGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGaa                          <  1:3338224/63‑1 (MQ=255)
                                                     ccGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGtttt               <  1:1395218/71‑1 (MQ=255)
                                                      cGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGt                                         >  1:3481566/1‑44 (MQ=255)
                                                      cGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGt                                         >  1:398567/1‑44 (MQ=255)
                                                      cGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTa              >  1:649091/1‑71 (MQ=255)
                                                        tCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTAc             <  1:1303366/70‑1 (MQ=255)
                                                         cATCGTCCCAGCTGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCg           >  1:531925/1‑71 (MQ=255)
                                                            cGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCg           >  1:1920615/1‑68 (MQ=255)
                                                             gTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGtttt               <  1:343695/63‑1 (MQ=255)
                                                             gTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCGGTg        <  1:2104026/70‑1 (MQ=255)
                                                              tCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCGGTGcc      >  1:903739/1‑71 (MQ=255)
                                                               cccAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCGGTg        <  1:3589536/68‑1 (MQ=255)
                                                               cccAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCGGTGCCg     >  1:1522498/1‑71 (MQ=255)
                                                               cccAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCGGTGCCg     >  1:3462469/1‑71 (MQ=255)
                                                               cccAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCGGTGCCg     >  1:3138598/1‑71 (MQ=255)
                                                                 cAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCGGTGCCGGa   <  1:2886926/71‑1 (MQ=255)
                                                                 cAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCGGTGCCGGa   <  1:1128329/71‑1 (MQ=255)
                                                                  aGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCGGTGCCGGaa  <  1:1604196/71‑1 (MQ=255)
                                                                   gCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGttt                <  1:492865/56‑1 (MQ=255)
                                                                   gCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGttt                <  1:100536/56‑1 (MQ=255)
                                                                   gCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTAc             >  1:392555/1‑59 (MQ=255)
                                                                   gCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTAc             >  1:1519290/1‑59 (MQ=255)
                                                                     |                                                                   
CTTTCGACAGCTTGATAGTTTTTGCGTAGCAGCCGCCTTCGAAGTTAAACACGCCGTCATCGTCCCAGCCGTGTTCGTCATCGCCAATCAGGCGACGTTTCGGGTCGGTGGAAAGGGTGGTTTTACCGGTGCCGGAA  >  minE/2687396‑2687532

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: