Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,606,276 A→G 47.5% R289R (CGA→CGG amiA → N‑acetylmuramoyl‑l‑alanine amidase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,606,2760AG47.5% 2.6 / 71.4 81R289R (CGA→CGGamiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base A (2/40);  new base G (22/16);  total (25/56)
Fisher's exact test for biased strand distribution p-value = 2.29e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

CGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAT  >  minE/1606209‑1606338
                                                                   |                                                              
cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt                                                               <  1:2802082/69‑1 (MQ=255)
cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt                                                               <  1:313405/69‑1 (MQ=255)
cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt                                                               <  1:1638799/69‑1 (MQ=255)
cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt                                                               <  1:3359583/69‑1 (MQ=255)
cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt                                                               <  1:1910305/69‑1 (MQ=255)
cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt                                                               <  1:2582861/69‑1 (MQ=255)
cGATTGCTGAACGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt                                                               <  1:3076937/69‑1 (MQ=255)
 gATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt                                                               <  1:3187376/68‑1 (MQ=255)
    tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt                                                           <  1:708536/69‑1 (MQ=255)
    tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt                                                           <  1:1366595/69‑1 (MQ=255)
    tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt                                                           <  1:574779/69‑1 (MQ=255)
    tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt                                                           <  1:147396/69‑1 (MQ=255)
    tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt                                                           <  1:2498404/69‑1 (MQ=255)
    tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt                                                           <  1:2974145/69‑1 (MQ=255)
    tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt                                                           <  1:740420/69‑1 (MQ=255)
    tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt                                                           <  1:944122/69‑1 (MQ=255)
    tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt                                                           <  1:1796640/69‑1 (MQ=255)
    tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt                                                           <  1:235721/69‑1 (MQ=255)
           ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                    <  1:2051070/69‑1 (MQ=255)
           ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                    <  1:1881221/69‑1 (MQ=255)
           ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                    <  1:2134287/69‑1 (MQ=255)
           ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                    <  1:1748171/69‑1 (MQ=255)
           ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                    <  1:1533701/69‑1 (MQ=255)
           ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                    <  1:2135826/69‑1 (MQ=255)
           ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                    <  1:2358139/69‑1 (MQ=255)
             cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                    <  1:1301254/67‑1 (MQ=255)
              gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                 <  1:2373057/69‑1 (MQ=255)
              gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                 <  1:2807089/69‑1 (MQ=255)
              gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                 <  1:1072763/69‑1 (MQ=255)
              gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                 <  1:1104405/69‑1 (MQ=255)
               tGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                 <  1:2874272/68‑1 (MQ=255)
                    aGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGca                                           <  1:2061367/69‑1 (MQ=255)
                    aGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGca                                           <  1:1524868/69‑1 (MQ=255)
                      ttATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGca                                           <  1:494087/67‑1 (MQ=255)
                       tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACg                                             <  1:1976084/64‑1 (MQ=255)
                       tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACg                                             <  1:2049228/64‑1 (MQ=255)
                       tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACg                                             <  1:673727/64‑1 (MQ=255)
                       tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCacac                                        <  1:2186453/69‑1 (MQ=255)
                       tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCacac                                        <  1:7288/69‑1 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGtgatga                                                      >  1:965430/1‑46 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAGACCCGACGCACACCAGGTTaa                                >  1:1631841/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGcaca                                         >  1:1276619/1‑59 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:1175700/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:666148/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:1529633/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:3428131/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:3365461/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:1611271/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:3350073/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:3328248/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:1765700/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:283026/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:1989972/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:2815214/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:2401693/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:2272388/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:249986/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:2203143/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:2703380/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:2734313/1‑68 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTa                                 >  1:2378217/1‑67 (MQ=255)
                                tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGCTAAGTGATGAAACCCGACGCACACCAGGTTaa                                >  1:878926/1‑68 (MQ=255)
                                 ggTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAAc                              <  1:2904586/69‑1 (MQ=255)
                                     cGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCg                                                >  1:2788395/1‑47 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:1951301/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:1178808/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:932091/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:784282/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:124013/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:1315863/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:135215/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:1844939/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:612957/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:1413607/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:3000247/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:364641/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:3433238/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:2521342/68‑1 (MQ=255)
                                        cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:1764234/68‑1 (MQ=255)
                                         aaCCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                        <  1:2040196/67‑1 (MQ=255)
                                                             aaaGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAt  >  1:2735172/1‑69 (MQ=255)
                                                                   |                                                              
CGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAT  >  minE/1606209‑1606338

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: