Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,782 A→T 13.8% G272G (GGT→GGA arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,7820AT13.8% 223.3 / 25.9 87G272G (GGT→GGAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (42/33);  new base T (9/3);  total (51/36)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATG  >  minE/2595713‑2595848
                                                                     |                                                                  
gCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAc                                                                   >  1:3936944/1‑71 (MQ=255)
   aTAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTccgg                                                                <  1:906572/71‑3 (MQ=255)
    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:4378702/1‑57 (MQ=255)
                  cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGtg                                                              >  1:5689818/1‑57 (MQ=255)
                  cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAc                                                 >  1:1060550/1‑71 (MQ=255)
                  cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAc                                                 >  1:180481/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)
                             agaaaaTAACGATCTGCCATGGCGCACCGCGCAGGACTTTTccg                                                                 >  1:1081720/1‑43 (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:4033741/71‑1 (MQ=255)
                                  aTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCgacgac                                 <  1:1089310/71‑1 (MQ=255)
                                   tGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGa                                                  >  1:3044607/1‑53 (MQ=255)
                                   tGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGa                                                  >  1:2307392/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:4115735/71‑1 (MQ=255)
                                     aCAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                              <  1:1027949/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:533624/1‑71 (MQ=255)
                                                         gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg          >  1:4817201/1‑71 (MQ=255)
                                                         gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg          >  1:4191834/1‑71 (MQ=255)
                                                         gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAg          >  1:2575954/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:4868977/1‑71 (MQ=255)
                                                            cAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTg       >  1:4706815/1‑71 (MQ=255)
                                                            cAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTg       >  1:2949618/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:5552354/37‑1 (MQ=255)
                                                                     aCCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                <  1:2451467/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)
                                                                     |                                                                  
GCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATG  >  minE/2595713‑2595848

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: