Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,943,994 A→G 45.5% P355P (CCA→CCG bglA → 6‑phospho‑beta‑glucosidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,943,9940AG45.5% 4.5 / 44.0 55P355P (CCA→CCGbglA6‑phospho‑beta‑glucosidase A
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base G (12/13);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGTATGAGCGTTATCAGAGGCCGCTGTTTATTG  >  minE/1943938‑1944055
                                                        |                                                             
ggCAGCGTACCAAACCCGTTTGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTgcgc                                                   <  1:3302572/69‑1 (MQ=255)
ggCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTgcgc                                                   <  1:2804944/69‑1 (MQ=255)
ggCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTgcgc                                                   <  1:716507/69‑1 (MQ=255)
ggCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTgcgc                                                   <  1:1928075/69‑1 (MQ=255)
ggCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTgcgc                                                   <  1:394261/69‑1 (MQ=255)
ggCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTgcgc                                                   <  1:1526891/69‑1 (MQ=255)
ggCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGATCTgcgc                                                   <  1:436210/69‑1 (MQ=255)
 gCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTgcgc                                                   <  1:2158182/68‑1 (MQ=255)
 gCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTgcgc                                                   <  1:2586047/68‑1 (MQ=255)
 gCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTgcgc                                                   <  1:3410632/68‑1 (MQ=255)
 gCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTgcgc                                                   <  1:496811/68‑1 (MQ=255)
        aCCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACt                                           <  1:1366614/69‑1 (MQ=255)
        aCCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACt                                           <  1:569488/69‑1 (MQ=255)
        aCCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACt                                           <  1:1074441/69‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:438135/68‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:3239119/68‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:847254/68‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:2677372/68‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:2187372/68‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:992903/68‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:1199612/68‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:1449904/68‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:2104768/68‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:2054752/68‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:2073786/68‑1 (MQ=255)
              tccGTATGTTAAAGCGTCTGACTGGGGCTGGCAGATTGATCCGGTTGGGCTGCGCTACGCACTTTGCGa                                     <  1:2981268/68‑1 (MQ=255)
              cccGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACt                                           <  1:2100564/63‑1 (MQ=255)
                    tGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGt                                >  1:2524906/1‑68 (MQ=255)
                    tGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGt                                >  1:3002375/1‑68 (MQ=255)
                    tGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGt                                >  1:1339520/1‑68 (MQ=255)
                    tGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGt                                >  1:44690/1‑68 (MQ=255)
                    tGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGt                                >  1:2624762/1‑68 (MQ=255)
                    tGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGt                                >  1:483129/1‑68 (MQ=255)
                    tGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGt                                >  1:1910316/1‑68 (MQ=255)
                      ttAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGTATg                             >  1:83558/1‑69 (MQ=255)
                      ttAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGTATg                             >  1:3253312/1‑69 (MQ=255)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGc                                              >  1:2173487/1‑43 (MQ=37)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGAGCTGTATGAGCGtt                       >  1:1336362/1‑66 (MQ=255)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGAGCTGTATGAGCGtt                       >  1:579077/1‑66 (MQ=255)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGAGCTGTATGAGCGtt                       >  1:1084922/1‑66 (MQ=255)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGAGCTGTATGAGCGtt                       >  1:3036607/1‑66 (MQ=255)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGAGCTGTATGAGCGtt                       >  1:1603640/1‑66 (MQ=255)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGAGCTGTATGAGCGtt                       >  1:1734114/1‑66 (MQ=255)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGAGCTGTATGAGCGtt                       >  1:2576541/1‑66 (MQ=255)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGAGCTGTATGAGCGtt                       >  1:2379569/1‑66 (MQ=255)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGAGCTGTATGAGCGtt                       >  1:2188459/1‑66 (MQ=255)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGAGCTGTATGAGCGt                        >  1:21974/1‑65 (MQ=255)
                               cTGACTGGGGCTGGCAGATTGATCCGGTTGGTCTGCGCTACGCAATTTGCGAGCTGTATGAGCGtt                       >  1:2189412/1‑66 (MQ=255)
                                        gCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGTATGAGCGTTATCAGAGGCCg            >  1:201986/1‑68 (MQ=255)
                                        gCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGTATGAGCGTTATCAGAGGCCGc           >  1:3416042/1‑69 (MQ=255)
                                        gCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGTATGAGCGTTATCAGAGGCCGc           >  1:1923780/1‑69 (MQ=255)
                                        gCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGTATGAGCGTTATCAGAGGCCGc           >  1:2699655/1‑69 (MQ=255)
                                        gCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGTATGAGCGTTATCAGAGGCCGc           >  1:1244490/1‑69 (MQ=255)
                                            gCAGATTGATCCGGTTGGTCTGCGCTACGCACTTTGCGa                                     <  1:1036783/39‑1 (MQ=37)
                                                 ttGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGTATGAGCGTTATCAGAGGCCGCTGTTTATTg  >  1:2374846/1‑69 (MQ=255)
                                                        |                                                             
GGCAGCGTACCAAACCCGTATGTTAAAGCATCTGACTGGGGCTGGCAGATTGATCCAGTAGGTCTGCGCTATGCACTTTGCGAACTGTATGAGCGTTATCAGAGGCCGCTGTTTATTG  >  minE/1943938‑1944055

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: