Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,380 C→G 7.1% G406G (GGG→GGC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,3800CG 7.1% 259.4 / 6.5 84G406G (GGG→GGCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (36/42);  new base G (2/4);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGC  >  minE/2595310‑2595449
                                                                      |                                                                     
aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATc                                                                       >  1:2465341/1‑71 (MQ=255)
aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATc                                                                       >  1:1923899/1‑71 (MQ=255)
  aaaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATccc                                                                     <  1:1333661/71‑1 (MQ=255)
  aaaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATccc                                                                     <  1:3610673/71‑1 (MQ=255)
    aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                                                   <  1:3023071/71‑1 (MQ=255)
     gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                                                   >  1:2408914/1‑70 (MQ=255)
     gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                                                   >  1:2602856/1‑70 (MQ=255)
     gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                                                   >  1:1961298/1‑70 (MQ=255)
        aaaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACg                                                               >  1:2924354/1‑71 (MQ=255)
         aaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                              <  1:1728026/71‑1 (MQ=255)
         aaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                              <  1:2684986/71‑1 (MQ=255)
         aaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                              <  1:3705634/71‑1 (MQ=255)
          agagagACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                              >  1:1833044/1‑70 (MQ=255)
          agagagACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                                              >  1:3004337/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:3399759/1‑71 (MQ=255)
               gaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                                                        <  1:778476/71‑1 (MQ=255)
                aCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                                                        >  1:537398/1‑70 (MQ=255)
                   tAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAata                                                      >  1:2637145/1‑69 (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:345183/1‑70 (MQ=255)
                             gcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                           >  1:1051765/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:84994/44‑1 (MQ=255)
                                       tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                                                           <  1:2691768/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:2165363/1‑71 (MQ=255)
                                         acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                              >  1:2506286/1‑71 (MQ=255)
                                          caAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttcttc                             >  1:2525046/1‑71 (MQ=255)
                                          caAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttcttc                             <  1:1889340/71‑1 (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:3072840/71‑1 (MQ=255)
                                                              aTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCa         >  1:1819255/1‑71 (MQ=255)
                                                               tGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCag        <  1:2091029/71‑1 (MQ=255)
                                                                gATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCat                                    >  1:1226182/1‑42 (MQ=255)
                                                                gATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa         >  1:2359636/1‑69 (MQ=255)
                                                                gATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa         >  1:2284123/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)
                                                                    aTCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta    <  1:2374533/70‑1 (MQ=255)
                                                                     tCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagc  <  1:3396641/71‑1 (MQ=255)
                                                                      |                                                                     
ACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGC  >  minE/2595310‑2595449

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: