Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,943,226 A→T 7.6% P99P (CCA→CCT bglA → 6‑phospho‑beta‑glucosidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,943,2260AT 7.6% 365.4 / 14.2 129P99P (CCA→CCTbglA6‑phospho‑beta‑glucosidase A
Reads supporting (aligned to +/- strand):  ref base A (66/53);  new base T (4/6);  total (70/59)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 137 total aligned reads displayed.

CATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAACTGCTG  >  minE/1943157‑1943295
                                                                     |                                                                     
cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                                                                      <  1:3915062/71‑1 (MQ=255)
cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                                                                      <  1:2737300/71‑1 (MQ=255)
cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                                                                      <  1:1660478/71‑1 (MQ=255)
  tCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaaaa                                                                    <  1:2251707/71‑1 (MQ=255)
   cAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaaa                                                                     >  1:5078748/1‑69 (MQ=255)
    aaGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAgg                                                                  >  1:3382894/1‑71 (MQ=255)
         aTTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGAt                                                              >  1:2868506/1‑70 (MQ=255)
          tttGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGacg                                                             >  1:2557695/1‑68 (MQ=255)
          tttGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGac                                                              >  1:5424863/1‑68 (MQ=255)
          tttGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGac                                                              >  1:5418485/1‑68 (MQ=255)
          tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATg                                                             <  1:782863/70‑1 (MQ=255)
          tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATg                                                             <  1:3887882/70‑1 (MQ=255)
          tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGa                                                            <  1:2378942/71‑1 (MQ=255)
            tGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGACGaa                                                           <  1:1182698/70‑1 (MQ=255)
            tGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAg                                                          >  1:5497366/1‑71 (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:4237026/65‑1 (MQ=255)
              ccGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGAt                                                              <  1:4395735/65‑1 (MQ=255)
              ccGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGc                                                         <  1:3439558/70‑1 (MQ=255)
              ccGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt                                                        <  1:132510/71‑1 (MQ=255)
                gAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt                                                        >  1:1439941/1‑69 (MQ=255)
                  aaTGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGc                                                    >  1:3092252/1‑71 (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:4422594/57‑1 (MQ=255)
                       gCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCaa                                                 <  1:5023357/69‑1 (MQ=255)
                       gCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCACATTTTTCCAAAAGGCGATg                                                             <  1:1586141/57‑1 (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)
                              aTGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTc                                                       >  1:757194/1‑56 (MQ=255)
                              aTGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGaagaag                                        >  1:2320890/1‑71 (MQ=255)
                              aTGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGaagaag                                        >  1:54746/1‑71 (MQ=255)
                               tGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAgg                                       >  1:5491372/1‑71 (MQ=255)
                                gTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGa                                                               >  1:3968923/1‑46 (MQ=255)
                                gTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGa                                                               >  1:4443691/1‑46 (MQ=255)
                                gTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAggg                                      >  1:828286/1‑71 (MQ=255)
                                gTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAgg                                       >  1:5163031/1‑70 (MQ=255)
                                  tttCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCt                                    >  1:5520228/1‑71 (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: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)
                                         cATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTc                             >  1:5297251/1‑71 (MQ=255)
                                            ccATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGc                                                         <  1:5413390/40‑1 (MQ=255)
                                            ccATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAggg                                      >  1:1128385/1‑59 (MQ=255)
                                            ccATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTAc                          >  1:1567920/1‑71 (MQ=255)
                                             cATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt                                                        >  1:1939601/1‑40 (MQ=255)
                                             cATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTAc                          >  1:5328887/1‑70 (MQ=255)
                                             cATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTAc                          <  1:4206074/70‑1 (MQ=255)
                                             cATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTAc                          <  1:1387380/70‑1 (MQ=255)
                                             cATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACg                         <  1:520970/71‑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)
                                                     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)
                                                       aCGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGtt               <  1:5131945/71‑1 (MQ=255)
                                                        cgcgCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTc              >  1:1742485/1‑71 (MQ=255)
                                                        cgcgCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTc              >  1:5013462/1‑71 (MQ=255)
                                                        cgcgCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTc              <  1:3287107/71‑1 (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:5697141/1‑49 (MQ=255)
                                                           gcATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGaa                                 >  1:720485/1‑49 (MQ=255)
                                                             aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACgaagaa                                         >  1:3942095/1‑39 (MQ=255)
                                                             aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAg                                <  1:4797649/48‑1 (MQ=255)
                                                             aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAg                                <  1:4597890/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:3942348/43‑1 (MQ=255)
                                                              tttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCt                                    <  1:1462416/43‑1 (MQ=255)
                                                              tttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGata                   <  1:361333/60‑1 (MQ=255)
                                                               ttttCCAAAAGGCGATGAAGCTTAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGa            >  1:282214/1‑66 (MQ=255)
                                                               ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGtt                              >  1:4741764/1‑48 (MQ=255)
                                                               ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGtt                              >  1:2644877/1‑48 (MQ=255)
                                                               ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGtt                              >  1:80023/1‑48 (MQ=255)
                                                               ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAAc       >  1:1181475/1‑71 (MQ=255)
                                                                tttCCAAAAGGCGATGATGCTCAGCCAAACGAAGAAGGGCTg                                   <  1:2978085/42‑1 (MQ=255)
                                                                tttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTa                           >  1:3207495/1‑50 (MQ=255)
                                                                 ttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTc                             <  1:749697/47‑1 (MQ=255)
                                                                 ttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGa         <  1:1905496/67‑1 (MQ=255)
                                                                   ccAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTa                           >  1:948559/1‑47 (MQ=255)
                                                                    cAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAActgctg  <  1:3199517/71‑1 (MQ=255)
                                                                     |                                                                     
CATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAACTGCTG  >  minE/1943157‑1943295

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: