Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,943,229 A→G 7.7% K100K (AAA→AAG) ‡ bglA → 6‑phospho‑beta‑glucosidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,943,2290AG 7.7% 360.9 / 14.9 129K100K (AAA→AAG) ‡bglA6‑phospho‑beta‑glucosidase A
Reads supporting (aligned to +/- strand):  ref base A (66/53);  new base G (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 = 5.94e-01

Only 100 of 139 total aligned reads displayed.

TCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAACTGCTGAAA  >  minE/1943159‑1943298
                                                                      |                                                                     
tCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaaaa                                                                       <  1:2251707/71‑1 (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:3887882/70‑1 (MQ=255)
        tttGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGa                                                               <  1:2378942/71‑1 (MQ=255)
          tGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGACGaa                                                              <  1:1182698/70‑1 (MQ=255)
          tGCCGAGATGGGCTTCACCTGTCTGCGTATTTCCATTGCCTGGGCGCGAATTTTCCCTCAGGGCGACGaa                                                              <  1:373058/70‑1 (MQ=255)
          tGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGa                                                                  >  1:2888843/1‑66 (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:4237026/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:40758/1‑70 (MQ=255)
                    ggCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCaa                                                    <  1:3538224/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)
                            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)
                              gTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAggg                                         >  1:828286/1‑71 (MQ=255)
                              gTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAgg                                          >  1:5163031/1‑70 (MQ=255)
                               ttttCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTc                                                          <  1:4460066/53‑1 (MQ=255)
                                  tCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCt                                                           <  1:5515062/49‑1 (MQ=255)
                                  tCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTg                                      <  1:4410441/70‑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:3188849/37‑1 (MQ=255)
                                      aCATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCg                                                                   <  1:5414125/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)
                                       cATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGttg                                <  1:3918682/71‑2 (MQ=255)
                                       cATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTc                                >  1:5297251/1‑71 (MQ=255)
                                       cATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTc                                >  1:4727373/1‑71 (MQ=255)
                                         tCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACgaaga                                             >  1:4877372/1‑56 (MQ=255)
                                         tCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACgaaga                                             >  1:5644705/1‑56 (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:4206074/70‑1 (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:406420/1‑39 (MQ=255)
                                                    gACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCaaa                                                   >  1:44252/1‑39 (MQ=255)
                                                    gACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGatat                     >  1:5517412/1‑69 (MQ=255)
                                                     aCGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGtt                  <  1:5131945/71‑1 (MQ=255)
                                                      cgcgCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTc                 >  1:5013462/1‑71 (MQ=255)
                                                      cgcgCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTc                 <  1:3287107/71‑1 (MQ=255)
                                                      cgcgCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTc                 <  1:1179960/71‑1 (MQ=255)
                                                      cgcgCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTc                 >  1:1742485/1‑71 (MQ=255)
                                                        cgcATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGa               <  1:3036843/71‑1 (MQ=255)
                                                         gcATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGaa                                    >  1:5697141/1‑49 (MQ=255)
                                                         gcATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGaa                                    >  1:720485/1‑49 (MQ=255)
                                                          cATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATg             <  1:4988778/71‑1 (MQ=255)
                                                           aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAg                                   <  1:4597890/48‑1 (MQ=255)
                                                           aTTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAg                                   <  1:4797649/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: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:80023/1‑48 (MQ=255)
                                                             ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGtt                                 >  1:2644877/1‑48 (MQ=255)
                                                             ttttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAAc          >  1:1181475/1‑71 (MQ=255)
                                                              tttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTa                              >  1:3207495/1‑50 (MQ=255)
                                                               ttCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGa            <  1:1905496/67‑1 (MQ=255)
                                                                 ccAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTa                              >  1:948559/1‑47 (MQ=255)
                                                                  cAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAActgctg     <  1:3199517/71‑1 (MQ=255)
                                                                  cAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAActgctg     <  1:764870/71‑1 (MQ=255)
                                                                    aaaGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAACTGCTGaa   <  1:1726144/71‑1 (MQ=255)
                                                                    aaaGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAACTGCTGaa   <  1:806231/71‑1 (MQ=255)
                                                                    aaaGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTATACga                           >  1:3859978/1‑47 (MQ=255)
                                                                     aaGGCTATGAAGCTCAGCCAAACGAATAAGGGCTGAAGTTCTAc                             <  1:378493/44‑1 (MQ=38)
                                                                     aaGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGaa           >  1:1524891/1‑62 (MQ=255)
                                                                     aaGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGaa           >  1:5476393/1‑62 (MQ=255)
                                                                     aaGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAACTGCTGaaa  <  1:3512501/71‑1 (MQ=255)
                                                                     aaGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAACTGCTGaaa  >  1:4092482/1‑71 (MQ=255)
                                                                      |                                                                     
TCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGGCGATGAAGCTCAGCCAAACGAAGAAGGGCTGAAGTTCTACGATGATATGTTCGATGAACTGCTGAAA  >  minE/1943159‑1943298

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: