Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,932,006 C→T 30.7% G208G (GGG→GGA yqeF ← predicted acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,932,0060CT30.7% 46.8 / 22.1 52G208G (GGG→GGAyqeFpredicted acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/10);  new base T (7/9);  total (33/19)
Fisher's exact test for biased strand distribution p-value = 6.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGC  >  minE/1931940‑1932068
                                                                  |                                                              
aCGGGCTAAGCCTTCTGCGCTTGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGcc                                                               <  1:464588/68‑1 (MQ=255)
aCGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGcc                                                               <  1:3096912/68‑1 (MQ=255)
aCGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGcc                                                               <  1:2114209/68‑1 (MQ=255)
aCGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGcc                                                               <  1:2069546/68‑1 (MQ=255)
aCGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGcc                                                               <  1:2016785/68‑1 (MQ=255)
aCGGGCTAAGCCTTCTGCCCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGcc                                                               <  1:2793794/68‑1 (MQ=255)
     cTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTAc                                                         <  1:621088/69‑1 (MQ=255)
     cTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTAc                                                         <  1:3350887/69‑1 (MQ=255)
                              gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                <  1:3199079/69‑1 (MQ=255)
                              gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                >  1:2349006/1‑69 (MQ=255)
                              gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                >  1:1722490/1‑69 (MQ=255)
                              gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                >  1:250401/1‑69 (MQ=255)
                              gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                >  1:151746/1‑69 (MQ=255)
                              gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                >  1:39248/1‑69 (MQ=255)
                              gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                >  1:427459/1‑69 (MQ=255)
                              gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                <  1:747990/69‑1 (MQ=255)
                              gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtct                                 >  1:134904/1‑68 (MQ=255)
                              gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtct                                 >  1:2178196/1‑68 (MQ=255)
                                 tGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATc                             >  1:2154103/1‑69 (MQ=255)
                                         cATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATc                     >  1:3071522/1‑69 (MQ=255)
                                         cATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATc                     >  1:2762128/1‑69 (MQ=255)
                                         cATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATc                     >  1:2300573/1‑69 (MQ=255)
                                         cATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATc                     >  1:2211021/1‑69 (MQ=255)
                                         cATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATc                     >  1:3418983/1‑69 (MQ=255)
                                         cATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATc                     >  1:923720/1‑69 (MQ=255)
                                           tCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCg                    <  1:2246271/68‑1 (MQ=255)
                                            cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                     <  1:3308660/50‑1 (MQ=39)
                                            cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                     <  1:3364621/50‑1 (MQ=39)
                                            cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                     <  1:2866277/50‑1 (MQ=39)
                                            cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                     <  1:262510/50‑1 (MQ=39)
                                            cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                     <  1:2682106/50‑1 (MQ=39)
                                            cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                     <  1:2551178/50‑1 (MQ=39)
                                            cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                     <  1:2316980/50‑1 (MQ=39)
                                            cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                     <  1:1364265/50‑1 (MQ=39)
                                            cGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTc                  >  1:746652/1‑69 (MQ=255)
                                            cGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTc                  >  1:3373029/1‑69 (MQ=255)
                                              gTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCg                >  1:1343158/1‑69 (MQ=255)
                                              gTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCg                >  1:568748/1‑69 (MQ=255)
                                                   aacaacCAACGTCTGCCCGTTACTTTGGGTCATTa                                             >  1:2326321/1‑35 (MQ=255)
                                                   aacaacCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCTTc           >  1:1357229/1‑69 (MQ=255)
                                                   aacaacCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGt            >  1:769913/1‑68 (MQ=255)
                                                   aacaacCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTc           >  1:2390645/1‑69 (MQ=255)
                                                   aacaacCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTc           >  1:482295/1‑69 (MQ=255)
                                                   aacaacCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTc           >  1:1832502/1‑69 (MQ=255)
                                                   aacaacCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTc           >  1:748529/1‑69 (MQ=255)
                                                     caaccaacGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCa          <  1:2684635/68‑1 (MQ=255)
                                                     caaccaacGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCa          <  1:454802/68‑1 (MQ=255)
                                                     caaccaacGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCa          <  1:2114381/68‑1 (MQ=255)
                                                     caaccaacGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCa          <  1:2123779/68‑1 (MQ=255)
                                                     caaccaacGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCa          <  1:860766/68‑1 (MQ=255)
                                                          aacGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                     <  1:1298125/36‑1 (MQ=37)
                                                              tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc  >  1:778704/1‑67 (MQ=255)
                                                              tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc  >  1:1233664/1‑67 (MQ=255)
                                                              tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc  >  1:1478944/1‑67 (MQ=255)
                                                              tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc  >  1:1599303/1‑67 (MQ=255)
                                                              tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc  >  1:3246543/1‑67 (MQ=255)
                                                              tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc  >  1:257162/1‑67 (MQ=255)
                                                              tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc  >  1:1151132/1‑67 (MQ=255)
                                                                  |                                                              
ACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGC  >  minE/1931940‑1932068

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: