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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5480GC 8.5% 277.5 / 12.9 94G350G (GGC→GGGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (40/45);  new base C (6/2);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGT  >  minE/2595479‑2595616
                                                                     |                                                                    
cGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctg                                                                      <  1:2801501/70‑1 (MQ=255)
 ggTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctgc                                                                     <  1:1404519/70‑1 (MQ=255)
     ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                                 <  1:4034034/70‑1 (MQ=255)
     ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                                 >  1:4091261/1‑70 (MQ=255)
       aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                              >  1:3418224/1‑71 (MQ=255)
       aaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                              >  1:4046826/1‑71 (MQ=255)
        aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                             <  1:2535050/71‑1 (MQ=255)
        aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                             >  1:4115897/1‑71 (MQ=255)
         aTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                             >  1:360019/1‑70 (MQ=255)
          tcacAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                                                           >  1:381027/4‑71 (MQ=255)
          tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                             <  1:3067505/69‑1 (MQ=255)
          tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                             <  1:766684/69‑1 (MQ=255)
          tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                             <  1:2573050/69‑1 (MQ=255)
            gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                          <  1:643250/70‑1 (MQ=255)
            gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                          >  1:1839690/1‑70 (MQ=255)
               gCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcg                                                      <  1:3030406/71‑1 (MQ=255)
                ccAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                                      >  1:3492090/1‑70 (MQ=255)
                  aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                   >  1:250142/1‑71 (MQ=255)
                  aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                   >  1:489177/1‑71 (MQ=255)
                  aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                   >  1:1982458/1‑71 (MQ=255)
                   aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                                                  <  1:1189527/71‑1 (MQ=255)
                    tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCa                                                 <  1:3039604/71‑1 (MQ=255)
                    tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCa                                                 <  1:3010525/71‑1 (MQ=255)
                     cacaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                              >  1:578681/1‑57 (MQ=255)
                      acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                               >  1:800071/1‑55 (MQ=255)
                      acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                               >  1:402603/1‑55 (MQ=255)
                      acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                               >  1:3450223/1‑55 (MQ=255)
                      acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                               >  1:1472393/1‑71 (MQ=255)
                      acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                               >  1:3792116/1‑71 (MQ=255)
                       caTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                               <  1:2885897/70‑1 (MQ=255)
                        aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                                                  >  1:3042308/1‑66 (MQ=255)
                        aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                             <  1:2423512/71‑1 (MQ=255)
                          tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAg                                            <  1:1339796/70‑1 (MQ=255)
                           ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                                          >  1:843990/1‑71 (MQ=255)
                           ggCATAAACCATCGCTTCTTTGATAACGCCAGAAGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                                          >  1:2676502/1‑71 (MQ=255)
                             cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                               <  1:920206/64‑1 (MQ=255)
                             cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                        <  1:2098857/71‑1 (MQ=255)
                             cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                        <  1:668832/71‑1 (MQ=255)
                             cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                        <  1:2447130/71‑1 (MQ=255)
                             cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                        <  1:2578237/71‑1 (MQ=255)
                               tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                             <  1:103480/48‑1 (MQ=255)
                               tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                             <  1:495801/48‑1 (MQ=255)
                               tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                               <  1:190327/62‑1 (MQ=255)
                               tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                        <  1:3149685/69‑1 (MQ=255)
                               tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                                      >  1:2218321/1‑71 (MQ=255)
                               tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                                      >  1:3616433/1‑71 (MQ=255)
                                aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                               >  1:3523999/1‑45 (MQ=255)
                                aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                     >  1:12138/1‑71 (MQ=255)
                                 aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                             >  1:357494/1‑62 (MQ=255)
                                 aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                    >  1:3014557/1‑71 (MQ=255)
                                 aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                    >  1:3678812/1‑71 (MQ=255)
                                 aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                    >  1:81717/1‑71 (MQ=255)
                                 aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                    >  1:3892744/1‑71 (MQ=255)
                                 aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                    >  1:3349748/1‑71 (MQ=255)
                                 aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                    <  1:310282/71‑1 (MQ=255)
                                 aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGAACAACGCGCCAACCAGTACCGTCg                                    <  1:2289317/71‑1 (MQ=255)
                                  aCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                          <  1:837545/48‑1 (MQ=255)
                                  aCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                          <  1:3931714/48‑1 (MQ=255)
                                    cATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCa                                 <  1:2187272/71‑1 (MQ=255)
                                     aTCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                          <  1:2311672/45‑1 (MQ=255)
                                     aTCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCCGCat                                <  1:1576686/71‑1 (MQ=255)
                                      tCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                             >  1:2786869/1‑57 (MQ=255)
                                      tCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCata                               <  1:3855249/71‑1 (MQ=255)
                                       cGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCatat                              >  1:1477518/1‑71 (MQ=255)
                                        gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCa                                                                >  1:2743605/1‑36 (MQ=255)
                                        gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                              >  1:3742028/1‑38 (MQ=255)
                                        gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                                                        >  1:2183550/1‑44 (MQ=255)
                                        gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt                             >  1:3478741/1‑71 (MQ=255)
                                        gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt                             >  1:2203847/1‑71 (MQ=255)
                                         cttcttTGATAACGCCAGATGCCGTGCTTCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTg                            <  1:1520909/71‑1 (MQ=255)
                                         cttcttTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                          <  1:2122617/41‑1 (MQ=255)
                                         cttcttTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTg                            <  1:1238127/71‑1 (MQ=255)
                                              ttGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCa                        >  1:1295209/1‑70 (MQ=255)
                                               tGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCat                       <  1:1604755/70‑1 (MQ=255)
                                                 aTAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataa                     <  1:976142/70‑1 (MQ=255)
                                                 aTAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                      <  1:2527304/69‑1 (MQ=255)
                                                   aaCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGt                           >  1:2260972/1‑62 (MQ=255)
                                                    aCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCatcac                    <  1:3165685/68‑4 (MQ=255)
                                                    aCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCatcac                    <  1:2500424/68‑4 (MQ=255)
                                                     cGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaa                <  1:477469/71‑1 (MQ=255)
                                                      gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCat                                <  1:424490/54‑1 (MQ=255)
                                                       ccAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga              <  1:1733267/71‑1 (MQ=255)
                                                        cAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga              <  1:57482/70‑1 (MQ=255)
                                                          gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag           <  1:902455/71‑1 (MQ=255)
                                                          gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag           <  1:897302/71‑1 (MQ=255)
                                                          gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag           <  1:273273/71‑1 (MQ=255)
                                                           tcgCCGTACTCCCGTCAATCGACAGCGCACCAATCATCACCGTCGGCATATTGTTAATAACTGAtgaga            >  1:1648825/3‑69 (MQ=255)
                                                           aTGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga              <  1:2344369/67‑1 (MQ=255)
                                                            tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                    >  1:2814061/1‑44 (MQ=255)
                                                            tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagaga            >  1:265722/1‑68 (MQ=255)
                                                            tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGa         <  1:3914127/71‑1 (MQ=255)
                                                            tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGCCATATTGTTCATAATAGAAGAGAgg          >  1:3063473/1‑70 (MQ=255)
                                                            cgCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGAtgagaga          >  1:301409/2‑69 (MQ=255)
                                                             gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAgg          >  1:2012762/1‑69 (MQ=255)
                                                             gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGa         <  1:2279910/70‑1 (MQ=255)
                                                                gTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCata                               >  1:1538798/1‑45 (MQ=255)
                                                                 tgctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataata                   >  1:1029990/1‑56 (MQ=255)
                                                                   ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                      >  1:3794525/1‑51 (MQ=255)
                                                                   ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt  <  1:2412370/71‑1 (MQ=255)
                                                                   ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt  >  1:1287844/1‑71 (MQ=255)
                                                                     |                                                                    
CGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGT  >  minE/2595479‑2595616

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: