Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,002,444 C→A 7.6% intergenic (+91/+32) yghU → / ← hybG predicted S‑transferase/hydrogenase 2 accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,002,4440CA 7.6% 286.5 / 5.6 92intergenic (+91/+32)yghU/hybGpredicted S‑transferase/hydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base C (44/40);  new base A (3/4);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

AGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGT  >  minE/2002375‑2002512
                                                                       |                                                                    
aGCCATCCAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGcc                                                                     >  1:300041/1‑71 (MQ=255)
aGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGc                                                                      >  1:5213317/1‑70 (MQ=255)
   cATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCcgcg                                                                  >  1:3412365/1‑71 (MQ=255)
       taGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCt                                                               >  1:1747449/2‑71 (MQ=11)
       cAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCt                                                               >  1:547717/1‑70 (MQ=255)
        aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                             <  1:2424158/71‑1 (MQ=255)
         ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCtt                                                              <  1:2369869/70‑1 (MQ=1)
         ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTa                                                             >  1:5575257/1‑71 (MQ=1)
         ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                                                            <  1:398235/71‑1 (MQ=255)
         ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                                                            <  1:1976584/71‑1 (MQ=255)
          gCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                                                            <  1:366548/71‑1 (MQ=1)
           ccGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCa                                                          <  1:2557997/71‑1 (MQ=255)
              gATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAgg                                                        >  1:103501/1‑71 (MQ=1)
               aTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAgac                                                       <  1:2140569/71‑3 (MQ=1)
               aTAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGcc                                                      >  1:299414/1‑71 (MQ=255)
                 aaGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTa                                                    >  1:502717/1‑71 (MQ=255)
                  aGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCt                                                     <  1:143691/70‑1 (MQ=11)
                   ggCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTAc                                                   <  1:4056090/71‑1 (MQ=1)
                    gCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGa                                                 <  1:5425052/71‑1 (MQ=255)
                     cGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaa                                                <  1:4972006/71‑1 (MQ=255)
                     cGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaa                                                <  1:2040050/71‑1 (MQ=255)
                        tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                             >  1:2818161/1‑71 (MQ=255)
                        tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                             >  1:445070/1‑71 (MQ=255)
                        tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                             >  1:2566937/1‑71 (MQ=255)
                        tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                             >  1:5073097/1‑71 (MQ=255)
                        tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                             >  1:1811793/1‑71 (MQ=255)
                         cACGCCGCATCCGGCAATCGTTCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                             <  1:2863666/70‑1 (MQ=255)
                         cACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGa                                                 >  1:448828/1‑66 (MQ=255)
                         cACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                             >  1:1488499/1‑70 (MQ=255)
                          aCGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCt                                            <  1:270106/70‑1 (MQ=255)
                          aCGCCGCATCCGGCAATCGAGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTa                                           <  1:3990030/71‑1 (MQ=255)
                             ccgcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                                          >  1:3214259/1‑69 (MQ=255)
                             ccgcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                                        <  1:5619053/71‑1 (MQ=255)
                               gcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                                      <  1:1283253/71‑1 (MQ=255)
                               gcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                                      <  1:1862813/71‑1 (MQ=255)
                                   ccGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcg                                  >  1:853808/1‑71 (MQ=255)
                                    cGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                                          >  1:1066588/1‑62 (MQ=255)
                                    cGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgc                                 >  1:3899269/1‑71 (MQ=255)
                                    cGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgc                                 >  1:2264226/1‑71 (MQ=255)
                                      gCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACgcgc                               >  1:313665/1‑71 (MQ=255)
                                       cAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCt                              >  1:3221371/1‑71 (MQ=255)
                                       cAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCt                              >  1:4051214/1‑71 (MQ=255)
                                        aaTCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaa                                                <  1:4997144/52‑1 (MQ=255)
                                        aaTCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTg                             >  1:4606671/1‑71 (MQ=255)
                                         aTCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg                            <  1:1021144/71‑1 (MQ=255)
                                          tCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                                        >  1:5196056/1‑58 (MQ=255)
                                          tCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg                            >  1:4487556/1‑70 (MQ=255)
                                          tCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGt                           <  1:5235391/71‑1 (MQ=255)
                                           cGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                                      >  1:4082724/1‑59 (MQ=255)
                                           cGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgcg                                >  1:3916511/1‑65 (MQ=255)
                                           cGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGt                           <  1:5005017/70‑1 (MQ=255)
                                            gTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatca                                       <  1:307170/57‑1 (MQ=255)
                                            gTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                         >  1:2209611/1‑71 (MQ=255)
                                             tGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAAt                        <  1:3547729/71‑1 (MQ=255)
                                               cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg                       <  1:3538600/70‑1 (MQ=255)
                                               cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg                       <  1:3390203/70‑1 (MQ=255)
                                               cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg                       <  1:2037908/70‑1 (MQ=255)
                                               cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                      <  1:887715/71‑1 (MQ=255)
                                               cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                      <  1:3471460/71‑1 (MQ=255)
                                               cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                      <  1:2162253/71‑1 (MQ=255)
                                               cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                      >  1:4828888/1‑71 (MQ=255)
                                                acaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                      <  1:5081828/70‑1 (MQ=255)
                                                 caATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACgcgc                               <  1:1787944/60‑1 (MQ=255)
                                                 caATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTg                             >  1:4726473/1‑62 (MQ=255)
                                                 caATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg                    <  1:615478/71‑1 (MQ=255)
                                                  aaTGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg                    >  1:3729234/1‑70 (MQ=255)
                                                    aTGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg                            <  1:1094859/61‑1 (MQ=255)
                                                    aTGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg                            <  1:3056775/61‑1 (MQ=255)
                                                      gCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                                      <  1:1928746/49‑1 (MQ=255)
                                                      gCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                         >  1:3842880/1‑62 (MQ=255)
                                                      gCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAAt                >  1:3180678/1‑71 (MQ=255)
                                                       ccTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAAt                <  1:4275079/70‑1 (MQ=255)
                                                       ccTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATc               >  1:499980/1‑71 (MQ=255)
                                                       ccTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATc               >  1:506836/1‑71 (MQ=255)
                                                           aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAg                                              >  1:5421369/1‑36 (MQ=255)
                                                           aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                                          <  1:16249/40‑1 (MQ=255)
                                                           aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCa             <  1:5367826/69‑1 (MQ=255)
                                                           aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCa             <  1:522617/69‑1 (MQ=255)
                                                           aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCa             >  1:4843897/1‑69 (MQ=255)
                                                           aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAtt           >  1:179636/1‑71 (MQ=255)
                                                           aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAt            <  1:3104497/70‑1 (MQ=255)
                                                            tGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                                      <  1:2751855/43‑1 (MQ=255)
                                                            tGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAttt          >  1:1786443/1‑71 (MQ=255)
                                                             gCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgc                                 >  1:5508818/1‑47 (MQ=255)
                                                             gCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                         <  1:5073057/55‑1 (MQ=255)
                                                             gCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                         <  1:369893/55‑1 (MQ=255)
                                                             gCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTg         >  1:1637329/1‑71 (MQ=255)
                                                             gCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTg         >  1:2480610/1‑71 (MQ=255)
                                                              cGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAAtcatca                                    >  1:1548345/1‑43 (MQ=255)
                                                                 cGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                                        >  1:1670158/1‑36 (MQ=255)
                                                                 cGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                                        >  1:4582596/1‑36 (MQ=255)
                                                                 cGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                         >  1:4973019/1‑51 (MQ=255)
                                                                 cGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGgcgc     <  1:4883610/71‑1 (MQ=255)
                                                                 cGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGgcg      >  1:3603018/1‑70 (MQ=255)
                                                                       ccGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCa             <  1:1994584/58‑1 (MQ=255)
                                                                       ccGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGt  <  1:2823075/69‑1 (MQ=255)
                                                                       |                                                                    
AGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGT  >  minE/2002375‑2002512

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: