Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 866,934 A→G 62.7% L126L (CTT→CTC btuR ← cob(I)alamin adenolsyltransferase/cobinamide ATP‑dependent adenolsyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE866,9340AG62.7% 55.7 / 55.5 91L126L (CTT→CTCbtuRcob(I)alamin adenolsyltransferase/cobinamide ATP‑dependent adenolsyltransferase
Reads supporting (aligned to +/- strand):  ref base A (16/18);  new base G (50/7);  total (66/25)
Fisher's exact test for biased strand distribution p-value = 5.76e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

TAACGCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAGACTT  >  minE/866871‑866992
                                                               |                                                          
tAACGCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCaa                                                       >  1:2282401/1‑69 (MQ=255)
tAACGCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCaa                                                       >  1:17593/1‑69 (MQ=255)
tAACGCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCaa                                                       >  1:916833/1‑69 (MQ=255)
    gCCTGCACCACTTCTTCCAGTGTCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAAc                                                    <  1:2303478/68‑1 (MQ=255)
    gCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAAc                                                    <  1:311990/68‑1 (MQ=255)
    gCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAAc                                                    <  1:306122/68‑1 (MQ=255)
    gCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAAc                                                    <  1:2414063/68‑1 (MQ=255)
    gCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAAc                                                    <  1:2284739/68‑1 (MQ=255)
    gCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAAc                                                    <  1:947811/68‑1 (MQ=255)
    gCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAAc                                                    <  1:1929552/68‑1 (MQ=255)
    gCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAAc                                                    <  1:784892/68‑1 (MQ=255)
    gCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAAc                                                    <  1:154131/68‑1 (MQ=255)
    gCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAAc                                                    <  1:2044320/68‑1 (MQ=255)
           ccacTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAACCATATcc                                             <  1:489697/68‑1 (MQ=255)
           ccacTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAACCATATcc                                             <  1:1962987/68‑1 (MQ=255)
           ccacTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAACCATATcc                                             <  1:2294648/68‑1 (MQ=255)
           ccacTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAACCATATcc                                             <  1:1099887/68‑1 (MQ=255)
           ccacTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAACCATATcc                                             <  1:2841626/68‑1 (MQ=255)
           ccacTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAACCATATcc                                             <  1:3232822/68‑1 (MQ=255)
           ccacTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAACCATATcc                                             <  1:1204884/68‑1 (MQ=255)
           ccacTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGGCAGTTCATCAAGCAAAACCATATcc                                             <  1:3299430/68‑1 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:3430688/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:3405306/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:2328917/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:2699596/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:110861/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:329658/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:2791907/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:3290819/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:543577/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:2841248/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:2928682/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:3126250/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:1539237/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:2072989/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:850082/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:1867489/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:1170071/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:1863634/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:964228/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:1429816/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:1740115/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAcc                                                   >  1:1473511/1‑60 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAc                                                    >  1:1317939/1‑59 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAc                                                    >  1:3234861/1‑59 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAAc                                                    >  1:1747798/1‑59 (MQ=255)
             acTTCTTCCAGTGGCAAATAGTCATACGCCACCATATAAGTCAGTTCATCGAGCAAAAcc                                                   >  1:914360/1‑60 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:979762/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:586077/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:314932/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:845414/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:3377045/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:842195/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:81866/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:813911/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:536532/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:793278/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:7838/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:115876/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:2126246/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:3127620/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:1990398/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:2298704/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:2384398/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:2419517/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:1714026/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:2677294/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:2773614/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:2901066/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGa                                     >  1:2986317/1‑67 (MQ=255)
                    ccAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGAAt                                   >  1:3433280/1‑69 (MQ=255)
                                       cGCCACCATATACGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTg                 >  1:1771199/1‑68 (MQ=255)
                                            ccaTATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGAATCAGCCAGCATCCGCTTTGCATGtt           <  1:909641/69‑1 (MQ=255)
                                            ccaTATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGAATCAGCCAGCATCCGCTTTGCATGtt           <  1:1212032/69‑1 (MQ=255)
                                            ccaTATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGAATCAGCCAGCATCCGCTTTGCATGtt           <  1:2047815/69‑1 (MQ=255)
                                            ccaTATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGAATCAGCCAGCATCCGCTTTGCATGtt           <  1:660748/69‑1 (MQ=255)
                                            ccaTATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGAATCAGCCAGCATCCGCTTTGCATGtt           <  1:773743/69‑1 (MQ=255)
                                            ccaTATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGAATCAGCCAGCATCCGCTTTGCATGtt           <  1:2793432/69‑1 (MQ=255)
                                             caTATACGTCAGTTCATCGAGCAAAACCATATCCAGTGAGGAATCAGCCAGCATCCGCTTTGCATGtt           <  1:1803343/68‑1 (MQ=255)
                                                 taCGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAg      >  1:571593/1‑69 (MQ=255)
                                                 taCGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAg      >  1:80809/1‑69 (MQ=255)
                                                 taCGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAg      >  1:2546920/1‑69 (MQ=255)
                                                 taCGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAg      >  1:2606428/1‑69 (MQ=255)
                                                 taCGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAg      >  1:168575/1‑69 (MQ=255)
                                                 taCGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAg      >  1:1530893/1‑69 (MQ=255)
                                                 taCGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAg      >  1:1484145/1‑69 (MQ=255)
                                                 taCGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAg      >  1:3166288/1‑69 (MQ=255)
                                                   cGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAg      >  1:2717195/1‑67 (MQ=255)
                                                   cGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAGAc    >  1:2186743/1‑69 (MQ=255)
                                                   cGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAGAc    >  1:1797161/1‑69 (MQ=255)
                                                     tcggTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCAAGACtt  >  1:1465846/4‑69 (MQ=255)
                                                               |                                                          
TAACGCCTGCACCACTTCTTCCAGTGGCAAATAGTCATACGCCACCATATACGTCAGTTCATCAAGCAAAACCATATCCAGTGAGGAATCAGCAAGCATCCGCTTTGCATGTTGCCAGACTT  >  minE/866871‑866992

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: