Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,595,608 G→A 32.5% F330F (TTC→TTT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,595,6080GA32.5% 46.3 / 31.6 43F330F (TTC→TTTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (18/11);  new base A (4/10);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 5.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTA  >  Exported/2595542‑2595675
                                                                  |                                                                   
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)
            aaTCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGaagc                                                     <  1:3490046/71‑3 (MQ=255)
            aaTCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGaagc                                                     <  1:4929070/71‑3 (MQ=255)
             aTGGACAACCCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATc                                                     >  1:3660895/1‑70 (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:4148179/70‑1 (MQ=255)
                    gcgcaCCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCg                                             <  1:3499369/70‑1 (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)
                              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)
                                                           aagaGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGtt       <  1:4530913/70‑1 (MQ=255)
                                                            tgagagaAATGCCGTCAAGAAGCCGGTGCCGAACGTCGCTGCCCATAACCCCTTGTCTGCCAGCa           <  1:580369/64‑1 (MQ=255)
                                                            tgagagaAATGCCGTCAAGAAGCCGGTGCCGAACGTCGCTGCCCATAACCCCTTGTCTGCCAGCa           <  1:4123460/64‑1 (MQ=255)
                                                             gagagGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAg    <  1:1046953/71‑1 (MQ=255)
                                                               gagGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTa  >  1:993494/1‑71 (MQ=255)
                                                                  |                                                                   
CGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTA  >  Exported/2595542‑2595675

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: