Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 799,144 C→G 39.3% P149P (CCG→CCC treA ← periplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE799,1440CG39.3% 32.4 / 51.7 56P149P (CCG→CCCtreAperiplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base C (32/2);  new base G (7/15);  total (39/17)
Fisher's exact test for biased strand distribution p-value = 1.00e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAAC  >  minE/799080‑799208
                                                                |                                                                
tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATGCCTCGCGAAAACGTCCGCCGGGc                                                               >  1:81341/1‑68 (MQ=255)
tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:1369795/1‑68 (MQ=255)
tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:223241/1‑68 (MQ=255)
tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:3381263/1‑68 (MQ=255)
tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:1477765/1‑68 (MQ=255)
tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:2396153/1‑68 (MQ=255)
tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:2501869/1‑68 (MQ=255)
           ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGt                                                   <  1:2721777/69‑1 (MQ=255)
           ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGt                                                   <  1:204758/69‑1 (MQ=255)
                    gTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGc                                          >  1:1163460/1‑69 (MQ=255)
                          tAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAg                                                     >  1:1653604/1‑52 (MQ=255)
                          tAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGt                                       >  1:2637191/1‑66 (MQ=255)
                          tAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGt                                       >  1:2887917/1‑66 (MQ=255)
                              tGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACaga                                 >  1:2449044/1‑68 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:2368546/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:2615958/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:2791089/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:2263411/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:222910/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:1981825/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:1972956/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:3005558/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:1814340/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:1745785/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:1313035/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:615884/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:80524/69‑1 (MQ=255)
                               gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAcaga                               <  1:3110854/69‑1 (MQ=255)
                                     gTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCa                         >  1:358886/1‑69 (MQ=255)
                                     gTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCa                         >  1:119942/1‑69 (MQ=255)
                                      tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCa                         >  1:3063869/1‑68 (MQ=255)
                                      tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt                        >  1:347214/1‑69 (MQ=255)
                                      tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt                        >  1:957497/1‑69 (MQ=255)
                                      tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt                        >  1:1991619/1‑69 (MQ=255)
                                      tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt                        >  1:1070586/1‑69 (MQ=255)
                                      tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt                        >  1:185862/1‑69 (MQ=255)
                                      tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt                        >  1:1265937/1‑69 (MQ=255)
                                      tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAAAAGAGAATCCCAt                        >  1:1185827/1‑69 (MQ=255)
                                       aaTATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAtt                       >  1:2965852/1‑69 (MQ=255)
                                       aaTATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:2281198/61‑1 (MQ=255)
                                        atatACCTCGCGAAAGCGTCCGCCCGGCACGACATaa                                                      >  1:2859876/1‑37 (MQ=255)
                                        atatACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAttt                      >  1:1013861/1‑69 (MQ=255)
                                        atatACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAttt                      >  1:2291147/1‑69 (MQ=255)
                                        atatACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAttt                      >  1:773300/1‑69 (MQ=255)
                                        atatACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAttt                      >  1:367176/1‑69 (MQ=255)
                                             ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:2527826/1‑69 (MQ=255)
                                             ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:416574/1‑69 (MQ=255)
                                             ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:1318125/1‑69 (MQ=255)
                                             ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:342748/1‑69 (MQ=255)
                                             ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:1556474/1‑69 (MQ=255)
                                             ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:1863319/1‑69 (MQ=255)
                                             ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:2992244/1‑69 (MQ=255)
                                             ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:2921012/1‑69 (MQ=255)
                                                            cGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAAc  >  1:1582394/1‑69 (MQ=255)
                                                            cGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAAc  >  1:1112481/1‑69 (MQ=255)
                                                            cGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAAc  >  1:2531539/1‑69 (MQ=255)
                                                                |                                                                
TTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAAC  >  minE/799080‑799208

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: