Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,596,163 G→A 24.7% F145F (TTC→TTT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,596,1630GA24.7% 116.4 / 37.1 69F145F (TTC→TTTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (19/32);  new base A (9/8);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGA  >  Exported/2596094‑2596233
                                                                     |                                                                      
aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGa                                                                       <  1:2147117/71‑1 (MQ=255)
 aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa                                                                      <  1:1532620/71‑1 (MQ=255)
  aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAt                                                                     >  1:2210402/1‑71 (MQ=255)
  aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAt                                                                     >  1:343780/1‑71 (MQ=255)
   gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                                    <  1:1231880/71‑1 (MQ=255)
   gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                                    >  1:8113/1‑71 (MQ=255)
   gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                                                    <  1:1474192/71‑1 (MQ=255)
      aTCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCgg                                                                 >  1:25208/1‑71 (MQ=255)
             gagaCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAt                                                          >  1:618636/4‑71 (MQ=255)
             gagaCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAt                                                          >  1:5038161/4‑71 (MQ=255)
             gAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                                          <  1:5361754/71‑1 (MQ=255)
             gAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                                          <  1:332899/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)
                     atTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATc                                                                    >  1:2517566/1‑53 (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:3866034/1‑56 (MQ=255)
                            aGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                                                          >  1:3513260/1‑56 (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)
                              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:3078189/71‑1 (MQ=255)
                                    gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                   <  1:5668941/71‑1 (MQ=255)
                                    gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                   <  1:1011863/71‑1 (MQ=255)
                                    gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt                                   <  1:2713506/71‑1 (MQ=255)
                                     aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATcct                                                                  <  1:3834837/39‑2 (MQ=21)
                                     aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                                  <  1:1373634/71‑4 (MQ=255)
                                     aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg                                  >  1:2704228/1‑68 (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:5329713/1‑71 (MQ=255)
                                                                     gAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGa  >  1:4338142/1‑71 (MQ=255)
                                                                     |                                                                      
AAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACTGAACCCTAAAGCCAGCAGCATGGCGATGA  >  Exported/2596094‑2596233

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: