Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,595,785 C→G 11.3% T271T (ACG→ACC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,595,7850CG11.3% 237.6 / 15.3 88T271T (ACG→ACCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (43/35);  new base G (8/2);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

TAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTC  >  Exported/2595717‑2595854
                                                                    |                                                                     
tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                                                                       <  1:299542/69‑1 (MQ=255)
   accaccAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTAtt                                                                  <  1:719264/71‑1 (MQ=255)
     caccaGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                                                                <  1:257140/71‑1 (MQ=255)
       ccaGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATc                                                              <  1:2879642/71‑1 (MQ=255)
        caGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                             <  1:3722434/71‑1 (MQ=255)
         aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                             >  1:4186020/1‑70 (MQ=255)
           aTACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                                                                >  1:2385229/1‑65 (MQ=255)
           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)
              cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGtg                                                                    >  1:5689818/1‑57 (MQ=255)
              cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGtg                                                                    >  1:4378702/1‑57 (MQ=255)
              cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAc                                                       >  1:180481/1‑71 (MQ=255)
              cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAc                                                       >  1:1060550/1‑71 (MQ=255)
                tGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                                                                 <  1:2975736/59‑1 (MQ=255)
                  ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                   <  1:2361138/71‑1 (MQ=255)
                   cgagcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGt                                                  <  1:4160915/71‑1 (MQ=255)
                     agcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGt                                                  >  1:3451039/1‑69 (MQ=255)
                     agcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGt                                                  >  1:1319482/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)
                              aTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                                                                 <  1:4941705/45‑1 (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)
                                    aTCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                                                                 <  1:5520304/39‑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:3221159/1‑71 (MQ=255)
                                                     gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg                >  1:4191834/1‑71 (MQ=255)
                                                     gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg                >  1:4817201/1‑71 (MQ=255)
                                                     gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg                >  1:2575954/1‑71 (MQ=255)
                                                     gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg                >  1:533624/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:2949618/1‑71 (MQ=255)
                                                        cAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTg             >  1:4868977/1‑71 (MQ=255)
                                                        cAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTg             >  1:4706815/1‑71 (MQ=255)
                                                          ggACTTTACCCGTATTAATCGCATGACCGCGTTTAGcgac                                          <  1:2011494/40‑1 (MQ=255)
                                                          ggACTTTACCCGTATTAATCGCATGACCGCGTTTAGcgac                                          <  1:1290140/40‑1 (MQ=255)
                                                          ggACTTTACCCGTATTAATCGCATGACCGCGTTTAGcgac                                          <  1:5274742/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:2451467/37‑1 (MQ=255)
                                                                 aCCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                      <  1:5552354/37‑1 (MQ=255)
                                                                 aCCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                      <  1:3747462/37‑1 (MQ=255)
                                                                 aCCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCt              >  1:2392865/1‑61 (MQ=255)
                                                                 aCCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCt              >  1:3693062/1‑61 (MQ=255)
                                                                   ccGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTc  <  1:4287425/71‑1 (MQ=255)
                                                                    |                                                                     
TAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTC  >  Exported/2595717‑2595854

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: