Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,127 G→A 9.5% S157S (TCC→TCT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1270GA 9.5% 221.4 / 10.5 74S157S (TCC→TCTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (34/32);  new base A (4/3);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTG  >  minE/2596057‑2596196
                                                                      |                                                                     
gCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                                                       <  1:1051795/71‑1 (MQ=255)
gCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                                                       <  1:200746/71‑1 (MQ=255)
 cACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTgg                                                                      <  1:310067/71‑1 (MQ=255)
  aCCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                                                     <  1:2253877/71‑1 (MQ=255)
     aTCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                                                  <  1:1835696/71‑1 (MQ=255)
     aTCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACAAGGTTGgaga                                                                   >  1:441133/1‑70 (MQ=255)
      tCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                                                  >  1:286727/1‑70 (MQ=255)
      tCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                                                  >  1:2976417/1‑70 (MQ=255)
      tCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                                                  <  1:3843012/70‑1 (MQ=255)
      tCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                                                  >  1:1716669/1‑70 (MQ=255)
        aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACaa                                                                <  1:3558223/70‑1 (MQ=255)
        aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACaa                                                                >  1:2273499/1‑70 (MQ=255)
        aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACaa                                                                >  1:2293903/1‑70 (MQ=255)
        aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACaa                                                                >  1:321323/1‑70 (MQ=255)
            aGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGgaga                                                                   <  1:1119755/63‑1 (MQ=255)
              gCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                                                       >  1:1275315/1‑57 (MQ=255)
                gTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCa                                                       <  1:2965322/71‑1 (MQ=255)
                 tATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCgg                                                         <  1:1864799/68‑1 (MQ=255)
                   ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                    >  1:3602550/1‑71 (MQ=255)
                   ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                    >  1:425572/1‑71 (MQ=255)
                   ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                    >  1:2903001/1‑71 (MQ=255)
                   ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                    >  1:481599/1‑71 (MQ=255)
                   ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                    >  1:1105007/1‑71 (MQ=255)
                    tCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATaa                                                             <  1:3453981/61‑1 (MQ=255)
                     cgcgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTg                                                  >  1:3162516/1‑71 (MQ=255)
                      gcgAATGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTgg                                                 <  1:80226/71‑1 (MQ=255)
                      cgcaaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTgg                                                 >  1:3307063/4‑71 (MQ=255)
                              cGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTCGAGACAATAAGCGGCAGGCTGGc                                                <  1:2813100/64‑1 (MQ=255)
                                aGGCCAATGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                       <  1:3031067/71‑1 (MQ=255)
                                aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                       <  1:1867257/71‑1 (MQ=255)
                                aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                       <  1:660342/71‑1 (MQ=255)
                                  gCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACaa                                                                >  1:3831951/1‑44 (MQ=255)
                                   ccAAAGAAATCAGCGGAAACGATATTCTCCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                    >  1:2133069/1‑71 (MQ=255)
                                   ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                    <  1:3797684/71‑1 (MQ=255)
                                   ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                    <  1:3334549/71‑1 (MQ=255)
                                    cAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                                                     >  1:3707010/1‑37 (MQ=255)
                                    cAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCgg                                        >  1:3883718/1‑66 (MQ=255)
                                     aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGa                                  <  1:2125149/71‑1 (MQ=255)
                                     aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGa                                  <  1:3242990/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)
                                        gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAg                                                            >  1:515603/1‑42 (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)
                                                cGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAg                                                      >  1:1870491/1‑40 (MQ=255)
                                                cGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAg                                                      >  1:761210/1‑40 (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)
                                                     aCGATATTCACCAGGTTGGAGACAATAAGCGGCAg                                                      <  1:580537/35‑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:806710/1‑69 (MQ=255)
                                                             cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc            >  1:3254166/1‑69 (MQ=255)
                                                             cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc            >  1:591246/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:3551316/1‑70 (MQ=255)
                                                                      ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg  <  1:2050154/70‑1 (MQ=255)
                                                                      |                                                                     
GCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTG  >  minE/2596057‑2596196

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: