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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1690GT 6.7% 299.7 / 4.4 89A143A (GCC→GCAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (42/41);  new base T (3/3);  total (45/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 = 7.60e-01

GCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATC  >  minE/2596104‑2596238
                                                                 |                                                                     
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:3354604/71‑3 (MQ=255)
        gATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                                          <  1:3350762/71‑3 (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)
                       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:1714446/1‑71 (MQ=255)
                            cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTa                                      <  1:3885421/71‑1 (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: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:1178081/1‑71 (MQ=255)
                                                     ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATgg             >  1:1195497/1‑71 (MQ=255)
                                                     ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATg              >  1:901798/1‑70 (MQ=255)
                                                      gCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGc            <  1:305547/71‑1 (MQ=255)
                                                      gCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGc            <  1:3502706/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:2740918/71‑1 (MQ=255)
                                                          tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATg        <  1:1123222/71‑1 (MQ=255)
                                                          tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATg        <  1:753872/71‑1 (MQ=255)
                                                          tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATg        <  1:381507/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:3225410/1‑41 (MQ=255)
                                                            aaTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                                    >  1:245192/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:2541705/1‑71 (MQ=255)
                                                                cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATc  <  1:3738856/71‑1 (MQ=255)
                                                                cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATc  >  1:3355953/1‑71 (MQ=255)
                                                                 |                                                                     
GCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATC  >  minE/2596104‑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: