Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,097,286 T→C 39.5% S174S (AGT→AGC ydiK → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,097,2860TC39.5% 14.8 / 27.3 43S174S (AGT→AGCydiKpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (18/8);  new base C (9/8);  total (27/16)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

CGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCT  >  minE/1097223‑1097348
                                                               |                                                              
cGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGc                                                              <  1:563959/66‑1 (MQ=255)
                    tCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAAc                                       >  1:3326169/1‑69 (MQ=255)
                    tCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAAc                                       >  1:2174991/1‑69 (MQ=255)
                    tCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAAc                                       >  1:925934/1‑69 (MQ=255)
                                  gTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAgg                          <  1:2945745/68‑1 (MQ=255)
                                  gTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAgg                          <  1:649473/68‑1 (MQ=255)
                                    gCATTGTGCGCTGATGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAt                       <  1:1324718/69‑1 (MQ=255)
                                    gCATTGTGCGCTGATGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAt                       <  1:1303878/69‑1 (MQ=255)
                                    gCATTGTGCGCTGATGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAt                       <  1:1592297/69‑1 (MQ=255)
                                    gCATTGTGCGCTGATGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAt                       <  1:1713416/69‑1 (MQ=255)
                                    gCATTGTGCGCTGATGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAt                       <  1:2692114/69‑1 (MQ=255)
                                    gCATTGTGCGCTGATGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAt                       <  1:1925863/69‑1 (MQ=255)
                                    gCATTGTGCGCTGATGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAt                       <  1:2616449/69‑1 (MQ=255)
                                    gCATTGTGCGCTGATGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAt                       <  1:1177205/69‑1 (MQ=255)
                                     cATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAtt                      >  1:1328571/1‑69 (MQ=255)
                                     cATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAtt                      >  1:3143593/1‑69 (MQ=255)
                                     cATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAtt                      >  1:1228944/1‑69 (MQ=255)
                                     cATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAtt                      >  1:3403025/1‑69 (MQ=255)
                                     cATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAt                       >  1:17685/1‑68 (MQ=255)
                                     cATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTG‑‑GCGGTGAACAGGTGGCACAAGGCATTc                     >  1:408406/1‑68 (MQ=255)
                                        tgtgCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGc                   <  1:3191867/69‑1 (MQ=255)
                                        tgtgCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGc                   <  1:1037979/69‑1 (MQ=255)
                                        tgtgCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGc                   <  1:2038285/69‑1 (MQ=255)
                                        tgtgCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGc                   <  1:63620/69‑1 (MQ=255)
                                        tgtgCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGc                   <  1:1751973/69‑1 (MQ=255)
                                         gtgCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGcc                  >  1:3030050/1‑69 (MQ=255)
                                                 aTGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAt                       >  1:1657175/1‑56 (MQ=255)
                                                 aTGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCa          >  1:319347/1‑69 (MQ=255)
                                                 aTGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCa          >  1:1040135/1‑69 (MQ=255)
                                                 aTGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCa          >  1:3082336/1‑69 (MQ=255)
                                                 aTGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCa          >  1:287161/1‑69 (MQ=255)
                                                 aTGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCa          >  1:2180242/1‑69 (MQ=255)
                                                 aTGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCa          >  1:1962839/1‑69 (MQ=255)
                                                 aTGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCa          >  1:1576368/1‑69 (MQ=255)
                                                 aTGCTTCTCTTCAGCGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCa          >  1:1524719/1‑69 (MQ=255)
                                                          ttCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCt  >  1:3014954/1‑68 (MQ=255)
                                                          ttCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCt  >  1:2895912/1‑68 (MQ=255)
                                                          ttCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCt  >  1:2805775/1‑68 (MQ=255)
                                                          ttCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCt  >  1:3421623/1‑68 (MQ=255)
                                                          ttCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCt  >  1:1963470/1‑68 (MQ=255)
                                                          ttCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCt  >  1:607977/1‑68 (MQ=255)
                                                          ttCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCt  >  1:1914891/1‑68 (MQ=255)
                                                          ttCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGAAACCCGTCt  >  1:2232796/1‑68 (MQ=255)
                                                               |                                                              
CGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCT  >  minE/1097223‑1097348

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: