Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,814,290 C→G 59.7% intergenic (‑375/‑55) 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,2900CG59.7% 29.0 / 54.1 57intergenic (‑375/‑55)yjcB/yjcCpredicted inner membrane protein/predicted signal transduction protein
Reads supporting (aligned to +/- strand):  ref base C (14/9);  new base G (27/7);  total (41/16)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

CAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTG  >  minE/2814241‑2814357
                                                 |                                                                   
cAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCa                                                                  >  1:939433/1‑53 (MQ=255)
cAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGtttt                                                           >  1:961534/1‑60 (MQ=255)
cAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGtttt                                                           >  1:3380283/1‑60 (MQ=255)
cAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGtttt                                                           >  1:1871432/1‑60 (MQ=255)
cAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGtttt                                                           >  1:419139/1‑60 (MQ=255)
cAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGtttt                                                           >  1:1295836/1‑60 (MQ=255)
cAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGtttt                                                           >  1:1158801/1‑60 (MQ=255)
cAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCa                                                  >  1:1552646/1‑69 (MQ=255)
cAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCa                                                  >  1:257376/1‑69 (MQ=255)
cAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCa                                                  >  1:3232746/1‑69 (MQ=255)
cAGGCAAAAGTGATATTTATACCGTAAGATTTATCTAAAGACGTCGGTAc                                                                     >  1:3189635/1‑50 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:1415400/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:3035270/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:270162/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:3095417/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:2459305/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:2311706/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:2239490/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:1997933/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:1995520/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:324184/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:1227565/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:532623/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:621900/1‑56 (MQ=255)
 aGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGt                                                              >  1:1342838/1‑56 (MQ=255)
        aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                         <  1:3238308/69‑1 (MQ=255)
        aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                         <  1:2987433/69‑1 (MQ=255)
        aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                         <  1:578851/69‑1 (MQ=255)
        aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                         <  1:2724986/69‑1 (MQ=255)
        aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                         <  1:1905066/69‑1 (MQ=255)
        aGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                         <  1:2693249/69‑1 (MQ=255)
                 tataCCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGtata                                     <  1:403885/65‑1 (MQ=255)
                 tataCCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGtata                                     <  1:2711607/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)
                              ttATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCa                                                  >  1:1515769/1‑39 (MQ=255)
                              ttATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCa                                                  <  1:1978052/39‑1 (MQ=255)
                              ttATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCa                                                  <  1:465300/39‑1 (MQ=255)
                               tATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCa                                                  <  1:1644111/38‑1 (MQ=255)
                                        aCGTCGGTAGCCAGGGTTTCACCTTGCAATAACCggg                                         <  1:461719/37‑1 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCATTGAGTCATc     >  1:3080792/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:1478700/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:656043/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:488769/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:1021446/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:1760979/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:355684/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:3369078/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:2281768/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:322369/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:2594754/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:3187634/1‑69 (MQ=255)
                                            cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc     >  1:3108455/1‑69 (MQ=255)
                                                aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg  <  1:1976477/69‑1 (MQ=255)
                                                aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg  <  1:3097213/69‑1 (MQ=255)
                                                aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg  <  1:3093138/69‑1 (MQ=255)
                                                aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg  <  1:1119463/69‑1 (MQ=255)
                                                 |                                                                   
CAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTG  >  minE/2814241‑2814357

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: