Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,595,598 T→C 26.5% I334V (ATT→GTT) ‡ arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,595,5980TC26.5% 86.4 / 32.9 59I334V (ATT→GTT) ‡arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base T (18/25);  new base C (5/10);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

ACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCAC  >  Exported/2595531‑2595667
                                                                   |                                                                     
aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATaga                                                                    <  1:1441744/71‑1 (MQ=255)
aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATaga                                                                    <  1:1448564/71‑1 (MQ=255)
aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                                                     <  1:5685742/70‑1 (MQ=255)
aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                                                     <  1:364533/70‑1 (MQ=255)
aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                                                     <  1:1061604/70‑1 (MQ=255)
aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                                                     <  1:4980414/70‑1 (MQ=255)
 cGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaa                                                                   <  1:464555/71‑1 (MQ=255)
 cGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaa                                                                   <  1:3890566/71‑1 (MQ=255)
  gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaag                                                                  <  1:1003764/71‑1 (MQ=255)
   ccAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga                                                                 <  1:1625058/71‑1 (MQ=255)
      gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag                                                              <  1:5283863/71‑1 (MQ=255)
           cGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGaaatg                                                         >  1:4332174/1‑69 (MQ=255)
             tgctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCg                                                       <  1:2098620/71‑1 (MQ=255)
                  acATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAAAAGAAGAGAGGAAGGCGGTGAg                                                  >  1:2566778/2‑71 (MQ=255)
                     tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGaa                                               <  1:4671686/71‑1 (MQ=255)
                      cAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAAt                                              >  1:606189/1‑71 (MQ=255)
                       aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagag                                                                >  1:2789313/1‑52 (MQ=255)
                       aaTCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGaagc                                             <  1:4929070/71‑3 (MQ=255)
                       aaTCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGaagc                                             <  1:3490046/71‑3 (MQ=255)
                        aTGGACAACCCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATc                                             >  1:3660895/1‑70 (MQ=255)
                         tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                                                     <  1:2407353/45‑1 (MQ=255)
                         tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                                                     <  1:4363712/45‑1 (MQ=255)
                         tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                                                     <  1:3637119/45‑1 (MQ=255)
                         tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                                                     <  1:3300775/45‑1 (MQ=255)
                         tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt                                                     <  1:218872/61‑1 (MQ=255)
                           gACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGt                                         <  1:3523388/71‑1 (MQ=255)
                           gACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGaa                                                           >  1:1057806/1‑53 (MQ=255)
                             cAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGc                                       >  1:3578976/1‑71 (MQ=255)
                               gcgcaCCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCg                                     <  1:3499369/70‑1 (MQ=255)
                               gcgcaCCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCg                                     <  1:4148179/70‑1 (MQ=255)
                                cgcgCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga                                                                 >  1:800238/1‑42 (MQ=255)
                                cgcgCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga                                                                 >  1:5598692/1‑42 (MQ=255)
                                cgcgCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga                                                                 >  1:2020014/1‑42 (MQ=255)
                                cgcgCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCgg                                          >  1:1755250/1‑65 (MQ=255)
                                   accAATCAGCACCGTCGGCATATTGTTCATGACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCGaac                                  <  1:4954810/69‑3 (MQ=255)
                                   accAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCGAACg                                 >  1:5332894/2‑71 (MQ=255)
                                    ccAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGt                                >  1:1209201/1‑71 (MQ=255)
                                     cAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCGAACGTc                               <  1:5586629/71‑1 (MQ=255)
                                     cAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGt                                <  1:4575829/70‑1 (MQ=255)
                                     cAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTc                               >  1:5083682/1‑71 (MQ=255)
                                      aatcAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCGAACGTcg                              >  1:1122560/4‑71 (MQ=255)
                                      aaCCAGTACCGTCGGCATATTGTTCATAATagaaga                                                                 <  1:579353/36‑1 (MQ=255)
                                         cAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCgg                           >  1:4912015/1‑71 (MQ=255)
                                             aCCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGccc                        <  1:3247073/70‑1 (MQ=255)
                                              ccGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCca                       <  1:2492122/70‑1 (MQ=255)
                                               cGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCGAACGTCGCTGCCcata                     <  1:86341/71‑3 (MQ=255)
                                                gTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCGAACGTCGCTGCCcata                     >  1:1935347/1‑68 (MQ=255)
                                                gTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAg                    <  1:1573231/71‑1 (MQ=255)
                                                 tCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCGAACGTCGCTGCCcataa                    >  1:4628111/1‑67 (MQ=255)
                                                  cGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGaa                                               <  1:1006725/42‑1 (MQ=255)
                                                   ggCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCGAACGt                                <  1:836104/56‑1 (MQ=255)
                                                   ggCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGcc                 >  1:1587922/1‑71 (MQ=255)
                                                    gCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCGAACGTCGCTGCCcat                      <  1:3142045/65‑2 (MQ=255)
                                                       tatTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGccc                        >  1:2083290/1‑60 (MQ=255)
                                                         ttGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCgg                                          >  1:3211892/1‑40 (MQ=255)
                                                         ttGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCgg                                          >  1:5231212/1‑40 (MQ=255)
                                                         ttGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCgg                                          >  1:5356225/1‑40 (MQ=255)
                                                         ttGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATc           >  1:4249838/1‑71 (MQ=255)
                                                            ttCATAATAGAAGAGGGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTa             <  1:2715517/66‑1 (MQ=255)
                                                                 aataGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCa   >  1:1830576/1‑71 (MQ=255)
                                                                  ataGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCAc  >  1:1467928/1‑71 (MQ=255)
                                                                   |                                                                     
ACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCAC  >  Exported/2595531‑2595667

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: