Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,595,796 C→A 8.4% A268S (GCG→TCG)  arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,595,7960CA 8.4% 301.6 / 8.2 94A268S (GCG→TCG) arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (49/37);  new base A (6/2);  total (55/39)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGCCG  >  Exported/2595728‑2595866
                                                                    |                                                                      
aTACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGc                                                                        <  1:5510319/69‑1 (MQ=255)
aTACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCAt                                                                      <  1:1668657/71‑1 (MQ=255)
 tACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAAt                                                                      <  1:1845789/70‑1 (MQ=255)
  aCATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGa                                                                    >  1:5107652/1‑71 (MQ=255)
   cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAc                                                                   >  1:1060550/1‑71 (MQ=255)
   cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAc                                                                   >  1:180481/1‑71 (MQ=255)
       ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                               <  1:2361138/71‑1 (MQ=255)
        cgagcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGt                                                              <  1:4160915/71‑1 (MQ=255)
          agcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGt                                                              >  1:1319482/1‑69 (MQ=255)
          agcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGt                                                              >  1:3451039/1‑69 (MQ=255)
               gaagaTGACAATCTTCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGtt                                                             >  1:5332342/1‑65 (MQ=255)
               gaagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGtt                                                             >  1:2317579/1‑65 (MQ=255)
               gaagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGtt                                                             >  1:189617/1‑65 (MQ=255)
               gaagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcga                                                       >  1:713254/1‑71 (MQ=255)
                aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgc                                                                >  1:1017289/1‑61 (MQ=255)
                   aTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcgacga                                                    >  1:4151440/1‑70 (MQ=255)
                   aTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCgacgac                                                   <  1:1089310/71‑1 (MQ=255)
                   aTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCgacgac                                                   <  1:4033741/71‑1 (MQ=255)
                    tGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGa                                                                    >  1:2307392/1‑53 (MQ=255)
                    tGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGa                                                                    >  1:3044607/1‑53 (MQ=255)
                    tGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                                  <  1:694332/71‑1 (MQ=255)
                     gACAATCTGCCAGGTGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGc                                                         >  1:2381533/1‑63 (MQ=255)
                     gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGc                                                         >  1:2823252/1‑63 (MQ=255)
                     gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGc                                                         >  1:1600690/1‑63 (MQ=255)
                     gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                                  >  1:1396365/1‑70 (MQ=255)
                      aCAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                                <  1:1027949/71‑1 (MQ=255)
                      aCAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                                <  1:4115735/71‑1 (MQ=255)
                       cAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                                >  1:2115101/1‑70 (MQ=255)
                       cAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaat                                               <  1:3595832/71‑1 (MQ=255)
                         aTCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGt                                                              <  1:3955823/54‑1 (MQ=255)
                         aTCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                                             <  1:5658511/71‑1 (MQ=255)
                          tCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCTTTTCGCCACCACaaacag                                             >  1:1109487/1‑67 (MQ=255)
                           cTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                                >  1:3463228/1‑66 (MQ=255)
                              ccATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCtttt                                                            >  1:334792/1‑51 (MQ=255)
                              ccATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCTTTTCGCCACCACaaacag                                             <  1:3102050/66‑4 (MQ=255)
                              ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                            <  1:5305533/67‑1 (MQ=255)
                              ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat                                          <  1:5287046/69‑1 (MQ=255)
                               cAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAg                                       >  1:4824907/1‑71 (MQ=255)
                                aGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGc                                                                        >  1:939516/1‑37 (MQ=255)
                                  ggggCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAgcgc                                    <  1:5058683/71‑1 (MQ=255)
                                      cACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGcc                                 <  1:3093418/70‑1 (MQ=255)
                                       aCCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCTTTTCGCCACCACaaa                                                >  1:4772379/1‑54 (MQ=255)
                                        ccGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCca                               >  1:317747/1‑70 (MQ=255)
                                         cgcgCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCcaca                             <  1:3079868/71‑1 (MQ=255)
                                         cgcgCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCca                               >  1:3932763/1‑69 (MQ=255)
                                          gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                                             >  1:3531616/1‑54 (MQ=255)
                                          gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAgcg                                     >  1:4576056/1‑62 (MQ=255)
                                          gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg                            >  1:4817201/1‑71 (MQ=255)
                                          gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg                            >  1:533624/1‑71 (MQ=255)
                                          gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg                            >  1:4191834/1‑71 (MQ=255)
                                          gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg                            >  1:2575954/1‑71 (MQ=255)
                                          gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg                            >  1:3221159/1‑71 (MQ=255)
                                           cgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGc                           >  1:3481100/1‑71 (MQ=255)
                                           cgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGc                           >  1:2619880/1‑71 (MQ=255)
                                            gcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGc                           >  1:4719319/1‑70 (MQ=255)
                                            gcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGc                           >  1:1332667/1‑70 (MQ=255)
                                            gcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCt                          <  1:649150/71‑1 (MQ=255)
                                             cAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTg                         >  1:4706815/1‑71 (MQ=255)
                                             cAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTg                         >  1:2949618/1‑71 (MQ=255)
                                             cAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTg                         >  1:4868977/1‑71 (MQ=255)
                                               ggACTTTACCCGTATTAATCGCATGACCGCGTTTAGcgac                                                      <  1:5274742/40‑1 (MQ=255)
                                               ggACTTTACCCGTATTAATCGCATGACCGCGTTTAGcgac                                                      <  1:1290140/40‑1 (MQ=255)
                                               ggACTTTACCCGTATTAATCGCATGACCGCGTTTAGcgac                                                      <  1:2011494/40‑1 (MQ=255)
                                               ggACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaat                                               <  1:484439/47‑1 (MQ=255)
                                               ggACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGc                        >  1:3915310/1‑70 (MQ=255)
                                                gACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCaa                      <  1:548972/71‑1 (MQ=255)
                                                gACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCaa                      <  1:1299457/71‑1 (MQ=255)
                                                 aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAataata                                           >  1:3168700/1‑49 (MQ=255)
                                                 aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAataata                                           >  1:5628449/1‑49 (MQ=255)
                                                    ttACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg                            <  1:2908000/61‑1 (MQ=255)
                                                     tACCCGTATTAATCGCATGACCGCGTTTAGcgacga                                                    >  1:4852765/1‑36 (MQ=255)
                                                     tACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATg                    >  1:1924365/1‑68 (MQ=255)
                                                      aCCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                                  <  1:5552354/37‑1 (MQ=255)
                                                      aCCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                                  <  1:3747462/37‑1 (MQ=255)
                                                      aCCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                                  <  1:2451467/37‑1 (MQ=255)
                                                      aCCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCt                          >  1:3693062/1‑61 (MQ=255)
                                                      aCCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCt                          >  1:2392865/1‑61 (MQ=255)
                                                        ccGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTc              <  1:4287425/71‑1 (MQ=255)
                                                           tATTAATCGCATGACCGCGTTTAGCGACGACAAataata                                           <  1:1426728/39‑1 (MQ=255)
                                                           tATTAATCGCATGACCGCGTTTAGCGACGACAAataata                                           <  1:1089890/39‑1 (MQ=255)
                                                           tATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat                                          <  1:2900300/40‑1 (MQ=255)
                                                           tATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCca                               >  1:1029097/1‑51 (MQ=255)
                                                           tATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCca                               >  1:462970/1‑51 (MQ=255)
                                                            aTTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCg          >  1:1817278/1‑71 (MQ=255)
                                                             ttAATCGCATGACCGCGTTTAGCGACGACAaataat                                            <  1:2693475/36‑1 (MQ=255)
                                                             ttAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCaa                      >  1:2240880/1‑58 (MQ=255)
                                                             ttAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCaa                      >  1:2064077/1‑58 (MQ=255)
                                                             ttAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCgg         <  1:2520740/71‑1 (MQ=255)
                                                             ttAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCgg         <  1:554643/71‑1 (MQ=255)
                                                             ttAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCg          >  1:4282271/1‑70 (MQ=255)
                                                              tAATCGCATGACCGCGTTTAGCGACGACAAataata                                           <  1:2746013/36‑1 (MQ=255)
                                                                 tCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCgg         <  1:1259079/67‑1 (MQ=255)
                                                                 tCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCgg         <  1:125589/67‑1 (MQ=255)
                                                                 tCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATg     >  1:2672228/1‑71 (MQ=255)
                                                                  cGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGaa       >  1:5333786/1‑68 (MQ=255)
                                                                   gCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGcc   >  1:5008642/1‑71 (MQ=255)
                                                                    cATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGcc   >  1:1683106/1‑70 (MQ=255)
                                                                    cATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGCCg  <  1:1921033/71‑1 (MQ=255)
                                                                    cATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGCCg  <  1:3792421/71‑1 (MQ=255)
                                                                    cATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGCCg  <  1:4518276/71‑1 (MQ=255)
                                                                    |                                                                      
ATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGCCG  >  Exported/2595728‑2595866

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: