Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,551 A→G 8.1% D349D (GAT→GAC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5510AG 8.1% 279.3 / 12.6 98D349D (GAT→GACarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (43/47);  new base G (6/2);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGT  >  minE/2595484‑2595616
                                                                   |                                                                 
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:2573050/69‑1 (MQ=255)
     tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                             <  1:3067505/69‑1 (MQ=255)
     tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                             <  1:766684/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:1982458/1‑71 (MQ=255)
             aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                   >  1:489177/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:3792116/1‑71 (MQ=255)
                 acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                               >  1:1472393/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:2578237/71‑1 (MQ=255)
                        cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                        <  1:668832/71‑1 (MQ=255)
                        cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                        <  1:2447130/71‑1 (MQ=255)
                        cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                        <  1:2098857/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:3616433/1‑71 (MQ=255)
                          tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                                      >  1:2218321/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:3678812/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)
                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                    >  1:81717/1‑71 (MQ=255)
                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                    >  1:3014557/1‑71 (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:897302/71‑1 (MQ=255)
                                                     gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag           <  1:902455/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)
                                                                   aTCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt  <  1:3900776/66‑1 (MQ=255)
                                                                   |                                                                 
CCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGT  >  minE/2595484‑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: