Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,157 G→T 10.1% A147A (GCC→GCA arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1570GT10.1% 224.9 / 7.4 69A147A (GCC→GCAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (26/36);  new base T (4/3);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

AGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGG  >  minE/2596089‑2596227
                                                                    |                                                                      
aGGCCAATGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                                                      <  1:3031067/71‑1 (MQ=255)
aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                                                      <  1:660342/71‑1 (MQ=255)
aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                                                      <  1:1867257/71‑1 (MQ=255)
   ccAAAGAAATCAGCGGAAACGATATTCTCCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                                                   >  1:2133069/1‑71 (MQ=255)
   ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                                                   <  1:3334549/71‑1 (MQ=255)
   ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                                                   <  1:3797684/71‑1 (MQ=255)
    cAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCgg                                                                       >  1:3883718/1‑66 (MQ=255)
     aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGa                                                                 <  1:2125149/71‑1 (MQ=255)
     aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGa                                                                 <  1:3242990/71‑1 (MQ=255)
      aagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATaaa                                                                >  1:750249/1‑71 (MQ=255)
      aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa                                                                <  1:3136089/71‑1 (MQ=255)
       aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAt                                                               <  1:2531797/71‑1 (MQ=255)
        gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                              >  1:2964234/1‑71 (MQ=255)
        gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                              >  1:1136039/1‑71 (MQ=255)
         aaaTCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                              <  1:2731646/70‑1 (MQ=255)
               gCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgccgc                                                       >  1:2606000/1‑71 (MQ=255)
                 ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCa                                                     <  1:2142609/71‑1 (MQ=255)
                 ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCa                                                     >  1:2265911/1‑71 (MQ=255)
                 ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCa                                                     >  1:437794/1‑71 (MQ=255)
                  gAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                                    <  1:1563665/71‑1 (MQ=255)
                       gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                               <  1:3350762/71‑3 (MQ=255)
                       gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                               <  1:3354604/71‑3 (MQ=255)
                        atatTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCg                                                                        <  1:2202926/45‑1 (MQ=255)
                           ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                                                      <  1:3358153/44‑1 (MQ=255)
                           ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                                    <  1:3587156/62‑1 (MQ=255)
                            tCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                           >  1:2080108/1‑70 (MQ=255)
                             cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                           >  1:3254166/1‑69 (MQ=255)
                             cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                           >  1:591246/1‑69 (MQ=255)
                             cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                           >  1:806710/1‑69 (MQ=255)
                              aCCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAg                                        >  1:338330/1‑71 (MQ=255)
                                cAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg                                      <  1:1680970/71‑1 (MQ=255)
                                cAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGACAGc                                       >  1:256378/1‑70 (MQ=255)
                                cAGGTTGGAGACAATAAGCCGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGc                                       >  1:3229683/1‑70 (MQ=255)
                                 aGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGt                                     <  1:411312/71‑1 (MQ=255)
                                 aGGTTAGAAACAATGCGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGc                                           >  1:958970/1‑65 (MQ=255)
                                 aGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGc                                           >  1:1915426/1‑65 (MQ=255)
                                 aGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgt                                     <  1:128093/71‑1 (MQ=255)
                                    ttGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg                                      <  1:3983386/67‑1 (MQ=255)
                                    ttGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGt                                  <  1:2078208/71‑1 (MQ=255)
                                     tGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg                                 >  1:3775205/1‑71 (MQ=255)
                                      ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                              >  1:32683/1‑41 (MQ=255)
                                      ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg                                 >  1:3551316/1‑70 (MQ=255)
                                      ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg                                 <  1:2050154/70‑1 (MQ=255)
                                       gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGcc                               <  1:2140427/71‑1 (MQ=255)
                                         gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                             <  1:309246/71‑1 (MQ=255)
                                          aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                            >  1:3790992/1‑68 (MQ=255)
                                           cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTa                           <  1:3885421/71‑1 (MQ=255)
                                           cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTa                           >  1:1714446/1‑71 (MQ=255)
                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                         <  1:2807375/70‑1 (MQ=255)
                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                         >  1:2121772/1‑70 (MQ=255)
                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTg                        <  1:334970/71‑1 (MQ=255)
                                              tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTg                        <  1:3601908/71‑1 (MQ=255)
                                               aaGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCATCGTAGTGCCTTTACTg                        >  1:2710512/1‑70 (MQ=255)
                                                aGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTa                                    >  1:1080377/1‑57 (MQ=255)
                                                aGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGaa                      >  1:2727881/1‑71 (MQ=255)
                                                 gCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAAc                     >  1:2215050/1‑71 (MQ=255)
                                                 gCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAAc                     >  1:1842847/1‑71 (MQ=255)
                                                     cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                 <  1:3277368/71‑1 (MQ=255)
                                                     cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                 <  1:92884/71‑1 (MQ=255)
                                                     cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                 <  1:2701181/71‑1 (MQ=255)
                                                        gCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAg              <  1:1033214/71‑1 (MQ=255)
                                                         cTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGc             <  1:470088/71‑1 (MQ=255)
                                                         cTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGc             <  1:3423796/71‑1 (MQ=255)
                                                              ggTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAg              <  1:2002415/65‑1 (MQ=255)
                                                               gTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcagcag       <  1:2371791/71‑1 (MQ=255)
                                                               gTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcagcag       <  1:1193840/71‑1 (MQ=255)
                                                                tATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTg                        >  1:2696969/1‑53 (MQ=255)
                                                                tATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCagcagc      <  1:3416614/71‑1 (MQ=255)
                                                                 aTCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTAc                          >  1:821160/1‑50 (MQ=255)
                                                                 aTCGGCAATGAATCCGACCGCCATCACGAACGCCa                                         >  1:413182/1‑35 (MQ=255)
                                                                  tCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg                                      <  1:2963690/37‑1 (MQ=255)
                                                                    ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGt                                     >  1:2521008/1‑36 (MQ=255)
                                                                    ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAg                                   >  1:1383372/1‑38 (MQ=255)
                                                                    ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCGTTTa                           >  1:2888205/1‑46 (MQ=255)
                                                                    ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATgg  >  1:1195497/1‑71 (MQ=255)
                                                                    ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATgg  >  1:81407/1‑71 (MQ=255)
                                                                    ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATgg  >  1:1178081/1‑71 (MQ=255)
                                                                    ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATgg  >  1:3601550/1‑71 (MQ=255)
                                                                    ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATg   >  1:901798/1‑70 (MQ=255)
                                                                    |                                                                      
AGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGG  >  minE/2596089‑2596227

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: