Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,609,767 C→G 35.7% L136L (CTG→CTC nikD ← nickel transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,609,7670CG35.7% 45.9 / 45.0 56L136L (CTG→CTCnikDnickel transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (16/20);  new base G (11/9);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAAGCGTAGC  >  minE/2609707‑2609829
                                                            |                                                              
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        <  1:468615/69‑1 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        <  1:225641/69‑1 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        >  1:94267/1‑69 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        <  1:1999453/69‑1 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        <  1:2304220/69‑1 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        >  1:2339526/1‑69 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        >  1:2355009/1‑69 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        >  1:2409328/1‑69 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        <  1:782109/69‑1 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        >  1:781389/1‑69 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        <  1:3429846/69‑1 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        <  1:127422/69‑1 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        <  1:1264568/69‑1 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        >  1:592859/1‑69 (MQ=255)
caGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        >  1:451537/1‑69 (MQ=255)
 aGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAg                                                        <  1:804650/68‑1 (MQ=255)
             gCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCgcggcg                                           >  1:1385272/1‑69 (MQ=255)
             gCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCgcggcg                                           >  1:1749076/1‑69 (MQ=255)
               aaTCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGtt                                         >  1:17862/1‑69 (MQ=255)
                   atcatGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCAcgcg                                                  >  1:1056128/1‑56 (MQ=255)
                   atcatGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTcc                                     >  1:2287920/1‑69 (MQ=255)
                   atcatGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTcc                                     >  1:511671/1‑69 (MQ=255)
                   atcatGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTcc                                     >  1:1832696/1‑69 (MQ=255)
                   atcatGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTcc                                     >  1:600127/1‑69 (MQ=255)
                   atcatGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTcc                                     >  1:2523096/1‑69 (MQ=255)
                         cgcTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAGCCCGCGCGCGGCGTTTTCCa                                    >  1:3310202/1‑64 (MQ=255)
                         cgcTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAGCACGCGCGCGGCGTTTTCCAGcccc                               >  1:2293539/1‑69 (MQ=255)
                         cgcTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAGCACGCGCGCGGCGTTTTCCAGcccc                               >  1:2170482/1‑69 (MQ=255)
                         cgcTGCAACATGCCGCCGCTCATCTCGAACGGGTAGAGCTTCAGCACGCGCGCGGCGTTTTCCAGcccc                               >  1:2779066/1‑69 (MQ=255)
                                 cATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTc                       <  1:516312/69‑1 (MQ=255)
                                 cATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTc                       <  1:1150459/69‑1 (MQ=255)
                                 cATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTc                       <  1:1309723/69‑1 (MQ=255)
                                 cATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTc                       <  1:634025/69‑1 (MQ=255)
                                 cATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTc                       <  1:1489144/69‑1 (MQ=255)
                                 cATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTc                       <  1:2551121/69‑1 (MQ=255)
                                 cATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTc                       <  1:247368/69‑1 (MQ=255)
                                 cATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTc                       <  1:2382911/69‑1 (MQ=255)
                                             cATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGt           <  1:2739045/69‑1 (MQ=255)
                                             cATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGt           <  1:2833895/69‑1 (MQ=255)
                                             cATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGt           <  1:1378773/69‑1 (MQ=255)
                                             cATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGt           <  1:1577483/69‑1 (MQ=255)
                                             cATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGt           <  1:1612181/69‑1 (MQ=255)
                                             cATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGt           <  1:1243921/69‑1 (MQ=255)
                                             cATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGt           <  1:1958856/69‑1 (MQ=255)
                                                ctcGAACGGGTACAGCTTCAACACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGt           >  1:944853/1‑66 (MQ=255)
                                                 tcGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAag        >  1:918578/1‑68 (MQ=255)
                                                 tcGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAag        >  1:1100577/1‑68 (MQ=255)
                                                 tcGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAag        >  1:2808952/1‑68 (MQ=255)
                                                 tcGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAag        >  1:2372009/1‑68 (MQ=255)
                                                 tcGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAag        >  1:2035411/1‑68 (MQ=255)
                                                 tcGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAag        >  1:199757/1‑68 (MQ=255)
                                                     aCGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAagcgt     <  1:3401189/67‑1 (MQ=255)
                                                      cGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAagcgtagc  <  1:347131/69‑1 (MQ=255)
                                                      cGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAagcgtagc  <  1:335903/69‑1 (MQ=255)
                                                      cGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAagcgtagc  <  1:2535770/69‑1 (MQ=255)
                                                      cGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAagcgtagc  <  1:94820/69‑1 (MQ=255)
                                                            cAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGt           <  1:2954426/54‑1 (MQ=255)
                                                            |                                                              
CAGCACCGCCATCGCAATCATCATGCGCTGCAACATGCCGCCGCTCATCTCGAACGGGTACAGCTTCAGCACGCGCGCGGCGTTTTCCAGCCCCACCGCTTCTATGGCAGCGGTAAGCGTAGC  >  minE/2609707‑2609829

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: