Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,376 T→C 6.7% I408V (ATC→GTC)  arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,3760TC 6.7% 258.2 / 7.9 89I408V (ATC→GTC) arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base T (41/42);  new base C (2/4);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGT  >  minE/2595307‑2595446
                                                                     |                                                                      
aTTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGAt                                                                        >  1:1450432/1‑70 (MQ=255)
aTTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATa                                                                       <  1:14396/71‑1 (MQ=255)
  tACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAAt                                                                     >  1:1778553/1‑71 (MQ=255)
   aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATc                                                                    >  1:2465341/1‑71 (MQ=255)
   aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATc                                                                    >  1:1923899/1‑71 (MQ=255)
     aaaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATccc                                                                  <  1:3610673/71‑1 (MQ=255)
     aaaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATccc                                                                  <  1:1333661/71‑1 (MQ=255)
       aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                                                <  1:3023071/71‑1 (MQ=255)
        gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                                                >  1:2408914/1‑70 (MQ=255)
        gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                                                >  1:1961298/1‑70 (MQ=255)
        gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                                                >  1:2602856/1‑70 (MQ=255)
          gAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTgatgat                                                                        <  1:1408265/60‑1 (MQ=255)
           aaaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACg                                                            >  1:2924354/1‑71 (MQ=255)
            aaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                           <  1:1728026/71‑1 (MQ=255)
            aaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                           <  1:3705634/71‑1 (MQ=255)
            aaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                           <  1:2684986/71‑1 (MQ=255)
             agagagACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                           >  1:3004337/1‑70 (MQ=255)
             agagagACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                           >  1:1833044/1‑70 (MQ=255)
             agagagACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                           >  1:2619540/1‑70 (MQ=255)
             agagagACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGa                                                          <  1:3008348/71‑1 (MQ=255)
                 agaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                                                        >  1:816757/1‑69 (MQ=255)
                 agaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                                                        >  1:821700/1‑69 (MQ=255)
                 agaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTa                                                      <  1:3709359/71‑1 (MQ=255)
                  gaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                                                     <  1:3628161/71‑1 (MQ=255)
                  gaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                                                     <  1:778476/71‑1 (MQ=255)
                  gaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                                                     >  1:3399759/1‑71 (MQ=255)
                   aCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                                                     >  1:537398/1‑70 (MQ=255)
                      tAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGAt                                                                        >  1:78461/1‑48 (MQ=255)
                      tAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAata                                                   >  1:2637145/1‑69 (MQ=255)
                       aGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATaa                                                                      >  1:3272811/1‑49 (MQ=255)
                       aGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATcc                                                                   >  1:2510901/1‑52 (MQ=255)
                        gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATccc                                               >  1:2700483/1‑71 (MQ=255)
                        gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATcc                                                <  1:75801/70‑1 (MQ=255)
                           ccAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                                                     <  1:631034/62‑1 (MQ=255)
                             aGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCt                                          <  1:3509221/71‑1 (MQ=255)
                             aGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATACCCAGCt                                          <  1:2398680/71‑1 (MQ=255)
                              gcgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTAc                                                             <  1:4074871/51‑1 (MQ=255)
                               cgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                                                     <  1:2804827/58‑1 (MQ=255)
                                gcggccagaGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAt                                       >  1:304701/4‑71 (MQ=255)
                                gcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                        >  1:1051765/1‑70 (MQ=255)
                                gcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                        >  1:345183/1‑70 (MQ=255)
                                   gccagaGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATcccc                                              <  1:849417/61‑1 (MQ=255)
                                    ccagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTg                                         <  1:3610481/65‑1 (MQ=255)
                                     cagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa                                  >  1:542016/1‑71 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                           <  1:2152524/41‑1 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                           <  1:394732/41‑1 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                           <  1:1036500/41‑1 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGaa                                                         >  1:3783589/1‑43 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                                                        <  1:2691768/44‑1 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                                                        <  1:84994/44‑1 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGt                                                       <  1:1568847/45‑1 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGt                                                       <  1:3844983/45‑1 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                                                     <  1:3121005/47‑1 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc                             <  1:561201/71‑1 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc                             <  1:1822834/71‑1 (MQ=255)
                                          tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc                             >  1:2031137/1‑71 (MQ=255)
                                           cacaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGt                                                       <  1:849909/44‑1 (MQ=255)
                                            acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACg                                                            >  1:1060898/1‑38 (MQ=255)
                                            acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAat                                                    >  1:3633580/1‑46 (MQ=255)
                                            acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCa                                             >  1:3681957/1‑53 (MQ=255)
                                            acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTg                                         <  1:1752835/57‑1 (MQ=255)
                                            acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTg                                         <  1:484389/57‑1 (MQ=255)
                                            acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                           >  1:2506286/1‑71 (MQ=255)
                                            acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                           >  1:2165363/1‑71 (MQ=255)
                                             caAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttcttc                          <  1:1889340/71‑1 (MQ=255)
                                             caAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttcttc                          >  1:2525046/1‑71 (MQ=255)
                                             caAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttcttc                          <  1:2508316/71‑1 (MQ=255)
                                                   gCACAGGCAGGGTCATGTTAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGaa                    <  1:3516601/71‑1 (MQ=255)
                                                     acaGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAg                  <  1:1201508/71‑1 (MQ=255)
                                                      caggcaggGTCATGATAATCCCTGTACGGAAGTAata                                                   <  1:348770/37‑1 (MQ=255)
                                                      caggcaggGTCATGATAATCCCTGTACGGAAGTAata                                                   <  1:32898/37‑1 (MQ=255)
                                                        ggcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTc                                   <  1:848417/51‑1 (MQ=255)
                                                        ggcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTc                                   <  1:3288447/51‑1 (MQ=255)
                                                         gcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACg              >  1:330065/1‑71 (MQ=255)
                                                         gcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACg              >  1:51357/1‑71 (MQ=255)
                                                            gggTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttcttc                          <  1:974816/56‑1 (MQ=255)
                                                             ggTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTg            >  1:76264/1‑69 (MQ=255)
                                                               tCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGc        <  1:770669/71‑1 (MQ=255)
                                                                cATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGc        >  1:3696913/1‑70 (MQ=255)
                                                                cATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGc        >  1:2458076/1‑70 (MQ=255)
                                                                cATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttttt                           <  1:315477/51‑1 (MQ=255)
                                                                 aTGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGc        >  1:2529828/1‑69 (MQ=255)
                                                                 aTGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa      <  1:3659215/71‑1 (MQ=255)
                                                                 aTGATAATCACTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTa                >  1:4099395/1‑61 (MQ=255)
                                                                 aTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCa      >  1:1819255/1‑71 (MQ=255)
                                                                 aTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCa      <  1:3072840/71‑1 (MQ=255)
                                                                  tGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCag     <  1:2091029/71‑1 (MQ=255)
                                                                   gATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCat                                 >  1:1226182/1‑42 (MQ=255)
                                                                   gATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa      >  1:2284123/1‑69 (MQ=255)
                                                                   gATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa      >  1:2359636/1‑69 (MQ=255)
                                                                    aTAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCagag   >  1:1817348/1‑71 (MQ=255)
                                                                    aTAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa      >  1:2737528/1‑68 (MQ=255)
                                                                     tAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgt  <  1:528847/71‑1 (MQ=255)
                                                                     |                                                                      
ATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGT  >  minE/2595307‑2595446

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: