Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,943,223 T→C 8.4% F98F (TTT→TTC bglA → 6‑phospho‑beta‑glucosidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,943,2230TC 8.4% 333.4 / 15.3 119F98F (TTT→TTCbglA6‑phospho‑beta‑glucosidase A
Reads supporting (aligned to +/- strand):  ref base T (58/51);  new base C (4/6);  total (62/57)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 136 total aligned reads displayed.

AAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAAC  >  minE/1943153‑1943290
                                                                      |                                                                   
aaGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCAttttt                                                                     >  1:2585227/1‑71 (MQ=255)
    cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                                                                 <  1:2737300/71‑1 (MQ=255)
      tCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaaaa                                                               <  1:2251707/71‑1 (MQ=255)
       cAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaaa                                                                >  1:1561382/1‑69 (MQ=255)
        aaGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAgg                                                             >  1:3382894/1‑71 (MQ=255)
        aaGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAgg                                                             >  1:5079192/1‑71 (MQ=255)
             aTTTGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGac                                                         >  1:2223538/1‑69 (MQ=255)
             aTTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGAt                                                         >  1:2868506/1‑70 (MQ=255)
              tttGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGac                                                         >  1:5418485/1‑68 (MQ=255)
              tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATg                                                        <  1:3887882/70‑1 (MQ=255)
              tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATg                                                        <  1:782863/70‑1 (MQ=255)
              tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGa                                                       <  1:2378942/71‑1 (MQ=255)
                tGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGACGaa                                                      <  1:1182698/70‑1 (MQ=255)
                tGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGa                                                          >  1:2888843/1‑66 (MQ=255)
                  ccGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGACGaa                                                      <  1:3346741/68‑1 (MQ=255)
                  ccGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGACGaa                                                      <  1:5681630/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)
                  ccGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt                                                   <  1:132510/71‑1 (MQ=255)
                    gAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt                                                   >  1:1439941/1‑69 (MQ=255)
                    gAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt                                                   >  1:3975968/1‑69 (MQ=255)
                      aaTGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGc                                               >  1:3092252/1‑71 (MQ=255)
                        tGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGcgacga                                                       <  1:2414136/61‑1 (MQ=255)
                         gggCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCa                                             >  1:4804842/1‑70 (MQ=255)
                         gggCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCa                                             >  1:40758/1‑70 (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)
                              tCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTcc                                                                   <  1:5015099/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)
                                 aaTGTTCTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAg                                                              >  1:12434/1‑45 (MQ=255)
                                  aTGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTc                                                  >  1:757194/1‑56 (MQ=255)
                                  aTGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGaagaag                                   >  1:54746/1‑71 (MQ=255)
                                   tGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAgg                                  >  1:5491372/1‑71 (MQ=255)
                                    gTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGa                                                          >  1:4443691/1‑46 (MQ=255)
                                    gTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAgg                                  >  1:5163031/1‑70 (MQ=255)
                                        tCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt                                                   <  1:5515062/49‑1 (MQ=255)
                                        tCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTg                              >  1:452740/1‑70 (MQ=255)
                                          gTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAg                           >  1:889944/1‑71 (MQ=255)
                                           tACATCCATTGCCTGGACGCGCATTTTTACAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGt                          <  1:5561714/71‑1 (MQ=255)
                                            aCATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCg                                                           <  1:5414125/37‑1 (MQ=255)
                                            aCATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCg                                                           <  1:3188849/37‑1 (MQ=255)
                                            aCATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGaa                            >  1:3119779/1‑68 (MQ=255)
                                            aCATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGtt                         >  1:35992/1‑71 (MQ=255)
                                            aCATCCATTGCCTGGACGCGCATTTATCCAAAAGGCg                                                           <  1:3012323/37‑1 (MQ=255)
                                             cATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACgaaga                                     <  1:3792148/58‑1 (MQ=255)
                                             cATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGttg                        <  1:3918682/71‑2 (MQ=255)
                                             cATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTc                        >  1:5297251/1‑71 (MQ=255)
                                               tCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACgaaga                                     >  1:5644705/1‑56 (MQ=255)
                                                ccATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGc                                                    <  1:5413390/40‑1 (MQ=255)
                                                ccATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAggg                                 >  1:1128385/1‑59 (MQ=255)
                                                 cATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt                                                   >  1:1939601/1‑40 (MQ=255)
                                                 cATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTAc                     <  1:1387380/70‑1 (MQ=255)
                                                 cATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTAc                     >  1:5328887/1‑70 (MQ=255)
                                                 cATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTAc                     >  1:87697/1‑70 (MQ=255)
                                                 cATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTAc                     <  1:4206074/70‑1 (MQ=255)
                                                 cATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAACGGCTGAAGTTCTa                      >  1:4928409/1‑69 (MQ=255)
                                                   ttGCCTGGACGCGCATTTTTCCAAAAGGCTATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTAc                     <  1:2714317/68‑1 (MQ=255)
                                                    tGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACgat                  >  1:3278914/1‑70 (MQ=255)
                                                     gCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACgatga                >  1:1582035/1‑71 (MQ=255)
                                                      ccTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACgatga                <  1:3718369/70‑1 (MQ=255)
                                                        tGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTAc                     <  1:4498266/63‑1 (MQ=255)
                                                        tGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTAc                     <  1:2602028/63‑1 (MQ=255)
                                                         ggACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTa                      >  1:527416/1‑61 (MQ=255)
                                                         ggACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATg            >  1:489374/1‑71 (MQ=255)
                                                          gACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCaaa                                           >  1:44252/1‑39 (MQ=255)
                                                          gACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCaaa                                           >  1:406420/1‑39 (MQ=255)
                                                          gACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGatat             >  1:5517412/1‑69 (MQ=255)
                                                          gACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGt           >  1:3306118/1‑71 (MQ=255)
                                                           aCGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGtt          <  1:5131945/71‑1 (MQ=255)
                                                            cgcgCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTc         >  1:1742485/1‑71 (MQ=255)
                                                            cgcgCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTc         <  1:1179960/71‑1 (MQ=255)
                                                              cgcATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGa       <  1:3036843/71‑1 (MQ=255)
                                                              cgcATTTTTCCAAAAGGCGATGAACCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGa       <  1:4649584/71‑1 (MQ=255)
                                                               gcATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGaa                            >  1:720485/1‑49 (MQ=255)
                                                               gcATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGaa                            >  1:5697141/1‑49 (MQ=255)
                                                                cATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATg     <  1:4988778/71‑1 (MQ=255)
                                                                 aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACgaagaa                                    >  1:3942095/1‑39 (MQ=255)
                                                                 aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAg                           <  1:4797649/48‑1 (MQ=255)
                                                                 aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAg                           <  1:4597890/48‑1 (MQ=255)
                                                                 aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAg                           <  1:652724/48‑1 (MQ=255)
                                                                 aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATg     >  1:1125801/1‑70 (MQ=255)
                                                                 aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATg     >  1:250038/1‑70 (MQ=255)
                                                                  tttttCCAAAAGGCGATGAAGCTCAGCCAAACgaag                                      <  1:1103294/36‑1 (MQ=255)
                                                                  tttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCt                               <  1:1462416/43‑1 (MQ=255)
                                                                  tttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCt                               <  1:3942348/43‑1 (MQ=255)
                                                                  tttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGata              <  1:361333/60‑1 (MQ=255)
                                                                   ttttCCAAAAGGCGATGAAGCTTAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGa       >  1:282214/1‑66 (MQ=255)
                                                                   ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGtt                         >  1:2644877/1‑48 (MQ=255)
                                                                   ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGtt                         >  1:4741764/1‑48 (MQ=255)
                                                                   ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGtt                         >  1:80023/1‑48 (MQ=255)
                                                                   ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGaa   >  1:3445961/1‑70 (MQ=255)
                                                                   ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAAc  >  1:1181475/1‑71 (MQ=255)
                                                                    tttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTa                      >  1:3207495/1‑50 (MQ=255)
                                                                    tttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTa                      >  1:1850605/1‑50 (MQ=255)
                                                                      |                                                                   
AAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAAC  >  minE/1943153‑1943290

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: