Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,521 T→C 26.2% E359E (GAA→GAG arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5210TC26.2% 299.1 / 123.0 191E359E (GAA→GAGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base T (78/63);  new base C (36/14);  total (114/77)
Fisher's exact test for biased strand distribution p-value = 4.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 201 total aligned reads displayed.

AGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATT  >  minE/2595452‑2595589
                                                                     |                                                                    
agcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGct                                                                      >  1:117541/1‑70 (MQ=255)
agcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAACCGCAGCCAATCACATTGGCATAAACCATCGctt                                                                     >  1:1749143/1‑71 (MQ=255)
 gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctc                                                                    >  1:138197/1‑71 (MQ=255)
    aGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                                                                 <  1:1136169/71‑1 (MQ=255)
      gCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGa                                                               >  1:3076353/1‑71 (MQ=255)
        tACCAATTGGGGTAATTTTAGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGAt                                                              >  1:3327418/1‑70 (MQ=255)
          ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCATTCACATTGGCATAAACCATCGCTTCTTTGa                                                               >  1:3344975/1‑67 (MQ=255)
          ccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGacg                                                             <  1:530882/69‑3 (MQ=255)
               tGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                                                         >  1:4094631/1‑68 (MQ=255)
                ggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgt                                                     >  1:1923256/1‑70 (MQ=255)
                ggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgt                                                     >  1:4442622/1‑70 (MQ=255)
                ggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                                                      <  1:3745095/70‑1 (MQ=255)
                ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGcttctt                                                                  >  1:2215367/1‑58 (MQ=255)
                 gggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                                                    >  1:5121328/1‑69 (MQ=255)
                 gggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                     <  1:2410889/70‑1 (MQ=255)
                  ggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGAt                                                              <  1:1616517/60‑1 (MQ=255)
                    tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                 <  1:848510/71‑1 (MQ=255)
                     aaTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg                                                >  1:1316500/1‑71 (MQ=255)
                      aTTTTCGGTCCCAAATCTCAGCCAATCACATTGGCATAAACAATCTATTCTTTGATAACGCCAGATGcc                                                 >  1:5660729/1‑69 (MQ=255)
                      aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                                                         <  1:3192559/61‑1 (MQ=255)
                       ttttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGta                                              >  1:3302211/1‑70 (MQ=255)
                        tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATa                                                             <  1:1283530/55‑1 (MQ=255)
                        tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                             >  1:2993080/1‑71 (MQ=255)
                        tttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtac                                             >  1:947643/1‑69 (MQ=255)
                        tttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtac                                             >  1:4252284/1‑69 (MQ=255)
                          tCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                            >  1:803359/1‑70 (MQ=255)
                             gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                                                                 >  1:2205122/1‑46 (MQ=255)
                             gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                                                                 <  1:1582934/46‑1 (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:5393428/1‑60 (MQ=255)
                               cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                 >  1:50638/1‑60 (MQ=255)
                               cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt                                       >  1:763433/1‑70 (MQ=255)
                                ctaaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTc                                      >  1:2891717/3‑70 (MQ=255)
                                ccAAATCGCAGCCAATCACATTGGCATAAACCATCGcttct                                                                   <  1:3713674/41‑1 (MQ=255)
                                ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                                                                 <  1:5252944/43‑1 (MQ=255)
                                 taaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCaa                                    >  1:3565030/2‑71 (MQ=255)
                                  aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                   <  1:4455755/71‑1 (MQ=255)
                                  aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                   <  1:833011/71‑1 (MQ=255)
                                   aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                                                    >  1:269945/1‑51 (MQ=255)
                                   aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAtc                                  <  1:5437133/71‑2 (MQ=255)
                                   aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                   >  1:2903195/1‑70 (MQ=255)
                                     tcacAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                >  1:2392320/4‑71 (MQ=255)
                                     tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                >  1:4288879/1‑71 (MQ=255)
                                     tCGCAGCCAATCACATTGGCATAAACAATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                >  1:167266/1‑71 (MQ=255)
                                      cacAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                                                      >  1:5465245/3‑48 (MQ=255)
                                       gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                               <  1:3386210/70‑1 (MQ=255)
                                       gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                               <  1:871666/70‑1 (MQ=255)
                                       acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGAc                               <  1:3159032/69‑1 (MQ=255)
                                       acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGAc                               <  1:1526969/69‑1 (MQ=255)
                                       acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACa                              <  1:214923/70‑1 (MQ=255)
                                        cAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                             <  1:770904/71‑1 (MQ=255)
                                          gCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                           >  1:4631505/1‑71 (MQ=255)
                                          gCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcg                           <  1:2923227/71‑1 (MQ=255)
                                           ccAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgc                          <  1:3262228/71‑1 (MQ=255)
                                           ccAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                             >  1:3954659/1‑52 (MQ=255)
                                           ccAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                             >  1:416460/1‑52 (MQ=255)
                                           ccAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgc                          >  1:4701257/1‑71 (MQ=255)
                                              aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                >  1:3782538/1‑62 (MQ=255)
                                              aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                >  1:66522/1‑62 (MQ=255)
                                              aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                >  1:2777726/1‑62 (MQ=255)
                                              aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                >  1:2176540/1‑62 (MQ=255)
                                              aTCACATTGGCATAAACCATCGCTTCTTTGATAACGcc                                                        <  1:553610/38‑1 (MQ=255)
                                              aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                       <  1:2079645/71‑1 (MQ=255)
                                               tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                >  1:2594187/1‑61 (MQ=255)
                                               tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                     >  1:3407748/1‑40 (MQ=255)
                                                cacaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGta                                              >  1:2939914/1‑45 (MQ=255)
                                                cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                     >  1:5013667/1‑71 (MQ=255)
                                                cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                     >  1:464410/1‑71 (MQ=255)
                                                 acaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgca                         >  1:223661/1‑65 (MQ=255)
                                                 acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                    <  1:1037860/71‑1 (MQ=255)
                                                 acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                    >  1:2308503/1‑71 (MQ=255)
                                                 acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                    <  1:3254161/71‑1 (MQ=255)
                                                 acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                    >  1:3459158/1‑71 (MQ=255)
                                                  caTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                    <  1:1966156/70‑1 (MQ=255)
                                                   aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAc                            <  1:2983404/61‑1 (MQ=255)
                                                   aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                  <  1:4849085/71‑1 (MQ=255)
                                                    ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtactc                                           >  1:2718459/1‑44 (MQ=255)
                                                    ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                   >  1:3608310/1‑53 (MQ=255)
                                                     tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                    >  1:2678707/1‑51 (MQ=255)
                                                     tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                >  1:2541286/1‑55 (MQ=255)
                                                     tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                >  1:1689067/1‑55 (MQ=255)
                                                     tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                     >  1:3920017/1‑66 (MQ=255)
                                                      ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGt                <  1:1822723/70‑1 (MQ=255)
                                                       gCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa               <  1:3422995/70‑1 (MQ=255)
                                                          tAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                               <  1:1312794/35‑1 (MQ=255)
                                                          tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                  <  1:135212/48‑1 (MQ=255)
                                                           aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAc              >  1:1249439/1‑67 (MQ=255)
                                                           aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAc              >  1:1088319/1‑67 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                             >  1:2762815/1‑51 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                         >  1:611440/1‑55 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc          <  1:4005571/70‑1 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg         >  1:190643/1‑71 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg         >  1:4762948/1‑71 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg         >  1:4571849/1‑71 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGAGCTGCCATCAATGGACAACGCGCCAAcc                   >  1:4243777/1‑61 (MQ=255)
                                                               cATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                      <  1:4131968/39‑1 (MQ=255)
                                                                aTCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                    >  1:4003632/1‑40 (MQ=255)
                                                                   gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt  >  1:285528/1‑71 (MQ=255)
                                                                     |                                                                    
AGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATT  >  minE/2595452‑2595589

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: