Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,130 G→A 9.6% V156V (GTC→GTT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1300GA 9.6% 221.5 / 11.2 73V156V (GTC→GTTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (34/32);  new base A (4/3);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTT  >  minE/2596062‑2596200
                                                                    |                                                                      
aTCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                                                      <  1:1835696/71‑1 (MQ=255)
aTCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACAAGGTTGgaga                                                                       >  1:441133/1‑70 (MQ=255)
 tCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                                                      <  1:3843012/70‑1 (MQ=255)
 tCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                                                      >  1:2976417/1‑70 (MQ=255)
 tCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                                                      >  1:1716669/1‑70 (MQ=255)
 tCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                                                      >  1:286727/1‑70 (MQ=255)
   aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACaa                                                                    <  1:3558223/70‑1 (MQ=255)
   aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACaa                                                                    >  1:321323/1‑70 (MQ=255)
   aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACaa                                                                    >  1:2293903/1‑70 (MQ=255)
   aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACaa                                                                    >  1:2273499/1‑70 (MQ=255)
       aGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGgaga                                                                       <  1:1119755/63‑1 (MQ=255)
           gTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCa                                                           <  1:2965322/71‑1 (MQ=255)
            tATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCgg                                                             <  1:1864799/68‑1 (MQ=255)
              ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                        >  1:481599/1‑71 (MQ=255)
              ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                        >  1:1105007/1‑71 (MQ=255)
              ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                        >  1:425572/1‑71 (MQ=255)
              ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                        >  1:2903001/1‑71 (MQ=255)
              ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                        >  1:3602550/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: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)
                                   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:761210/1‑40 (MQ=255)
                                           cGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAg                                                          >  1:1870491/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:591246/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:3551316/1‑70 (MQ=255)
                                                                 ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg      <  1:2050154/70‑1 (MQ=255)
                                                                  gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGcc    <  1:2140427/71‑1 (MQ=255)
                                                                    gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt  <  1:309246/71‑1 (MQ=255)
                                                                    |                                                                      
ATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTT  >  minE/2596062‑2596200

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: