Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,562 A→G 16.5% L346L (TTG→CTG)  arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5620AG16.5% 269.3 / 41.3 121L346L (TTG→CTG) arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (53/48);  new base G (9/11);  total (62/59)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 127 total aligned reads displayed.

CAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGC  >  minE/2595492‑2595630
                                                                      |                                                                    
cAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                                                                      <  1:770904/71‑1 (MQ=255)
  gCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                                                    <  1:2046393/71‑1 (MQ=255)
  gCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                                                                    >  1:4631505/1‑71 (MQ=255)
      aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                                                                <  1:5071896/71‑1 (MQ=255)
       tCACATTGGCATAAATCATTGCCTCTTTGACGCCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCAcc                                                                >  1:2467144/1‑70 (MQ=255)
        cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                                              >  1:464410/1‑71 (MQ=255)
         acaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgca                                                                  >  1:3832121/1‑65 (MQ=255)
         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                                              <  1:466127/70‑1 (MQ=255)
         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                             <  1:1037860/71‑1 (MQ=255)
         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                             <  1:3254161/71‑1 (MQ=255)
         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                             >  1:3459158/1‑71 (MQ=255)
         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                             <  1:5472608/71‑1 (MQ=255)
         acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                             >  1:2308503/1‑71 (MQ=255)
          caTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAcc                                                            <  1:1634087/71‑1 (MQ=255)
          caTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc                                                             <  1:1966156/70‑1 (MQ=255)
           aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAc                                                                     <  1:2983404/61‑1 (MQ=255)
           aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                                           <  1:4849085/71‑1 (MQ=255)
             tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                                              >  1:3920017/1‑66 (MQ=255)
             tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa                                                              >  1:3055542/1‑66 (MQ=255)
              ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGt                                                         <  1:1472548/70‑1 (MQ=255)
              ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGt                                                         <  1:1822723/70‑1 (MQ=255)
              ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGt                                                         <  1:4262823/70‑1 (MQ=255)
               gCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                                                        <  1:3422995/70‑1 (MQ=255)
               gCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                                                        <  1:2396974/70‑1 (MQ=255)
                  tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                                                                      >  1:5394825/1‑53 (MQ=255)
                  tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                                                    >  1:1124572/1‑71 (MQ=255)
                   aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAc                                                       >  1:1249439/1‑67 (MQ=255)
                   aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAc                                                       >  1:1088319/1‑67 (MQ=255)
                   aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                                   >  1:2560031/1‑71 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                                                                      >  1:2762815/1‑51 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                                                  >  1:5035245/1‑55 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                                                  >  1:2718291/1‑55 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                                                  >  1:611440/1‑55 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                                                  >  1:3736814/1‑55 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                                   <  1:4005571/70‑1 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                  >  1:4571849/1‑71 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                  >  1:190643/1‑71 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                  >  1:4762948/1‑71 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                  >  1:4005558/1‑71 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                  >  1:1911825/1‑71 (MQ=255)
                    aaCCATCGCTTCTTTGATAACGCCAGATGCCGAGCTGCCATCAATGGACAACGCGCCAAcc                                                            >  1:4243777/1‑61 (MQ=255)
                    aaCCATAGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                                                                  >  1:3215243/1‑55 (MQ=255)
                           gcctcttTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCatat                                            >  1:300896/4‑70 (MQ=255)
                           gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt                                           >  1:285528/1‑71 (MQ=255)
                            cctcttTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTg                                          <  1:4234961/69‑1 (MQ=255)
                               cttTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgc                                                                   <  1:5218190/43‑1 (MQ=255)
                               cttTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGt                                         <  1:1743118/69‑1 (MQ=255)
                                tttGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGtt                                        <  1:1311838/69‑1 (MQ=255)
                                 ttGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                                 >  1:1388473/1‑43 (MQ=255)
                                 ttGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                                 >  1:2577139/1‑43 (MQ=255)
                                 ttGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                                                                 >  1:367765/1‑43 (MQ=255)
                                  tGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                                                    >  1:389031/1‑55 (MQ=255)
                                  tGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCat                                     <  1:278853/70‑1 (MQ=255)
                                   gATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataa                                   >  1:467426/1‑71 (MQ=255)
                                    aTAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGt                                         >  1:5610185/1‑64 (MQ=255)
                                     cgaCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTc                                       <  1:2852098/63‑1 (MQ=255)
                                       aCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCat                                     <  1:5301797/65‑1 (MQ=255)
                                       aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATaga                               <  1:1441744/71‑1 (MQ=255)
                                       aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                <  1:5685742/70‑1 (MQ=255)
                                       aCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                <  1:4980414/70‑1 (MQ=255)
                                        cGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaa                              <  1:3890566/71‑1 (MQ=255)
                                         gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                                           >  1:2004573/1‑41 (MQ=255)
                                         gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                                                           >  1:1560367/1‑41 (MQ=255)
                                         gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaag                             <  1:1003764/71‑1 (MQ=255)
                                          ccAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga                            <  1:1625058/71‑1 (MQ=255)
                                                gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                                   >  1:1026562/1‑42 (MQ=255)
                                                gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                                   >  1:3029323/1‑42 (MQ=255)
                                                gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                                   >  1:5313032/1‑42 (MQ=255)
                                                gCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTg                                          >  1:3789575/1‑51 (MQ=255)
                                                gCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGt                                         >  1:2675742/1‑52 (MQ=255)
                                                 ccGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGtt                                        >  1:4158798/1‑52 (MQ=255)
                                                 ccGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCat                                     <  1:1134433/55‑1 (MQ=255)
                                                 ccGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCat                                     <  1:4648604/55‑1 (MQ=255)
                                                  cGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGaaatg                    >  1:4332174/1‑69 (MQ=255)
                                                    tgctgcCATCAATGGACAACGCGCCAACCAGTACCGTCgg                                                 <  1:2187126/40‑1 (MQ=255)
                                                    tgctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCg                  <  1:2098620/71‑1 (MQ=255)
                                                      ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCat                                              >  1:569487/1‑41 (MQ=255)
                                                         acATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAAAAGAAGAGAGGAAGGCGGTGAg             >  1:2566778/2‑71 (MQ=255)
                                                          cgtCAATCGACAGCGCACCAATCAGTACCGTCGGCATATTGTTCat                                     >  1:1992271/3‑46 (MQ=25)
                                                          cgtCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCat                                     >  1:566624/3‑46 (MQ=255)
                                                            tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTg                                          <  1:2668579/39‑1 (MQ=255)
                                                            tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTg                                          <  1:1539955/39‑1 (MQ=255)
                                                            tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTg                                          <  1:3412685/39‑1 (MQ=255)
                                                            tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGaa          <  1:4671686/71‑1 (MQ=255)
                                                             cAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAAt         >  1:606189/1‑71 (MQ=255)
                                                              aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                                    >  1:1108406/1‑43 (MQ=255)
                                                              aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                                    >  1:2090024/1‑43 (MQ=255)
                                                              aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                                    >  1:1973001/1‑43 (MQ=255)
                                                              aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                                    >  1:2973408/1‑43 (MQ=255)
                                                              aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagag                           >  1:2789313/1‑52 (MQ=255)
                                                              aaTCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGaagc        <  1:4929070/71‑3 (MQ=255)
                                                              aaTCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGaagc        <  1:3490046/71‑3 (MQ=255)
                                                               aTGGACAACCCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATc        >  1:3660895/1‑70 (MQ=255)
                                                                tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                <  1:2407353/45‑1 (MQ=255)
                                                                tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                <  1:3637119/45‑1 (MQ=255)
                                                                tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATag                                <  1:3300775/45‑1 (MQ=255)
                                                                tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt                <  1:218872/61‑1 (MQ=255)
                                                                  gACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGt    <  1:3523388/71‑1 (MQ=255)
                                                                  gACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGaa                      >  1:1057806/1‑53 (MQ=255)
                                                                    cAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGc  >  1:3578976/1‑71 (MQ=255)
                                                                      |                                                                    
CAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGC  >  minE/2595492‑2595630

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: