Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,943,227 A→C 7.7% K100Q (AAA→CAA) ‡ bglA → 6‑phospho‑beta‑glucosidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,943,2270AC 7.7% 328.0 / 19.8 129K100Q (AAA→CAA) ‡bglA6‑phospho‑beta‑glucosidase A
Reads supporting (aligned to +/- strand):  ref base A (66/53);  new base C (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.88e-01

Only 100 of 140 total aligned reads displayed.

CATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAACTGCTGAA  >  minE/1943157‑1943297
                                                                      |                                                                      
cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                                                                        <  1:3915062/71‑1 (MQ=255)
cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                                                                        <  1:1202158/71‑1 (MQ=255)
cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                                                                        <  1:2737300/71‑1 (MQ=255)
cATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaa                                                                        <  1:1660478/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)
          tttGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGacg                                                               >  1:2557695/1‑68 (MQ=255)
          tttGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGac                                                                >  1:5418485/1‑68 (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)
            tGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAg                                                            >  1:5497366/1‑71 (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)
              ccGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt                                                          <  1:132510/71‑1 (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)
                      ggCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCaa                                                   <  1:3538224/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)
                           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)
                                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)
                                    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)
                                        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)
                                         cATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTc                               <  1:410016/71‑1 (MQ=255)
                                           tCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACgaaga                                            >  1:5644705/1‑56 (MQ=255)
                                           tCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACgaaga                                            >  1:4877372/1‑56 (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)
                                             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)
                                                    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:5013462/1‑71 (MQ=255)
                                                        cgcgCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTc                <  1:1179960/71‑1 (MQ=255)
                                                        cgcgCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTc                >  1:1742485/1‑71 (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)
                                                             aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAg                                  <  1:652724/48‑1 (MQ=255)
                                                             aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAg                                  <  1:4797649/48‑1 (MQ=255)
                                                             aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAg                                  <  1:4597890/48‑1 (MQ=255)
                                                             aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATg            >  1:1125801/1‑70 (MQ=255)
                                                              tttttCCAAAAGGCGATGAAGCTCAGCCAAACgaag                                             <  1:1103294/36‑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:80023/1‑48 (MQ=255)
                                                               ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGtt                                >  1:4741764/1‑48 (MQ=255)
                                                               ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGaa          >  1:3445961/1‑70 (MQ=255)
                                                               ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAAc         >  1:1181475/1‑71 (MQ=255)
                                                                tttCCAAAAGGCGATGATGCTCAGCCAAACGAAGAAGGGCTg                                     <  1:2978085/42‑1 (MQ=255)
                                                                tttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTa                             >  1:1850605/1‑50 (MQ=255)
                                                                tttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTa                             >  1:3207495/1‑50 (MQ=255)
                                                                   ccAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTa                             >  1:948559/1‑47 (MQ=255)
                                                                    cAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTg                                     <  1:316091/38‑1 (MQ=255)
                                                                    cAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAActgctg    <  1:3199517/71‑1 (MQ=255)
                                                                    cAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAActgctg    <  1:764870/71‑1 (MQ=255)
                                                                     aaaaGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAActgctg    <  1:1608034/70‑1 (MQ=255)
                                                                      aaaGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAACTGCTGaa  <  1:806231/71‑1 (MQ=255)
                                                                      |                                                                      
CATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAACTGCTGAA  >  minE/1943157‑1943297

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: