Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,529 T→C 27.5% I357V (ATC→GTC)  arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5290TC27.5% 280.8 / 130.4 189I357V (ATC→GTC) arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base T (80/57);  new base C (37/15);  total (117/72)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

Only 100 of 195 total aligned reads displayed.

TACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAA  >  minE/2595460‑2595597
                                                                     |                                                                    
tACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGAt                                                                      >  1:2367820/1‑70 (MQ=255)
     atcggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                                                                 <  1:97747/67‑1 (MQ=255)
      ttGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                                                                 >  1:3748577/1‑69 (MQ=255)
      ttGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                                                                 >  1:3177689/1‑69 (MQ=255)
        ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                                                                    >  1:1869377/1‑64 (MQ=255)
         gggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                                                            >  1:3301054/1‑69 (MQ=255)
            tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                         >  1:4639318/1‑71 (MQ=255)
            tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                         <  1:848510/71‑1 (MQ=255)
             aaTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg                                                        >  1:1316500/1‑71 (MQ=255)
              aTTTTCGGTCCCAAATCTCAGCCAATCACATTGGCATAAACAATCTATTCTTTGATAACGCCAGATGcc                                                         >  1:5660729/1‑69 (MQ=255)
               ttttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                       <  1:2214288/70‑1 (MQ=255)
               ttttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                       <  1:3674096/70‑1 (MQ=255)
                tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGACAGATGCCGtg                                                      <  1:1286167/70‑1 (MQ=255)
                tttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtac                                                     >  1:4252284/1‑69 (MQ=255)
                 ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                    >  1:1945590/1‑71 (MQ=255)
                  tCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                    >  1:803359/1‑70 (MQ=255)
                      tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt                                               >  1:5619693/1‑71 (MQ=255)
                       cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                         >  1:5393428/1‑60 (MQ=255)
                       cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                         >  1:50638/1‑60 (MQ=255)
                       cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                         >  1:3994363/1‑60 (MQ=255)
                       cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                         >  1:3866857/1‑60 (MQ=255)
                       cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt                                               >  1:763433/1‑70 (MQ=255)
                        ctaaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTc                                              >  1:2891717/3‑70 (MQ=255)
                        ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                              <  1:3949922/70‑1 (MQ=255)
                        ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                              <  1:1722350/70‑1 (MQ=255)
                        ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                              <  1:5299032/70‑1 (MQ=255)
                         taaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCaa                                            >  1:3565030/2‑71 (MQ=255)
                          aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                     >  1:1087065/1‑61 (MQ=255)
                          aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                           >  1:2927406/1‑71 (MQ=255)
                          aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                           <  1:833011/71‑1 (MQ=255)
                          aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                           <  1:4455755/71‑1 (MQ=255)
                           aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                          <  1:334155/71‑1 (MQ=255)
                           aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                                                            >  1:269945/1‑51 (MQ=255)
                           aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAtc                                          <  1:5437133/71‑2 (MQ=255)
                           aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAtc                                          <  1:911476/71‑2 (MQ=255)
                           aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                           >  1:2903195/1‑70 (MQ=255)
                             tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                        >  1:4288879/1‑71 (MQ=255)
                              cGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGAt                                                                      <  1:1812747/40‑1 (MQ=255)
                               gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                       <  1:3386210/70‑1 (MQ=255)
                               gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                       <  1:871666/70‑1 (MQ=255)
                               gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                       <  1:4587681/70‑1 (MQ=255)
                               gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                       <  1:3255432/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)
                                cAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                                     <  1:770904/71‑1 (MQ=255)
                                 aGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                       <  1:2764875/52‑1 (MQ=255)
                                  gCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                   <  1:2046393/71‑1 (MQ=255)
                                  gCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcg                                   <  1:2923227/71‑1 (MQ=255)
                                   ccAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                     >  1:1266942/1‑52 (MQ=255)
                                   ccAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                     >  1:416460/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:66522/1‑62 (MQ=255)
                                      aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                               <  1:5071896/71‑1 (MQ=255)
                                      aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                               <  1:2079645/71‑1 (MQ=255)
                                       tCACATTGGCATAAATCATTGCCTCTTTGACGCCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCAcc                               >  1:2467144/1‑70 (MQ=255)
                                       tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                        >  1:4010419/1‑61 (MQ=255)
                                       tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                        >  1:2594187/1‑61 (MQ=255)
                                       tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                             >  1:3407748/1‑40 (MQ=255)
                                       tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                             >  1:2980187/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)
                                         acaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgca                                 >  1:3832121/1‑65 (MQ=255)
                                         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                             <  1:466127/70‑1 (MQ=255)
                                         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                            <  1:1037860/71‑1 (MQ=255)
                                         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                            <  1:5472608/71‑1 (MQ=255)
                                         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                            >  1:3459158/1‑71 (MQ=255)
                                          caTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAcc                           <  1:1634087/71‑1 (MQ=255)
                                           aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                                        >  1:168534/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)
                                            ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                           >  1:2988091/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)
                                              ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                            >  1:3308327/1‑50 (MQ=255)
                                              ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                            >  1:3142306/1‑50 (MQ=255)
                                               gCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                       <  1:2396974/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)
                                                  tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                   >  1:1124572/1‑71 (MQ=255)
                                                   aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                  >  1:2560031/1‑71 (MQ=255)
                                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                 >  1:611440/1‑55 (MQ=255)
                                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                 >  1:3736814/1‑55 (MQ=255)
                                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                 >  1:4762948/1‑71 (MQ=255)
                                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                 >  1:1911825/1‑71 (MQ=255)
                                                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                 >  1:4005558/1‑71 (MQ=255)
                                                    aaCCATAGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                 >  1:3215243/1‑55 (MQ=255)
                                                       cATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                              <  1:4131968/39‑1 (MQ=255)
                                                           gcctcttTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCatat           >  1:300896/4‑70 (MQ=255)
                                                           gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt          >  1:285528/1‑71 (MQ=255)
                                                            cctcttTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTg         <  1:4234961/69‑1 (MQ=255)
                                                                tttGATAACGCCAGATGCCGTGCTGCCATCAATGGACa                                      >  1:2753541/1‑38 (MQ=255)
                                                                 ttGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                >  1:1388473/1‑43 (MQ=255)
                                                                     |                                                                    
TACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAA  >  minE/2595460‑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: