Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,596,184 C→G 12.1% A138A (GCG→GCC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,596,1840CG12.1% 195.2 / 15.6 74A138A (GCG→GCCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (37/28);  new base G (5/4);  total (42/32)
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.85e-01

TTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATAA  >  Exported/2596116‑2596251
                                                                    |                                                                   
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:1011863/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)
                     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:3868169/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)
                                                         gccATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGgtgt          >  1:2643244/1‑71 (MQ=255)
                                                         gccATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGgtgt          >  1:2634198/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)
                                                               aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTaaa                                       >  1:1523426/1‑36 (MQ=255)
                                                               aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTaaa                                       >  1:1143062/1‑36 (MQ=255)
                                                               aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCaaaa     >  1:4032289/1‑70 (MQ=255)
                                                               aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCaaaa     >  1:3223943/1‑70 (MQ=255)
                                                               aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCaaaa     >  1:4825592/1‑70 (MQ=255)
                                                               aCGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCaaaa     >  1:1783282/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)
                                                                    |                                                                   
TTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGACAATCGGTGTCAAAATAA  >  Exported/2596116‑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: