Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,087,622 G→A 43.3% L39L (CTG→CTA yhaH → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,087,6220GA43.3% 15.6 / 43.4 60L39L (CTG→CTAyhaHpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (23/11);  new base A (22/4);  total (45/15)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

AAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTA  >  minE/2087567‑2087686
                                                       |                                                                
aaGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGt                                                     <  1:3400366/69‑1 (MQ=255)
aaGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGt                                                     <  1:2964988/69‑1 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTgg                                                              >  1:3148102/1‑59 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATaaaa                                                        >  1:2711065/1‑65 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:3094440/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:1949418/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:2274625/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:2579286/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:2742936/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:2903084/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:3024045/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:2838256/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:496913/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:555134/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:695515/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:841651/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:904113/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:1528762/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:1276352/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:1317976/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:134443/1‑69 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGt                                                     >  1:2024654/1‑68 (MQ=255)
 agagTACTGGATGTTTATTCTGGTCAACAGCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:3110696/1‑69 (MQ=39)
 agagTACTGGATGTTTATTCTGCTCAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGtt                                                    >  1:611062/1‑69 (MQ=39)
     tACTGGATGTTTATTCTGGTCAACATTATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGc                                                >  1:1337240/1‑69 (MQ=255)
     tACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGc                                                >  1:3411998/1‑69 (MQ=255)
     tACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGc                                                >  1:1474854/1‑69 (MQ=255)
     tACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGc                                                >  1:2957935/1‑69 (MQ=255)
     tACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGc                                                >  1:3208072/1‑69 (MQ=255)
     tACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGc                                                >  1:3312598/1‑69 (MQ=255)
     tACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGc                                                >  1:2648139/1‑69 (MQ=255)
     tACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGc                                                >  1:454029/1‑69 (MQ=255)
     tACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGc                                                >  1:2097921/1‑69 (MQ=255)
     tACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGc                                                >  1:862028/1‑69 (MQ=255)
     tACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGc                                                >  1:1126365/1‑69 (MQ=255)
         ggATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTgg                                             >  1:1629040/1‑68 (MQ=255)
         ggATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGc                                            >  1:3033915/1‑69 (MQ=255)
         ggATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGc                                            >  1:3029844/1‑69 (MQ=255)
         ggATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGc                                            >  1:2888437/1‑69 (MQ=255)
         ggATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGc                                            >  1:2221961/1‑69 (MQ=255)
         ggATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGc                                            >  1:1611453/1‑69 (MQ=255)
                tATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAAcgcg                                      >  1:3287471/1‑68 (MQ=255)
                tATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACgcgc                                     >  1:2597761/1‑69 (MQ=255)
                  ttCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCt                                    <  1:1680653/68‑1 (MQ=255)
                  ttCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCt                                    <  1:1624723/68‑1 (MQ=255)
                       gTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCtggtgg                               >  1:2009448/1‑68 (MQ=255)
                         cAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGTTAGGCTGGCAACGTGCCGGTGGCGa                            <  1:1615437/69‑1 (MQ=255)
                         cAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGTTAGGCTGGCAACGTGCCGGTGGCGa                            <  1:1974608/69‑1 (MQ=255)
                         cAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGTTAGGCTGGCAACGTGCCGGTGGCGa                            <  1:1208343/69‑1 (MQ=255)
                         cAACATCATCTTTACGTTCGTGCTGGGGCTACTGGATAAAATGTTAGGCTGGCAACGTGCCGGTGGCGa                            <  1:2239594/69‑1 (MQ=255)
                              tcatcTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAAGGCGATGGTGGCGAAGGCa                       >  1:51440/1‑69 (MQ=255)
                              tcatcTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGGTGGTGGCGAAGGCa                       >  1:890512/1‑69 (MQ=255)
                              tcatcTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCa                       <  1:1027391/69‑1 (MQ=255)
                              tcatcTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCa                       >  1:2417298/1‑69 (MQ=255)
                                       cGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTgacg               <  1:2853171/68‑1 (MQ=255)
                                       cGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTgacg               <  1:651725/68‑1 (MQ=255)
                                       cGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTgacg               <  1:3366574/68‑1 (MQ=255)
                                          tCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACaa           <  1:790511/69‑1 (MQ=255)
                                          tCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACaa           <  1:1487476/69‑1 (MQ=255)
                                          tCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCCGACGACaa           <  1:1264850/69‑1 (MQ=255)
                                                    gctgctGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTa  <  1:1682146/68‑1 (MQ=255)
                                                    gctgctGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTa  <  1:1399122/68‑1 (MQ=255)
                                                    gctgctGGATAAAATGTTAGGCTAGCAACGCGCTGGTGTCGAAGGCACCCTGACGACAATCTACGGTa  <  1:3045237/68‑1 (MQ=255)
                                                     ctgctgGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACaa           <  1:3395913/58‑1 (MQ=255)
                                                       |                                                                
AAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTA  >  minE/2087567‑2087686

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: