Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,782 A→T 6.7% 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,7820AT 6.7% 209.3 / 6.5 75G272G (GGT→GGAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (40/30);  new base T (2/3);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

CATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCG  >  minE/2595715‑2595851
                                                                   |                                                                     
cATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                                                                    <  1:474058/71‑1 (MQ=255)
cATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                                                                    >  1:161659/1‑71 (MQ=255)
 aTAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                                                   <  1:351432/71‑1 (MQ=255)
  tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                                                   <  1:1543030/70‑1 (MQ=255)
  tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                                                   >  1:4079911/1‑70 (MQ=255)
   aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa                                                                 >  1:1050198/1‑71 (MQ=255)
     accaccAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTAtt                                                               <  1:3981654/71‑1 (MQ=255)
      ccaccaGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                                                              <  1:3118844/71‑1 (MQ=255)
       caccaGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                                                             >  1:1221211/1‑71 (MQ=255)
          caGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                          <  1:873797/71‑1 (MQ=255)
          caGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                          >  1:3218807/1‑71 (MQ=255)
           aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa                                                                 >  1:2292813/1‑63 (MQ=255)
            gATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                                                        >  1:479754/1‑71 (MQ=255)
              tACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                                                             >  1:1446015/1‑64 (MQ=255)
              tACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAAt                                                       <  1:3560186/70‑1 (MQ=255)
              tACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATg                                                      <  1:1272004/71‑1 (MQ=255)
               aCATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGa                                                     >  1:1143159/1‑71 (MQ=255)
                cATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                          <  1:1546797/65‑1 (MQ=255)
                    ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                >  1:1150191/1‑71 (MQ=255)
                    ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                <  1:3724475/71‑1 (MQ=255)
                    ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                >  1:3715574/1‑71 (MQ=255)
                    ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgc                                                 <  1:1239249/70‑1 (MQ=255)
                     cgagcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgc                                                 >  1:1707907/1‑69 (MQ=255)
                      gagcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGt                                                                  >  1:3120323/1‑51 (MQ=255)
                      gagcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGtt                                              >  1:2354371/1‑71 (MQ=255)
                       agcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcg                                                  <  1:542549/66‑1 (MQ=255)
                        gcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTa                                                                       >  1:575450/1‑44 (MQ=255)
                        gcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTAtt                                                               >  1:4146121/1‑52 (MQ=255)
                         tgagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGttg                                                              <  1:1255131/51‑2 (MQ=255)
                         cgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                                                                    >  1:310340/1‑46 (MQ=255)
                          gagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCTTTTCGc                                          >  1:2720152/1‑71 (MQ=255)
                               gaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                                                   <  1:2428047/41‑1 (MQ=255)
                                aTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                                         <  1:3595693/66‑1 (MQ=255)
                                  gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcgac                                       >  1:3410138/1‑66 (MQ=255)
                                  gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                   >  1:2972399/1‑70 (MQ=255)
                                   aCAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                 >  1:3746971/1‑71 (MQ=255)
                                    cAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                 >  1:3127947/1‑70 (MQ=255)
                                      aTCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                                                             <  1:2015137/40‑1 (MQ=255)
                                      aTCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                              <  1:649158/71‑1 (MQ=255)
                                       tCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                              >  1:2263615/1‑70 (MQ=255)
                                        cTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                                                              >  1:424177/1‑37 (MQ=255)
                                         tGCCAGGGGGCACCGCTCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat                           >  1:1346485/1‑71 (MQ=255)
                                         tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCAt                                                       >  1:2794700/1‑43 (MQ=255)
                                         tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                              >  1:339359/1‑68 (MQ=255)
                                         tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat                           >  1:1679962/1‑71 (MQ=255)
                                         tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat                           >  1:3684755/1‑71 (MQ=255)
                                          gCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcga                                        >  1:63436/1‑57 (MQ=255)
                                           ccAGTGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                                         <  1:2647464/55‑1 (MQ=255)
                                           ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                                         <  1:349721/55‑1 (MQ=255)
                                           ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                                         <  1:1266369/55‑1 (MQ=255)
                                           ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                             <  1:1589338/67‑1 (MQ=255)
                                           ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAataata                            <  1:4052833/68‑1 (MQ=255)
                                           ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCa                         <  1:959771/71‑1 (MQ=255)
                                            cAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAg                        <  1:936051/71‑1 (MQ=255)
                                              gggggCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAgc                       <  1:70213/70‑1 (MQ=255)
                                                  gCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGcc                  <  1:3954906/71‑1 (MQ=255)
                                                  gCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGcc                  <  1:472007/71‑1 (MQ=255)
                                                     ccGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCTTTTCGcca                                        >  1:16448/1‑44 (MQ=255)
                                                     ccGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaata                               <  1:902947/55‑1 (MQ=255)
                                                     ccGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaata                               <  1:937985/55‑1 (MQ=255)
                                                      cgcgCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcga                                        >  1:1090519/1‑45 (MQ=255)
                                                      cgcgCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                             >  1:3025188/1‑56 (MQ=255)
                                                      cgcgCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                             >  1:2578243/1‑56 (MQ=255)
                                                       gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGc                                          <  1:2715114/42‑1 (MQ=255)
                                                       gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                   >  1:2985067/1‑49 (MQ=255)
                                                           aGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGc         <  1:2579417/71‑1 (MQ=255)
                                                             gACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCaa       <  1:787428/71‑1 (MQ=255)
                                                               cTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATg     >  1:2107911/1‑71 (MQ=255)
                                                                tttACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                             >  1:2183796/1‑46 (MQ=255)
                                                                tttACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCaa       <  1:2987129/68‑1 (MQ=255)
                                                                 ttACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGgc   >  1:1996529/1‑71 (MQ=255)
                                                                  tACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                 >  1:3192892/1‑40 (MQ=255)
                                                                  tACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                 >  1:3010115/1‑40 (MQ=255)
                                                                  tACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa                                 >  1:666368/1‑40 (MQ=255)
                                                                  tACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGgcg  >  1:3088745/1‑71 (MQ=255)
                                                                  tACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGgcg  >  1:1553068/1‑71 (MQ=255)
                                                                   |                                                                     
CATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCG  >  minE/2595715‑2595851

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: