Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,159,315 C→T 49.8% T171T (ACC→ACT ispB → octaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,159,3150CT49.8% ‑5.2 / 79.7 87T171T (ACC→ACTispBoctaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (22/21);  new base T (23/20);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

AACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACGCCGGAGGAG  >  minE/2159262‑2159378
                                                     |                                                               
aaCGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGt                                                       >  1:1028377/1‑64 (MQ=255)
aaCGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAg                                                  >  1:2725862/1‑69 (MQ=255)
aaCGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAg                                                  >  1:2773112/1‑69 (MQ=255)
aaCGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAg                                                  >  1:1423846/1‑69 (MQ=255)
aaCGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAg                                                  >  1:1002856/1‑69 (MQ=255)
     tCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTg                                              >  1:364071/1‑68 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATTCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:3010681/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGcc                                            >  1:912414/1‑69 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGcc                                            >  1:53109/1‑69 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:3311624/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:259015/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:2534863/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:3265777/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:2970614/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:2146381/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:2128005/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:311310/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:1833970/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:1803859/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:335575/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:731390/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:1401893/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:3198495/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:1780156/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:3264618/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:3410466/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                                      <  1:723548/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGAGCGTCTGtt                                                      <  1:1396449/59‑1 (MQ=255)
      ccGGACATAACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAAcc                                                                 >  1:321682/1‑48 (MQ=255)
            aTCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAg                                      <  1:3055218/69‑1 (MQ=255)
            aTCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAg                                      <  1:3089133/69‑1 (MQ=255)
            aTCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAg                                      <  1:1590495/69‑1 (MQ=255)
              cACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTg                                                        <  1:1833737/49‑1 (MQ=255)
              cACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTg                                                        <  1:3429543/49‑1 (MQ=255)
               aCTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAg                                      >  1:622369/1‑66 (MQ=255)
                  gaagaaAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTcc                                >  1:2588689/1‑69 (MQ=255)
                  gaagaaAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTcc                                >  1:1775251/1‑69 (MQ=255)
                     gaaAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCggg                             >  1:2197288/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:3413234/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:2937490/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:3136097/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:3285346/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:402494/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:1124293/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:1151930/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:1232388/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:1308805/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:1333437/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:1422769/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:1433622/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:154484/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:710042/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:698735/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:1965989/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:2112398/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:2309855/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:2353837/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:469184/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:2675892/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTAAAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:613303/1‑69 (MQ=255)
                          cTACATACGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct                        >  1:441200/1‑69 (MQ=255)
                            aCATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGc                                             <  1:938797/46‑1 (MQ=255)
                               tGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTctggct                    <  1:3121982/68‑1 (MQ=255)
                               tGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTctggct                    <  1:3427795/68‑1 (MQ=255)
                               tGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTctggct                    <  1:681942/68‑1 (MQ=255)
                               tGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTctggct                    <  1:2389233/68‑1 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCggg                             >  1:2743183/1‑51 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGAt                           >  1:1226861/1‑53 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGAt                           >  1:1354442/1‑53 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGAt                           >  1:1719751/1‑53 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGAt                           >  1:2497647/1‑53 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGAt                           >  1:1926673/1‑53 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGgctggctg               <  1:2258562/65‑1 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGgctggctg               <  1:616624/65‑1 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGgctggctg               <  1:241508/65‑1 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACg           >  1:3252233/1‑69 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACg           >  1:461385/1‑69 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACg           >  1:343077/1‑69 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACg           >  1:3124867/1‑69 (MQ=255)
                                       aTCTATAGCAAAACAGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGAt                           >  1:2919305/1‑53 (MQ=255)
                                                aaaaCCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACGCCggaggag  <  1:2422050/69‑1 (MQ=255)
                                                aaaaCCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACGCCggaggag  <  1:692931/69‑1 (MQ=255)
                                                aaaaCCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACGCCggaggag  <  1:1915757/69‑1 (MQ=255)
                                                aaaaCCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACGCCggaggag  <  1:700840/69‑1 (MQ=255)
                                                aaaaCCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACGCCggaggag  <  1:1898724/69‑1 (MQ=255)
                                                aaaaCCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACGCCggaggag  <  1:2771319/69‑1 (MQ=255)
                                                aaaaCCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACGCCggaggag  <  1:73951/69‑1 (MQ=255)
                                                aaaaCCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACGCCggaggag  <  1:1284277/69‑1 (MQ=255)
                                                     |                                                               
AACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACGCCGGAGGAG  >  minE/2159262‑2159378

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