Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,169 G→T 22.0% 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,1690GT22.0% 156.8 / 27.1 68A143A (GCC→GCAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (23/30);  new base T (8/7);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

ATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCG  >  minE/2596100‑2596239
                                                                     |                                                                      
aTCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCgg                                                                       >  1:25208/1‑71 (MQ=255)
       gagaCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAt                                                                >  1:618636/4‑71 (MQ=255)
       gagaCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAt                                                                >  1:5038161/4‑71 (MQ=255)
       gAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                                                <  1:332899/71‑1 (MQ=255)
       gAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                                                <  1:5361754/71‑1 (MQ=255)
           cGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACg                                                            <  1:5620776/71‑1 (MQ=255)
           cGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACg                                                            <  1:1992297/71‑1 (MQ=255)
           cGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAca                                                            <  1:2088555/71‑2 (MQ=255)
            gATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgccgcc                                                                  <  1:1043670/64‑1 (MQ=255)
            gATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGa                                                           <  1:1032259/71‑1 (MQ=255)
            gATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGa                                                           <  1:1039141/71‑1 (MQ=255)
                ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                                       <  1:2205322/71‑1 (MQ=255)
                ttCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGc                                                       <  1:2890870/71‑1 (MQ=255)
                ttCACCAGCTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAt                                                                <  1:2379629/62‑1 (MQ=255)
                  cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGcc                                                      <  1:5543577/70‑1 (MQ=255)
                  cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGcc                                                      <  1:4629170/70‑1 (MQ=255)
                  cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                                                       >  1:3401340/1‑69 (MQ=255)
                  cACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGcc                                                      <  1:2777360/70‑1 (MQ=255)
                     cAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg                                                  <  1:1548326/71‑1 (MQ=255)
                     cAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtg                                                  >  1:3791224/1‑69 (MQ=255)
                      aGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                                                >  1:3513260/1‑56 (MQ=255)
                      aGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                                                >  1:3866034/1‑56 (MQ=255)
                      aGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCa                                                     >  1:2594433/1‑67 (MQ=255)
                       ggTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg                                                  >  1:4563082/1‑69 (MQ=255)
                       ggTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgtc                                                >  1:1878442/1‑70 (MQ=255)
                       ggTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgtc                                                >  1:2786185/1‑70 (MQ=255)
                        gTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                                                           <  1:4699307/59‑3 (MQ=255)
                        gTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtg                                                  <  1:204203/68‑3 (MQ=255)
                            gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCg                                                                        <  1:4551296/42‑1 (MQ=255)
                            gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGcc                                           <  1:2455605/71‑1 (MQ=255)
                              gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                         <  1:1011863/71‑1 (MQ=255)
                              gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                         <  1:3078189/71‑1 (MQ=255)
                              gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                         <  1:5668941/71‑1 (MQ=255)
                              gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                         <  1:2713506/71‑1 (MQ=255)
                               aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                                        >  1:2704228/1‑68 (MQ=255)
                               aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                                        <  1:1373634/71‑4 (MQ=255)
                                cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCttt                                        <  1:2081361/70‑1 (MQ=255)
                                cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCttt                                        <  1:2206291/70‑1 (MQ=255)
                                   tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt                                     >  1:3766970/1‑70 (MQ=255)
                                     aGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGaa                                  <  1:679941/71‑1 (MQ=255)
                                     aGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGa                                   <  1:3919204/70‑1 (MQ=255)
                                     aGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATAACGaa                                                          >  1:2443821/1‑47 (MQ=255)
                                     aGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACaaa                                                          >  1:4412467/1‑47 (MQ=255)
                                     aGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGCCCTGGCTGaa                                  >  1:1293074/1‑71 (MQ=255)
                                             gCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAg                          <  1:2946016/71‑1 (MQ=255)
                                                 gCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcag                      >  1:1390934/1‑71 (MQ=255)
                                                     tATCGGCAATTAATCCGGCCGCCATCACGAACGCCAGCGt                                                 <  1:1241997/40‑1 (MQ=255)
                                                      aTCGGCAATGAATCCGGCCGCCATCACGAACGCCa                                                     >  1:3517808/1‑35 (MQ=255)
                                                      aTCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAg                                               >  1:788957/1‑41 (MQ=255)
                                                      aTCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGaa                                  >  1:4666977/1‑54 (MQ=255)
                                                       tCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAgcagca                 <  1:925310/70‑1 (MQ=255)
                                                         ggCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGaa                                  >  1:2865199/1‑51 (MQ=255)
                                                           cAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGc             >  1:94738/1‑70 (MQ=255)
                                                           cAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGc             >  1:1820521/1‑70 (MQ=255)
                                                            aaTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGa           >  1:2228028/1‑71 (MQ=255)
                                                             aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                         <  1:299745/40‑1 (MQ=255)
                                                             aTGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGAt          <  1:5395251/71‑1 (MQ=255)
                                                              tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGcc                                           <  1:5668350/37‑1 (MQ=255)
                                                              tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                         <  1:5542945/39‑1 (MQ=255)
                                                              tGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCca                       <  1:2866961/57‑1 (MQ=255)
                                                               gAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAAcc                                >  1:2135383/1‑47 (MQ=255)
                                                               gAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGa        >  1:4338142/1‑71 (MQ=255)
                                                               gAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGa        >  1:5329713/1‑71 (MQ=255)
                                                                 aTCCGGCCGCCATCACGAACGCCAGCGTAGTGCCttt                                        <  1:5331328/37‑1 (MQ=255)
                                                                  tCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACaa     >  1:3221777/1‑71 (MQ=255)
                                                                    cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGAc       >  1:1018803/1‑67 (MQ=255)
                                                                    cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGAc       >  1:1849414/1‑67 (MQ=255)
                                                                    cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGAc       >  1:3868169/1‑67 (MQ=255)
                                                                    cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATc   >  1:2272012/1‑71 (MQ=255)
                                                                    cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATc   <  1:5071510/71‑1 (MQ=255)
                                                                     ggCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCg  <  1:2678734/71‑1 (MQ=255)
                                                                     |                                                                      
ATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCG  >  minE/2596100‑2596239

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