Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,283,867 C→T 100% A1882A (GCC→GCT yeeJ → adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,283,8670CT83.5% 34.0 / 0.4 24A1882A (GCC→GCTyeeJadhesin
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base T (0/20);  total (4/20)
Fisher's exact test for biased strand distribution p-value = 9.41e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACAACTGGAGGCTATTGATGAAAAACTGACACTCACGGCGTCCAGTCCGCTTATCGGTGTCTATGCCCCTACAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAATGGCACTCCAGTGGAGGGTCAGGTC  >  minE/1283801‑1283933
                                                                  |                                                                  
acaacTGGAGGCTATTGATGAAAAACTGACACTCACGGCGTCCAGTCCGCTTATCGGTGTCTATGcccc                                                                  >  1:1059451/1‑69 (MQ=255)
acaacTGGAGGCTATTGATGAAAAACTGACACTCACGGCGTCCAGTCCGCTTATCGGTGTCTATGcccc                                                                  >  1:956656/1‑69 (MQ=255)
acaacTGGAGGCTATTGATGAAAAACTGACACTCACGGCGTCCAGTCCGCTTATCGGTGTCTATGcccc                                                                  >  1:497613/1‑69 (MQ=255)
acaacTGGAGGCTATTGATGAAAAACTGACACTCACGGCGTCCAGTCCGCTTATCGGTGTCTATGcccc                                                                  >  1:2394109/1‑69 (MQ=255)
acaacTGGAGGCTATTGATGAAAAACTGACACTCACGGCGTCCAGTCCGCTTATCGGTGTCTATGcccc                                                                  >  1:1025128/1‑69 (MQ=255)
acaacTGGAGGCTATTGATGAAAAACTGACACTCACGGCGTCCAGTCCGCTTATCGGTGTCTATGccc                                                                   >  1:3060845/1‑68 (MQ=255)
                         cTGACACTCACGGCGTCCAGTCCGCTTATCGGTGTCTATGCCCCTACAGGCGCTACTCTGACGGCAAc                                          >  1:78385/1‑68 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAATGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:3033954/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:2687106/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:543986/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:1122416/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:45325/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:383116/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:2937748/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:2915612/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:2871590/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:2618463/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:1825674/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:129814/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:1580979/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:1667608/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:1733708/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:2567730/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:1839301/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:1916408/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:2345439/69‑1 (MQ=255)
                                        tCCAGTCCGCTTATCGGTGTCAATGCTCCTAAAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAAt                          <  1:1144247/69‑1 (MQ=255)
                                                             tATGCCCCTACAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAATGGCACTCCAGTGGAGggtcag     >  1:1994599/1‑69 (MQ=255)
                                                             tATGCCCCTACAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAATGGCACTCCAGTGGAGggtcag     >  1:1766839/1‑69 (MQ=255)
                                                             tATGCCCCTACAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAATGGCACTCCAGTGGAGggtcag     >  1:419018/1‑69 (MQ=255)
                                                                 ccccTACAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAATGGCACt                    >  1:3255198/1‑50 (MQ=255)
                                                                 ccccTACAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAATGGCACTCCAGTGGAGggtcaggtc  >  1:3016655/1‑68 (MQ=255)
                                                                 ccccTACAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAATGGCACTCCAGTGGAGggtcaggtc  >  1:1161320/1‑68 (MQ=255)
                                                                 ccccTACAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAATGGCACTCCAGTGGAGggtcaggtc  >  1:1516141/1‑68 (MQ=255)
                                                                 ccccTACAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAATGGCACTCCAGTGGAGggtcaggt   >  1:2818802/1‑67 (MQ=255)
                                                                  |                                                                  
ACAACTGGAGGCTATTGATGAAAAACTGACACTCACGGCGTCCAGTCCGCTTATCGGTGTCTATGCCCCTACAGGCGCTACTCTGACGGCAACGCTAACCTCTGCAAATGGCACTCCAGTGGAGGGTCAGGTC  >  minE/1283801‑1283933

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: