Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,266,786 G→A 45.6% C207C (TGC→TGT yecA ← conserved metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,266,7860GA45.6% 6.1 / 34.4 46C207C (TGC→TGTyecAconserved metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (19/6);  new base A (13/8);  total (32/14)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCC  >  minE/1266723‑1266846
                                                               |                                                            
gACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAgg                                                         <  1:2827329/69‑1 (MQ=255)
                                 cTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGtt                                        >  1:1043721/1‑53 (MQ=255)
                                 cTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGAt                        >  1:85748/1‑69 (MQ=255)
                                      tAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGt                    <  1:373346/68‑1 (MQ=255)
                                      tAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGt                    <  1:3052364/68‑1 (MQ=255)
                                      tAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGt                    <  1:1477379/68‑1 (MQ=255)
                                      tAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGt                    <  1:217249/68‑1 (MQ=255)
                                      tAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGt                    <  1:213190/68‑1 (MQ=255)
                                           ttCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTg               >  1:2622521/1‑68 (MQ=255)
                                           ttCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTg               >  1:812292/1‑68 (MQ=255)
                                           ttCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTg               >  1:60339/1‑68 (MQ=255)
                                           ttCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTg               >  1:3342923/1‑68 (MQ=255)
                                           ttCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTg               >  1:2430014/1‑68 (MQ=255)
                                           ttCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTg               >  1:1428097/1‑68 (MQ=255)
                                            tCTTACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGaa             <  1:2078400/69‑1 (MQ=255)
                                            tCTTACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGaa             <  1:1950647/69‑1 (MQ=255)
                                            tCTTACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGaa             <  1:3381995/69‑1 (MQ=255)
                                            tCTTACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGaa             <  1:146434/69‑1 (MQ=255)
                                            tCTTACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGaa             <  1:2657612/69‑1 (MQ=255)
                                            tCTTACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGaa             <  1:2659363/69‑1 (MQ=255)
                                            tCTTACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGaa             <  1:2932500/69‑1 (MQ=255)
                                            tCTTACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGaa             <  1:1435000/69‑1 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCg                                           >  1:585807/1‑36 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTctct                                   >  1:2857294/1‑44 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg          >  1:1084652/1‑69 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg          >  1:788656/1‑69 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg          >  1:120530/1‑69 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg          >  1:1418349/1‑69 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg          >  1:305623/1‑69 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg          >  1:2712925/1‑69 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg          >  1:2735903/1‑69 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg          >  1:2590865/1‑69 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg          >  1:1687286/1‑69 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg          >  1:240443/1‑69 (MQ=255)
                                               tACCACTACCGCACGGACAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg          >  1:2317415/1‑69 (MQ=255)
                                                       ccGCACGGGCAAGGATCGTTACGCCCCTGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc  >  1:1433906/1‑69 (MQ=255)
                                                       ccGCACGGGCAAGGATCGTTACGCCCCGTTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc  >  1:1093105/1‑69 (MQ=255)
                                                       ccGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc  >  1:1989521/1‑69 (MQ=255)
                                                       ccGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc  >  1:746254/1‑69 (MQ=255)
                                                       ccGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc  >  1:726693/1‑69 (MQ=255)
                                                       ccGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc  >  1:2317736/1‑69 (MQ=255)
                                                       ccGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc  >  1:2385853/1‑69 (MQ=255)
                                                       ccGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc  >  1:572641/1‑69 (MQ=255)
                                                       ccGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc  >  1:434311/1‑69 (MQ=255)
                                                       ccGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc  >  1:296886/1‑69 (MQ=255)
                                                       ccGCACGGGCAAGGATCGTTACGCCACGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc  >  1:2614480/1‑69 (MQ=255)
                                                               |                                                            
GACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCC  >  minE/1266723‑1266846

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