Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,595,634 G→A 29.3% L322F (CTC→TTC)  arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,595,6340GA29.3% 54.4 / 25.4 41L322F (CTC→TTC) arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (14/15);  new base A (6/6);  total (20/21)
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.80e-01

GCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCC  >  Exported/2595566‑2595700
                                                                    |                                                                  
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)
                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)
                      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)
                                                ggTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGaaa                  <  1:2801304/71‑1 (MQ=255)
                                                 gTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAg                 <  1:3897867/71‑1 (MQ=255)
                                                 gTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAg                 <  1:3026881/71‑1 (MQ=255)
                                                    aagAAGCCGGTGCCGAACGTCGCTGCCCATAACCCCTTGTCTGCCAGCAGATTCAGCACGCCAGACAGgta              >  1:3969711/3‑71 (MQ=255)
                                                    aagAAGCCGGTGCCGAACGTCGCTGCCCATAACCCCTTGTCTGCCAGCAGATTCAGCACGCCAGACAGgta              >  1:373104/3‑71 (MQ=255)
                                                    aGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACg                                  <  1:4723700/51‑1 (MQ=255)
                                                        aaccGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTcc          <  1:3970417/69‑1 (MQ=255)
                                                        aTCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTcc          <  1:2700788/71‑1 (MQ=255)
                                                         tCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGaaa                  >  1:4175024/1‑62 (MQ=255)
                                                         tCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGata              >  1:4961856/1‑66 (MQ=255)
                                                          ccGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATc                                            >  1:3994860/1‑35 (MQ=255)
                                                          ccGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATc                                            >  1:3724830/1‑35 (MQ=255)
                                                          ccGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATc                                            >  1:3653529/1‑35 (MQ=255)
                                                             gTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTaa     <  1:153207/71‑1 (MQ=255)
                                                                ccGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATcc  <  1:5288169/71‑1 (MQ=255)
                                                                    |                                                                  
GCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCC  >  Exported/2595566‑2595700

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: