Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,113,349 T→C 42.5% K485K (AAA→AAG pps ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,113,3490TC42.5% 18.4 / 55.4 73K485K (AAA→AAGppsphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base T (18/24);  new base C (15/16);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GGCGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATA  >  minE/1113304‑1113417
                                             |                                                                    
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                               <  1:1156858/69‑1 (MQ=255)
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                               <  1:3342444/69‑1 (MQ=255)
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                               <  1:2112250/69‑1 (MQ=255)
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                               <  1:3165183/69‑1 (MQ=255)
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                               <  1:2797426/69‑1 (MQ=255)
ggcGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGt                                               <  1:1882391/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)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACgg                                                              >  1:745622/1‑51 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACgg                                                              >  1:2379823/1‑51 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACgg                                                              >  1:738631/1‑51 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACgg                                                              >  1:587152/1‑51 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACgg                                                              >  1:936695/1‑51 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTc                                            <  1:651237/69‑1 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTc                                            <  1:452961/69‑1 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTc                                            <  1:218451/69‑1 (MQ=255)
   gAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTc                                            <  1:3390048/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:585011/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:1133534/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                       >  1:1352898/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                       >  1:203200/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                       >  1:1988479/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                       >  1:1414268/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                       >  1:1483526/1‑68 (MQ=255)
         gAAAGCACGGTCCGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                       >  1:1751046/1‑68 (MQ=255)
         gAAAGCACGGTACGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGc                                       >  1:3280761/1‑68 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCtt                                     <  1:488280/59‑2 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                     <  1:2092723/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                     <  1:1206272/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                     <  1:3328126/59‑1 (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:91073/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                     <  1:502102/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                     <  1:1833843/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                     <  1:2803214/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                     <  1:1953006/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                     <  1:2618240/59‑1 (MQ=255)
                    ccGGGTTACCGACGTTCATCATCACCTTCAACGGCAGGTCCGGCATCGTTTCTACGCTg                                     <  1:781099/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:910285/1‑69 (MQ=255)
                                    catcatcaCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCaa           >  1:35854/1‑69 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      <  1:894471/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      <  1:654221/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      <  1:947814/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      <  1:3225511/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      <  1:3330957/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      <  1:2980504/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      <  1:1557132/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      <  1:1665190/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      <  1:2688601/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      <  1:2674673/68‑1 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      >  1:1931398/1‑68 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      >  1:2317721/1‑68 (MQ=255)
                                          caCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      >  1:231583/1‑68 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCtt                               >  1:1114822/1‑42 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCtt                               >  1:1957813/1‑42 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCtt                               >  1:775271/1‑42 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      >  1:866856/1‑67 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      >  1:245195/1‑67 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      >  1:2630235/1‑67 (MQ=255)
                                           aCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCg      >  1:560677/1‑67 (MQ=255)
                                             tttCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATa  <  1:2917393/69‑1 (MQ=255)
                                             tttCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATa  <  1:3098133/69‑1 (MQ=255)
                                             |                                                                    
GGCGAAGTCGAAAGCACGGTCCGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATA  >  minE/1113304‑1113417

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: