Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,595,656 A→G 6.8% D314D (GAT→GAC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,595,6560AG 6.8% 166.2 / 2.4 59D314D (GAT→GACarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (26/29);  new base G (2/2);  total (28/31)
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.61e-01

TTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCAC  >  Exported/2595588‑2595724
                                                                    |                                                                    
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:373104/3‑71 (MQ=255)
                              aagAAGCCGGTGCCGAACGTCGCTGCCCATAACCCCTTGTCTGCCAGCAGATTCAGCACGCCAGACAGgta                                      >  1:3969711/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:3653529/1‑35 (MQ=255)
                                    ccGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATc                                                                    >  1:3724830/1‑35 (MQ=255)
                                    ccGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATc                                                                    >  1:3994860/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:1999490/1‑71 (MQ=255)
                                                   gcgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTgcg                 >  1:3402723/1‑71 (MQ=255)
                                                    cgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGaaa                                          <  1:3236981/45‑1 (MQ=255)
                                                    cgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTa                              <  1:2347561/57‑1 (MQ=255)
                                                    cgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTa                              <  1:34383/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:2875794/1‑71 (MQ=255)
                                                          acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCa          >  1:1896900/1‑71 (MQ=255)
                                                          acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCa          >  1:2682705/1‑71 (MQ=255)
                                                          acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCa          >  1:4029242/1‑71 (MQ=255)
                                                                cGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGa                                        <  1:1618024/35‑1 (MQ=255)
                                                                cGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCg                                 >  1:1189166/1‑42 (MQ=255)
                                                                cGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCg                                 >  1:4684057/1‑42 (MQ=255)
                                                                cGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCg                                 >  1:3618126/1‑42 (MQ=255)
                                                                cGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAacc    <  1:5526257/71‑1 (MQ=255)
                                                                  ttATCCGCCAGCACGTTGAGTACGCCAGAAAGata                                      <  1:3502368/35‑1 (MQ=255)
                                                                  ttATCCGCCAGCACGTTGAGTACGCCAGAAAGata                                      <  1:2678556/35‑1 (MQ=255)
                                                                  ttATCCGCCAGCACGTTGAGTACGCCAGAAAGata                                      <  1:1857714/35‑1 (MQ=255)
                                                                  ttATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCgg                        <  1:127529/49‑1 (MQ=255)
                                                                  ttATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCgg                        <  1:5618776/49‑1 (MQ=255)
                                                                  ttATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCgg                        <  1:5383454/49‑1 (MQ=255)
                                                                  ttATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCgg                        <  1:2278032/49‑1 (MQ=255)
                                                                  ttATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATaaa      <  1:4828598/67‑1 (MQ=255)
                                                                   tATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAaccac  <  1:2842828/70‑1 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCa                      >  1:1257627/1‑49 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCa                      >  1:4861178/1‑49 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAgg             >  1:2411165/1‑58 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAaccac  >  1:1762316/1‑69 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAaccac  >  1:4305840/1‑69 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAaccac  >  1:2871876/1‑69 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAacca   >  1:5442086/1‑68 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAacca   >  1:2537540/1‑68 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAacca   >  1:1326877/1‑68 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAacc    >  1:2400230/1‑67 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAacc    >  1:4164990/1‑67 (MQ=255)
                                                                    aTCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAac     >  1:1798524/1‑66 (MQ=255)
                                                                    |                                                                    
TTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCAC  >  Exported/2595588‑2595724

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: