Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,596,136 A→G 30.5% L154L (CTT→CTC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,596,1360AG30.5% 101.7 / 61.6 82L154L (CTT→CTCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (24/33);  new base G (10/15);  total (34/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

GCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACT  >  Exported/2596071‑2596204
                                                                 |                                                                    
gCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCgg                                                                 >  1:1004157/1‑71 (MQ=255)
 cGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGc                                                                <  1:4355585/71‑1 (MQ=255)
    aTTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                                                                      <  1:1196447/62‑1 (MQ=255)
     ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGAAATTCACAAGGTTGGAGACAATAAGCGGCAgg                                                             >  1:2493236/1‑70 (MQ=255)
      tCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAg                                                                    >  1:5291843/1‑62 (MQ=255)
        cgcaaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTgg                                                         >  1:3393996/4‑71 (MQ=255)
        cgcaaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTgg                                                         >  1:3387305/4‑71 (MQ=255)
         cgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                                                        >  1:351590/1‑71 (MQ=255)
         cgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                                                        <  1:3513256/71‑1 (MQ=255)
         cgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                                                        >  1:1196194/1‑71 (MQ=255)
         cgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                                                        <  1:1877612/71‑1 (MQ=255)
          gAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCg                                                       <  1:3270281/71‑1 (MQ=255)
           aaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCgg                                                                 <  1:3444343/60‑1 (MQ=255)
           aaaGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCgg                                                                 <  1:11626/60‑1 (MQ=255)
           aaaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTTGAGACAATAAGCGGCAGGCTGGCgg                                                      >  1:2759706/1‑71 (MQ=255)
             aGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTa                                                    >  1:4325507/1‑71 (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:2210402/1‑71 (MQ=255)
                         aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAt                                        >  1:343780/1‑71 (MQ=255)
                          gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                                                                      >  1:627901/1‑40 (MQ=255)
                          gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                                                                      >  1:4757351/1‑40 (MQ=255)
                          gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTgg                                                         >  1:5434049/1‑53 (MQ=255)
                          gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                       <  1:1231880/71‑1 (MQ=255)
                          gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                       >  1:8113/1‑71 (MQ=255)
                          gAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATc                                       <  1:1474192/71‑1 (MQ=255)
                             aTCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCgg                                    >  1:25208/1‑71 (MQ=255)
                               cAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                                                      >  1:3068564/1‑51 (MQ=255)
                                  aggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                                               <  1:920563/54‑1 (MQ=255)
                                    gagaCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAt                             >  1:5038161/4‑71 (MQ=255)
                                    gagaCGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAt                             >  1:618636/4‑71 (MQ=255)
                                    gAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                             <  1:5361754/71‑1 (MQ=255)
                                    gAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                             <  1:332899/71‑1 (MQ=255)
                                       aCGATATTCACCAGGTTGGAGACAATAAGCGTCAGCCt                                                           >  1:1199707/1‑38 (MQ=38)
                                        cGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACg                         <  1:1992297/71‑1 (MQ=255)
                                        cGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACg                         <  1:5620776/71‑1 (MQ=255)
                                        cGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAca                         <  1:2088555/71‑2 (MQ=255)
                                         gATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGc                                                        >  1:580845/1‑39 (MQ=255)
                                         gATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgccgcc                               <  1:1043670/64‑1 (MQ=255)
                                         gATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGa                        <  1:1039141/71‑1 (MQ=255)
                                         gATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGa                        <  1:1032259/71‑1 (MQ=255)
                                          atatTCACCAGGTTGGAGACAATAAGCGGCAGGCTgg                                                         >  1:3356018/1‑37 (MQ=255)
                                            atTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATc                                       >  1:2517566/1‑53 (MQ=255)
                                             ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTAt                                                   >  1:1198090/1‑40 (MQ=255)
                                             ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                            <  1:86279/47‑1 (MQ=255)
                                             ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                                            <  1:5580356/47‑1 (MQ=255)
                                             ttCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                    <  1:2205322/71‑1 (MQ=255)
                                             ttCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGc                    <  1:2890870/71‑1 (MQ=255)
                                             ttCACCAGCTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCAt                             <  1:2379629/62‑1 (MQ=255)
                                               cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGcc                   <  1:4629170/70‑1 (MQ=255)
                                               cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGcc                   <  1:5543577/70‑1 (MQ=255)
                                               cACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGc                    >  1:3401340/1‑69 (MQ=255)
                                               cACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGcc                   <  1:2777360/70‑1 (MQ=255)
                                                  cAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg               <  1:1548326/71‑1 (MQ=255)
                                                  cAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtg               >  1:3791224/1‑69 (MQ=255)
                                                   aGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                             >  1:3866034/1‑56 (MQ=255)
                                                   aGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCAt                             >  1:3513260/1‑56 (MQ=255)
                                                   aGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCa                  >  1:2594433/1‑67 (MQ=255)
                                                    ggTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCg               >  1:4563082/1‑69 (MQ=255)
                                                    ggTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgtc             >  1:2786185/1‑70 (MQ=255)
                                                    ggTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtgtc             >  1:1878442/1‑70 (MQ=255)
                                                     gTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGAcaa                        <  1:4699307/59‑3 (MQ=255)
                                                     gTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAgtg               <  1:204203/68‑3 (MQ=255)
                                                        ggAGACAATAAGCGGCAGGCTGGCGGTATCGGaaa                                             >  1:995273/1‑35 (MQ=255)
                                                         gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCg                                     <  1:4551296/42‑1 (MQ=255)
                                                         gagaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGcc        <  1:2455605/71‑1 (MQ=255)
                                                           gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt      <  1:5668941/71‑1 (MQ=255)
                                                           gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt      <  1:3078189/71‑1 (MQ=255)
                                                           gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt      <  1:2713506/71‑1 (MQ=255)
                                                           gaCAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCtt      <  1:1011863/71‑1 (MQ=255)
                                                            aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATcct                                     <  1:3834837/39‑2 (MQ=21)
                                                            aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg     >  1:2704228/1‑68 (MQ=255)
                                                            aCAATGAGCGGCAGGCTGGCCGTATCTGCAATAAATCCTGCAGCCATGACAAAGGCCAGTGTCGTGccctg     <  1:1373634/71‑4 (MQ=255)
                                                             cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCttt     <  1:2206291/70‑1 (MQ=255)
                                                             cAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCttt     <  1:2081361/70‑1 (MQ=255)
                                                                tAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACt  >  1:3766970/1‑70 (MQ=255)
                                                                 |                                                                    
GCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCATCACGAACGCCAGCGTAGTGCCTTTACT  >  Exported/2596071‑2596204

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