Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,538 A→T 24.9% S354T (TCT→ACT)  arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5380AT24.9% 328.1 / 122.3 181S354T (TCT→ACT) arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (71/65);  new base T (32/13);  total (103/78)
Fisher's exact test for biased strand distribution p-value = 3.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 185 total aligned reads displayed.

GGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAG  >  minE/2595469‑2595607
                                                                     |                                                                     
gggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                                       <  1:720008/70‑1 (MQ=255)
gggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                                       <  1:2410889/70‑1 (MQ=255)
   tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                                   <  1:848510/71‑1 (MQ=255)
      ttttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                 <  1:3674096/70‑1 (MQ=255)
       tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGACAGATGCCGtg                                                                <  1:1286167/70‑1 (MQ=255)
        ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                              >  1:1945590/1‑71 (MQ=255)
         tCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                              >  1:803359/1‑70 (MQ=255)
              cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                                   >  1:3866857/1‑60 (MQ=255)
              cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                                   >  1:50638/1‑60 (MQ=255)
              cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                                   >  1:5393428/1‑60 (MQ=255)
              cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt                                                         >  1:3396011/1‑70 (MQ=255)
               ctaaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTc                                                        >  1:2891717/3‑70 (MQ=255)
               ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                        <  1:3367170/70‑1 (MQ=255)
               ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                        <  1:1722350/70‑1 (MQ=255)
                 aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                               >  1:1087065/1‑61 (MQ=255)
                 aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                     <  1:833011/71‑1 (MQ=255)
                 aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                     >  1:2927406/1‑71 (MQ=255)
                 aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                     <  1:4455755/71‑1 (MQ=255)
                  aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg                                                                  <  1:131323/57‑1 (MQ=255)
                  aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg                                                                  <  1:2150286/57‑1 (MQ=255)
                  aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAtc                                                    <  1:911476/71‑2 (MQ=255)
                  aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAtc                                                    <  1:5437133/71‑2 (MQ=255)
                  aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                     >  1:2903195/1‑70 (MQ=255)
                    tcacAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                  >  1:2392320/4‑71 (MQ=255)
                      gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                 <  1:4587681/70‑1 (MQ=255)
                      gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                 <  1:871666/70‑1 (MQ=255)
                      acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                  >  1:2101087/2‑69 (MQ=255)
                      acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGAc                                                 <  1:1526969/69‑1 (MQ=255)
                      acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGAc                                                 <  1:3159032/69‑1 (MQ=255)
                       cAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                                               <  1:770904/71‑1 (MQ=255)
                        aGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                 <  1:4159363/52‑1 (MQ=255)
                         gCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                             <  1:2046393/71‑1 (MQ=255)
                         gCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                             >  1:4631505/1‑71 (MQ=255)
                          ccAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgc                                            <  1:3262228/71‑1 (MQ=255)
                          ccAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                               >  1:3954659/1‑52 (MQ=255)
                          ccAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgc                                            >  1:4701257/1‑71 (MQ=255)
                             aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                  >  1:2181547/1‑62 (MQ=255)
                             aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                  >  1:2777726/1‑62 (MQ=255)
                              tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                  >  1:2221376/1‑61 (MQ=255)
                              tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                  >  1:2594187/1‑61 (MQ=255)
                              tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                  >  1:4010419/1‑61 (MQ=255)
                              tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                                       >  1:3407748/1‑40 (MQ=255)
                              tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                                       >  1:2980187/1‑40 (MQ=255)
                               cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctgc                                                            >  1:1351659/1‑50 (MQ=255)
                               cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                       >  1:464410/1‑71 (MQ=255)
                                acaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgca                                           >  1:3832121/1‑65 (MQ=255)
                                acaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgca                                           >  1:223661/1‑65 (MQ=255)
                                acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                       <  1:466127/70‑1 (MQ=255)
                                acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                       >  1:3466013/1‑70 (MQ=255)
                                acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                      >  1:157761/1‑71 (MQ=255)
                                acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                      <  1:5472608/71‑1 (MQ=255)
                                acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                      <  1:1037860/71‑1 (MQ=255)
                                acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                      <  1:3254161/71‑1 (MQ=255)
                                 caTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAcc                                     <  1:1634087/71‑1 (MQ=255)
                                  aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                                  >  1:5175147/1‑57 (MQ=255)
                                  aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                                  >  1:168534/1‑57 (MQ=255)
                                  aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                    <  1:4849085/71‑1 (MQ=255)
                                   ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                     >  1:3608310/1‑53 (MQ=255)
                                   ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                     >  1:2988091/1‑53 (MQ=255)
                                   ttGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAg                                   <  1:1979696/71‑1 (MQ=255)
                                    tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                      >  1:2678707/1‑51 (MQ=255)
                                    tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                       >  1:3920017/1‑66 (MQ=255)
                                    tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                       >  1:3055542/1‑66 (MQ=255)
                                     ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                      >  1:3308327/1‑50 (MQ=255)
                                     ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGt                                  <  1:1822723/70‑1 (MQ=255)
                                     ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGt                                  <  1:4262823/70‑1 (MQ=255)
                                      gCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                                 <  1:3422995/70‑1 (MQ=255)
                                         tAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                 <  1:1312794/35‑1 (MQ=255)
                                         tAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                 <  1:2747557/35‑1 (MQ=255)
                                         tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                             >  1:1124572/1‑71 (MQ=255)
                                          aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAc                                >  1:1249439/1‑67 (MQ=255)
                                          aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                            >  1:2560031/1‑71 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                                               >  1:2762815/1‑51 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                           >  1:2718291/1‑55 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                           >  1:3736814/1‑55 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                           >  1:611440/1‑55 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                           >  1:5035245/1‑55 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                            <  1:4005571/70‑1 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                           >  1:190643/1‑71 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGAGCTGCCATCAATGGACAACGCGCCAAcc                                     >  1:4243777/1‑61 (MQ=255)
                                              cATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                        <  1:769114/39‑1 (MQ=255)
                                               aTCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                      >  1:4003632/1‑40 (MQ=255)
                                                  gcctcttTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCatat                     >  1:300896/4‑70 (MQ=255)
                                                  gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                     >  1:5639749/1‑38 (MQ=255)
                                                  gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt                    >  1:285528/1‑71 (MQ=255)
                                                  gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt                    <  1:4176290/71‑1 (MQ=255)
                                                        ttGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                          >  1:1388473/1‑43 (MQ=255)
                                                         tGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCat              <  1:278853/70‑1 (MQ=255)
                                                           aTAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGt                  >  1:5610185/1‑64 (MQ=255)
                                                              aCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCat              <  1:5301797/65‑1 (MQ=255)
                                                              aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATaga        <  1:1448564/71‑1 (MQ=255)
                                                              aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag         <  1:4980414/70‑1 (MQ=255)
                                                              aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag         <  1:5685742/70‑1 (MQ=255)
                                                               cGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaa       <  1:3890566/71‑1 (MQ=255)
                                                                gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                    >  1:2674055/1‑41 (MQ=255)
                                                                gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                    >  1:2004573/1‑41 (MQ=255)
                                                                gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaag      <  1:1003764/71‑1 (MQ=255)
                                                                 ccAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                               <  1:2483012/45‑1 (MQ=255)
                                                                 ccAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga     <  1:1625058/71‑1 (MQ=255)
                                                                    gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag  <  1:5283863/71‑1 (MQ=255)
                                                                     |                                                                     
GGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAG  >  minE/2595469‑2595607

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: