Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,650 G→C 8.5% G316G (GGC→GGG arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,6500GC 8.5% 131.9 / 3.2 47G316G (GGC→GGGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (24/19);  new base C (2/2);  total (26/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.39e-01

GGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCA  >  minE/2595582‑2595716
                                                                    |                                                                  
ggCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGcc                                                                  >  1:1587922/1‑71 (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:3026881/71‑1 (MQ=255)
                                 gTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAg                                 <  1:3897867/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:3402723/1‑71 (MQ=255)
                                                         gcgGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTgcg         >  1:981857/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:2526325/1‑71 (MQ=255)
                                                              ccACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGc    >  1:5602778/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:2682705/1‑71 (MQ=255)
                                                                acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCa  >  1:1896900/1‑71 (MQ=255)
                                                                acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCa  >  1:4029242/1‑71 (MQ=255)
                                                                acaGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCa  >  1:2875794/1‑71 (MQ=255)
                                                                    |                                                                  
GGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCA  >  minE/2595582‑2595716

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: