Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,530 A→G 28.1% V356V (GTT→GTC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5300AG28.1% 268.6 / 129.8 185V356V (GTT→GTCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (77/56);  new base G (37/15);  total (114/71)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 189 total aligned reads displayed.

CCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAA  >  minE/2595462‑2595597
                                                                    |                                                                   
  aaTCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                                                                 <  1:2905130/71‑1 (MQ=255)
    ttGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                                                                 >  1:3177689/1‑69 (MQ=255)
     tGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                                                                 >  1:4094631/1‑68 (MQ=255)
      ggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgt                                                             >  1:4442622/1‑70 (MQ=255)
      ggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                                                              <  1:3745095/70‑1 (MQ=255)
      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                                                                    >  1:1869377/1‑64 (MQ=255)
      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                                                                    >  1:2880212/1‑64 (MQ=255)
      ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                                                                    >  1:623043/1‑64 (MQ=255)
       gggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                                                            >  1:5121328/1‑69 (MQ=255)
       gggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                             <  1:720008/70‑1 (MQ=255)
       gggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                             <  1:2410889/70‑1 (MQ=255)
          tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                         >  1:4639318/1‑71 (MQ=255)
           aaTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg                                                        >  1:1316500/1‑71 (MQ=255)
            aTTTTCGGTCCCAAATCTCAGCCAATCACATTGGCATAAACAATCTATTCTTTGATAACGCCAGATGcc                                                         >  1:5660729/1‑69 (MQ=255)
            aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                                                                 <  1:3192559/61‑1 (MQ=255)
             ttttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                       <  1:2214288/70‑1 (MQ=255)
             ttttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGta                                                      >  1:870299/1‑70 (MQ=255)
             ttttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGta                                                      >  1:3302211/1‑70 (MQ=255)
              tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGACAGATGCCGtg                                                      <  1:1286167/70‑1 (MQ=255)
              tttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtac                                                     >  1:947643/1‑69 (MQ=255)
               ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                    >  1:1945590/1‑71 (MQ=255)
                    tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt                                               >  1:5619693/1‑71 (MQ=255)
                     cctaaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTc                                              <  1:2033291/68‑1 (MQ=255)
                     cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                         >  1:3866857/1‑60 (MQ=255)
                     cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                         >  1:3994363/1‑60 (MQ=255)
                     cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                         >  1:50638/1‑60 (MQ=255)
                     cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                         >  1:5393428/1‑60 (MQ=255)
                      ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                              <  1:3949922/70‑1 (MQ=255)
                      ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                              <  1:5299032/70‑1 (MQ=255)
                        aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                     >  1:1087065/1‑61 (MQ=255)
                        aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                           >  1:2927406/1‑71 (MQ=255)
                         aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                          <  1:334155/71‑1 (MQ=255)
                         aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                                                            >  1:269945/1‑51 (MQ=255)
                         aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAtc                                          <  1:5437133/71‑2 (MQ=255)
                           tcacAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                        >  1:2392320/4‑71 (MQ=255)
                           tCGCAGCCAATCACATTGGCATAAACAATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                        >  1:167266/1‑71 (MQ=255)
                             gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                       <  1:4587681/70‑1 (MQ=255)
                             gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                       <  1:3386210/70‑1 (MQ=255)
                             acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGAc                                       <  1:1526969/69‑1 (MQ=255)
                             acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACa                                      <  1:214923/70‑1 (MQ=255)
                              cAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACag                                     <  1:1179477/71‑2 (MQ=255)
                              cAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                           <  1:84031/65‑1 (MQ=255)
                               aGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                       <  1:4159363/52‑1 (MQ=255)
                               aGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                       <  1:2764875/52‑1 (MQ=255)
                                 ccAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgc                                  <  1:3262228/71‑1 (MQ=255)
                                 ccAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                     >  1:1266942/1‑52 (MQ=255)
                                 ccAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgc                                  >  1:4701257/1‑71 (MQ=255)
                                    aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                        >  1:1608917/1‑62 (MQ=255)
                                    aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                        >  1:2176540/1‑62 (MQ=255)
                                    aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                        >  1:3782538/1‑62 (MQ=255)
                                    aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                               <  1:2079645/71‑1 (MQ=255)
                                     tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                        >  1:2594187/1‑61 (MQ=255)
                                     tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                             >  1:2980187/1‑40 (MQ=255)
                                     tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                             >  1:3407748/1‑40 (MQ=255)
                                     tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCa                              <  1:2781269/71‑1 (MQ=255)
                                      cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATgg                                         >  1:5564043/1‑59 (MQ=255)
                                      cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                             >  1:5013667/1‑71 (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:1037860/71‑1 (MQ=255)
                                       acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                            <  1:5472608/71‑1 (MQ=255)
                                       acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                            <  1:3254161/71‑1 (MQ=255)
                                       acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                            >  1:3459158/1‑71 (MQ=255)
                                         aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                        >  1:823567/1‑57 (MQ=255)
                                         aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                        >  1:5175147/1‑57 (MQ=255)
                                         aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                          <  1:4849085/71‑1 (MQ=255)
                                          ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtactc                                                   >  1:2718459/1‑44 (MQ=255)
                                          ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCaa                                            >  1:79368/1‑52 (MQ=255)
                                          ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                           >  1:2988091/1‑53 (MQ=255)
                                          ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                           >  1:3608310/1‑53 (MQ=255)
                                          ttGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAg                         <  1:1979696/71‑1 (MQ=255)
                                          ttGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAg                         <  1:3909527/71‑1 (MQ=255)
                                           tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                            >  1:2678707/1‑51 (MQ=255)
                                           tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                        >  1:1689067/1‑55 (MQ=255)
                                           tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                        >  1:2541286/1‑55 (MQ=255)
                                            ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                            >  1:3142306/1‑50 (MQ=255)
                                            ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGt                        <  1:4262823/70‑1 (MQ=255)
                                            ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGt                        <  1:1472548/70‑1 (MQ=255)
                                             gCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                       <  1:3422995/70‑1 (MQ=255)
                                                tAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                       <  1:1312794/35‑1 (MQ=255)
                                                tAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                       <  1:2900094/35‑1 (MQ=255)
                                                tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                          <  1:135212/48‑1 (MQ=255)
                                                tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                   >  1:1124572/1‑71 (MQ=255)
                                                 aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                  >  1:2560031/1‑71 (MQ=255)
                                                  aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                 >  1:611440/1‑55 (MQ=255)
                                                  aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                  <  1:4005571/70‑1 (MQ=255)
                                                  aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                 >  1:4571849/1‑71 (MQ=255)
                                                  aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                 >  1:1911825/1‑71 (MQ=255)
                                                  aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                 >  1:4762948/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)
                                                         gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt          <  1:4176290/71‑1 (MQ=255)
                                                         gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt          >  1:285528/1‑71 (MQ=255)
                                                             cttTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGt        <  1:1743118/69‑1 (MQ=255)
                                                              tttGATAACGCCAGATGCCGTGCTGCCATCAATGGACa                                      >  1:2753541/1‑38 (MQ=255)
                                                               ttGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                >  1:367765/1‑43 (MQ=255)
                                                                 gATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataa  >  1:900912/1‑71 (MQ=255)
                                                                    |                                                                   
CCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAA  >  minE/2595462‑2595597

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: