Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,256,681 G→A 7.1% E22K (GAA→AAA)  yecO → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,256,6810GA 7.1% 178.0 / 3.3 56E22K (GAA→AAA) yecOpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (27/25);  new base A (2/2);  total (29/27)
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.34e-01

TATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTC  >  minE/1256617‑1256748
                                                                |                                                                   
tATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACggg                                                               >  1:1180772/1‑71 (MQ=255)
 aTGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACggg                                                               >  1:1985181/1‑70 (MQ=255)
      ccACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTg                                                         <  1:39778/70‑1 (MQ=255)
       cACCGCGACAGGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACgg                                                                <  1:2443449/63‑1 (MQ=255)
       cACCGCGACACGCTATTTTCTGCCCATATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTg                                                         >  1:2103497/1‑70 (MQ=255)
        aCCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTg                                                         <  1:4081437/69‑1 (MQ=255)
        aCCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGa                                                        <  1:938640/70‑1 (MQ=255)
        aCCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATAAACGGGTAGCTg                                                         >  1:1878899/1‑69 (MQ=255)
             gACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATAAACGGGTAGCTGAAGTCt                                                   <  1:2428511/70‑1 (MQ=255)
                acGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTccc                                               >  1:970242/1‑71 (MQ=255)
                      ttttCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATAAACGGGTAGCTGAAGTCTTCCCGGat                                           <  1:3961218/69‑1 (MQ=255)
                        ttCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGaa                                                                   >  1:4130394/1‑43 (MQ=255)
                        ttCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGCCTTCCCGGATATg                                        <  1:1113993/70‑1 (MQ=255)
                         tCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGa                                       <  1:3780241/70‑1 (MQ=255)
                          cTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATc                                     >  1:1430035/1‑71 (MQ=255)
                           tGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACgg                                                                <  1:1987322/43‑1 (MQ=255)
                           tGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATAAACGGGTAGCTGAAGTCTTCCCGGATATGATc                                     >  1:3461507/1‑70 (MQ=255)
                                cTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGt                               >  1:3021824/1‑71 (MQ=255)
                                cTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGt                               >  1:3893368/1‑71 (MQ=255)
                                 tATCGCCAGACTGGGCGACTGGACCTTTGATGAAc                                                                  <  1:885526/35‑1 (MQ=255)
                                 tATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGt                                                              <  1:2435302/39‑1 (MQ=255)
                                 tATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGt                                                              <  1:2484062/39‑1 (MQ=255)
                                  aTCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGaa                                                       <  1:1272815/45‑1 (MQ=255)
                                  aTCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGaa                                                       <  1:1165227/45‑1 (MQ=255)
                                   tCGCCAGACTGTGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAg                                  <  1:2555172/65‑1 (MQ=255)
                                   tCGCCAGACTGGGCGACTGGACCTTTGATGAACggg                                                               >  1:3547177/1‑36 (MQ=255)
                                   tCGCCAGACTGGGCGACTGGACCTTTGATGAACggg                                                               >  1:1580958/1‑36 (MQ=255)
                                   tCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGa                                                        >  1:596854/1‑43 (MQ=255)
                                   tCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTc                                                 >  1:1063469/1‑50 (MQ=255)
                                   tCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTcc                            >  1:420519/1‑71 (MQ=255)
                                    cGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCg                           >  1:3822151/1‑71 (MQ=255)
                                     gCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCg                                              <  1:2193810/51‑1 (MQ=255)
                                     gCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCg                           >  1:2650434/1‑70 (MQ=255)
                                     gCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCg                           >  1:2991810/1‑70 (MQ=255)
                                       cAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGat                                           <  1:1983480/52‑1 (MQ=255)
                                         gACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTAccc                      >  1:2881014/1‑71 (MQ=255)
                                          aCTGGGCGACTGGACCTTTGATGAACGGGTAGCTGa                                                        >  1:2582094/1‑36 (MQ=255)
                                          aCTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAg                                                      <  1:3799885/38‑1 (MQ=255)
                                          aCTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCg                                              <  1:475614/46‑1 (MQ=255)
                                            tGGGCGACTGGACCTTTGATGAACGGGTAGCTGaa                                                       >  1:1966309/1‑35 (MQ=255)
                                              ggCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGc                   <  1:2615551/69‑1 (MQ=255)
                                              ggCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCt                  <  1:1189594/70‑1 (MQ=255)
                                                cGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGc                   <  1:2295989/67‑1 (MQ=255)
                                                  aCTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTcc             <  1:238252/71‑1 (MQ=255)
                                                  aCTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTc              >  1:1045371/1‑70 (MQ=255)
                                                    tGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGc                   <  1:3075019/63‑1 (MQ=255)
                                                      gACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGt                          >  1:73987/1‑54 (MQ=255)
                                                      gACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCaa           <  1:3320229/69‑1 (MQ=255)
                                                       aCCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCTTACCCGGCTATTCCAatat        <  1:3623449/71‑1 (MQ=255)
                                                       aCCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAata         >  1:3846403/1‑70 (MQ=255)
                                                        ccTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAAtatt       <  1:637733/71‑1 (MQ=255)
                                                         cTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAAtatta      >  1:4133576/1‑71 (MQ=255)
                                                          tttGATGAACGGGTAGCTGAAGTCTTCCCGGatat                                         >  1:168812/1‑35 (MQ=255)
                                                           ttGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATattatt    >  1:338623/1‑71 (MQ=255)
                                                           ttGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATattatt    >  1:4114666/1‑71 (MQ=255)
                                                           ttGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATattatt    >  1:1043597/1‑71 (MQ=255)
                                                             gatgaACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTc  >  1:2348756/1‑71 (MQ=255)
                                                                |                                                                   
TATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATGATCCAGCGTTCCGTACCCGGCTATTCCAATATTATTTC  >  minE/1256617‑1256748

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: