Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,178 C→G 17.7% 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,1780CG17.7% 146.2 / 23.2 68V140V (GTG→GTCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base G (6/6);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAA  >  minE/2596111‑2596247
                                                                   |                                                                     
cGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACg                                                                    <  1:5620776/71‑1 (MQ=255)
cGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACg                                                                    <  1:1992297/71‑1 (MQ=255)
cGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAca                                                                    <  1:2088555/71‑2 (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: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)
           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)
                 gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGcc                                                   <  1:2455605/71‑1 (MQ=255)
                   gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                                 <  1:1011863/71‑1 (MQ=255)
                   gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                                 <  1:2713506/71‑1 (MQ=255)
                   gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                                 <  1:3078189/71‑1 (MQ=255)
                   gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                                 <  1:5668941/71‑1 (MQ=255)
                    aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                                                >  1:2704228/1‑68 (MQ=255)
                    aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                                                <  1:1373634/71‑4 (MQ=255)
                     cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCttt                                                <  1:2206291/70‑1 (MQ=255)
                     cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCttt                                                <  1:2081361/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:5071510/71‑1 (MQ=255)
                                                         cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATc           >  1:2272012/1‑71 (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:2634198/1‑71 (MQ=255)
                                                              gccATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGgtgt      >  1:2643244/1‑71 (MQ=255)
                                                                  tCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcag                              >  1:5582416/1‑43 (MQ=255)
                                                                  tCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCcag                              >  1:4216609/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)
                                                                   |                                                                     
CGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAA  >  minE/2596111‑2596247

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: