Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,107,580 Δ1 bp 8.6% coding (518/1155 nt) agaS → tagatose‑6‑phosphate ketose/aldose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,107,5800C. 8.6% 228.2 / 12.1 70coding (518/1155 nt)agaStagatose‑6‑phosphate ketose/aldose isomerase
Reads supporting (aligned to +/- strand):  ref base C (28/36);  new base . (2/4);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

CAACAGCGATAACGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTGAGACGATCAACAGCCAAAC  >  minE/2107515‑2107649
                                                                 |                                                                     
cAACAGCGATAACGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACcag                                                                    >  1:2743016/1‑69 (MQ=255)
 aaCAGCGATAACGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACcagcag                                                                 <  1:3902207/71‑1 (MQ=255)
  aCAGCGATAACGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCagcagc                                                                <  1:3885543/71‑1 (MQ=255)
  aCAGCGATAACGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCagcagc                                                                <  1:1816684/71‑1 (MQ=255)
   cAGCGATAACGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAgcagca                                                               >  1:2112192/1‑71 (MQ=255)
   cAGCGATAACGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAgcagca                                                               >  1:352058/1‑71 (MQ=255)
       gATAACGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTac                                                           <  1:866249/71‑1 (MQ=255)
         tAACGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTacca                                                         >  1:1087611/1‑71 (MQ=255)
         tAACGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACAGCAGCATTaccac                                                        <  1:3712311/71‑1 (MQ=255)
           aCGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTaccacc                                                       <  1:1757361/71‑1 (MQ=255)
           aCGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTaccacc                                                       <  1:1538669/71‑1 (MQ=255)
             gcgTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCat                                                     <  1:4024413/71‑1 (MQ=255)
             gcgTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCat                                                     <  1:1911042/71‑1 (MQ=255)
               gTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCat                                                     <  1:3964107/69‑1 (MQ=255)
                tttGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCatga                                                   >  1:2889463/1‑70 (MQ=255)
                  tGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATgg                                                >  1:1358065/1‑71 (MQ=255)
                  tGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACAgcagca                                                               >  1:3203737/1‑55 (MQ=255)
                     ccTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTaccac                                                        >  1:542018/1‑60 (MQ=255)
                     ccTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCat                                                     <  1:3729958/63‑1 (MQ=255)
                       tgctgATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCat                                                     <  1:841718/61‑1 (MQ=255)
                       tgctgATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGc                                           >  1:2354246/1‑71 (MQ=255)
                        gctgATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCt                                          >  1:282291/1‑71 (MQ=255)
                        gctgATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCt                                          >  1:2029299/1‑71 (MQ=255)
                         ctgATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGAc                                                                       >  1:2840326/1‑41 (MQ=255)
                         ctgATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCatga                                                   >  1:1582890/1‑61 (MQ=255)
                         ctgATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCt                                          >  1:4061356/1‑70 (MQ=255)
                             tGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCat                                                     <  1:2080994/55‑1 (MQ=255)
                              gCCCGCAGAAACGGACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCatga                                                   >  1:548341/1‑56 (MQ=255)
                              gCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACAGCAGCATTACCACCATGATGGCCAGCTGCCTcgc                                   >  1:793527/1‑71 (MQ=255)
                               cccGCAGAAACGCACGATCGCGGCTTTGCGATGACCagcagc                                                                <  1:1027772/42‑1 (MQ=255)
                                 cGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCat                                                     <  1:396581/51‑1 (MQ=255)
                                  gCAGAAACGCACGATCGCGGCTTTGCGATGACAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGtt                               <  1:2611147/71‑1 (MQ=255)
                                   cAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTacca                                                         >  1:3614039/1‑45 (MQ=255)
                                   cAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTacca                                                         >  1:3749513/1‑45 (MQ=255)
                                   cAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGtt                               <  1:2723828/71‑1 (MQ=255)
                                   cAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGtt                               >  1:3776482/1‑71 (MQ=255)
                                    aGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGttt                              >  1:2779330/1‑71 (MQ=255)
                                     gAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCatg                                                    <  1:4098732/48‑1 (MQ=255)
                                      aaaCGCACGATCGCGGCTTTGCGATGACAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCg                           <  1:3206515/71‑1 (MQ=255)
                                       aaCGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGc                                               <  1:3211708/51‑1 (MQ=255)
                                       aaCGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCg                           >  1:1086344/1‑71 (MQ=255)
                                       aaCGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCg                           <  1:1739897/71‑1 (MQ=255)
                                        aCGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCatga                                                   >  1:829153/1‑46 (MQ=255)
                                        aCGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGc                                           >  1:298750/1‑54 (MQ=255)
                                        aCGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGc                                           >  1:3830217/1‑54 (MQ=255)
                                        aCGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGtt                               <  1:224720/66‑1 (MQ=255)
                                        aCGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGtt                               >  1:1176014/1‑66 (MQ=255)
                                        aCGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGtt                               >  1:3389671/1‑66 (MQ=255)
                                        aCGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGc                          <  1:1983767/71‑1 (MQ=255)
                                         cGCACGATCGCGGCTTTGCGATGACCAGCAGCAtt                                                             <  1:812596/35‑1 (MQ=255)
                                           cACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCat                                                     <  1:2484184/41‑1 (MQ=255)
                                            aCGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCt                      <  1:2867371/71‑1 (MQ=255)
                                            aCGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCt                      <  1:1272960/71‑1 (MQ=255)
                                            aCGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCt                      <  1:1276449/71‑1 (MQ=255)
                                             cGATCGCGGCTTTGCGATGACAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTga                    <  1:2511651/71‑1 (MQ=255)
                                              gATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCt                      >  1:2020251/1‑69 (MQ=255)
                                              gATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTga                    <  1:867055/71‑1 (MQ=255)
                                               aTCGCGGCTTTGCGATGACCAGCAGCATTACCACCat                                                     <  1:3308578/37‑1 (MQ=255)
                                               aTCGCGGCTTTGCGATGACCAGCAGCATTACCACCat                                                     <  1:2749906/37‑1 (MQ=255)
                                               aTCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTgag                   <  1:2539405/71‑1 (MQ=255)
                                               aTCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTgag                   <  1:2095942/71‑1 (MQ=255)
                                                tCGCGGCTTTGCGATGACCAGCAGCATTACCACCatga                                                   >  1:2917036/1‑38 (MQ=255)
                                                tCGCGGCTTTGCGATGACCAGCAGCATTACCACCatga                                                   >  1:1586726/1‑38 (MQ=255)
                                                tCGCGGCTTTGCGATGACCAGCAGCATTACCACCatga                                                   >  1:767564/1‑38 (MQ=255)
                                                    ggCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTGAGACGAt              <  1:2382148/71‑1 (MQ=255)
                                                     gCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTGAGACGATc             <  1:649323/71‑1 (MQ=255)
                                                         tGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTGAGACGATCaa           <  1:3159680/69‑1 (MQ=255)
                                                            gATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTgag                   <  1:2406790/58‑1 (MQ=255)
                                                              tGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTGAGACGATCAACAGCCaa    <  1:3489366/71‑1 (MQ=255)
                                                              tGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTGAGACGATCAACAGCCaa    <  1:3484750/71‑1 (MQ=255)
                                                                aCCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTGAGACGATCAACAGCCAAAc  >  1:2245116/1‑71 (MQ=255)
                                                                 |                                                                     
CAACAGCGATAACGCGTTTGCCCTGCTGATGCCCGCAGAAACGCACGATCGCGGCTTTGCGATGACCAGCAGCATTACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTGAGACGATCAACAGCCAAAC  >  minE/2107515‑2107649

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