Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,539 T→C 9.5% A353A (GCA→GCG arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5390TC 9.5% 294.0 / 18.7 106A353A (GCA→GCGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base T (43/52);  new base C (8/2);  total (52/54)
Fisher's exact test for biased strand distribution p-value = 4.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 113 total aligned reads displayed.

GGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGA  >  minE/2595469‑2595609
                                                                      |                                                                      
     aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                   <  1:651935/71‑1 (MQ=255)
     aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                   <  1:3253416/71‑1 (MQ=255)
     aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                                                   <  1:3358893/71‑1 (MQ=255)
       tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATg                                                                       <  1:3513434/65‑1 (MQ=255)
       tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                                                 >  1:2636073/1‑71 (MQ=255)
        ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                                >  1:2552317/1‑71 (MQ=255)
        ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                                >  1:1271173/1‑71 (MQ=255)
        ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                                                >  1:1971017/1‑71 (MQ=255)
        ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGACAGATGCCGtgc                                                                 >  1:2625970/1‑70 (MQ=255)
          cGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctg                                                               <  1:2801501/70‑1 (MQ=255)
           ggTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctgc                                                              <  1:1404519/70‑1 (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:2535050/71‑1 (MQ=255)
                  aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                                      >  1:4115897/1‑71 (MQ=255)
                  aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGta                                                                  >  1:2271685/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)
                      gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                   <  1:643250/70‑1 (MQ=255)
                      gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                   >  1:1839690/1‑70 (MQ=255)
                          ccAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                               >  1:3492090/1‑70 (MQ=255)
                            aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                            >  1:1982458/1‑71 (MQ=255)
                            aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                            >  1:250142/1‑71 (MQ=255)
                            aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                            >  1:489177/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)
                                acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                                        >  1:402603/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)
                                     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:2578237/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: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:3892744/1‑71 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                             >  1:81717/1‑71 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                             >  1:3678812/1‑71 (MQ=255)
                                           aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                             >  1:3014557/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)
                                                 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)
                                                   cttcttTGATAACGCCAGATGCCGTGCTTCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTg                     <  1:1520909/71‑1 (MQ=255)
                                                   cttcttTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                                                   <  1:2122617/41‑1 (MQ=255)
                                                   cttcttTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTg                     <  1:1238127/71‑1 (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:3165685/68‑4 (MQ=255)
                                                              aCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCatcac             <  1:2500424/68‑4 (MQ=255)
                                                               cGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaa         <  1:477469/71‑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)
                                                                      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)
                                                                      |                                                                      
GGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGA  >  minE/2595469‑2595609

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: