Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,596,172 C→A 22.0% A142A (GCG→GCT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,596,1720CA22.0% 155.6 / 26.0 68A142A (GCG→GCTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (23/30);  new base A (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 = 9.51e-01

GAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGG  >  Exported/2596107‑2596240
                                                                 |                                                                    
gagaCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAt                                                                 >  1:5038161/4‑71 (MQ=255)
gagaCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAt                                                                 >  1:618636/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: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)
               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)
                     gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGcc                                            <  1:2455605/71‑1 (MQ=255)
                       gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                          <  1:5668941/71‑1 (MQ=255)
                       gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                          <  1:2713506/71‑1 (MQ=255)
                       gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                          <  1:1011863/71‑1 (MQ=255)
                       gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                          <  1:3078189/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:1820521/1‑70 (MQ=255)
                                                    cAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGc              >  1:94738/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:3868169/1‑67 (MQ=255)
                                                             cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGAc        >  1:1018803/1‑67 (MQ=255)
                                                             cGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGAc        >  1:1849414/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)
                                                                 |                                                                    
GAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGG  >  Exported/2596107‑2596240

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: