Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,112,851 T→A 68.3% V651V (GTA→GTT pps ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,112,8510TA68.3% 92.7 / 59.5 85V651V (GTA→GTTppsphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base T (19/8);  new base A (16/42);  total (35/50)
Fisher's exact test for biased strand distribution p-value = 3.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CCCGTTCTCGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTT  >  minE/1112788‑1112916
                                                               |                                                                 
cccGTTCTCGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACgg                                                               <  1:3418793/68‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:317995/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:1925809/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:3328127/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:2288944/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:544195/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:2628138/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:1034492/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:1570227/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:3227292/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:2737193/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:1291284/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:2817975/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:2820836/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:2916267/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:3108989/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:929303/69‑1 (MQ=255)
cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACCGt                                                              <  1:3163005/69‑1 (MQ=255)
 ccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:493064/68‑1 (MQ=255)
 ccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt                                                              <  1:188708/68‑1 (MQ=255)
      ctcGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCa                                                         >  1:2785149/1‑68 (MQ=255)
      ctcGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCa                                                         >  1:3197036/1‑68 (MQ=255)
      ctcGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCa                                                         >  1:3227671/1‑68 (MQ=255)
      ctcGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCa                                                         >  1:2248657/1‑68 (MQ=255)
      ctcGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCa                                                         >  1:731676/1‑68 (MQ=255)
      ctcGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCa                                                         >  1:1453614/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:774937/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:340341/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:1543061/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:272334/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:2728363/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:2742419/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:824543/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:1233903/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:179038/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:3013331/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:3049002/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:1147270/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:1874741/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:340276/1‑68 (MQ=255)
              gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:1879779/1‑68 (MQ=255)
              gTTTCAGCCCCTGATACGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg                                                 >  1:1614958/1‑68 (MQ=255)
                                     tcttcAACCACCGCTTTCGCCTGATCTACGGTACGCAc                                                        >  1:2634094/1‑38 (MQ=255)
                                     tcttcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGt                         >  1:2679807/1‑69 (MQ=255)
                                     tcttcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGt                         >  1:1448764/1‑69 (MQ=255)
                                      cttcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATGTCAACGTTGGTCAGt                         >  1:2844170/1‑68 (MQ=255)
                                      cttcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAg                          >  1:954172/1‑67 (MQ=255)
                                      cttcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGt                         >  1:2565099/1‑68 (MQ=255)
                                      cttcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGt                         >  1:3243025/1‑68 (MQ=255)
                                      cttcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGt                         >  1:3297703/1‑68 (MQ=255)
                                      cttcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGt                         >  1:1740035/1‑68 (MQ=255)
                                      cttcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGt                         >  1:2267364/1‑68 (MQ=255)
                                      cttcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGt                         >  1:2147378/1‑68 (MQ=255)
                                      cttcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGt                         >  1:3381236/1‑68 (MQ=255)
                                        tcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTcc                       <  1:2502539/68‑1 (MQ=255)
                                        tcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTcc                       <  1:2768308/68‑1 (MQ=255)
                                        tcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTcc                       <  1:321662/68‑1 (MQ=255)
                                        tcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTcc                       <  1:2838268/68‑1 (MQ=255)
                                        tcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTcc                       <  1:123030/68‑1 (MQ=255)
                                        tcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTcc                       <  1:197099/68‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:2647247/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:850119/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:3180785/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:52477/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:2814617/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:1155375/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:1206488/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:1252576/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:1665811/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:1736428/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:1839312/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:188405/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:1910734/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:2183962/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:2344671/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:2515422/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:2623971/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:2664366/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:2720124/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:3077796/63‑1 (MQ=255)
                                         cAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTAAACGTTAGTCa                           <  1:1369199/63‑1 (MQ=255)
                                               ccGCTTTCGCCTGATCAACGGTACGCACGAACGGGATCATGATCTCAACGTTAGTCa                           <  1:1661876/57‑1 (MQ=255)
                                                  cTTTCGCCTGATCAACGGTACGCACGAACGGTATCATGATCTCAACGTTAGTCa                           <  1:1628676/54‑1 (MQ=255)
                                                              cTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGttt  <  1:2140540/67‑1 (MQ=255)
                                                              cTACGGTACGCACGAACGGGATCATGATCTAACGTTGGTCAGTCCCATGTCGtt              >  1:113492/1‑54 (MQ=255)
                                                               |                                                                 
CCCGTTCTCGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTT  >  minE/1112788‑1112916

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: