Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,159,306 T→C 53.4% Y168Y (TAT→TAC ispB → octaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,159,3060TC53.4% 10.1 / 65.0 80Y168Y (TAT→TACispBoctaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (24/13);  new base C (22/20);  total (47/33)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACG  >  minE/2159262‑2159369
                                            |                                                               
aaCGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGt                                              >  1:1028377/1‑64 (MQ=255)
aaCGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAg                                         >  1:2773112/1‑69 (MQ=255)
aaCGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAg                                         >  1:1002856/1‑69 (MQ=255)
aaCGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAg                                         >  1:2725862/1‑69 (MQ=255)
aaCGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAg                                         >  1:1423846/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:2128005/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:2970614/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:311310/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:3198495/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:3264618/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:3311624/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:335575/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:3410466/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:2146381/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:731390/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:723548/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:1803859/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:1401893/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:1780156/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGtt                                             <  1:1833970/59‑1 (MQ=255)
      ccGGACATCACTGAAGAAAACTACATGCGCGTTATCTACAGCAAGACTGAGCGTCTGtt                                             <  1:1396449/59‑1 (MQ=255)
      ccGGACATAACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAAcc                                                        >  1:321682/1‑48 (MQ=255)
            aTCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAg                             <  1:1590495/69‑1 (MQ=255)
            aTCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAg                             <  1:3055218/69‑1 (MQ=255)
            aTCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAg                             <  1:3089133/69‑1 (MQ=255)
              cACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTg                                               <  1:3429543/49‑1 (MQ=255)
              cACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTg                                               <  1:1833737/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:3285346/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:469184/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:698735/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:710042/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:3413234/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:3136097/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:1308805/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:402494/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:1422769/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:1433622/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:1232388/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:1124293/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:1151930/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:1333437/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:2675892/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:154484/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:1965989/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:2353837/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:2309855/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:2937490/1‑69 (MQ=255)
                          cTACATGCGCGTTATCTACAGCAAGACTGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTct               >  1:2112398/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:3427795/68‑1 (MQ=255)
                               tGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTctggct           <  1:2389233/68‑1 (MQ=255)
                               tGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTctggct           <  1:3121982/68‑1 (MQ=255)
                               tGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTctggct           <  1:681942/68‑1 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCggg                    >  1:2743183/1‑51 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGAt                  >  1:1926673/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:1226861/1‑53 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGgctggctg      <  1:2258562/65‑1 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGgctggctg      <  1:241508/65‑1 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGgctggctg      <  1:616624/65‑1 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACg  >  1:343077/1‑69 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACg  >  1:461385/1‑69 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACg  >  1:3252233/1‑69 (MQ=255)
                                       aTCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACg  >  1:3124867/1‑69 (MQ=255)
                                       aTCTATAGCAAAACAGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGAt                  >  1:2919305/1‑53 (MQ=255)
                                            |                                                               
AACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGCTGGCTGTACG  >  minE/2159262‑2159369

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: