Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,112,816 C→T 81.3% R663H (CGT→CAT)  pps ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,112,8160CT81.3% 76.2 / 9.9 53R663H (CGT→CAT) ppsphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base C (6/4);  new base T (17/26);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CACATCATGATGATTTTCAGCCCGTTCTCGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGG  >  minE/1112768‑1112869
                                                |                                                     
cacaTCATGATGATTTTCAGCCCGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGtt                                              <  1:970334/58‑1 (MQ=39)
cacaTCATGATGATTTTCAGCCCGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGtt                                              <  1:3049525/58‑1 (MQ=39)
cacaTCATGATGATTTTCAGCCCGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGtt                                              <  1:2469558/58‑1 (MQ=39)
cacaTCATGATGATTTTCAGCCCGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGtt                                              <  1:2292960/58‑1 (MQ=39)
cacaTCATGATGATTTTCAGCCCGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGtt                                              <  1:1742099/58‑1 (MQ=39)
cacaTCATGATGATTTTCAGCCCGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGtt                                              <  1:1366189/58‑1 (MQ=39)
cacaTCATGATGATTTTCAGCCCGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGtt                                              <  1:1362386/58‑1 (MQ=39)
           gatTTTCAGCCCGTTCTCGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTg                        <  1:1735074/69‑1 (MQ=255)
           gatTTTCAGCCCGTTCTCGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTg                        <  1:3089460/69‑1 (MQ=255)
           gatTTTCAGCCCGTTCTCGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTg                        <  1:2060894/69‑1 (MQ=255)
                    cccGTTCTCGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACgg                <  1:3418793/68‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:2916267/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:2820836/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:2817975/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:3328127/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:1034492/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:2737193/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:3108989/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:2628138/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:929303/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:1291284/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:2288944/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:1925809/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:544195/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:317995/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:1570227/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:3227292/69‑1 (MQ=255)
                    cccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACCGt               <  1:3163005/69‑1 (MQ=255)
                     ccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:188708/68‑1 (MQ=255)
                     ccGTTCTCACCACGTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGt               <  1:493064/68‑1 (MQ=255)
                          ctcGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCa          >  1:3197036/1‑68 (MQ=255)
                          ctcGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCa          >  1:3227671/1‑68 (MQ=255)
                          ctcGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCa          >  1:2785149/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)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCg                            >  1:2743188/1‑42 (MQ=37)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:340276/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:2742419/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:340341/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:774937/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:824543/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:3049002/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:3013331/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:2728363/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:272334/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:1879779/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:1874741/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:179038/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:1543061/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:1233903/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATGCGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:1147270/1‑68 (MQ=255)
                                  gTTTCAGCCCCTGATACGCCAGTTCCTCAACCACCGCTTTCGCCTGATCAACGGTACGCACGAACggg  >  1:1614958/1‑68 (MQ=255)
                                                |                                                     
CACATCATGATGATTTTCAGCCCGTTCTCGCCACGTTTCAGCCCCTGACGCGCCAGTTCTTCAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGG  >  minE/1112768‑1112869

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: