Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,113,361 A→G 41.9% D481D (GAT→GAC pps ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,113,3610AG41.9% 20.1 / 54.9 74D481D (GAT→GACppsphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base A (13/30);  new base G (15/16);  total (28/46)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

GGCGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGG  >  minE/1113304‑1113428
                                                         |                                                                   
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                                          <  1:1156858/69‑1 (MQ=255)
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                                          <  1:2112250/69‑1 (MQ=255)
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                                          <  1:3342444/69‑1 (MQ=255)
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                                          <  1:1882391/69‑1 (MQ=255)
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                                          <  1:2797426/69‑1 (MQ=255)
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                                          <  1:3165183/69‑1 (MQ=255)
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                                          <  1:1742403/69‑1 (MQ=255)
 gcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGtt                                                         <  1:475191/69‑1 (MQ=255)
 gcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGtt                                                         <  1:153027/69‑1 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTc                                                       <  1:651237/69‑1 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTc                                                       <  1:452961/69‑1 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTc                                                       <  1:3390048/69‑1 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTc                                                       <  1:218451/69‑1 (MQ=255)
    aaGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTc                                                       <  1:1601628/68‑1 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGATGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACg                                                   >  1:2662468/1‑67 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATGGTTTCTACGc                                                  >  1:3104207/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:3078085/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:1133534/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:3189165/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:2621115/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:3306947/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:585011/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:203200/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:1988479/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:1483526/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:1751046/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:1414268/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:1352898/1‑68 (MQ=255)
         gAAAGCACGGTACGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                                  >  1:3280761/1‑68 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCtt                                                <  1:488280/59‑2 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:3306065/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:3282898/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:1489128/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:3328126/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:2803214/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:1833843/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:502102/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:1206272/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:2618240/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:781099/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:1953006/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:91073/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:2092723/59‑1 (MQ=255)
                     cGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGTCATCGTTTCTACGCTg                                                <  1:148864/58‑1 (MQ=255)
                     cGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                                <  1:3005303/58‑1 (MQ=255)
                                 gTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCa                          >  1:306744/1‑68 (MQ=255)
                                    catcatcaCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCaa                      >  1:35854/1‑69 (MQ=255)
                                    catcatcaCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCaa                      >  1:910285/1‑69 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 <  1:654221/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 <  1:894471/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 <  1:3330957/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 <  1:947814/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 <  1:1557132/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 <  1:1665190/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 >  1:231583/1‑68 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 <  1:2980504/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 >  1:2317721/1‑68 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 <  1:2688601/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 <  1:2674673/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 <  1:3225511/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 >  1:1931398/1‑68 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCtt                                          >  1:1957813/1‑42 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCtt                                          >  1:775271/1‑42 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCtt                                          >  1:1114822/1‑42 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 >  1:560677/1‑67 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 >  1:866856/1‑67 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 >  1:245195/1‑67 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 >  1:2630235/1‑67 (MQ=255)
                                             tttCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATa             <  1:2917393/69‑1 (MQ=255)
                                             tttCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATa             <  1:3098133/69‑1 (MQ=255)
                                              ttCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATa             <  1:1684667/68‑1 (MQ=255)
                                                      cAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg                 <  1:1834251/56‑1 (MQ=255)
                                                        gATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCgg  <  1:380918/69‑1 (MQ=255)
                                                        gATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCgg  <  1:2031060/69‑1 (MQ=255)
                                                         |                                                                   
GGCGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGG  >  minE/1113304‑1113428

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: