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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1810GA11.8% 199.6 / 13.4 76F139F (TTC→TTTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (38/29);  new base A (5/4);  total (43/33)
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.99e-01

CGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATAA  >  minE/2596111‑2596251
                                                                      |                                                                      
cGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACg                                                                        <  1:1992297/71‑1 (MQ=255)
cGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACg                                                                        <  1:5620776/71‑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)
       cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGcc                                                                  <  1:4629170/70‑1 (MQ=255)
       cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGcc                                                                  <  1:5543577/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)
           aGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCa                                                                 >  1:2594433/1‑67 (MQ=255)
            ggTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg                                                              >  1:4563082/1‑69 (MQ=255)
            ggTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgtc                                                            >  1:2786185/1‑70 (MQ=255)
            ggTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgtc                                                            >  1:1878442/1‑70 (MQ=255)
             gTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtg                                                              <  1:204203/68‑3 (MQ=255)
                 gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGcc                                                       <  1:2455605/71‑1 (MQ=255)
                   gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                                     <  1:3078189/71‑1 (MQ=255)
                   gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                                     <  1:2713506/71‑1 (MQ=255)
                   gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                                     <  1:5668941/71‑1 (MQ=255)
                   gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                                     <  1:1011863/71‑1 (MQ=255)
                    aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                                                    <  1:1373634/71‑4 (MQ=255)
                    aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                                                    >  1:2704228/1‑68 (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:5329713/1‑71 (MQ=255)
                                                    gAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGa                    >  1:4338142/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)
                                                           gccgccATCACGAACGCCAGCGTAGTGCCTTTAATGAACCCTAAAGCCAGCAGCATGGCGATGACAATCgg             >  1:5388546/1‑71 (MQ=255)
                                                            ccgccATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCagcagc                              >  1:4858470/1‑53 (MQ=255)
                                                              gccATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGgtgt          >  1:2643244/1‑71 (MQ=255)
                                                              gccATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGgtgt          >  1:2634198/1‑71 (MQ=255)
                                                                  tCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcag                                  >  1:4216609/1‑43 (MQ=255)
                                                                  tCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcag                                  >  1:5582416/1‑43 (MQ=255)
                                                                  tCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcag                                  >  1:4140311/1‑43 (MQ=255)
                                                                   cACGAACGCCAGCGTAGTGCCTTTACTGAACCCTaaa                                       >  1:5196959/1‑37 (MQ=255)
                                                                   cACGAACGCCAGCGTAGTGCCTTTACTGAACCCTaaa                                       >  1:2194918/1‑37 (MQ=255)
                                                                   cACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAAt                >  1:2524078/1‑60 (MQ=255)
                                                                   cACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCaaa      <  1:162405/70‑1 (MQ=255)
                                                                    aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTaaa                                       >  1:1523426/1‑36 (MQ=255)
                                                                    aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTaaa                                       >  1:1143062/1‑36 (MQ=255)
                                                                    aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCaaaa     >  1:4825592/1‑70 (MQ=255)
                                                                    aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCaaaa     >  1:4032289/1‑70 (MQ=255)
                                                                    aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCaaaa     >  1:1783282/1‑70 (MQ=255)
                                                                    aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCaaaa     >  1:3223943/1‑70 (MQ=255)
                                                                    aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAAt    <  1:2323789/71‑1 (MQ=255)
                                                                     cGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCAt                            >  1:3217842/1‑46 (MQ=255)
                                                                     cGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATa   <  1:1441707/71‑1 (MQ=255)
                                                                      gAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATaa  <  1:2795679/71‑1 (MQ=255)
                                                                      |                                                                      
CGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATAA  >  minE/2596111‑2596251

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: