Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,943,186 A→C 24.7% K86T (AAA→ACA) ‡ bglA → 6‑phospho‑beta‑glucosidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,943,1860AC24.7% 133.0 / 65.1 97K86T (AAA→ACA) ‡bglA6‑phospho‑beta‑glucosidase A
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base C (9/15);  total (43/54)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

AACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAG  >  minE/1943116‑1943254
                                                                      |                                                                    
aaCCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCaa                                                                      <  1:4409780/71‑1 (MQ=255)
  ccATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAAt                                                                    >  1:2176368/1‑71 (MQ=255)
  ccATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAAt                                                                    >  1:5057066/1‑71 (MQ=255)
       aaGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGttt                                                                <  1:1281601/70‑1 (MQ=255)
       aaGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGtt                                                                 >  1:4994179/1‑69 (MQ=255)
         gCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCg                                                             >  1:5194858/1‑71 (MQ=255)
         gCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCg                                                             >  1:5188412/1‑71 (MQ=255)
            gTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTAc                                                          <  1:1646768/71‑1 (MQ=255)
             ttGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGtttt                                                               <  1:2610467/65‑1 (MQ=255)
               gATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGtt                                                                 >  1:2394812/1‑61 (MQ=255)
               gATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTa                                                           >  1:1758445/1‑67 (MQ=255)
               gATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATc                                                       >  1:4981092/1‑71 (MQ=255)
                  ttttATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCAt                                                    >  1:2646213/1‑71 (MQ=255)
                         gTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTg                                              <  1:543105/70‑1 (MQ=255)
                              tataAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACAt                                                        <  1:495589/55‑1 (MQ=255)
                               ataAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCa                                                     >  1:3619433/1‑57 (MQ=255)
                               ataAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGAcgc                                         >  1:2506636/1‑69 (MQ=255)
                                  aGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTc                                                              >  1:4256618/1‑45 (MQ=255)
                                  aGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCAtt                                    >  1:253814/1‑71 (MQ=255)
                                  aGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCAtt                                    >  1:3264055/1‑71 (MQ=255)
                                   ggAAGATATCGCGTTATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAAttt                                   <  1:1477950/71‑1 (MQ=255)
                                   ggAAGATATCGCGTTATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAAttt                                   >  1:764619/1‑71 (MQ=255)
                                   ggAAGATATCGCGTTATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGACGCGAAttt                                   <  1:2554187/71‑1 (MQ=255)
                                   ggAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATc                                                       <  1:2639763/51‑1 (MQ=255)
                                   ggAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATc                                                       <  1:4451693/51‑1 (MQ=255)
                                    gAAGATATCGCGTTATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAAtttt                                  >  1:532169/1‑71 (MQ=255)
                                    gAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGtttt                                                               >  1:3709977/1‑42 (MQ=255)
                                     aaGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACAt                                                        <  1:1327729/48‑1 (MQ=255)
                                     aaGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCAttttt                                 >  1:2585227/1‑71 (MQ=255)
                                        aCATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACgcgc                                       <  1:3753628/62‑1 (MQ=255)
                                         cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                             <  1:3915062/71‑1 (MQ=255)
                                         cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                             <  1:1660478/71‑1 (MQ=255)
                                         cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                             <  1:1202158/71‑1 (MQ=255)
                                         cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                             <  1:2737300/71‑1 (MQ=255)
                                           tCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaaaa                           <  1:2251707/71‑1 (MQ=255)
                                            cAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATc                                                       <  1:3208253/42‑1 (MQ=255)
                                            cAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaaa                            >  1:1561382/1‑69 (MQ=255)
                                            cAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaaa                            >  1:5078748/1‑69 (MQ=255)
                                             aaGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGa                                            >  1:1166436/1‑52 (MQ=255)
                                             aaGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGa                                            >  1:1521333/1‑52 (MQ=255)
                                             aaGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGa                                            >  1:2398283/1‑52 (MQ=255)
                                             aaGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAgg                         >  1:5079192/1‑71 (MQ=255)
                                             aaGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAgg                         >  1:3382894/1‑71 (MQ=255)
                                             aaGCTATTTGCCAAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGa                                            >  1:2704324/1‑52 (MQ=255)
                                                ttATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGc                                                 <  1:2787904/43‑1 (MQ=255)
                                                ttATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCt                                               <  1:1014966/45‑1 (MQ=255)
                                                ttATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCt                                               <  1:1751890/45‑1 (MQ=255)
                                                ttATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCt                                               <  1:2433945/45‑1 (MQ=255)
                                                ttATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCt                                               <  1:2951430/45‑1 (MQ=255)
                                                ttATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCt                                               <  1:2860158/45‑1 (MQ=255)
                                                ttATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTg                                              >  1:2727003/2‑47 (MQ=255)
                                                ttATTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTg                                              >  1:2378438/2‑47 (MQ=255)
                                                 tATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATc                                                       <  1:2929240/37‑1 (MQ=255)
                                                  aTTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGac                     >  1:2223538/1‑69 (MQ=255)
                                                  aTTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGAt                     >  1:2868506/1‑70 (MQ=255)
                                                   tttGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAAtttt                                  <  1:4656174/56‑1 (MQ=255)
                                                   tttGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGacg                    >  1:2557695/1‑68 (MQ=255)
                                                   tttGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGac                     >  1:5424863/1‑68 (MQ=255)
                                                   tttGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGac                     >  1:5418485/1‑68 (MQ=255)
                                                   tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCAttt                                   <  1:5500840/55‑1 (MQ=255)
                                                   tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCAttt                                   <  1:5487623/55‑1 (MQ=255)
                                                   tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTcc                               <  1:5689610/59‑1 (MQ=255)
                                                   tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATg                    <  1:2761226/70‑1 (MQ=255)
                                                   tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATg                    <  1:3887882/70‑1 (MQ=255)
                                                   tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGa                   <  1:2378942/71‑1 (MQ=255)
                                                   tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGa                   <  1:1408563/71‑1 (MQ=255)
                                                     tGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGACGaa                  <  1:373058/70‑1 (MQ=255)
                                                     tGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGACGaa                  <  1:1182698/70‑1 (MQ=255)
                                                     tGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGa                      >  1:2888843/1‑66 (MQ=255)
                                                     tGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAg                 >  1:5497366/1‑71 (MQ=255)
                                                      gCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAAtt                                    >  1:1556532/1‑51 (MQ=255)
                                                       ccGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGACGaa                  <  1:5681630/68‑1 (MQ=255)
                                                       ccGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGACGaa                  <  1:3346741/68‑1 (MQ=255)
                                                       ccGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGACGAAgtc               <  1:1774428/71‑3 (MQ=255)
                                                       ccGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGAt                     <  1:4395735/65‑1 (MQ=255)
                                                       ccGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGAt                     <  1:4237026/65‑1 (MQ=255)
                                                       ccGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGc                <  1:3439558/70‑1 (MQ=255)
                                                       ccGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt               <  1:132510/71‑1 (MQ=255)
                                                         gAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt               >  1:3975968/1‑69 (MQ=255)
                                                         gAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt               >  1:1439941/1‑69 (MQ=255)
                                                           aaTGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGc           >  1:3092252/1‑71 (MQ=255)
                                                             tGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGcgacga                   <  1:2414136/61‑1 (MQ=255)
                                                              gggCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCa         >  1:40758/1‑70 (MQ=255)
                                                              gggCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCa         >  1:4804842/1‑70 (MQ=255)
                                                               ggCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCaa        <  1:3538224/70‑1 (MQ=255)
                                                               ggCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCaa        <  1:3531599/70‑1 (MQ=255)
                                                                gCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATg                    <  1:5446878/57‑1 (MQ=255)
                                                                gCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATg                    <  1:4422594/57‑1 (MQ=255)
                                                                gCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCaa        <  1:5023357/69‑1 (MQ=255)
                                                                gCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCACATTTTTCCAAAAGGCGATg                    <  1:1586141/57‑1 (MQ=255)
                                                                 cTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCAttt                                   <  1:4405322/41‑1 (MQ=255)
                                                                 cTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCAttt                                   <  1:2040899/41‑1 (MQ=255)
                                                                 cTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCAttt                                   <  1:293134/41‑1 (MQ=255)
                                                                   tCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTcc                               <  1:5015099/43‑1 (MQ=255)
                                                                   tCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTcc                               <  1:5401820/43‑1 (MQ=255)
                                                                    cAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAgg                         >  1:3078846/1‑48 (MQ=255)
                                                                    cAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACg     >  1:4576113/1‑68 (MQ=255)
                                                                     aaaTGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACgaag  <  1:3468139/70‑1 (MQ=255)
                                                                      aaTGTTTTCGTACATCCATTGCCTGGAGGCGCATTTTTCCAAAAGGCGATGAAGCTCAGcc          >  1:3087242/1‑61 (MQ=255)
                                                                      aaTGTTCTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAg                          >  1:12434/1‑45 (MQ=255)
                                                                      |                                                                    
AACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAG  >  minE/1943116‑1943254

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: