Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,178 C→G 7.0% V140V (GTG→GTC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1780CG 7.0% 261.1 / 7.5 86V140V (GTG→GTCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (40/40);  new base G (3/3);  total (43/43)
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.18e-01

GATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATC  >  minE/2596112‑2596238
                                                                  |                                                            
gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                                          <  1:3354604/71‑3 (MQ=255)
gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                                          <  1:3350762/71‑3 (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)
               ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg                                            <  1:2050154/70‑1 (MQ=255)
               ggAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTg                                            >  1:3551316/1‑70 (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:2121772/1‑70 (MQ=255)
                       tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                                    <  1:2807375/70‑1 (MQ=255)
                       tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTg                                   <  1:3601908/71‑1 (MQ=255)
                       tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTg                                   <  1:334970/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:92884/71‑1 (MQ=255)
                              cAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTa                            <  1:3277368/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)
                                              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:970924/40‑1 (MQ=255)
                                                 aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                        <  1:1090316/40‑1 (MQ=255)
                                                 aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt         <  1:3310172/71‑1 (MQ=255)
                                                 aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt         <  1:3874024/71‑1 (MQ=255)
                                                  tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATg        <  1:753872/71‑1 (MQ=255)
                                                  tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATg        <  1:1123222/71‑1 (MQ=255)
                                                  tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATg        <  1:381507/71‑1 (MQ=255)
                                                  tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATg        <  1:2740918/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)
                                                    aaTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTa                                      >  1:891372/1‑39 (MQ=255)
                                                    aaTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                                    >  1:245192/1‑41 (MQ=255)
                                                    aaTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                                    >  1:3225410/1‑41 (MQ=255)
                                                    aaTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGa          >  1:177270/1‑67 (MQ=255)
                                                    aaTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt         <  1:2752507/68‑1 (MQ=255)
                                                    aaTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGa       >  1:1024029/1‑70 (MQ=255)
                                                    aaTCCGGCCGCCATCACGAACGCCAGAGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt         <  1:3357043/68‑1 (MQ=255)
                                                     aTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt         <  1:44709/67‑1 (MQ=255)
                                                     aTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACa     <  1:2354083/71‑1 (MQ=255)
                                                       ccGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAAt   >  1:3158369/1‑71 (MQ=255)
                                                        cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATc  >  1:3355953/1‑71 (MQ=255)
                                                        cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATc  >  1:2541705/1‑71 (MQ=255)
                                                        cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATc  <  1:3738856/71‑1 (MQ=255)
                                                          gccgccATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGa          >  1:2381085/1‑61 (MQ=255)
                                                             gccATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt         <  1:1760618/59‑1 (MQ=255)
                                                                 tCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCagcagc                 <  1:3386815/47‑1 (MQ=255)
                                                                  cACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcag                     >  1:2107312/1‑42 (MQ=255)
                                                                  |                                                            
GATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATC  >  minE/2596112‑2596238

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: