Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,485 C→T 24.8% L371L (TTG→TTA arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,4850CT24.8% 227.4 / 73.9 129L371L (TTG→TTAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (43/54);  new base T (21/11);  total (64/65)
Fisher's exact test for biased strand distribution p-value = 4.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

Only 100 of 144 total aligned reads displayed.

TTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAA  >  minE/2595417‑2595555
                                                                    |                                                                      
 tcttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTccc                                                                        >  1:1699749/1‑68 (MQ=255)
 tcttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTccc                                                                        >  1:1469863/1‑68 (MQ=255)
  cttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCaaa                                                                     <  1:4285372/70‑1 (MQ=255)
   ttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATc                                                                   >  1:696216/1‑71 (MQ=255)
    tctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCg                                                                  >  1:855357/1‑71 (MQ=255)
       gCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAg                                                               >  1:2687148/1‑71 (MQ=255)
       gCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAg                                                               >  1:4940678/1‑71 (MQ=255)
       gCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAg                                                               >  1:3120286/1‑71 (MQ=255)
        cGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGc                                                                 <  1:2793191/68‑1 (MQ=255)
         gAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATc                                                                   <  1:3178037/65‑1 (MQ=255)
         gAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGc                                                              >  1:3121949/1‑70 (MQ=255)
          aaaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCa                                                            <  1:1218542/71‑1 (MQ=255)
          aaaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCa                                                            <  1:1807966/71‑1 (MQ=255)
           aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCaa                                                           <  1:1582911/71‑1 (MQ=255)
           aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCaa                                                           <  1:3251813/71‑1 (MQ=255)
            aGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAAt                                                          >  1:4120768/1‑71 (MQ=255)
            aGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAAt                                                          >  1:3615097/1‑71 (MQ=255)
              tACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATca                                                        <  1:2751941/71‑1 (MQ=255)
               acatGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                                                        >  1:1514912/4‑70 (MQ=255)
                 gTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACAt                                                     <  1:4778729/71‑1 (MQ=255)
                 gTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACAt                                                     <  1:3790476/71‑1 (MQ=255)
                 atGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACAt                                                     <  1:3814260/70‑1 (MQ=255)
                      gCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACAtt                                                    <  1:3803179/67‑1 (MQ=255)
                        cAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                              <  1:1972875/71‑1 (MQ=255)
                         agagTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATaa                                             <  1:4839785/71‑1 (MQ=255)
                         agagTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                              >  1:2716251/1‑70 (MQ=255)
                           agTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAc                                           <  1:3884190/71‑1 (MQ=255)
                           agTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAc                                           <  1:5198531/71‑1 (MQ=255)
                             tagcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCa                                                                       <  1:507659/41‑1 (MQ=255)
                             tagcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCa                                                                       <  1:892936/41‑1 (MQ=255)
                             tagcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCa                                                                       <  1:698790/41‑1 (MQ=255)
                                  ttgcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATcac                                                       <  1:5693837/50‑1 (MQ=255)
                                  ttgcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCAt                                        >  1:308917/3‑67 (MQ=255)
                                  ttgcCAGACTGCCAACCGGGGTGATTTTCGGGCCTAAATCACAGCCAATcac                                                       <  1:5169098/50‑1 (MQ=255)
                                  tagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATcac                                                       <  1:4060700/52‑1 (MQ=255)
                                  tagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCAt                                        <  1:2479251/67‑1 (MQ=255)
                                  tagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCAt                                        <  1:5261264/67‑1 (MQ=255)
                                   agcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAACCGCAGCCAATCACATTGGCATAAACCATCGctt                                   >  1:1749143/1‑71 (MQ=255)
                                    gcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                              <  1:4551350/59‑1 (MQ=255)
                                    gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctc                                  >  1:138197/1‑71 (MQ=255)
                                    gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctc                                  >  1:4111743/1‑71 (MQ=255)
                                     ccAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAcc                                          <  1:1278617/62‑1 (MQ=255)
                                       aGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                               <  1:1136169/71‑1 (MQ=255)
                                         gCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                              <  1:5621723/54‑1 (MQ=255)
                                         gCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                              <  1:5274453/54‑1 (MQ=255)
                                         gCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGa                             >  1:3076353/1‑71 (MQ=255)
                                          ctgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctct                                 >  1:3569147/4‑66 (MQ=255)
                                          ctgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctct                                 >  1:2284202/4‑66 (MQ=255)
                                           tgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                                                        <  1:5516054/40‑1 (MQ=16)
                                           tgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                                                        <  1:5408309/40‑1 (MQ=16)
                                           tACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGAt                            >  1:2367820/1‑70 (MQ=255)
                                           tACCAATTGGGGTAATTTTAGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGAt                            >  1:3327418/1‑70 (MQ=255)
                                             ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCATTCACATTGGCATAAACCATCGCTTCTTTGa                             >  1:3344975/1‑67 (MQ=255)
                                             ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATc                                                         >  1:1156533/1‑39 (MQ=255)
                                             ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATc                                       <  1:4385994/57‑1 (MQ=255)
                                             ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAAc                         >  1:3451052/1‑71 (MQ=255)
                                             ccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGacg                           <  1:530882/69‑3 (MQ=255)
                                              cAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATcaca                                                      <  1:5446397/41‑1 (MQ=255)
                                              cAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATcaca                                                      <  1:2562623/41‑1 (MQ=255)
                                               aaTCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                       <  1:2905130/71‑1 (MQ=255)
                                                 ttGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                       >  1:3177689/1‑69 (MQ=255)
                                                  tGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                       >  1:4094631/1‑68 (MQ=255)
                                                   ggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgt                   >  1:1923256/1‑70 (MQ=255)
                                                   ggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgt                   >  1:4442622/1‑70 (MQ=255)
                                                   ggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                    <  1:3745095/70‑1 (MQ=255)
                                                   ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                          >  1:2880212/1‑64 (MQ=255)
                                                   ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                          >  1:623043/1‑64 (MQ=255)
                                                   ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                       >  1:4763680/1‑67 (MQ=255)
                                                   ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                       >  1:4738066/1‑67 (MQ=255)
                                                    gggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                  >  1:3301054/1‑69 (MQ=255)
                                                    gggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                  >  1:5121328/1‑69 (MQ=255)
                                                    gggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                   <  1:720008/70‑1 (MQ=255)
                                                    gggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                   <  1:2410889/70‑1 (MQ=255)
                                                      gTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCa                                         <  1:5212239/46‑1 (MQ=255)
                                                      gTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTg                              <  1:4444327/57‑1 (MQ=255)
                                                       tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc               >  1:4639318/1‑71 (MQ=255)
                                                       tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc               <  1:848510/71‑1 (MQ=255)
                                                        aaTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                              <  1:5461221/39‑1 (MQ=255)
                                                        aaTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg              >  1:1316500/1‑71 (MQ=255)
                                                         aTTTTCGGTCCCAAATCTCAGCCAATCACATTGGCATAAACAATCTATTCTTTGATAACGCCAGATGcc               >  1:5660729/1‑69 (MQ=255)
                                                          ttttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                               <  1:4251608/52‑1 (MQ=255)
                                                          ttttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt             <  1:2214288/70‑1 (MQ=255)
                                                           tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc           >  1:2993080/1‑71 (MQ=255)
                                                           tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGACAGATGCCGtg            <  1:1286167/70‑1 (MQ=255)
                                                           tttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtac           >  1:947643/1‑69 (MQ=255)
                                                           tttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtac           >  1:4252284/1‑69 (MQ=255)
                                                            ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct          >  1:1945590/1‑71 (MQ=255)
                                                             tCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct          >  1:803359/1‑70 (MQ=255)
                                                                gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                               >  1:2205122/1‑46 (MQ=255)
                                                                gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                               <  1:1582934/46‑1 (MQ=255)
                                                                gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGcc                      >  1:1162773/1‑55 (MQ=255)
                                                                  cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc               >  1:5393428/1‑60 (MQ=255)
                                                                  cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc               >  1:3994363/1‑60 (MQ=255)
                                                                  cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc               >  1:3866857/1‑60 (MQ=255)
                                                                  cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt     >  1:763433/1‑70 (MQ=255)
                                                                   ccAAATCGCAGCCAATCACATTGGCATAAACCATCGcttct                                 <  1:3713674/41‑1 (MQ=255)
                                                                   ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc    <  1:5299032/70‑1 (MQ=255)
                                                                   ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc    <  1:1722350/70‑1 (MQ=255)
                                                                   ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc    <  1:3367170/70‑1 (MQ=255)
                                                                    cAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa  <  1:3422265/71‑1 (MQ=255)
                                                                    |                                                                      
TTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAA  >  minE/2595417‑2595555

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: