Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,602,887 A→G 45.2% T239T (ACT→ACC yfeX ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,602,8870AG45.2% 7.0 / 50.0 62T239T (ACT→ACCyfeXconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/6);  new base G (18/10);  total (46/16)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCG  >  minE/1602823‑1602952
                                                                |                                                                 
gctgctCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGc                                                               >  1:3173321/1‑69 (MQ=255)
      cAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTAc                                                          <  1:2729112/68‑1 (MQ=255)
      cAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTAc                                                          >  1:1394353/1‑68 (MQ=255)
      cAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTAc                                                          <  1:1517069/68‑1 (MQ=255)
             aTGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCGGGct                                                  >  1:362646/1‑69 (MQ=255)
             aTGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGct                                                  >  1:1155532/1‑69 (MQ=255)
             aTGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGct                                                  >  1:1033607/1‑69 (MQ=255)
             aTGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGct                                                  >  1:744452/1‑69 (MQ=255)
             aTGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGct                                                  >  1:627063/1‑69 (MQ=255)
             aTGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGct                                                  >  1:197181/1‑69 (MQ=255)
                    cgcgcgCAGTAGGCGCAGAAGTACAGTCCGTGAGTGCCACTGGCa                                                                   <  1:509966/45‑1 (MQ=255)
                    cgcgcgCAGTAGGCGCAGAAGTACAGACCGTGATTGCCACTGGCa                                                                   <  1:2672782/45‑1 (MQ=255)
                    cgcgcgCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCa                                                                   <  1:1954916/45‑1 (MQ=255)
                    cgcgcgCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCa                                                                   <  1:2813501/45‑1 (MQ=255)
                    cgcgcgCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCa                                                                   <  1:3315556/45‑1 (MQ=255)
                    cgcgcgCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCa                                                                   <  1:1148765/45‑1 (MQ=255)
                    cgcgcgCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCa                                                                   <  1:1536687/45‑1 (MQ=255)
                    cgcgcgCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCa                                                                   <  1:1482213/45‑1 (MQ=255)
                    cgcgcgCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCg                                            >  1:604395/1‑68 (MQ=255)
                    cgcgcgCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCg                                            >  1:608602/1‑68 (MQ=255)
                        cgcAGTAGGCGCAGAAGTACAGCCCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACaa                                       >  1:2259127/1‑69 (MQ=255)
                        cgcAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACaa                                       >  1:1051991/1‑69 (MQ=255)
                        cgcAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACaa                                       >  1:2528255/1‑69 (MQ=255)
                        cgcAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACaa                                       >  1:2524372/1‑69 (MQ=255)
                        cgcAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACaa                                       >  1:2263172/1‑69 (MQ=255)
                         gcAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAAt                                      <  1:3186270/69‑1 (MQ=255)
                         gcAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAAt                                      <  1:2961535/69‑1 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGCCGAACAATCTTCAGccc                             >  1:569035/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:920619/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:3244002/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:3314284/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:962921/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:354087/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:362729/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:970749/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:2142427/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:1908525/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:1552398/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:2044281/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:2079806/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:213127/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:1415997/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:1215879/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:2518006/1‑69 (MQ=255)
                                  gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                             >  1:1182025/1‑69 (MQ=255)
                                         tACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACaa                                       >  1:240294/1‑52 (MQ=255)
                                         tACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCa                      >  1:673127/1‑69 (MQ=255)
                                         tACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCa                      >  1:643064/1‑69 (MQ=255)
                                         tACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCa                      >  1:935524/1‑69 (MQ=255)
                                         tACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCa                      >  1:2744181/1‑69 (MQ=255)
                                           cAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCa                      >  1:1441776/1‑67 (MQ=255)
                                           cAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCAtc                    >  1:1966076/1‑69 (MQ=255)
                                           cAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCAt                     >  1:714254/1‑68 (MQ=255)
                                           cAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCAt                     >  1:443523/1‑68 (MQ=255)
                                           cAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCAt                     >  1:2052418/1‑68 (MQ=255)
                                           cAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCAt                     >  1:2201295/1‑68 (MQ=255)
                                           cAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCAt                     >  1:3120329/1‑68 (MQ=255)
                                           cAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCAt                     >  1:2863794/1‑68 (MQ=255)
                                                       gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc        <  1:1734659/69‑1 (MQ=255)
                                                       gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc        <  1:1633577/69‑1 (MQ=255)
                                                       gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc        <  1:609600/69‑1 (MQ=255)
                                                       gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc        <  1:2065976/69‑1 (MQ=255)
                                                       gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc        <  1:3291963/69‑1 (MQ=255)
                                                       gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc        <  1:713357/69‑1 (MQ=255)
                                                       gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc        <  1:855566/69‑1 (MQ=255)
                                                       gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc        <  1:2758014/69‑1 (MQ=255)
                                                       gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc        <  1:119743/69‑1 (MQ=255)
                                                       gCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc        <  1:2696782/69‑1 (MQ=255)
                                                        ccACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc        <  1:760614/68‑1 (MQ=255)
                                                             ggCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAAcgcg  <  1:2577051/69‑1 (MQ=255)
                                                                |                                                                 
GCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCG  >  minE/1602823‑1602952

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: