Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,091,278 C→T 33.2% A141A (GCG→GCA)  sufB ← component of SufBCD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,091,2780CT33.2% 20.4 / 12.1 30A141A (GCG→GCA) sufBcomponent of SufBCD complex
Reads supporting (aligned to +/- strand):  ref base C (16/4);  new base T (5/5);  total (21/9)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

TGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTAC  >  minE/1091216‑1091342
                                                              |                                                                
tGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCAc                                                             <  1:190907/68‑1 (MQ=255)
             tttCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTcc                                                >  1:493202/1‑68 (MQ=255)
                            tGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCa                                >  1:49816/1‑69 (MQ=255)
                            tGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCa                                >  1:2142978/1‑69 (MQ=255)
                            tGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCa                                >  1:2242112/1‑69 (MQ=255)
                            tGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCa                                >  1:3190711/1‑69 (MQ=255)
                            tGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCa                                >  1:2866422/1‑69 (MQ=255)
                              gCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAAc                              <  1:3198548/69‑1 (MQ=255)
                              gCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAAc                              <  1:3052424/69‑1 (MQ=255)
                              gCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAAc                              <  1:2496527/69‑1 (MQ=255)
                                aaCCGAAACAGAGTCGAAAATGGCATCCACTGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACt                             <  1:405307/68‑1 (MQ=255)
                                aaCCGAAACAGAGTCGAAAATGGCATCCACTGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACt                             <  1:3290466/68‑1 (MQ=255)
                                aaCCGAAACAGAGTCGAAAATGGCATCCACTGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACt                             <  1:3251249/68‑1 (MQ=255)
                                aaCCGAAACAGAGTCGAAAATGGCATCCACTGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACt                             <  1:3213570/68‑1 (MQ=255)
                                              cGAAAATGGCATCCACTGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACgccgcc               >  1:523728/1‑68 (MQ=255)
                                              cGAAAATGGCATCCACTGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACgccgcc               >  1:2535669/1‑68 (MQ=255)
                                              cGAAAATGGCATCCACTGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACgccgcc               >  1:2457152/1‑68 (MQ=255)
                                              cGAAAATGGCATCCACTGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACgccgcc               >  1:131528/1‑68 (MQ=255)
                                              cGAAAATGGCATCCACTGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACgccgcc               >  1:1155284/1‑68 (MQ=255)
                                                     ggCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACctc       >  1:1148982/1‑69 (MQ=255)
                                                     ggCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACctc       >  1:1069043/1‑69 (MQ=255)
                                                      gCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACctc       >  1:2209829/1‑68 (MQ=255)
                                                          ccACTGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACt                             <  1:2448674/42‑1 (MQ=255)
                                                           cACCGCCACCTCTTTGCCTTCCCGCGCGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTag  >  1:2135415/1‑67 (MQ=255)
                                                           cACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAAc                              >  1:3203443/1‑40 (MQ=255)
                                                           cACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTAc  >  1:3107774/1‑68 (MQ=255)
                                                           cACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTAc  >  1:3234617/1‑68 (MQ=255)
                                                           cACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTAc  >  1:2124270/1‑68 (MQ=255)
                                                           cACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTAc  >  1:1611716/1‑68 (MQ=255)
                                                           cACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTAc  >  1:419226/1‑68 (MQ=255)
                                                              |                                                                
TGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTAC  >  minE/1091216‑1091342

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: ‑