Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 303,058 G→A 32.0% D191D (GAC→GAT tesB ← acyl‑CoA thioesterase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE303,0580GA32.0% 74.3 / 40.1 84D191D (GAC→GATtesBacyl‑CoA thioesterase II
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base A (12/15);  total (36/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAC  >  minE/303007‑303119
                                                   |                                                             
cAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTAcc                                              <  1:2867206/69‑1 (MQ=255)
        tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:968107/1‑68 (MQ=255)
        tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:861970/1‑68 (MQ=255)
        tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:491556/1‑68 (MQ=255)
        tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:419406/1‑68 (MQ=255)
        tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:369033/1‑68 (MQ=255)
        tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:3228233/1‑68 (MQ=255)
        tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:1792018/1‑68 (MQ=255)
        tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:2088219/1‑68 (MQ=255)
        tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:2696236/1‑68 (MQ=255)
        tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:1506718/1‑68 (MQ=255)
        tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:2280196/1‑68 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCAt                                            >  1:746174/1‑42 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:2035747/1‑67 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:1093473/1‑67 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:2185183/1‑67 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:2402553/1‑67 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:2540670/1‑67 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:2993874/1‑67 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:3118205/1‑67 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:3237213/1‑67 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:1101636/1‑67 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:1420078/1‑67 (MQ=255)
                             gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:1622812/1‑67 (MQ=255)
                               aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:2942130/69‑1 (MQ=255)
                               aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:1068201/69‑1 (MQ=255)
                               aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:1112356/69‑1 (MQ=255)
                               aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:3147161/69‑1 (MQ=255)
                               aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:2945936/69‑1 (MQ=255)
                               aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:2906028/69‑1 (MQ=255)
                               aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:2713862/69‑1 (MQ=255)
                               aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:2177785/69‑1 (MQ=255)
                                   aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:2223021/1‑68 (MQ=255)
                                   aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:1119662/1‑68 (MQ=255)
                                   aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:1240142/1‑68 (MQ=255)
                                   aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:3105058/1‑68 (MQ=255)
                                   aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:2878620/1‑68 (MQ=255)
                                   aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:2061708/1‑68 (MQ=255)
                                   aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:150081/1‑68 (MQ=255)
                                   aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:2606718/1‑68 (MQ=255)
                                   aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:2631960/1‑68 (MQ=255)
                                   aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:949379/1‑68 (MQ=255)
                                         gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGt     <  1:2720363/69‑1 (MQ=255)
                                         gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGt     <  1:894715/69‑1 (MQ=255)
                                         gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGt     <  1:367495/69‑1 (MQ=255)
                                         gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGt     <  1:433660/69‑1 (MQ=255)
                                         gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGt     <  1:2137572/69‑1 (MQ=255)
                                         gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:1130791/69‑1 (MQ=255)
                                         gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:1662781/69‑1 (MQ=255)
                                         gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:43773/69‑1 (MQ=255)
                                         gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:94334/69‑1 (MQ=255)
                                         gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:61664/69‑1 (MQ=255)
                                         gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:2920777/69‑1 (MQ=255)
                                         gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:220760/69‑1 (MQ=255)
                                         gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:2408064/69‑1 (MQ=255)
                                         gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:2760041/69‑1 (MQ=255)
                                         gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:2754100/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:496739/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:1115361/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:771355/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:483842/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:836672/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    >  1:858129/1‑69 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:2232916/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:884585/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:1415039/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:2293220/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    >  1:2948861/1‑69 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    >  1:1708104/1‑69 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:189847/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:2911532/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:1470764/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:2737515/69‑1 (MQ=255)
                                          aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:2213932/69‑1 (MQ=255)
                                          aaCGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:3316344/68‑1 (MQ=255)
                                          aaCGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:806696/68‑1 (MQ=255)
                                            cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAc  <  1:1882526/69‑1 (MQ=255)
                                            cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAc  <  1:1624051/69‑1 (MQ=255)
                                            cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAc  <  1:903075/69‑1 (MQ=255)
                                            cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAc  <  1:1735277/69‑1 (MQ=255)
                                            cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAc  <  1:2998196/69‑1 (MQ=255)
                                            cgcgCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGTgcc  <  1:2212479/69‑3 (MQ=255)
                                            cgcgCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGTgcc  <  1:2801069/69‑3 (MQ=255)
                                            cgcgCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGTgcc  <  1:1361838/69‑3 (MQ=255)
                                                   |                                                             
CAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAC  >  minE/303007‑303119

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