Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,136 A→G 11.3% L154L (CTT→CTC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1360AG11.3% 192.0 / 14.8 71L154L (CTT→CTCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base G (5/3);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

GTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTG  >  minE/2596073‑2596205
                                                               |                                                                     
gTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCa                                                                <  1:2965322/71‑1 (MQ=255)
 tATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCgg                                                                  <  1:1864799/68‑1 (MQ=255)
   ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                             >  1:1105007/1‑71 (MQ=255)
   ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                             >  1:2903001/1‑71 (MQ=255)
   ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                             >  1:3602550/1‑71 (MQ=255)
   ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                             >  1:425572/1‑71 (MQ=255)
   ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGc                                                             >  1:481599/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)
                   ccAAAGAAATCAGCGGAAACGATATTCTCCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                             >  1:2133069/1‑71 (MQ=255)
                   ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                             <  1:3797684/71‑1 (MQ=255)
                   ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                             <  1:3334549/71‑1 (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:2964234/1‑71 (MQ=255)
                        gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                        >  1:1136039/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: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: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:1714446/1‑71 (MQ=255)
                                                           cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTa     <  1:3885421/71‑1 (MQ=255)
                                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt   <  1:2807375/70‑1 (MQ=255)
                                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt   >  1:2121772/1‑70 (MQ=255)
                                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTg  <  1:3601908/71‑1 (MQ=255)
                                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTg  <  1:334970/71‑1 (MQ=255)
                                                               aaGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCATCGTAGTGCCTTTACTg  >  1:2710512/1‑70 (MQ=255)
                                                               |                                                                     
GTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTG  >  minE/2596073‑2596205

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: