Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,932,021 C→T 37.2% M203I (ATG→ATA yqeF ← predicted acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,932,0210CT37.2% 35.1 / 33.8 59M203I (ATG→ATAyqeFpredicted acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/11);  new base T (7/15);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 6.42e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTTTGTTGCGAGCT  >  minE/1931970‑1932089
                                                   |                                                                    
gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                                     <  1:3199079/69‑1 (MQ=255)
gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                                     >  1:2349006/1‑69 (MQ=255)
gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                                     >  1:250401/1‑69 (MQ=255)
gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                                     >  1:427459/1‑69 (MQ=255)
gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                                     <  1:747990/69‑1 (MQ=255)
gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                                     >  1:151746/1‑69 (MQ=255)
gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                                     >  1:39248/1‑69 (MQ=255)
gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtctc                                                     >  1:1722490/1‑69 (MQ=255)
gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtct                                                      >  1:2178196/1‑68 (MQ=255)
gCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtct                                                      >  1:134904/1‑68 (MQ=255)
   tGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATc                                                  >  1:2154103/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:2762128/1‑69 (MQ=255)
           cATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATc                                          >  1:3071522/1‑69 (MQ=255)
           cATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATc                                          >  1:923720/1‑69 (MQ=255)
             tCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCg                                         <  1:2246271/68‑1 (MQ=255)
              cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                                          <  1:2682106/50‑1 (MQ=39)
              cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                                          <  1:262510/50‑1 (MQ=39)
              cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                                          <  1:2551178/50‑1 (MQ=39)
              cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                                          <  1:2866277/50‑1 (MQ=39)
              cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                                          <  1:3364621/50‑1 (MQ=39)
              cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                                          <  1:2316980/50‑1 (MQ=39)
              cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                                          <  1:3308660/50‑1 (MQ=39)
              cGGTATCAACAACCAACGTCTGTCCGTTACTTTGGGTTATTACCGGGACg                                                          <  1:1364265/50‑1 (MQ=39)
              cGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTc                                       >  1:3373029/1‑69 (MQ=255)
              cGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTc                                       >  1:746652/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:482295/1‑69 (MQ=255)
                     aacaacCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTc                                >  1:748529/1‑69 (MQ=255)
                     aacaacCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTc                                >  1:1832502/1‑69 (MQ=255)
                     aacaacCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTc                                >  1:2390645/1‑69 (MQ=255)
                       caaccaacGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCa                               <  1:2684635/68‑1 (MQ=255)
                       caaccaacGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCa                               <  1:2114381/68‑1 (MQ=255)
                       caaccaacGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCa                               <  1:2123779/68‑1 (MQ=255)
                       caaccaacGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCa                               <  1:454802/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:1151132/1‑67 (MQ=255)
                                tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc                       >  1:3246543/1‑67 (MQ=255)
                                tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc                       >  1:1233664/1‑67 (MQ=255)
                                tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc                       >  1:257162/1‑67 (MQ=255)
                                tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc                       >  1:1478944/1‑67 (MQ=255)
                                tCTGTCCGTTACTTTGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATcgccgc                       >  1:1599303/1‑67 (MQ=255)
                                      cGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCt                <  1:2731148/68‑1 (MQ=255)
                                      cGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCt                <  1:1773729/68‑1 (MQ=255)
                                            tttGGGTTATTACCGGGACGATCTCATCTTTAAATGGTCCGGCGTCAATCGCCGCTCGCGCTTTttgt          <  1:2757249/68‑1 (MQ=255)
                                            tttGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTttgt          <  1:3138924/68‑1 (MQ=255)
                                            tttGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTttgt          <  1:731928/68‑1 (MQ=255)
                                            tttGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTttgt          <  1:652886/68‑1 (MQ=255)
                                            tttGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTttgt          <  1:204381/68‑1 (MQ=255)
                                            tttGGGTTATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTttgt          <  1:2671829/68‑1 (MQ=255)
                                                ggTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTTTGTTGCg      <  1:1043784/68‑1 (MQ=255)
                                                   cATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTTTGTTGCGAGCt  <  1:1763255/69‑1 (MQ=255)
                                                   cATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTTTGTTGCGAGCt  <  1:222273/69‑1 (MQ=255)
                                                   cATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTTTGTTGCGAGCt  <  1:2272201/69‑1 (MQ=255)
                                                   cATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTTTGTTGCGAGCt  <  1:141321/69‑1 (MQ=255)
                                                   cATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTTTGTTGCGAGCt  <  1:963143/69‑1 (MQ=255)
                                                   |                                                                    
GCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTCGCGCTTTTTGTTGCGAGCT  >  minE/1931970‑1932089

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: