Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,538 A→T 9.3% S354T (TCT→ACT)  arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5380AT 9.3% 301.0 / 19.4 108S354T (TCT→ACT) arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (46/52);  new base T (8/2);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 9.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 110 total aligned reads displayed.

GGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAG  >  minE/2595469‑2595607
                                                                     |                                                                     
gggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                                                       <  1:1527354/70‑1 (MQ=255)
    aaTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg                                                                  >  1:2777550/1‑71 (MQ=255)
     aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                 <  1:3358893/71‑1 (MQ=255)
     aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                 <  1:3253416/71‑1 (MQ=255)
     aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                 <  1:651935/71‑1 (MQ=255)
     aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCg                                                                  >  1:2604821/1‑70 (MQ=255)
     aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCg                                                                  >  1:1822071/1‑70 (MQ=255)
       tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATg                                                                     <  1:3513434/65‑1 (MQ=255)
       tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                               >  1:2636073/1‑71 (MQ=255)
        ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                              >  1:1271173/1‑71 (MQ=255)
        ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                              >  1:2552317/1‑71 (MQ=255)
        ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGACAGATGCCGtgc                                                               >  1:2625970/1‑70 (MQ=255)
               ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGAt                                                                      >  1:2458431/1‑56 (MQ=255)
               ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGAt                                                                      >  1:265329/1‑56 (MQ=255)
               ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                                                   >  1:3030404/1‑59 (MQ=255)
               ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                        >  1:4091261/1‑70 (MQ=255)
               ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                                        <  1:4034034/70‑1 (MQ=255)
                 aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                                     >  1:3418224/1‑71 (MQ=255)
                 aaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                     >  1:4046826/1‑71 (MQ=255)
                  aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                    >  1:4115897/1‑71 (MQ=255)
                  aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                    <  1:2535050/71‑1 (MQ=255)
                  aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGta                                                                >  1:2271685/1‑58 (MQ=255)
                  aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGta                                                                >  1:3257359/1‑58 (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:1839690/1‑70 (MQ=255)
                         gCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcg                                             <  1:3030406/71‑1 (MQ=255)
                          ccAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                             >  1:3492090/1‑70 (MQ=255)
                            aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                          >  1:1982458/1‑71 (MQ=255)
                            aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                          >  1:489177/1‑71 (MQ=255)
                            aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                          >  1:250142/1‑71 (MQ=255)
                             aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                                         <  1:1189527/71‑1 (MQ=255)
                              tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCa                                        <  1:3010525/71‑1 (MQ=255)
                              tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCa                                        <  1:3039604/71‑1 (MQ=255)
                               cacaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                                     >  1:578681/1‑57 (MQ=255)
                               cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                               >  1:1219385/1‑47 (MQ=255)
                                acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                      >  1:3450223/1‑55 (MQ=255)
                                acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                      >  1:800071/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)
                                     ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                               >  1:1094884/1‑41 (MQ=255)
                                     ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                 <  1:2651554/39‑1 (MQ=255)
                                     ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                 <  1:1120083/39‑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: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:3892744/1‑71 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                           >  1:3678812/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)
                                               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:2203847/1‑71 (MQ=255)
                                                  gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt                    >  1:3478741/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)
                                                           aTAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata             <  1:2527304/69‑1 (MQ=255)
                                                             aaCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGt                  >  1:2260972/1‑62 (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:273273/71‑1 (MQ=255)
                                                                    gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag  <  1:897302/71‑1 (MQ=255)
                                                                    gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag  <  1:902455/71‑1 (MQ=255)
                                                                     aTGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga     <  1:2344369/67‑1 (MQ=255)
                                                                     |                                                                     
GGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAG  >  minE/2595469‑2595607

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: