Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,595,653 G→C 7.8% N315K (AAC→AAG arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,595,6530GC 7.8% 147.8 / 3.1 51N315K (AAC→AAGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base C (2/2);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

TTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACC  >  Exported/2595588‑2595722
                                                                 |                                                                     
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)
                                                 tCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCgg                      <  1:2943725/66‑1 (MQ=255)
                                                   gcgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTgcg               >  1:981857/1‑71 (MQ=255)
                                                   gcgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTgcg               >  1:3402723/1‑71 (MQ=255)
                                                   gcgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTgcg               >  1:1999490/1‑71 (MQ=255)
                                                    cgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGaaa                                        <  1:3236981/45‑1 (MQ=255)
                                                    cgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTa                            <  1:34383/57‑1 (MQ=255)
                                                    cgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTa                            <  1:2347561/57‑1 (MQ=255)
                                                    cgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTgcgc              <  1:3485629/71‑1 (MQ=255)
                                                     ggCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTgcg               >  1:3567205/1‑69 (MQ=255)
                                                      gCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCAtt                  <  1:2292966/65‑1 (MQ=255)
                                                        ccACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGatat                                   <  1:1465895/46‑1 (MQ=255)
                                                        ccACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTgc                >  1:3961938/1‑65 (MQ=255)
                                                        ccACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGc          >  1:5602778/1‑71 (MQ=255)
                                                        ccACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGc          >  1:2526325/1‑71 (MQ=255)
                                                        ccACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGc          >  1:3780043/1‑71 (MQ=255)
                                                         cacaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGc          >  1:4903128/1‑70 (MQ=255)
                                                          acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTaa                           <  1:5096946/52‑1 (MQ=255)
                                                          acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCa        >  1:4126379/1‑71 (MQ=255)
                                                          acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCa        >  1:1896900/1‑71 (MQ=255)
                                                          acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCa        >  1:4029242/1‑71 (MQ=255)
                                                          acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCa        >  1:2682705/1‑71 (MQ=255)
                                                          acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCa        >  1:2875794/1‑71 (MQ=255)
                                                                cGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGa                                      <  1:1618024/35‑1 (MQ=255)
                                                                cGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCg                               >  1:4684057/1‑42 (MQ=255)
                                                                cGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCg                               >  1:3618126/1‑42 (MQ=255)
                                                                cGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCg                               >  1:1189166/1‑42 (MQ=255)
                                                                cGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAacc  <  1:5526257/71‑1 (MQ=255)
                                                                 |                                                                     
TTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACC  >  Exported/2595588‑2595722

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: