Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,649 G→A 8.5% L317L (CTG→TTG) ‡ arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,6490GA 8.5% 131.2 / 2.1 47L317L (CTG→TTG) ‡arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (24/19);  new base A (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 = 1.00e+00

GTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCA  >  minE/2595579‑2595716
                                                                      |                                                                   
gTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAg                                                                     <  1:1573231/71‑1 (MQ=255)
   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)
                                                                      |                                                                   
GTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCA  >  minE/2595579‑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: