Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,305,959 A→G 6.7% Y34H (TAT→CAT)  yeeE ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,305,9590AG 6.7% 222.1 / 5.3 74Y34H (TAT→CAT) yeeEpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (26/43);  new base G (2/3);  total (28/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

AAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACA  >  minE/1305889‑1306028
                                                                      |                                                                     
aaaaCCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGAtata                                                                       <  1:3291725/71‑1 (MQ=255)
 aaaCCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATAc                                                                      <  1:2086725/71‑1 (MQ=255)
   aCCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACat                                                                    <  1:4810692/71‑1 (MQ=255)
    cccccACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACata                                                                   <  1:4326728/71‑1 (MQ=255)
    cccccACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACata                                                                   <  1:4951621/71‑1 (MQ=255)
    cccccACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACata                                                                   <  1:3931386/71‑1 (MQ=255)
     ccccACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACata                                                                   <  1:1246210/70‑1 (MQ=255)
           gCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:3746792/71‑1 (MQ=255)
           gCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:3586935/71‑1 (MQ=255)
           gCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:2814528/71‑1 (MQ=255)
           gCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:2760997/71‑1 (MQ=255)
           gCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:5160418/71‑1 (MQ=255)
            cTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            >  1:439366/1‑70 (MQ=255)
            cTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            >  1:4457464/1‑70 (MQ=255)
            cTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            >  1:4604796/1‑70 (MQ=255)
            cTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            >  1:3187694/1‑70 (MQ=255)
            cTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            >  1:3186079/1‑70 (MQ=255)
               tGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAac                                                         >  1:572557/1‑70 (MQ=255)
                gTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAacca                                                       >  1:2479079/1‑71 (MQ=255)
                           gCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGc                                             <  1:902378/70‑1 (MQ=255)
                           gCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGc                                             <  1:4621231/70‑1 (MQ=255)
                               tCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCaaaaa                                        >  1:3420142/1‑71 (MQ=255)
                               tCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCaaaaa                                        >  1:2040608/1‑71 (MQ=255)
                                 aGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCa                                                >  1:1101156/1‑61 (MQ=255)
                                 aGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCa                                                >  1:1382139/1‑61 (MQ=255)
                                 aGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCa                                                >  1:5169781/1‑61 (MQ=255)
                                 aGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCa                                                >  1:2644936/1‑61 (MQ=255)
                                 aGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAAc                                       <  1:838385/70‑1 (MQ=255)
                                 aGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAAc                                       <  1:2771751/70‑1 (MQ=255)
                                 aGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACg                                      >  1:80739/1‑71 (MQ=255)
                                 aGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACg                                      >  1:1588494/1‑71 (MQ=255)
                                 aGCAAGGCGTAAAACATGAGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCa                                                >  1:1265188/1‑61 (MQ=255)
                                       gCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACg                                >  1:1353956/1‑71 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:3347921/40‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:4246192/40‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:692412/40‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:177111/40‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:1390294/40‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:860585/40‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:2393029/40‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:4872831/40‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                 <  1:3235217/51‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                 <  1:4630912/51‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                 <  1:1009198/51‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                 <  1:2061096/51‑1 (MQ=255)
                                          tAAAACATGCGATTATTTTTCACGATATAAATATCGCGaa                                                            <  1:2262307/40‑1 (MQ=255)
                                           aaaaCATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:827242/39‑1 (MQ=255)
                                           aaaaCATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:4799582/39‑1 (MQ=255)
                                           aaaaCATGCGATTATTTTTCACGATATACATATCGCGaa                                                            <  1:3041435/39‑1 (MQ=255)
                                            aaaCATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                 >  1:3214995/1‑49 (MQ=255)
                                            aaaCATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                 >  1:1871662/1‑49 (MQ=255)
                                            aaaCATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                 >  1:1232850/1‑49 (MQ=255)
                                            aaaCATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                 >  1:618475/1‑49 (MQ=255)
                                                      ttattTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCa                                            <  1:918220/44‑1 (MQ=255)
                                                      ttattTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCa                                            <  1:3919602/44‑1 (MQ=255)
                                                      ttattTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCa                                            <  1:1823206/44‑1 (MQ=255)
                                                      ttattTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATcc                 <  1:3587608/71‑1 (MQ=255)
                                                       tattTTTCACGATATACATATCGCGAAAACCACCTGTCAg                                               >  1:3889553/1‑40 (MQ=255)
                                                       tattTTTCACGATATACATATCGCGAAAACCACCTGTCAg                                               >  1:1527246/1‑40 (MQ=255)
                                                       tattTTTCACGATATACATATCGCGAAAACCACCTGTCAg                                               >  1:2346320/1‑40 (MQ=255)
                                                       tattTTTCACGATATACATATCGCGAAAACCACCTGTCAg                                               >  1:404673/1‑40 (MQ=255)
                                                       tattTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATcc                 >  1:1476000/1‑70 (MQ=255)
                                                         tttttCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcag              <  1:2490811/71‑1 (MQ=255)
                                                         tttttCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcag              <  1:1906756/71‑1 (MQ=255)
                                                         tttttCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcag              <  1:2357071/71‑1 (MQ=255)
                                                         tttttCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcag              <  1:2380298/71‑1 (MQ=255)
                                                         tttttCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcag              <  1:4823649/71‑1 (MQ=255)
                                                         tttttCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcag              <  1:3518613/71‑1 (MQ=255)
                                                          ttttCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcag              <  1:1140211/70‑1 (MQ=255)
                                                           tttCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcag              <  1:3785861/69‑1 (MQ=255)
                                                              cgcGATATGTATATCGTGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCagag         >  1:4712559/3‑71 (MQ=38)
                                                              cgcGATATGTATATCGTGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCagag         >  1:1532021/3‑71 (MQ=38)
                                                               gcGATATGTATATCGTGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCagag         <  1:2801422/69‑1 (MQ=38)
                                                               gcGATATGTATATCGTGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCagag         <  1:3318746/69‑1 (MQ=38)
                                                               gcGATATGTATATCGTGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCagag         <  1:1183334/69‑1 (MQ=38)
                                                                     taCATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGcaccac   >  1:3588689/1‑70 (MQ=255)
                                                                     taCATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGcaccac   >  1:1704623/1‑70 (MQ=255)
                                                                     taCATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGcaccac   >  1:3839646/1‑70 (MQ=255)
                                                                     taCATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACcaca  <  1:859913/71‑1 (MQ=255)
                                                                     taCATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACcaca  <  1:3850982/71‑1 (MQ=255)
                                                                      |                                                                     
AAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACA  >  minE/1305889‑1306028

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: