Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 900,712 C→T 44.1% T186M (ACG→ATG)  ycjM → predicted glucosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE900,7120CT44.1% 34.8 / 94.0 120T186M (ACG→ATG) ycjMpredicted glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (51/16);  new base T (32/21);  total (83/37)
Fisher's exact test for biased strand distribution p-value = 7.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 129 total aligned reads displayed.

CTCAGCGCCGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTAGCCCTGAA  >  minE/900645‑900773
                                                                   |                                                             
ctcAGCGCCGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCTATCATTCAAcg                                                              >  1:2822660/1‑69 (MQ=255)
ctcAGCGCCGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAAcg                                                              >  1:2391261/1‑69 (MQ=255)
ctcAGCGCCGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAAcg                                                              >  1:1529983/1‑69 (MQ=255)
ctcAGCGCCGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAAcg                                                              >  1:2205857/1‑69 (MQ=255)
ctcAGCGCCGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAAcg                                                              >  1:3223453/1‑69 (MQ=255)
ctcAGCGCCGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAAcg                                                              >  1:2516188/1‑69 (MQ=255)
ctcAGCGCCGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAAcg                                                              >  1:2854994/1‑69 (MQ=255)
ctcAGCGCCGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAAcg                                                              >  1:1040106/1‑69 (MQ=255)
       ccGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCAtt                                                       <  1:710031/69‑1 (MQ=255)
       ccGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCAtt                                                       <  1:2841593/69‑1 (MQ=255)
        cGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCAtt                                                       <  1:1446217/68‑1 (MQ=255)
        cGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCAtt                                                       <  1:3392722/68‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:2706836/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:2504765/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:2071873/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:3344350/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:3424166/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:500922/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:690329/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:686581/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:1415179/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:1283986/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:1619616/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:854243/69‑1 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACg                                  >  1:1331254/1‑68 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACg                                  >  1:259106/1‑68 (MQ=255)
                             aCCGTTATTAACGCCATTCCAGAT‑‑GCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGAt                                >  1:526271/1‑68 (MQ=255)
                             aCCGTTATTAACGCCATCCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACg                                  >  1:2471502/1‑68 (MQ=255)
                              ccGTTATTAACGCCATTCCAGATTCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:1399192/68‑1 (MQ=255)
                                gTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTATTGACGATCa                              >  1:2408430/1‑69 (MQ=255)
                                gTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCa                              >  1:1229569/1‑69 (MQ=255)
                                gTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCa                              >  1:2463628/1‑69 (MQ=255)
                                gTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCa                              >  1:2973594/1‑69 (MQ=255)
                                gTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCa                              >  1:1589496/1‑69 (MQ=255)
                                gTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCa                              >  1:2093027/1‑69 (MQ=255)
                                   attaACGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCaaa                            >  1:3266761/1‑68 (MQ=255)
                                      aaCGCCATTCCAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:1718519/60‑1 (MQ=255)
                                      aaCGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGac                                                >  1:3429198/1‑45 (MQ=255)
                                      aaCGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGa                        >  1:62959/1‑69 (MQ=255)
                                      aaCGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGa                        >  1:3374152/1‑69 (MQ=255)
                                      aaCGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGa                        >  1:2884984/1‑69 (MQ=255)
                                      aaCGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGa                        >  1:3250941/1‑69 (MQ=255)
                                      aaCGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGa                        >  1:2135918/1‑69 (MQ=255)
                                      aaCGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGa                        >  1:625113/1‑69 (MQ=255)
                                      aaCGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGa                        >  1:2829154/1‑69 (MQ=255)
                                          ccATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTg                    <  1:358193/69‑1 (MQ=255)
                                          ccATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTg                    <  1:1999516/69‑1 (MQ=255)
                                          ccATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTg                    <  1:2983192/69‑1 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGacc                                               >  1:3291764/1‑37 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:322179/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:609469/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:2891939/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:330581/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:3374321/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:899037/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:458881/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:48015/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:506013/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:865417/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:1852182/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:217509/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:2184543/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:2007800/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:2242801/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:2374294/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:1806082/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:2467332/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:1674123/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:1644342/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:1127841/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGATCa                              >  1:1079317/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAATGCGCCATTTATGGACAACCTTTAGCGACGATCa                              >  1:2326041/1‑54 (MQ=255)
                                               ccAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACTATCAAATTGACCTGAAtt                >  1:1007462/1‑68 (MQ=255)
                                               ccAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTa               >  1:606658/1‑69 (MQ=255)
                                               ccAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTa               >  1:2092595/1‑69 (MQ=255)
                                               ccAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTa               >  1:1909634/1‑69 (MQ=255)
                                               ccAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTa               >  1:151115/1‑69 (MQ=255)
                                               ccAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTa               >  1:457838/1‑69 (MQ=255)
                                               ccAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTa               >  1:2020836/1‑69 (MQ=255)
                                               ccAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTa               >  1:2749640/1‑69 (MQ=255)
                                               ccAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTa               >  1:2950390/1‑69 (MQ=255)
                                               ccACATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGCGACGAta                               >  1:2569513/1‑52 (MQ=255)
                                                  gATGCGCGATCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:1901337/48‑1 (MQ=255)
                                                    tGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTg                                     >  1:2734518/1‑42 (MQ=255)
                                                    tGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTa          >  1:2053553/1‑69 (MQ=255)
                                                    tGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTa          >  1:1013490/1‑69 (MQ=255)
                                                    tGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTa          >  1:2717592/1‑69 (MQ=255)
                                                    tGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTa          >  1:1411426/1‑69 (MQ=255)
                                                    tGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTa          >  1:1495363/1‑69 (MQ=255)
                                                    tGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTa          >  1:2449776/1‑69 (MQ=255)
                                                    tGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTa          >  1:1874460/1‑69 (MQ=255)
                                                    tGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTa          >  1:2206332/1‑69 (MQ=255)
                                                          aTCATTCAATGCGCCATTTATGGACCACCTTTAGcgacga                                 <  1:3303070/40‑1 (MQ=38)
                                                           tCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTAGCCCTGa   <  1:2902777/69‑1 (MQ=255)
                                                           tCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTAGCCCTGa   <  1:2259481/69‑1 (MQ=255)
                                                           tCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTAGCCCTGa   <  1:2220960/69‑1 (MQ=255)
                                                           tCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTAGCCCTGa   <  1:2737136/69‑1 (MQ=255)
                                                           tCATTCAACGCGCCATTTATGGAACACCTTTAGTGACGATCAAATTGACCTGAATTACCGTAGCCCTGa   <  1:2391991/69‑1 (MQ=255)
                                                            cATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTAGCCCTGaa  >  1:435948/1‑69 (MQ=255)
                                                            cATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTAGCCCTGa   <  1:2626269/68‑1 (MQ=255)
                                                                   |                                                             
CTCAGCGCCGTCACTCGCCCGCGTGCGTTACCGTTATTAACGCCATTCCAGATGCGCGATCATTCAACGCGCCATTTATGGACCACCTTTAGTGACGATCAAATTGACCTGAATTACCGTAGCCCTGAA  >  minE/900645‑900773

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: