Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,404 G→A 29.4% S398S (AGC→AGT) ‡ arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,4040GA29.4% 196.8 / 97.1 136S398S (AGC→AGT) ‡arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (51/45);  new base A (18/22);  total (69/67)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

CCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAA  >  minE/2595334‑2595474
                                                                      |                                                                      
ccAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAg                                                                        >  1:4473558/1‑71 (MQ=255)
  aGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCt                                                                      <  1:130726/71‑1 (MQ=255)
   gcgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCt                                                                      >  1:2839713/1‑70 (MQ=255)
    cgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                                                    <  1:1243161/71‑1 (MQ=255)
    cgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                                                    <  1:29933/71‑1 (MQ=255)
    cgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                                                    <  1:513441/71‑1 (MQ=255)
       ggccagaGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAttg                                                                 >  1:3870623/2‑69 (MQ=255)
          cagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa                                                              <  1:214073/71‑1 (MQ=255)
           agcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa                                                              >  1:5646237/1‑70 (MQ=255)
           agcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa                                                              <  1:1674147/70‑1 (MQ=255)
             cGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCatat                                                           <  1:614030/71‑1 (MQ=255)
             cGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCatat                                                           <  1:4148128/71‑1 (MQ=255)
              gTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCatat                                                           <  1:2985680/70‑1 (MQ=255)
               tgacaAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttt                                                         <  1:1203956/69‑2 (MQ=255)
               tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                                                    <  1:4270512/60‑1 (MQ=255)
               tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc                                                         <  1:1791542/71‑1 (MQ=255)
                cacaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTc                                                               >  1:5575767/1‑64 (MQ=255)
                 acaAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttttt                                                       <  1:3852905/71‑1 (MQ=255)
                 acaAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAtttt                                                        <  1:2152288/70‑3 (MQ=255)
                 acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAg                                                                        >  1:2919375/1‑54 (MQ=255)
                 acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                                                       >  1:4352398/1‑71 (MQ=255)
                 acaAACAGCAAAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                                                       >  1:2334813/1‑71 (MQ=255)
                  caAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAtttttc                                                      <  1:3068514/71‑1 (MQ=255)
                   aaaCAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATttttct                                                     >  1:1942030/1‑71 (MQ=255)
                   aaaCAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATttttct                                                     >  1:760969/1‑71 (MQ=255)
                    aaCAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTg                                                    >  1:1783083/1‑71 (MQ=255)
                     aCAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGc                                                   >  1:4643370/1‑71 (MQ=255)
                      cagcagAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGc                                                   >  1:5555800/1‑70 (MQ=255)
                            aggcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGt                                                                <  1:2549641/51‑1 (MQ=255)
                             ggcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTAc                                           <  1:1130225/71‑1 (MQ=255)
                             ggcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTAc                                           <  1:4430395/71‑1 (MQ=255)
                               caggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGt                                         <  1:1486229/71‑1 (MQ=255)
                               caggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGt                                         <  1:4120927/71‑1 (MQ=255)
                                aggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCata                                                            >  1:3604778/1‑51 (MQ=255)
                                acagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg                                        <  1:5285263/68‑1 (MQ=255)
                                acagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg                                        <  1:555069/68‑1 (MQ=255)
                                acagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg                                        >  1:1389504/4‑71 (MQ=255)
                                 gggTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGc                                       <  1:1524422/71‑1 (MQ=255)
                                 cagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCa                                                              <  1:3979752/46‑1 (MQ=16)
                                 cagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg                                        <  1:2483560/68‑1 (MQ=255)
                                 cagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGc                                       <  1:725182/69‑1 (MQ=255)
                                  ggTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc                                                         <  1:3415756/52‑1 (MQ=255)
                                    tCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGaa                                                <  1:5615276/59‑1 (MQ=255)
                                    tCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGc                                    >  1:4912920/1‑71 (MQ=255)
                                    tCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttt                                                         <  1:4769993/50‑2 (MQ=255)
                                    tCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAg                                              >  1:2554648/1‑61 (MQ=255)
                                    tCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAg                                              >  1:774505/1‑61 (MQ=255)
                                     cATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAtt                                                          >  1:5272059/1‑48 (MQ=19)
                                     cATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGc                                    >  1:3535977/1‑70 (MQ=255)
                                      aTGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTc                                                               >  1:431553/1‑42 (MQ=255)
                                      aTGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTc                                                               >  1:4134779/1‑42 (MQ=255)
                                      aTGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa                                  <  1:5321996/71‑1 (MQ=255)
                                       tGATAATCCCTGTACGGAAGTAATATCCCCAGCTg                                                                     <  1:3026791/35‑1 (MQ=255)
                                         aTAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                                                       >  1:4746272/1‑47 (MQ=255)
                                         aTAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCagag                               <  1:5564881/71‑1 (MQ=255)
                                          tAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGt                                                                >  1:700295/1‑37 (MQ=255)
                                          tAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAg                                              >  1:470386/1‑55 (MQ=255)
                                          tAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAg                                              >  1:4243826/1‑55 (MQ=255)
                                          tAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAg                                              >  1:5424356/1‑55 (MQ=255)
                                           aaTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTg                                                    >  1:2184508/1‑48 (MQ=19)
                                             tCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagc                           >  1:4444805/1‑71 (MQ=255)
                                              cccTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta                             >  1:4924942/1‑68 (MQ=255)
                                              cccTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta                             >  1:516239/1‑68 (MQ=255)
                                              cccTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta                             >  1:5274802/1‑68 (MQ=255)
                                              cccTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta                             >  1:2801405/1‑68 (MQ=255)
                                               ccTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcg                          >  1:4551306/1‑70 (MQ=255)
                                                cggTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAgcagagtt                        <  1:5530447/69‑2 (MQ=255)
                                                 tGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgtag                       >  1:3436204/1‑71 (MQ=255)
                                                  gTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgtagc                      >  1:4700516/1‑71 (MQ=255)
                                                  gTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCAcat                                         <  1:2890320/52‑3 (MQ=255)
                                                   tACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTa                                            >  1:4353467/1‑48 (MQ=255)
                                                   tACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTa                                            >  1:1377180/1‑48 (MQ=255)
                                                   tACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa                                  <  1:1988895/58‑1 (MQ=255)
                                                   tACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgtag                       <  1:846704/69‑1 (MQ=255)
                                                    aCGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAgcagagttgc                      >  1:4250600/1‑69 (MQ=255)
                                                     cGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAg                   >  1:1359458/1‑71 (MQ=255)
                                                      ggAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAgg                  >  1:2068162/1‑71 (MQ=255)
                                                        aaGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGc                 >  1:4607899/1‑70 (MQ=255)
                                                        aaGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGc                 >  1:1061447/1‑70 (MQ=255)
                                                        aaGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCac                                 <  1:4735523/54‑2 (MQ=255)
                                                         aGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACtg               >  1:2813062/1‑70 (MQ=255)
                                                         aGTAATATCCCCAAGTGATTGTCATATTTTTCAGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACt                <  1:878483/70‑1 (MQ=255)
                                                          gTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg                                        <  1:457075/45‑1 (MQ=255)
                                                          gTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACtgc              <  1:1763288/71‑3 (MQ=255)
                                                           tAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACtgc              >  1:3296041/1‑68 (MQ=255)
                                                            aaTATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCa                                      >  1:3256388/1‑45 (MQ=255)
                                                             atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcg                          >  1:2705721/1‑56 (MQ=255)
                                                             atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCaa           >  1:2033974/1‑71 (MQ=255)
                                                             atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCa            <  1:5022952/70‑1 (MQ=255)
                                                              tatCCCCAGCTGATAGTCATATTCTTCTGCGaaat                                              <  1:367818/35‑2 (MQ=255)
                                                              tatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCATAGTAGCGTAGCCAGGCTACCaa           <  1:4347204/70‑1 (MQ=255)
                                                              tatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAAt          <  1:2615888/71‑1 (MQ=255)
                                                               atCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAAtt         >  1:2975719/1‑71 (MQ=255)
                                                                tCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTg        >  1:3338200/1‑71 (MQ=255)
                                                                  cccAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTggg      >  1:4068552/1‑71 (MQ=255)
                                                                  cccAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTggg      >  1:1531737/1‑71 (MQ=255)
                                                                    cAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCaga                                <  1:1632740/43‑1 (MQ=255)
                                                                      gCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagc                           >  1:1134345/1‑46 (MQ=255)
                                                                      gCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTaa  <  1:1272417/71‑1 (MQ=255)
                                                                      gCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTaa  <  1:2839072/71‑1 (MQ=255)
                                                                      |                                                                      
CCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAA  >  minE/2595334‑2595474

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: