Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,163 G→A 7.4% F145F (TTC→TTT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1630GA 7.4% 259.6 / 7.8 81F145F (TTC→TTTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (34/41);  new base A (3/3);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

AAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGA  >  minE/2596094‑2596233
                                                                     |                                                                      
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:1136039/1‑71 (MQ=255)
   gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                                    >  1:2964234/1‑71 (MQ=255)
    aaaTCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                                    <  1:2731646/70‑1 (MQ=255)
          gCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgccgc                                                             >  1:2606000/1‑71 (MQ=255)
            ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCa                                                           >  1:437794/1‑71 (MQ=255)
            ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCa                                                           >  1:2265911/1‑71 (MQ=255)
            ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCa                                                           <  1:2142609/71‑1 (MQ=255)
             gAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                                          <  1:1563665/71‑1 (MQ=255)
                  gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                                     <  1:3350762/71‑3 (MQ=255)
                  gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                                     <  1:3354604/71‑3 (MQ=255)
                      ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                                          <  1:3587156/62‑1 (MQ=255)
                       tCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                                 >  1:2080108/1‑70 (MQ=255)
                        cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                                 >  1:591246/1‑69 (MQ=255)
                        cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                                 >  1:3254166/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:2701181/71‑1 (MQ=255)
                                                cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                       <  1:3277368/71‑1 (MQ=255)
                                                cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                       <  1:92884/71‑1 (MQ=255)
                                                   gCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAg                    <  1:1033214/71‑1 (MQ=255)
                                                    cTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGc                   <  1:3423796/71‑1 (MQ=255)
                                                    cTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGc                   <  1:470088/71‑1 (MQ=255)
                                                         ggTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAg                    <  1:2002415/65‑1 (MQ=255)
                                                          gTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcagcag             <  1:1193840/71‑1 (MQ=255)
                                                          gTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcagcag             <  1:2371791/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:3601550/1‑71 (MQ=255)
                                                               ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATgg        >  1:81407/1‑71 (MQ=255)
                                                               ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATgg        >  1:1195497/1‑71 (MQ=255)
                                                               ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATgg        >  1:1178081/1‑71 (MQ=255)
                                                               ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATg         >  1:901798/1‑70 (MQ=255)
                                                                gCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGc       <  1:3502706/71‑1 (MQ=255)
                                                                gCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGc       <  1:305547/71‑1 (MQ=255)
                                                                  aaTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                   <  1:4118040/41‑1 (MQ=255)
                                                                  aaTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGa     >  1:2992834/1‑71 (MQ=255)
                                                                   aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                   <  1:1090316/40‑1 (MQ=255)
                                                                   aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                   <  1:970924/40‑1 (MQ=255)
                                                                   aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt    <  1:3874024/71‑1 (MQ=255)
                                                                   aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt    <  1:3310172/71‑1 (MQ=255)
                                                                    tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATg   <  1:381507/71‑1 (MQ=255)
                                                                    tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATg   <  1:753872/71‑1 (MQ=255)
                                                                    tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATg   <  1:2740918/71‑1 (MQ=255)
                                                                    tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATg   <  1:1123222/71‑1 (MQ=255)
                                                                     gAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTa                                 >  1:630507/1‑40 (MQ=255)
                                                                     gAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTAc                                >  1:1799294/1‑41 (MQ=255)
                                                                     gAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGa  >  1:2622781/1‑71 (MQ=255)
                                                                     |                                                                      
AAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGA  >  minE/2596094‑2596233

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: