Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,943,161 A→G 5.9% K78E (AAG→GAG) ‡ bglA → 6‑phospho‑beta‑glucosidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,943,1610AG 5.9% 193.3 / 2.5 67K78E (AAG→GAG) ‡bglA6‑phospho‑beta‑glucosidase A
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base G (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GCCAGGAAAATACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGG  >  minE/1943097‑1943231
                                                                |                                                                      
gCCAGGAAAATACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCt                                                                    <  1:2063208/69‑1 (MQ=255)
gCCAGGAAAATACTATCCAAACCATGAAGCCATTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTa                                                                   >  1:4559184/1‑70 (MQ=255)
 ccAGGAAAATACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTa                                                                   <  1:916333/69‑1 (MQ=255)
 ccAGGAAAATACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTAt                                                                  <  1:2964851/70‑1 (MQ=255)
  cAGGAAAATACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTAttt                                                                <  1:5597187/71‑1 (MQ=255)
   aGGAAAATACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTg                                                               >  1:666073/1‑71 (MQ=255)
    ggAAAATACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGc                                                              <  1:5199738/71‑1 (MQ=255)
       aaaTACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTAtt                                                                 >  1:4103314/1‑65 (MQ=255)
       aaaTACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGa                                                           <  1:1215014/71‑1 (MQ=255)
       aaaTACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGa                                                           <  1:2714342/71‑1 (MQ=255)
       aaaTACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGa                                                           <  1:5507287/71‑1 (MQ=255)
         aTACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGaaa                                                         >  1:4811708/1‑71 (MQ=255)
           aCTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATg                                                       >  1:5543989/1‑71 (MQ=255)
            cTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATgg                                                      <  1:21432/71‑1 (MQ=255)
            cTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATg                                                       <  1:2554015/70‑1 (MQ=255)
                 cAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCtt                                                  <  1:1933464/70‑1 (MQ=255)
                 cAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCtt                                                  >  1:589269/1‑70 (MQ=255)
                   aaCCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCtt                                                  <  1:3115914/68‑1 (MQ=255)
                   aaCCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCaa                                               <  1:4409780/71‑1 (MQ=255)
                   aaCCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCa                                                >  1:869867/1‑70 (MQ=255)
                     ccATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAAt                                             >  1:2176368/1‑71 (MQ=255)
                     ccATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAAt                                             >  1:5057066/1‑71 (MQ=255)
                        tGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCt                                                                    <  1:438776/45‑1 (MQ=255)
                        tGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTAt                                                                  >  1:967911/1‑47 (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)
                                ttGATTTTTATGGTCACTATAAGGAAGACATCAAGCTAt                                                                  <  1:15386/39‑1 (MQ=255)
                                ttGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGa                                                           <  1:472836/46‑1 (MQ=255)
                                ttGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGa                                                           <  1:3593195/46‑1 (MQ=255)
                                ttGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGtttt                                        <  1:2610467/65‑1 (MQ=255)
                                 tGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGa                                                           <  1:1574393/45‑1 (MQ=255)
                                  gATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGcc                                                             >  1:2196286/1‑42 (MQ=255)
                                  gATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGtt                                          >  1:2394812/1‑61 (MQ=255)
                                  gATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTa                                    >  1:1758445/1‑67 (MQ=255)
                                  gATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATc                                >  1:4981092/1‑71 (MQ=255)
                                   aTTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTg                                                               <  1:3607882/39‑1 (MQ=255)
                                   aTTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATggg                                                     >  1:3871363/1‑49 (MQ=255)
                                     ttttATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCAt                             >  1:2646213/1‑71 (MQ=255)
                                         aTGGTCACTATAAGGAAGACATCAAGCTATTTGCCGaaa                                                         >  1:529061/1‑39 (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:3264055/1‑71 (MQ=255)
                                                     aGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCAtt             >  1:253814/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:1202158/71‑1 (MQ=255)
                                                            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)
                                                               cAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATc                                <  1:3208253/42‑1 (MQ=255)
                                                               cAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaaa     >  1:5078748/1‑69 (MQ=255)
                                                               cAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCaaa     >  1:1561382/1‑69 (MQ=255)
                                                                aaGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGa                     >  1:2398283/1‑52 (MQ=255)
                                                                aaGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGa                     >  1:1521333/1‑52 (MQ=255)
                                                                aaGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGa                     >  1:1166436/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)
                                                                |                                                                      
GCCAGGAAAATACTATCCAAACCATGAAGCCGTTGATTTTTATGGTCACTATAAGGAAGACATCAAGCTATTTGCCGAAATGGGCTTCAAATGTTTTCGTACATCCATTGCCTGGACGCGCATTTTTCCAAAAGG  >  minE/1943097‑1943231

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: