Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 511,646 T→A 44.9% S563S (TCT→TCA)  ybhJ → predicted hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE511,6460TA44.9% 17.0 / 73.1 69S563S (TCT→TCA) ybhJpredicted hydratase
Reads supporting (aligned to +/- strand):  ref base T (10/28);  new base A (11/20);  total (21/48)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAG  >  minE/511591‑511711
                                                       |                                                                 
tCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAAt                                                       >  1:2343647/1‑68 (MQ=255)
tCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAAt                                                       >  1:3225434/1‑68 (MQ=255)
tCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAAt                                                       >  1:3330267/1‑68 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:3398199/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:2646150/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:2644016/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:2635210/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:2723086/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:2444513/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:2423195/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:237039/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:2240044/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:33236/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:3418584/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:369749/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:1181678/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:555801/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:1206955/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:980240/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:1137158/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:1331755/61‑1 (MQ=255)
 ccTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:1540167/61‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGAtt                                                 <  1:1690394/67‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:2610250/69‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:3189917/69‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:2417674/69‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:2549669/69‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:2369116/69‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:2188431/69‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:2165048/69‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:993353/69‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:605804/69‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:1949134/69‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:162313/69‑1 (MQ=255)
       cgaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:528054/69‑1 (MQ=255)
        gaAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:1970466/68‑1 (MQ=255)
                        acgaACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:689154/52‑1 (MQ=255)
                        acgaACTGATTCCTTCCGGTGAAACCTCTTCATATCGt                                                             <  1:705461/38‑1 (MQ=255)
                          gaACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTgg                                               <  1:471892/50‑1 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATcc                                                     <  1:1060342/39‑1 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATcc                                                     <  1:2792402/39‑1 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATcc                                                     <  1:3052929/39‑1 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATcc                                                     <  1:1193270/39‑1 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATcc                                                     <  1:3222281/39‑1 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATcc                                                     <  1:1244808/39‑1 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATcc                                                     <  1:3373777/39‑1 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATcc                                                     <  1:2147161/39‑1 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGtctc                        >  1:364335/1‑68 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTcg                       >  1:109350/1‑69 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTcg                       >  1:1554261/1‑69 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTcg                       >  1:2325798/1‑69 (MQ=255)
                               gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTcg                       >  1:2135118/1‑69 (MQ=255)
                                            gAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGt          >  1:1786213/1‑69 (MQ=255)
                                            gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCAcgccgat                  >  1:1412195/1‑59 (MQ=255)
                                            gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAggg          >  1:3176068/1‑68 (MQ=255)
                                            gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAgg           >  1:1107020/1‑68 (MQ=255)
                                            gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt          >  1:2762203/1‑69 (MQ=255)
                                            gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt          >  1:44516/1‑69 (MQ=255)
                                            gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt          >  1:1934504/1‑69 (MQ=255)
                                            gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt          >  1:2297477/1‑69 (MQ=255)
                                            gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt          >  1:2550664/1‑69 (MQ=255)
                                            gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt          >  1:2921522/1‑69 (MQ=255)
                                            gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt          >  1:24264/1‑69 (MQ=255)
                                            gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt          >  1:2438641/1‑69 (MQ=255)
                                               aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTa        <  1:2559533/68‑1 (MQ=255)
                                               aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTa        <  1:2224982/68‑1 (MQ=255)
                                               aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTa        >  1:1179/1‑68 (MQ=255)
                                               aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTa        <  1:491834/68‑1 (MQ=255)
                                               aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTa        <  1:340364/68‑1 (MQ=255)
                                                  tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt    <  1:386894/69‑1 (MQ=255)
                                                  tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt    <  1:209880/69‑1 (MQ=255)
                                                  tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt    <  1:533397/69‑1 (MQ=255)
                                                  tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt    <  1:1644619/69‑1 (MQ=255)
                                                  tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt    <  1:1184506/69‑1 (MQ=255)
                                                  tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt    <  1:2189004/69‑1 (MQ=255)
                                                  tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt    <  1:2220091/69‑1 (MQ=255)
                                                      cttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAg  <  1:2108237/67‑1 (MQ=255)
                                                      cttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAg  <  1:1677678/67‑1 (MQ=255)
                                                      cttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAg  <  1:2797675/67‑1 (MQ=255)
                                                       |                                                                 
TCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAG  >  minE/511591‑511711

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