Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,602,908 G→A 39.7% R232R (CGC→CGT yfeX ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,602,9080GA39.7% 32.0 / 44.3 69R232R (CGC→CGTyfeXconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/19);  new base A (17/10);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACG  >  minE/1602843‑1602976
                                                                 |                                                                    
cgcgcgCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCg                                                                    >  1:608602/1‑68 (MQ=255)
cgcgcgCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCg                                                                    >  1:604395/1‑68 (MQ=255)
    cgcAGTAGGCGCAGAAGTACAGCCCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACaa                                                               >  1:2259127/1‑69 (MQ=255)
    cgcAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACaa                                                               >  1:1051991/1‑69 (MQ=255)
    cgcAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACaa                                                               >  1:2263172/1‑69 (MQ=255)
    cgcAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACaa                                                               >  1:2528255/1‑69 (MQ=255)
    cgcAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACaa                                                               >  1:2524372/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:2079806/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:970749/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:213127/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:2142427/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:2518006/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:3244002/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:3314284/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:354087/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:362729/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:920619/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:962921/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:1182025/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:1552398/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:1908525/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:1415997/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:1215879/1‑69 (MQ=255)
              gcAGAAGTACAGACCGTGAGTGCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGccc                                                     >  1:2044281/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:3120329/1‑68 (MQ=255)
                       cAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCAt                                             >  1:2863794/1‑68 (MQ=255)
                       cAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCAt                                             >  1:2201295/1‑68 (MQ=255)
                                   gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc                                <  1:119743/69‑1 (MQ=255)
                                   gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc                                <  1:855566/69‑1 (MQ=255)
                                   gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc                                <  1:713357/69‑1 (MQ=255)
                                   gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc                                <  1:609600/69‑1 (MQ=255)
                                   gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc                                <  1:1633577/69‑1 (MQ=255)
                                   gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc                                <  1:3291963/69‑1 (MQ=255)
                                   gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc                                <  1:1734659/69‑1 (MQ=255)
                                   gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc                                <  1:2065976/69‑1 (MQ=255)
                                   gCCACTGGCGGTGCCGTACGGCAGGCTCTGACGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATc                                <  1:2758014/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)
                                             gTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGggtg                      <  1:1769768/69‑1 (MQ=255)
                                             gTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGggtg                      <  1:735695/69‑1 (MQ=255)
                                             gTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGggtg                      <  1:538159/69‑1 (MQ=255)
                                             gTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGggtg                      <  1:2691570/69‑1 (MQ=255)
                                             gTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGggtg                      <  1:2033776/69‑1 (MQ=255)
                                             gTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGggtg                      <  1:1808056/69‑1 (MQ=255)
                                             gTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGggtg                      <  1:3286268/69‑1 (MQ=255)
                                                       gCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGagag             <  1:440030/68‑1 (MQ=255)
                                                       gCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGagag             <  1:2477314/68‑1 (MQ=255)
                                                       gCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGagag             <  1:2401491/68‑1 (MQ=255)
                                                       gCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGagag             <  1:2579221/68‑1 (MQ=255)
                                                       gCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGagag             <  1:3102995/68‑1 (MQ=255)
                                                       gCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGagag             <  1:2534282/68‑1 (MQ=255)
                                                            ctctGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTc        >  1:2247804/1‑68 (MQ=255)
                                                            ctctGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTc        >  1:219563/1‑68 (MQ=255)
                                                               tGGCGAACAATCTTCATCCCTTTGCCATCTTCTTTCAGATCAACGCGggtg                      <  1:1947815/51‑1 (MQ=255)
                                                               tGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGggtg                      <  1:1943768/51‑1 (MQ=255)
                                                                 gCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACg  >  1:1478983/1‑69 (MQ=255)
                                                                 gCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACg  >  1:1827220/1‑69 (MQ=255)
                                                                 gCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACg  >  1:29954/1‑69 (MQ=255)
                                                                 gCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACg  >  1:1008209/1‑69 (MQ=255)
                                                                 |                                                                    
CGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACG  >  minE/1602843‑1602976

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: