Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,910,405 C→T 6.4% V117M (GTG→ATG)  ppdB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,910,4050CT 6.4% 303.5 / 8.2 94V117M (GTG→ATG) ppdBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (52/35);  new base T (3/3);  total (55/39)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

ATCCGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGT  >  minE/1910337‑1910470
                                                                    |                                                                 
aTCCGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACAt                                                                 >  1:671063/1‑71 (MQ=255)
aTCCGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACAt                                                                 >  1:476226/1‑71 (MQ=255)
aTCCGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACAt                                                                 <  1:324445/71‑1 (MQ=255)
aTCCGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACAt                                                                 >  1:1829563/1‑71 (MQ=255)
   cGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGc                                                               <  1:1063562/70‑1 (MQ=255)
   cGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGc                                                               >  1:3410093/1‑70 (MQ=255)
   cGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCt                                                              >  1:485404/1‑71 (MQ=255)
   cGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCt                                                              >  1:974548/1‑71 (MQ=255)
   cGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCt                                                              >  1:3207845/1‑71 (MQ=255)
   cGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCt                                                              >  1:3178319/1‑71 (MQ=255)
   cGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCt                                                              >  1:3095044/1‑71 (MQ=255)
   cGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCt                                                              >  1:2183761/1‑71 (MQ=255)
    ggATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGc                                                               <  1:2828128/69‑1 (MQ=255)
      aTTAGTGAGTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCt                                                           >  1:3399860/1‑71 (MQ=255)
      aTTAGTGACTTTATCCCAGCCCTTACCTTTACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCt                                                           <  1:1029452/71‑1 (MQ=255)
      aTTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCt                                                           <  1:346332/71‑1 (MQ=255)
      aTTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCt                                                           >  1:3063867/1‑71 (MQ=255)
       ttAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGcaca                                                                  >  1:2445197/1‑63 (MQ=255)
       ttAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCtt                                                          <  1:1717656/71‑1 (MQ=255)
       ttAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCtt                                                          <  1:2877395/71‑1 (MQ=255)
        tAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCtt                                                          >  1:1325013/1‑70 (MQ=255)
        tAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCtt                                                          <  1:3369286/70‑1 (MQ=255)
        tAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCtt                                                          >  1:1107418/1‑70 (MQ=255)
        tAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCtt                                                          >  1:1119564/1‑70 (MQ=255)
          gTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAg                                                       <  1:605827/71‑1 (MQ=255)
           tGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCa                                                        <  1:1041491/69‑1 (MQ=255)
            gACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCa                                                        <  1:3243085/68‑1 (MQ=255)
             aCTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGAc                                                     >  1:3457388/1‑70 (MQ=255)
              cTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGa                                                   <  1:316333/71‑1 (MQ=255)
              cTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGa                                                   <  1:2211531/71‑1 (MQ=255)
              cTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGa                                                   <  1:1120204/71‑1 (MQ=255)
              cTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGa                                                   <  1:1198705/71‑1 (MQ=255)
               tttATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGa                                                   <  1:3107262/70‑1 (MQ=255)
                 tATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGaaa                                                 >  1:694072/1‑70 (MQ=255)
                 tATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGaaa                                                 >  1:110916/1‑70 (MQ=255)
                 tATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAAt                                                >  1:1180294/1‑71 (MQ=255)
                   tCCCAGCCCTTCCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCa                                             <  1:1906998/71‑1 (MQ=255)
                     ccAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCa                                             <  1:2543326/70‑1 (MQ=255)
                        gCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCt                                         >  1:3136228/1‑71 (MQ=255)
                        gCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCt                                         >  1:3285023/1‑71 (MQ=255)
                              aCCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCaa                                            >  1:417861/1‑62 (MQ=255)
                              aCCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCaa                                            >  1:2981047/1‑62 (MQ=255)
                              aCCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCaa                                            >  1:1585026/1‑62 (MQ=255)
                              aCCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCaa                                            >  1:3419276/1‑62 (MQ=255)
                              aCCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTGCAGACGAAATCCaa                                            >  1:316787/1‑62 (MQ=255)
                              aCCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCATTCAGACGAAATCCaa                                            >  1:1738356/1‑62 (MQ=255)
                                cTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGAc                                 <  1:943535/71‑1 (MQ=255)
                                  tCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGActc                               >  1:2816209/1‑71 (MQ=255)
                                    acaGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGActct                              >  1:2420437/1‑70 (MQ=255)
                                    acaGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCtt                             >  1:1145429/1‑71 (MQ=255)
                                      aGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCttt                            >  1:2973923/1‑70 (MQ=255)
                                      aGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCttt                            >  1:1807333/1‑70 (MQ=255)
                                            tCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGt                     <  1:2025017/71‑1 (MQ=255)
                                               cACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTcg                  <  1:2323276/71‑1 (MQ=255)
                                               cACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTcg                  <  1:898901/71‑1 (MQ=255)
                                               cACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTcg                  <  1:1871126/71‑1 (MQ=255)
                                               cACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTcg                  <  1:1905807/71‑1 (MQ=255)
                                                aCCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATc                                          >  1:2673645/1‑46 (MQ=255)
                                                aCCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATc                                          >  1:2558999/1‑46 (MQ=255)
                                                aCCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATc                                          >  1:685148/1‑46 (MQ=255)
                                                aCCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATc                                          >  1:1927402/1‑46 (MQ=255)
                                                 ccGCGTAGCGTTTCCAGCATATGCTCCTTCAGACGa                                                   <  1:1242799/36‑1 (MQ=255)
                                                 ccGCGTAGCGTTTCCAGCATATGCTCCTTCAGACGa                                                   <  1:1185896/36‑1 (MQ=255)
                                                 ccGCGTAGCGTTTCCAGCATATGCTCCTTCAGACGa                                                   <  1:837445/36‑1 (MQ=255)
                                                 ccGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTcgcg                <  1:1590980/71‑1 (MQ=255)
                                                 ccGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTcgcg                <  1:1586846/71‑1 (MQ=255)
                                                 ccGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTcgc                 >  1:1455604/1‑70 (MQ=255)
                                                   gcgTAGCGTTTCCAGCATATGCTCCTTCAGACGaaa                                                 >  1:1948415/1‑36 (MQ=255)
                                                   gcgTAGCGTTTCCAGCATATGCTCCTTCAGACGaaa                                                 >  1:676077/1‑36 (MQ=255)
                                                   gcgTAGCGTTTCCAGCATATGCTCCTTCAGACGaaa                                                 >  1:745774/1‑36 (MQ=255)
                                                   gcgTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGa               >  1:3130481/1‑70 (MQ=255)
                                                   gcgTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGa               >  1:2701499/1‑70 (MQ=255)
                                                      tAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCa                                             <  1:636790/37‑1 (MQ=255)
                                                      tAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCa                                             <  1:1210415/37‑1 (MQ=255)
                                                      tAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATc                                          <  1:2330495/40‑1 (MQ=255)
                                                      tAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATc                                          <  1:715156/40‑1 (MQ=255)
                                                      tAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATc                                          <  1:2615510/40‑1 (MQ=255)
                                                        gCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTg                                        >  1:580392/1‑40 (MQ=255)
                                                        gCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTcgc                 <  1:1976123/63‑1 (MQ=255)
                                                        gCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGAt              >  1:752686/1‑66 (MQ=255)
                                                        gCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGAt              >  1:989743/1‑66 (MQ=255)
                                                        gCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGAt              >  1:3201235/1‑66 (MQ=255)
                                                        gCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGAt              >  1:2992895/1‑66 (MQ=255)
                                                        gCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGAt              >  1:152282/1‑66 (MQ=255)
                                                         cGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCgg                      <  1:1523648/57‑1 (MQ=255)
                                                         cGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGa               >  1:2655848/1‑64 (MQ=255)
                                                         cGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGa        <  1:3443916/71‑1 (MQ=255)
                                                         cGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGACGGACTCTTTTACCGGTTCGCGATCCCAGa        <  1:3118631/71‑1 (MQ=255)
                                                           tttCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATcc            >  1:731932/1‑65 (MQ=255)
                                                           tttCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATcc            >  1:981010/1‑65 (MQ=255)
                                                           tttCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATa      >  1:2418488/1‑71 (MQ=255)
                                                               cAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCgg                                   >  1:853434/1‑38 (MQ=255)
                                                               cAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCgg                                   >  1:1914761/1‑38 (MQ=255)
                                                               cAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGt  >  1:3154711/1‑71 (MQ=255)
                                                                 gCAGATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTcg                  <  1:2857445/53‑1 (MQ=255)
                                                                   acaTGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCa          <  1:2009034/59‑1 (MQ=255)
                                                                    |                                                                 
ATCCGGATTAGTGACTTTATCCCAGCCCTTACCTTCACAGGATGTCGCACCGCGTAGCGTTTCCAGCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGT  >  minE/1910337‑1910470

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: