Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,485,005 C→T 47.9% D279N (GAT→AAT)  menF ← isochorismate synthase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,485,0050CT47.9% 10.0 / 103.7 73D279N (GAT→AAT) menFisochorismate synthase 2
Reads supporting (aligned to +/- strand):  ref base C (21/17);  new base T (12/23);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGC  >  minE/1484941‑1485074
                                                                |                                                                     
gTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGc                                                                 <  1:3381455/71‑1 (MQ=255)
gTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGc                                                                 <  1:3794435/71‑1 (MQ=255)
 tATCGGCCTGTAATCGTTTACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGc                                                                 <  1:1369914/70‑1 (MQ=255)
 tATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGcc                                                                <  1:1936043/71‑1 (MQ=255)
 tATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGc                                                                 <  1:1517077/70‑1 (MQ=255)
 tATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGc                                                                 <  1:1401997/70‑1 (MQ=255)
 tATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGc                                                                 <  1:2983637/70‑1 (MQ=255)
     ggCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTatca                                                                      >  1:3053872/1‑61 (MQ=255)
     ggCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTatca                                                                      >  1:3629177/1‑61 (MQ=255)
     ggCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTatca                                                                      >  1:1004736/1‑61 (MQ=255)
     ggCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTatca                                                                      >  1:1069280/1‑61 (MQ=255)
     ggCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTatca                                                                      >  1:2510889/1‑61 (MQ=255)
       ccTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGc                                                          <  1:2933476/71‑1 (MQ=255)
        cTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGc                                                          <  1:128873/70‑1 (MQ=255)
        cTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGc                                                          <  1:812983/70‑1 (MQ=255)
           tAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCa                                                        <  1:2827448/69‑1 (MQ=255)
             aTCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCg                                                                  >  1:1561803/1‑57 (MQ=255)
             aTCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCActct                                                    >  1:2463468/1‑71 (MQ=255)
             aTCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCActct                                                    >  1:3374202/1‑71 (MQ=255)
             aTCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCActct                                                    >  1:117171/1‑71 (MQ=255)
             aTCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCActct                                                    >  1:3702273/1‑71 (MQ=255)
             aTCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCActc                                                     >  1:1925340/1‑70 (MQ=255)
                gTTGACATATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACTCTCCt                                                 <  1:765108/71‑1 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTatca                                                                      <  1:1868673/50‑1 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTatca                                                                      <  1:1045095/50‑1 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACTCTcc                                                  >  1:1887261/1‑70 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACTCTCCt                                                 <  1:3752482/71‑1 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACTCTCCt                                                 <  1:598735/71‑1 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACTCTCCt                                                 <  1:2621159/71‑1 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTcc                                                  >  1:28097/1‑70 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTcc                                                  >  1:2738685/1‑70 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTcc                                                  >  1:2352628/1‑70 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTcc                                                  >  1:1951373/1‑70 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTcc                                                  >  1:2137974/1‑70 (MQ=255)
                gTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCt                                                 <  1:1221603/71‑1 (MQ=255)
                   gACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAAc                                              >  1:3829121/1‑71 (MQ=255)
                   gACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAAc                                              >  1:2755584/1‑71 (MQ=255)
                   gACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAAc                                              >  1:2973216/1‑71 (MQ=255)
                    aCAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGc                                                                 <  1:2241802/51‑1 (MQ=255)
                        atatCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGcc                                                         <  1:3681918/55‑1 (MQ=255)
                        atatCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGcc                                                         <  1:3938574/55‑1 (MQ=255)
                        atatCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCt                                                 <  1:1508187/63‑1 (MQ=255)
                        atatCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCt                                                 <  1:3804101/63‑1 (MQ=255)
                        atatCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCt                                                 <  1:2608086/63‑1 (MQ=255)
                           tCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTatca                                                                      <  1:3689111/39‑1 (MQ=255)
                           tCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTgcg                                       <  1:1389395/70‑1 (MQ=255)
                           tCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTgcg                                       <  1:911966/70‑1 (MQ=255)
                             ttCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCa                                                            >  1:3074828/1‑47 (MQ=255)
                             ttCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCa                                                            >  1:3113838/1‑47 (MQ=255)
                               ccaccaccaGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCATCCa                                                        <  1:2528388/49‑1 (MQ=39)
                               ccaccaccaGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCa                                                        <  1:2149468/49‑1 (MQ=255)
                               ccaccaccaGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCa                                                        <  1:2186634/49‑1 (MQ=255)
                               ccaccaccaGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCa                                                        <  1:3288475/49‑1 (MQ=255)
                                 accaccaGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCAct                                                      >  1:1970820/1‑49 (MQ=255)
                                 accaccaGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCAct                                                      >  1:4000462/1‑49 (MQ=255)
                                 accaccaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCt                                 >  1:2255675/1‑70 (MQ=255)
                                   caccaGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTat                              <  1:3186377/71‑1 (MQ=255)
                                   caccaGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTat                              <  1:851227/71‑1 (MQ=255)
                                   caccaGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTat                              <  1:438398/71‑1 (MQ=255)
                                   caccaGCATGTTCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTat                              <  1:682394/71‑1 (MQ=255)
                                    accaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTat                              >  1:120158/1‑70 (MQ=255)
                                    accaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTat                              >  1:918256/1‑70 (MQ=255)
                                    accaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTat                              >  1:99914/1‑70 (MQ=255)
                                    accaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTat                              >  1:2718338/1‑70 (MQ=255)
                                    accaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTa                               >  1:1799227/1‑69 (MQ=255)
                                             ttCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACtctc                                                   >  1:377923/1‑40 (MQ=255)
                                             ttCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACtctc                                                   >  1:259173/1‑40 (MQ=255)
                                             ttCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACtctc                                                   >  1:3477983/1‑40 (MQ=255)
                                             ttCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACtctc                                                   >  1:1706102/1‑40 (MQ=255)
                                             ttCTCGCGCTGGTTTTTATTATCCGCCATCAGCCACtctc                                                   >  1:498489/1‑40 (MQ=255)
                                             ttCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGat                     <  1:3747726/70‑1 (MQ=255)
                                             ttCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGat                     >  1:2158148/1‑70 (MQ=255)
                                              tctcGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCt                                    >  1:807645/1‑54 (MQ=255)
                                               ctcGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGatta                   <  1:3144210/70‑1 (MQ=255)
                                                     cTGGTTTTTATTATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGatt                    <  1:1393439/63‑1 (MQ=255)
                                                     cTGGTTTTTATTATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGatt                    <  1:1219255/63‑1 (MQ=255)
                                                     cTGGTTTTTATTATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGatt                    <  1:2972961/63‑1 (MQ=255)
                                                     cTGGTTTTTATTATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGatt                    <  1:675046/63‑1 (MQ=255)
                                                     cTGGTTTTTATTATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGatt                    <  1:1918706/63‑1 (MQ=255)
                                                     cTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTa            >  1:1254456/1‑71 (MQ=255)
                                                        gTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTAc           <  1:248967/69‑1 (MQ=255)
                                                               tcatcCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTa            <  1:3618953/61‑1 (MQ=255)
                                                               tcatcCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGc  <  1:2565754/71‑1 (MQ=255)
                                                                |                                                                     
GTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGC  >  minE/1484941‑1485074

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: ‑