Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,596,106 G→C 15.4% S164S (TCC→TCG)  arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,596,1060GC15.4% 232.3 / 31.3 97S164S (TCC→TCG) arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (36/46);  new base C (7/8);  total (43/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

ATCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGG  >  Exported/2596038‑2596170
                                                                     |                                                                
aTCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGa                                                                <  1:350982/71‑1 (MQ=255)
   gcggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                              <  1:4626850/70‑1 (MQ=255)
   gcggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                              <  1:4620175/70‑1 (MQ=255)
     ggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGa                                                           <  1:1529366/71‑1 (MQ=255)
       cgATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                         >  1:5208722/1‑71 (MQ=255)
          tatCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGa                                                                <  1:2569275/61‑1 (MQ=255)
           atCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCa                                                     >  1:3750295/1‑71 (MQ=255)
            tCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                         >  1:4212338/1‑66 (MQ=255)
            tCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                         >  1:4246305/1‑66 (MQ=255)
            tCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                         >  1:3468187/1‑66 (MQ=255)
            tCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTc                                                      >  1:2542745/1‑69 (MQ=255)
              cACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAcc                                                   <  1:4238865/70‑1 (MQ=255)
              cACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCa                                                  <  1:2734764/71‑1 (MQ=255)
                ccGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGa                                                                <  1:245251/55‑1 (MQ=255)
                 cGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGt                                               <  1:3055166/71‑1 (MQ=255)
                  ggCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                              >  1:878055/1‑55 (MQ=255)
                  ggCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGtt                                              <  1:2563728/71‑1 (MQ=255)
                   gCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                             <  1:3655638/71‑1 (MQ=255)
                    gatcATAACGGAGGCGTACT‑GCGCAAAGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGt                                               >  1:1748718/4‑68 (MQ=255)
                      ccATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                             <  1:3571826/52‑1 (MQ=255)
                       cATAACGGAGGCGTACT‑GCGCAAAGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAgaaa                                         >  1:2358490/1‑69 (MQ=255)
                       cATAACGGAGGCGTACT‑GCGCAAAGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAgaaa                                         >  1:3207204/1‑69 (MQ=255)
                       cATAACGGAGGCGTACT‑GCGCAAAGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAgaa                                          >  1:4646796/1‑69 (MQ=255)
                         tCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                             <  1:2312640/49‑1 (MQ=255)
                           aCCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                             >  1:467677/1‑63 (MQ=255)
                           aCCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACaa                                      >  1:2641538/1‑70 (MQ=255)
                           aCCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAAt                                     <  1:1330378/71‑1 (MQ=255)
                           aCCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAAt                                     <  1:249723/71‑1 (MQ=255)
                             cGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGa                                                           <  1:1179556/47‑1 (MQ=255)
                             cGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAAt                                     >  1:3443319/1‑69 (MQ=255)
                                 gCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCgg                               >  1:1004157/1‑71 (MQ=255)
                                  cGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGc                              <  1:4355585/71‑1 (MQ=255)
                                    tATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                             <  1:2006210/38‑1 (MQ=255)
                                    tATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACg                                                            <  1:3665311/39‑1 (MQ=255)
                                    tATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACg                                                            <  1:5605021/39‑1 (MQ=255)
                                    tATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACg                                                            <  1:3508099/39‑1 (MQ=255)
                                    tATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAg                                                 <  1:1461946/50‑1 (MQ=255)
                                    tATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAg                                                 <  1:3006230/50‑1 (MQ=255)
                                    tATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGACATTCACCAg                                                 <  1:1082108/50‑1 (MQ=255)
                                     aTTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAgg                                                >  1:349147/1‑50 (MQ=255)
                                     aTTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                                    <  1:1196447/62‑1 (MQ=255)
                                      ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                             <  1:2061554/36‑1 (MQ=255)
                                      ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                                    <  1:223291/45‑1 (MQ=255)
                                      ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                             >  1:2173524/1‑52 (MQ=255)
                                      ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                             >  1:4434754/1‑52 (MQ=255)
                                      ttCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGAAATTCACAAGGTTGGAGACAATAAGCGGCAgg                           >  1:2493236/1‑70 (MQ=255)
                                       tCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                             <  1:3576795/35‑1 (MQ=255)
                                       tCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                         >  1:2091012/1‑39 (MQ=255)
                                       tCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAg                                  >  1:5291843/1‑62 (MQ=255)
                                        cgcg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                                    >  1:5292188/1‑43 (MQ=255)
                                         gcg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                         <  1:1122080/37‑1 (MQ=255)
                                          cgcaaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTgg                       >  1:3387305/4‑71 (MQ=255)
                                          cgcaaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTgg                       >  1:3393996/4‑71 (MQ=255)
                                          cg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                      >  1:1196194/1‑71 (MQ=255)
                                          cg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                      <  1:1877612/71‑1 (MQ=255)
                                          cg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                      >  1:351590/1‑71 (MQ=255)
                                          cg‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                      <  1:3513256/71‑1 (MQ=255)
                                           g‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGtt                                              <  1:519162/46‑1 (MQ=255)
                                           g‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGtt                                              <  1:4460923/46‑1 (MQ=255)
                                           g‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCg                     <  1:3270281/71‑1 (MQ=255)
                                            caaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACa                                       <  1:5531407/52‑1 (MQ=255)
                                             aaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCgg                               <  1:3444343/60‑1 (MQ=255)
                                             aaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCgg                               <  1:11626/60‑1 (MQ=255)
                                             aaaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTTGAGACAATAAGCGGCAGGCTGGCgg                    >  1:2759706/1‑71 (MQ=255)
                                               aGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCa                                                  <  1:5092252/39‑1 (MQ=255)
                                               aGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTa                  >  1:4325507/1‑71 (MQ=255)
                                                gCCGAGGCCAAATAAATCAGCGGAAACGATATTCACCAGGtt                                              <  1:4819783/42‑1 (MQ=255)
                                                     ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                           >  1:5066925/1‑40 (MQ=255)
                                                     ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                           >  1:4091978/1‑40 (MQ=255)
                                                     ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                           >  1:4005825/1‑40 (MQ=255)
                                                     ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                           >  1:3708164/1‑40 (MQ=255)
                                                     ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                           >  1:1260618/1‑40 (MQ=255)
                                                     ggCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                           >  1:3288686/1‑40 (MQ=255)
                                                      gCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                             <  1:15690/37‑1 (MQ=255)
                                                      gCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                                    <  1:303257/46‑1 (MQ=255)
                                                       ccgaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCaa           <  1:4561973/67‑1 (MQ=255)
                                                       ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                                    <  1:5329182/45‑1 (MQ=255)
                                                       ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt          <  1:61551/71‑1 (MQ=255)
                                                        cgaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAta         <  1:1027594/69‑2 (MQ=255)
                                                        cgaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAt          >  1:568668/3‑70 (MQ=255)
                                                         gaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAtaa        <  1:1074104/70‑3 (MQ=255)
                                                         gaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAtaa        <  1:2907473/70‑3 (MQ=255)
                                                         aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCa            <  1:254872/67‑1 (MQ=255)
                                                         aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGa        <  1:2147117/71‑1 (MQ=255)
                                                          aagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGc             <  1:2102362/65‑1 (MQ=255)
                                                          aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa       <  1:1532620/71‑1 (MQ=255)
                                                           aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCgg                               <  1:2897764/46‑1 (MQ=255)
                                                           aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAt      >  1:343780/1‑71 (MQ=255)
                                                           aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAt      >  1:2210402/1‑71 (MQ=255)
                                                            gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                                    >  1:4757351/1‑40 (MQ=255)
                                                            gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                                    >  1:627901/1‑40 (MQ=255)
                                                            gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTgg                       >  1:5434049/1‑53 (MQ=255)
                                                            gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc     <  1:1474192/71‑1 (MQ=255)
                                                            gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc     <  1:1231880/71‑1 (MQ=255)
                                                            gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc     >  1:8113/1‑71 (MQ=255)
                                                               aTCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCgg  >  1:25208/1‑71 (MQ=255)
                                                                 cAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                    >  1:3068564/1‑51 (MQ=255)
                                                                    aggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc             <  1:920563/54‑1 (MQ=255)
                                                                     |                                                                
ATCGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCG‑AAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGG  >  Exported/2596038‑2596170

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: ‑