Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,623 T→G 37.3% G325G (GGA→GGC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,6230TG37.3% 14.2 / 47.2 43G325G (GGA→GGCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base T (15/12);  new base G (6/10);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

CAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCC  >  minE/2595553‑2595692
                                                                      |                                                                     
cAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAAt                                                                       >  1:606189/1‑71 (MQ=255)
  aTGGACAACCCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATc                                                                      >  1:3660895/1‑70 (MQ=255)
     gACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGt                                                                  <  1:3523388/71‑1 (MQ=255)
       cAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGc                                                                >  1:3578976/1‑71 (MQ=255)
         gcgcaCCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCg                                                              <  1:4148179/70‑1 (MQ=255)
         gcgcaCCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCg                                                              <  1:3499369/70‑1 (MQ=255)
          cgcgCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCgg                                                                   >  1:1755250/1‑65 (MQ=255)
             accAATCAGCACCGTCGGCATATTGTTCATGACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCGaac                                                           <  1:4954810/69‑3 (MQ=255)
             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)
                                   ttGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCgg                                                                   >  1:3211892/1‑40 (MQ=255)
                                   ttGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCgg                                                                   >  1:5231212/1‑40 (MQ=255)
                                   ttGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCgg                                                                   >  1:5356225/1‑40 (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:4123460/64‑1 (MQ=255)
                                                 tgagagaAATGCCGTCAAGAAGCCGGTGCCGAACGTCGCTGCCCATAACCCCTTGTCTGCCAGCa                            <  1:580369/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)
                                                                     aTCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTcc  <  1:2700788/71‑1 (MQ=255)
                                                                      tCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGaaa          >  1:4175024/1‑62 (MQ=255)
                                                                      tCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGata      >  1:4961856/1‑66 (MQ=255)
                                                                      |                                                                     
CAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGCCCACAGGCCGTTATCCGCCAGCACGTTGAGTACGCCAGAAAGATATTCC  >  minE/2595553‑2595692

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: