Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,151 G→C 10.9% T149T (ACC→ACG arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1510GC10.9% 177.7 / 11.1 64T149T (ACC→ACGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base C (4/3);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

AGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGC  >  minE/2596089‑2596216
                                                              |                                                                 
aGGCCAATGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                                           <  1:3031067/71‑1 (MQ=255)
aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                                           <  1:660342/71‑1 (MQ=255)
aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                                           <  1:1867257/71‑1 (MQ=255)
   ccAAAGAAATCAGCGGAAACGATATTCTCCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                                        >  1:2133069/1‑71 (MQ=255)
   ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                                        <  1:3334549/71‑1 (MQ=255)
   ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                                        <  1:3797684/71‑1 (MQ=255)
    cAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCgg                                                            >  1:3883718/1‑66 (MQ=255)
     aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGa                                                      <  1:3242990/71‑1 (MQ=255)
     aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGa                                                      <  1:2125149/71‑1 (MQ=255)
      aagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATaaa                                                     >  1:750249/1‑71 (MQ=255)
      aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                                                                  <  1:920293/58‑1 (MQ=255)
      aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa                                                     <  1:3136089/71‑1 (MQ=255)
       aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAt                                                    <  1:2531797/71‑1 (MQ=255)
        gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGt                                                                 <  1:3708694/57‑1 (MQ=255)
        gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                   >  1:1136039/1‑71 (MQ=255)
        gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                   >  1:2964234/1‑71 (MQ=255)
         aaaTCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                   <  1:2731646/70‑1 (MQ=255)
              aGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                                                                  <  1:2174770/50‑1 (MQ=255)
               gCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                                                                  >  1:3634524/1‑49 (MQ=255)
               gCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgccgc                                            >  1:2606000/1‑71 (MQ=255)
                 ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCa                                          >  1:437794/1‑71 (MQ=255)
                 ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCa                                          >  1:2265911/1‑71 (MQ=255)
                 ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCa                                          <  1:2142609/71‑1 (MQ=255)
                  gAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                         <  1:1563665/71‑1 (MQ=255)
                       gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                    <  1:3354604/71‑3 (MQ=255)
                       gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                    <  1:3350762/71‑3 (MQ=255)
                        atatTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCg                                                             <  1:2202926/45‑1 (MQ=255)
                           ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                                           <  1:3358153/44‑1 (MQ=255)
                           ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                         <  1:3587156/62‑1 (MQ=255)
                            tCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                >  1:2080108/1‑70 (MQ=255)
                             cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                >  1:591246/1‑69 (MQ=255)
                             cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                >  1:806710/1‑69 (MQ=255)
                             cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                >  1:3254166/1‑69 (MQ=255)
                              aCCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAg                             >  1:338330/1‑71 (MQ=255)
                                cAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg                           <  1:1680970/71‑1 (MQ=255)
                                cAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGACAGc                            >  1:256378/1‑70 (MQ=255)
                                cAGGTTGGAGACAATAAGCCGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGc                            >  1:3229683/1‑70 (MQ=255)
                                 aGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGt                          <  1:411312/71‑1 (MQ=255)
                                 aGGTTAGAAACAATGCGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGc                                >  1:958970/1‑65 (MQ=255)
                                 aGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGc                                >  1:1915426/1‑65 (MQ=255)
                                 aGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgt                          <  1:128093/71‑1 (MQ=255)
                                    ttGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg                           <  1:3983386/67‑1 (MQ=255)
                                    ttGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGt                       <  1:2078208/71‑1 (MQ=255)
                                     tGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg                      >  1:3775205/1‑71 (MQ=255)
                                      ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                   >  1:32683/1‑41 (MQ=255)
                                      ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg                      <  1:2050154/70‑1 (MQ=255)
                                      ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg                      >  1:3551316/1‑70 (MQ=255)
                                       gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGcc                    <  1:2140427/71‑1 (MQ=255)
                                         gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                  <  1:309246/71‑1 (MQ=255)
                                          aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                 >  1:3790992/1‑68 (MQ=255)
                                           cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTa                <  1:3885421/71‑1 (MQ=255)
                                           cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTa                >  1:1714446/1‑71 (MQ=255)
                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt              >  1:2121772/1‑70 (MQ=255)
                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt              <  1:2807375/70‑1 (MQ=255)
                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTg             <  1:3601908/71‑1 (MQ=255)
                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTg             <  1:334970/71‑1 (MQ=255)
                                               aaGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCATCGTAGTGCCTTTACTg             >  1:2710512/1‑70 (MQ=255)
                                                aGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTa                         >  1:1080377/1‑57 (MQ=255)
                                                aGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGaa           >  1:2727881/1‑71 (MQ=255)
                                                 gCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAAc          >  1:2215050/1‑71 (MQ=255)
                                                 gCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAAc          >  1:1842847/1‑71 (MQ=255)
                                                     cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa      <  1:2701181/71‑1 (MQ=255)
                                                     cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa      <  1:3277368/71‑1 (MQ=255)
                                                     cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa      <  1:92884/71‑1 (MQ=255)
                                                        gCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAg   <  1:1033214/71‑1 (MQ=255)
                                                         cTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGc  <  1:470088/71‑1 (MQ=255)
                                                         cTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGc  <  1:3423796/71‑1 (MQ=255)
                                                              ggTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAg   <  1:2002415/65‑1 (MQ=255)
                                                              |                                                                 
AGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGC  >  minE/2596089‑2596216

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: