Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,151 G→C 28.8% T149T (ACC→ACG arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1510GC28.8% 86.3 / 40.5 59T149T (ACC→ACGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (17/25);  new base C (7/10);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

GAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAG  >  minE/2596081‑2596219
                                                                      |                                                                    
gAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCg                                                                      <  1:2789780/71‑1 (MQ=255)
  aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                                                                     >  1:2877939/1‑70 (MQ=255)
  aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                                                                     >  1:2715622/1‑70 (MQ=255)
  aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                                                                     >  1:635676/1‑70 (MQ=255)
  aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGt                                                                    <  1:203151/71‑1 (MQ=255)
    gCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTAt                                                                  <  1:2529413/71‑1 (MQ=255)
     ccGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATc                                                                 >  1:3677093/1‑71 (MQ=255)
      cGAGGCCAAAGAAATCAGCGTAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCg                                                                <  1:3030136/71‑1 (MQ=255)
        aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                                              <  1:1731717/71‑1 (MQ=255)
          gCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCaa                                                            >  1:2084086/1‑71 (MQ=255)
           ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                                           <  1:2061612/71‑1 (MQ=255)
            cgaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTa                                                                   <  1:2746529/60‑1 (MQ=255)
             gaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAtaa                                                         <  1:1599426/70‑3 (MQ=255)
              aagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATaaa                                                        <  1:2122001/71‑1 (MQ=255)
              aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa                                                        >  1:3449102/1‑71 (MQ=255)
              aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa                                                        >  1:333578/1‑71 (MQ=255)
               aGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATct                                                                >  1:675806/1‑61 (MQ=255)
               aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCg                                                                      >  1:3522283/1‑56 (MQ=255)
               aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAt                                                       >  1:3510425/1‑71 (MQ=255)
                        cGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgccgc                                               >  1:70191/1‑70 (MQ=255)
                         ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGTTATCGGCAATGAATCCGGCCGCCa                                             >  1:3568100/1‑71 (MQ=255)
                            gaCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAtga                                          <  1:2985791/70‑3 (MQ=255)
                             aCGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgcc                                                 >  1:56520/1‑63 (MQ=255)
                             aCGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCAc                                         <  1:1297938/71‑1 (MQ=255)
                             aCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAtaa                                                         >  1:2573727/1‑53 (MQ=255)
                               gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAt                                                           <  1:2702354/51‑1 (MQ=255)
                               gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                       <  1:2419478/71‑3 (MQ=255)
                               gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAca                                        <  1:1542783/70‑2 (MQ=255)
                                 tatTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAAc                                     >  1:2097649/1‑71 (MQ=255)
                                   ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGt                                                                    <  1:3217631/38‑1 (MQ=255)
                                   ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCAc                                         <  1:487693/65‑1 (MQ=255)
                                    tCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGaa                                      <  1:3465117/67‑1 (MQ=255)
                                    tCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGc                                   >  1:3485191/1‑70 (MQ=255)
                                      aCCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCg                                                    <  1:3889235/51‑1 (MQ=255)
                                      aCCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCa                                 <  1:819963/70‑1 (MQ=255)
                                       ccAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATcc                                                     >  1:991115/1‑49 (MQ=255)
                                       ccAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgt                               <  1:307056/71‑2 (MQ=255)
                                        cAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgccgcc                                              <  1:416266/55‑1 (MQ=255)
                                        cAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg                              <  1:3783559/71‑1 (MQ=255)
                                         aGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTgca                                                 >  1:3436126/1‑50 (MQ=255)
                                         aGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgt                             >  1:2497266/1‑71 (MQ=255)
                                          ggTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgtc                            >  1:3315800/1‑70 (MQ=255)
                                           gTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                                           <  1:1555969/39‑1 (MQ=255)
                                           gTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGc                                                              <  1:955960/36‑1 (MQ=12)
                                              ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAg                           >  1:1119177/1‑68 (MQ=255)
                                              ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg                         <  1:2402450/70‑1 (MQ=255)
                                               gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGc                        <  1:2267886/70‑1 (MQ=255)
                                                agaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCgg                                                   <  1:58925/42‑1 (MQ=255)
                                                agaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCt                      <  1:1856390/71‑1 (MQ=255)
                                                aaaCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACATAGGCCAGCGTCGTGccc                      >  1:548452/3‑70 (MQ=255)
                                                 gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGt                             <  1:3057455/63‑1 (MQ=255)
                                                  aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                    <  1:2091714/71‑4 (MQ=255)
                                                    aaTAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGcc                       <  1:2481318/66‑1 (MQ=255)
                                                     aTAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                 <  1:868664/71‑1 (MQ=255)
                                                         gCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAAc             <  1:3213076/71‑1 (MQ=255)
                                                            gCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                 >  1:1278145/1‑64 (MQ=255)
                                                            gCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAAccc           >  1:2921311/1‑70 (MQ=255)
                                                             cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAAcc            >  1:1064160/1‑68 (MQ=255)
                                                             cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa         <  1:857920/71‑1 (MQ=255)
                                                             cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa         <  1:988578/71‑1 (MQ=255)
                                                              aGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa         >  1:2518166/1‑70 (MQ=255)
                                                               ggCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa         >  1:2130549/1‑69 (MQ=255)
                                                                 cTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGc     <  1:3818098/71‑1 (MQ=255)
                                                                     cGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcag  >  1:2087774/1‑70 (MQ=255)
                                                                      |                                                                    
GAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAG  >  minE/2596081‑2596219

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: