Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,814,312 G→A 38.4% intergenic (‑397/‑33) yjcB ← / → yjcC predicted inner membrane protein/predicted signal transduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,814,3120GA38.4% 26.5 / 30.5 52intergenic (‑397/‑33)yjcB/yjcCpredicted inner membrane protein/predicted signal transduction protein
Reads supporting (aligned to +/- strand):  ref base G (12/20);  new base A (13/7);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 8.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

AGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATT  >  minE/2814249‑2814370
                                                               |                                                          
aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                                      <  1:3238308/69‑1 (MQ=255)
aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                                      <  1:2987433/69‑1 (MQ=255)
aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                                      <  1:1905066/69‑1 (MQ=255)
aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                                      <  1:2724986/69‑1 (MQ=255)
aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                                      <  1:2693249/69‑1 (MQ=255)
aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                                      <  1:578851/69‑1 (MQ=255)
         tataCCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGtata                                                  <  1:2711607/65‑1 (MQ=255)
         tataCCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGtata                                                  <  1:403885/65‑1 (MQ=255)
         tataCCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAAca                                              >  1:1488352/1‑69 (MQ=255)
                 aaGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGaa                                      >  1:659808/1‑69 (MQ=255)
                 aaGATTTATCTAAAGACGTCGGTACCAAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGaa                                      >  1:603210/1‑69 (MQ=255)
                                aCGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                                      <  1:461719/37‑1 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCATTGAGTCATc                  >  1:3080792/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:2281768/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:2594754/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:3108455/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:3187634/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:322369/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:3369078/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:355684/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:488769/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:1760979/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:1478700/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:656043/1‑69 (MQ=255)
                                    cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                  >  1:1021446/1‑69 (MQ=255)
                                        aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg               <  1:1119463/69‑1 (MQ=255)
                                        aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg               <  1:3097213/69‑1 (MQ=255)
                                        aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg               <  1:3093138/69‑1 (MQ=255)
                                        aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg               <  1:1976477/69‑1 (MQ=255)
                                                ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGc              >  1:2554444/1‑62 (MQ=255)
                                                ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca        >  1:908047/1‑68 (MQ=255)
                                                ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca        >  1:501819/1‑68 (MQ=255)
                                                ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca        >  1:464038/1‑68 (MQ=255)
                                                ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca        >  1:1969554/1‑68 (MQ=255)
                                                ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca        >  1:2074376/1‑68 (MQ=255)
                                                ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca        >  1:2289377/1‑68 (MQ=255)
                                                ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca        >  1:3316117/1‑68 (MQ=255)
                                                ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca        >  1:2432490/1‑68 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGATACCAAt   <  1:439765/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:3404393/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:3421042/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:458023/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:249436/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:2876235/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:1890428/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:2911339/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:1663858/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:1365494/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:1329331/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:1272227/68‑1 (MQ=255)
                                                     aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt   <  1:3002310/68‑1 (MQ=255)
                                                     aCCTTGAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAtt  <  1:2380075/68‑1 (MQ=255)
                                                               |                                                          
AGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATT  >  minE/2814249‑2814370

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: