Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,764 G→C 18.0% A278A (GCC→GCG arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,7640GC18.0% 163.2 / 26.4 72A278A (GCC→GCGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base C (6/7);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGC  >  minE/2595703‑2595830
                                                             |                                                                  
cATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCAccag                                                             <  1:677764/69‑3 (MQ=255)
  ttGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCgcgc                                                           <  1:70289/69‑1 (MQ=255)
   tGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAg                                                         >  1:3425185/1‑70 (MQ=255)
       cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                                    <  1:325150/71‑1 (MQ=255)
       cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                                    <  1:3630922/71‑1 (MQ=255)
       cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                                    <  1:2898541/71‑1 (MQ=255)
       cAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                                                    <  1:2720411/71‑1 (MQ=255)
       cAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                                                    <  1:3579096/71‑1 (MQ=255)
        aGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACttt                                                   <  1:2918027/71‑1 (MQ=255)
        aGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                                                    <  1:3526430/70‑1 (MQ=255)
         ggCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCgcgc                                                           <  1:3944700/62‑1 (MQ=255)
          gCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGa                                                       >  1:37375/1‑65 (MQ=255)
            cATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                                               >  1:161659/1‑71 (MQ=255)
            cATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                                               <  1:474058/71‑1 (MQ=255)
             aTAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                              <  1:351432/71‑1 (MQ=255)
              tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggggg                                                                   <  1:1820424/49‑1 (MQ=255)
              tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAgggg                                                                    <  1:299629/48‑1 (MQ=255)
              tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                              <  1:1543030/70‑1 (MQ=255)
              tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                              >  1:4079911/1‑70 (MQ=255)
               aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa                                            >  1:1050198/1‑71 (MQ=255)
                aaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCa                                                                 >  1:2956606/1‑49 (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:3218807/1‑71 (MQ=255)
                      caGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                     <  1:873797/71‑1 (MQ=255)
                       aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa                                            >  1:2292813/1‑63 (MQ=255)
                        gATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                                   >  1:479754/1‑71 (MQ=255)
                         aTACATGCCGAGCGAGAAGATGACAATCTGCCAggggg                                                                   <  1:19379/38‑1 (MQ=255)
                         aTACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCa                                                                 <  1:1026382/40‑1 (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)
                            cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAg                                                         >  1:1187780/1‑45 (MQ=255)
                            cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAg                                                         >  1:2411812/1‑45 (MQ=255)
                            cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                                    >  1:3468886/1‑50 (MQ=255)
                            cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                                    >  1:2056368/1‑50 (MQ=255)
                                ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                           >  1:3715574/1‑71 (MQ=255)
                                ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                           >  1:1150191/1‑71 (MQ=255)
                                ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                           <  1:3724475/71‑1 (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:3684755/1‑71 (MQ=255)
                                                     tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat      >  1:1679962/1‑71 (MQ=255)
                                                      gCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcga                   >  1:63436/1‑57 (MQ=255)
                                                       ccAGTGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                    <  1:2647464/55‑1 (MQ=255)
                                                       ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                    <  1:1266369/55‑1 (MQ=255)
                                                       ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                    <  1:349721/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)
                                                             |                                                                  
CATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGC  >  minE/2595703‑2595830

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: