Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 376,561 C→T 19.4% V176I (GTC→ATC)  fepD ← iron‑enterobactin transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE376,5610CT19.4% 75.6 / 9.3 51V176I (GTC→ATC) fepDiron‑enterobactin transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (26/15);  new base T (6/4);  total (32/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

TTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGC  >  minE/376496‑376623
                                                                 |                                                              
ttttAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTc                                                             <  1:1207137/69‑1 (MQ=255)
ttttAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTc                                                             <  1:2899981/69‑1 (MQ=255)
    aaGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAg                                                           <  1:277897/67‑1 (MQ=255)
    aaGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAg                                                           <  1:3297801/67‑1 (MQ=255)
    aaGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAg                                                           <  1:2918977/67‑1 (MQ=255)
    aaGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAg                                                           <  1:2741029/67‑1 (MQ=255)
    aaGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAg                                                           <  1:2682595/67‑1 (MQ=255)
    aaGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAg                                                           <  1:1521047/67‑1 (MQ=255)
    aaGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAg                                                           <  1:1579286/67‑1 (MQ=255)
            tAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGATGTCAGGATTAAGCAg                                                 >  1:1683720/1‑69 (MQ=255)
            tAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGATGTCAGGATTAAGCAg                                                 >  1:1545179/1‑69 (MQ=255)
            tAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGATGTCAGGATTAAGCAg                                                 >  1:252245/1‑69 (MQ=255)
            tAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGATGTCAGGATTAAGCAg                                                 >  1:697542/1‑69 (MQ=255)
            tAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGATGTCAGGATTAAGCAg                                                 >  1:194929/1‑69 (MQ=255)
            tAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGATGTCAGGATTAAGCAg                                                 >  1:1832563/1‑69 (MQ=255)
            tAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAg                                                 <  1:558617/69‑1 (MQ=255)
            tAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAg                                                 <  1:512408/69‑1 (MQ=255)
            tAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAg                                                 <  1:2437500/69‑1 (MQ=255)
            tAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAg                                                 <  1:1566767/69‑1 (MQ=255)
                ttGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGc                                              >  1:291442/1‑68 (MQ=255)
                ttGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGc                                              >  1:3408435/1‑68 (MQ=255)
                ttGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGc                                              >  1:2604807/1‑68 (MQ=255)
                ttGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGc                                              >  1:1679664/1‑68 (MQ=255)
                ttGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGc                                              >  1:1994884/1‑68 (MQ=255)
                ttGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGc                                              >  1:1304291/1‑68 (MQ=255)
                ttGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGc                                              >  1:28267/1‑68 (MQ=255)
                       tatCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAg                                                           <  1:3063126/48‑1 (MQ=255)
                       tatCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGc                                      >  1:3375597/1‑69 (MQ=255)
                       tatCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGc                                      >  1:825431/1‑69 (MQ=255)
                       tatCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGc                                      >  1:2617697/1‑69 (MQ=255)
                       tatCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGc                                      >  1:225576/1‑69 (MQ=255)
                         tCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTg                                    >  1:2870207/1‑69 (MQ=255)
                                  ccGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCt                           <  1:2898874/69‑1 (MQ=255)
                                                  gCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGt                              >  1:947342/1‑50 (MQ=255)
                                                  gCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGt                              >  1:839425/1‑50 (MQ=255)
                                                  gCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGt                              >  1:3302250/1‑50 (MQ=255)
                                                  gCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGt                              >  1:2726632/1‑50 (MQ=255)
                                                  gCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGt                              >  1:717975/1‑50 (MQ=255)
                                                  gCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGt                              >  1:2116292/1‑50 (MQ=255)
                                                  gCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGt                              >  1:603600/1‑50 (MQ=255)
                                                  gCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGt                              >  1:3014723/1‑50 (MQ=255)
                                                      cTGATCGTAGATGTCAGGATTAAGCAGGGCGATGCCGCTGGTCGGTCCTTCCAGCACCGCCGCCAgcgc       <  1:1016150/69‑1 (MQ=255)
                                                      cTGATCGTAGATGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAgcgc       <  1:2434572/69‑1 (MQ=255)
                                                      cTGATCGTAGATGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAgcgc       <  1:384169/69‑1 (MQ=255)
                                                      cTGATCGTAGATGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAgcgc       <  1:2469683/69‑1 (MQ=255)
                                                           cGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCAcgc  >  1:2237480/1‑69 (MQ=255)
                                                           cGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCAcgc  >  1:616433/1‑69 (MQ=255)
                                                           cGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCAcgc  >  1:2871486/1‑69 (MQ=255)
                                                           cGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCAcg   >  1:1902622/1‑68 (MQ=255)
                                                           cGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCAcg   >  1:3301283/1‑68 (MQ=255)
                                                           cGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCAAGCACCGCCGCCAGCGCCAcgc  >  1:385433/1‑69 (MQ=255)
                                                                 |                                                              
TTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGC  >  minE/376496‑376623

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: