Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,814,333 A→G 31.7% intergenic (‑418/‑12) 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,3330AG31.7% 52.1 / 29.6 63intergenic (‑418/‑12)yjcB/yjcCpredicted inner membrane protein/predicted signal transduction protein
Reads supporting (aligned to +/- strand):  ref base A (21/22);  new base G (13/7);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTT  >  minE/2814266‑2814401
                                                                   |                                                                    
aaGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGaa                                                                     >  1:659808/1‑69 (MQ=255)
aaGATTTATCTAAAGACGTCGGTACCAAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGaa                                                                     >  1:603210/1‑69 (MQ=255)
                   cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCATTGAGTCATc                                                 >  1:3080792/1‑69 (MQ=255)
                   cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                                                 >  1:1021446/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:1760979/1‑69 (MQ=255)
                   cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                                                 >  1:355684/1‑69 (MQ=255)
                   cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                                                 >  1:3187634/1‑69 (MQ=255)
                   cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                                                 >  1:3108455/1‑69 (MQ=255)
                   cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                                                 >  1:1478700/1‑69 (MQ=255)
                   cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                                                 >  1:488769/1‑69 (MQ=255)
                   cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                                                 >  1:2594754/1‑69 (MQ=255)
                   cGGTAGCCAGGGTTTCACCTTGCAATAACCGGGTATAAACAGGCAGGGAATTGATAGCAATGAGTCATc                                                 >  1:656043/1‑69 (MQ=255)
                       aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg                                              <  1:3097213/69‑1 (MQ=255)
                       aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg                                              <  1:3093138/69‑1 (MQ=255)
                       aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg                                              <  1:1976477/69‑1 (MQ=255)
                       aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTg                                              <  1:1119463/69‑1 (MQ=255)
                               ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGc                                             >  1:2554444/1‑62 (MQ=255)
                               ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca                                       >  1:2074376/1‑68 (MQ=255)
                               ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca                                       >  1:3316117/1‑68 (MQ=255)
                               ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca                                       >  1:1969554/1‑68 (MQ=255)
                               ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca                                       >  1:2289377/1‑68 (MQ=255)
                               ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca                                       >  1:2432490/1‑68 (MQ=255)
                               ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca                                       >  1:464038/1‑68 (MQ=255)
                               ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca                                       >  1:501819/1‑68 (MQ=255)
                               ttttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGaca                                       >  1:908047/1‑68 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGATACCAAt                                  <  1:439765/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:3404393/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:2876235/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:3421042/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:458023/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:3002310/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:2911339/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:1272227/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:249436/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:1329331/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:1365494/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:1663858/68‑1 (MQ=255)
                                    aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAt                                  <  1:1890428/68‑1 (MQ=255)
                                    aCCTTGAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAtt                                 <  1:2380075/68‑1 (MQ=255)
                                                       taAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCa                >  1:1013706/1‑67 (MQ=255)
                                                       taAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCa                >  1:3384393/1‑67 (MQ=255)
                                                       taAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCAtt              >  1:915771/1‑69 (MQ=255)
                                                       taAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCAtt              >  1:76484/1‑69 (MQ=255)
                                                       taAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCAtt              >  1:742937/1‑69 (MQ=255)
                                                       taAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCAtt              >  1:3234347/1‑69 (MQ=255)
                                                       taAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCAt               >  1:2005615/1‑68 (MQ=255)
                                                       taAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCAt               >  1:2624632/1‑68 (MQ=255)
                                                       taAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCAt               >  1:3434512/1‑68 (MQ=255)
                                                         aaCAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTa             <  1:245130/68‑1 (MQ=255)
                                                         aaCAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTa             <  1:1734054/68‑1 (MQ=255)
                                                           caggcaggAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCt          >  1:2242775/1‑69 (MQ=255)
                                                           caggcaggAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCt          >  1:516627/1‑69 (MQ=255)
                                                           caggcaggAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCt          >  1:1180577/1‑69 (MQ=255)
                                                             ggcaggGAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTTCCGGGCATTATCttt        <  1:1727869/69‑1 (MQ=255)
                                                             ggcaggGAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTTCCGGGCATTATCttt        <  1:1334811/69‑1 (MQ=255)
                                                             ggcaggGAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTTCCGGGCATTATCttt        <  1:2132565/69‑1 (MQ=255)
                                                             ggcaggGAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTTCCGGGCATTATCttt        <  1:2515618/69‑1 (MQ=255)
                                                             ggcaggGAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTTCCGGGCATTATCttt        <  1:1237195/69‑1 (MQ=255)
                                                             ggcaggGAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTTCCGGGCATTATCttt        <  1:1150389/69‑1 (MQ=255)
                                                             ggcaggGAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTTCCGGGCATTATCttt        <  1:3032852/69‑1 (MQ=255)
                                                              gcaggAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCtttt       <  1:724717/69‑1 (MQ=255)
                                                              gcaggAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCtttt       <  1:3353566/69‑1 (MQ=255)
                                                                   aaaTTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGtt  >  1:2810907/1‑69 (MQ=255)
                                                                   |                                                                    
AAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTT  >  minE/2814266‑2814401

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: