Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,575 T→C 10.5% V341V (GTA→GTG arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5750TC10.5% 260.8 / 19.0 95V341V (GTA→GTGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base T (50/35);  new base C (7/3);  total (57/38)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGC  >  minE/2595506‑2595643
                                                                     |                                                                    
ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                                                                     >  1:843990/1‑71 (MQ=255)
ggCATAAACCATCGCTTCTTTGATAACGCCAGAAGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa                                                                     >  1:2676502/1‑71 (MQ=255)
  cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                                                   <  1:668832/71‑1 (MQ=255)
  cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                                                   <  1:2578237/71‑1 (MQ=255)
  cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                                                   <  1:2447130/71‑1 (MQ=255)
  cATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                                                   <  1:2098857/71‑1 (MQ=255)
    tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAcc                                                                   <  1:3149685/69‑1 (MQ=255)
    tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                                                                 >  1:3616433/1‑71 (MQ=255)
    tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGt                                                                 >  1:2218321/1‑71 (MQ=255)
     aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc                                                                >  1:12138/1‑71 (MQ=255)
      aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                               >  1:81717/1‑71 (MQ=255)
      aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                               <  1:310282/71‑1 (MQ=255)
      aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                               >  1:3678812/1‑71 (MQ=255)
      aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                               >  1:3349748/1‑71 (MQ=255)
      aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                               >  1:3892744/1‑71 (MQ=255)
      aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                               >  1:3014557/1‑71 (MQ=255)
      aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGAACAACGCGCCAACCAGTACCGTCg                                                               <  1:2289317/71‑1 (MQ=255)
         cATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCa                                                            <  1:2187272/71‑1 (MQ=255)
          aTCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCCGCat                                                           <  1:1576686/71‑1 (MQ=255)
           tCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCata                                                          <  1:3855249/71‑1 (MQ=255)
            cGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCatat                                                         >  1:1477518/1‑71 (MQ=255)
             gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt                                                        >  1:3478741/1‑71 (MQ=255)
             gCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATAtt                                                        >  1:2203847/1‑71 (MQ=255)
              cttcttTGATAACGCCAGATGCCGTGCTTCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTg                                                       <  1:1520909/71‑1 (MQ=255)
              cttcttTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTg                                                       <  1:1238127/71‑1 (MQ=255)
                   ttGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCa                                                   >  1:1295209/1‑70 (MQ=255)
                    tGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCat                                                  <  1:1604755/70‑1 (MQ=255)
                      aTAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataa                                                <  1:976142/70‑1 (MQ=255)
                      aTAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                                                 <  1:2527304/69‑1 (MQ=255)
                        aaCGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGt                                                      >  1:2260972/1‑62 (MQ=255)
                         aCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCatcac                                               <  1:3165685/68‑4 (MQ=255)
                         aCGCCAGTCGCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCatcac                                               <  1:2500424/68‑4 (MQ=255)
                          cGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaa                                           <  1:477469/71‑1 (MQ=255)
                           gCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCat                                                           <  1:424490/54‑1 (MQ=255)
                            ccAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga                                         <  1:1733267/71‑1 (MQ=255)
                             cAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga                                         <  1:57482/70‑1 (MQ=255)
                               gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag                                      <  1:897302/71‑1 (MQ=255)
                               gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag                                      <  1:902455/71‑1 (MQ=255)
                               gATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag                                      <  1:273273/71‑1 (MQ=255)
                                tcgCCGTACTCCCGTCAATCGACAGCGCACCAATCATCACCGTCGGCATATTGTTAATAACTGAtgaga                                       >  1:1648825/3‑69 (MQ=255)
                                aTGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATagaaga                                         <  1:2344369/67‑1 (MQ=255)
                                 tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg                                                               >  1:2814061/1‑44 (MQ=255)
                                 tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagaga                                       >  1:265722/1‑68 (MQ=255)
                                 tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGa                                    <  1:3914127/71‑1 (MQ=255)
                                 tGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGCCATATTGTTCATAATAGAAGAGAgg                                     >  1:3063473/1‑70 (MQ=255)
                                 cgCCGTACTCCCGTCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGAtgagaga                                     >  1:301409/2‑69 (MQ=255)
                                  gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAgg                                     >  1:2012762/1‑69 (MQ=255)
                                  gCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGa                                    <  1:2279910/70‑1 (MQ=255)
                                     gTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCata                                                          >  1:1538798/1‑45 (MQ=255)
                                      tgctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataata                                              >  1:1029990/1‑56 (MQ=255)
                                        ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                                                 >  1:3794525/1‑51 (MQ=255)
                                        ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt                             <  1:2412370/71‑1 (MQ=255)
                                        ctgcCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt                             >  1:1287844/1‑71 (MQ=255)
                                             aTCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt                             <  1:3900776/66‑1 (MQ=255)
                                              tCAATGGACAACGCGCCAACCAGTACCGTCGGGATATTGTTCata                                                 >  1:2555150/1‑45 (MQ=255)
                                              tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCata                                                 >  1:393022/1‑45 (MQ=255)
                                              tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagagag                                      <  1:2739145/56‑1 (MQ=255)
                                              tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGa                           >  1:3784368/1‑67 (MQ=255)
                                              tCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAgg                         <  1:2350959/69‑1 (MQ=255)
                                              tCAATCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTc                                                    >  1:3653462/1‑42 (MQ=19)
                                               cAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGaa                                   >  1:4045008/1‑58 (MQ=255)
                                               cAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGaa                                   >  1:1597676/1‑58 (MQ=255)
                                               cAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGa                                    >  1:2348929/1‑57 (MQ=255)
                                               cAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGa                                    >  1:2060411/1‑57 (MQ=255)
                                                aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCataa                                                >  1:3300282/1‑44 (MQ=255)
                                                aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCa                                                   >  1:1483563/1‑41 (MQ=255)
                                                aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGa                           >  1:740986/1‑65 (MQ=255)
                                                aaTGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATc                     >  1:3072922/1‑71 (MQ=255)
                                                aaTCGACAGCGCACCAATCAGCACCGTCGGCATAtt                                                        >  1:1908294/1‑36 (MQ=16)
                                                aaTCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGAtgagagaa                                    >  1:698983/1‑54 (MQ=255)
                                                aaTCGACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGAtgagagaa                                    >  1:2820743/1‑54 (MQ=255)
                                                  tcgACAGCGCACCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCg                   >  1:4096321/3‑71 (MQ=255)
                                                  tGGATAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAgg                                 <  1:581612/57‑1 (MQ=255)
                                                  tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCg                   >  1:1849194/1‑71 (MQ=255)
                                                  tGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAAAAGAAGAGAGGAAGGCGGTGAGGAATCCg                   >  1:4035344/1‑71 (MQ=255)
                                                    gACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCg                               >  1:3341765/1‑57 (MQ=255)
                                                    gACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGt                             >  1:515988/1‑59 (MQ=255)
                                                    gACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCgg                  >  1:1366220/1‑70 (MQ=255)
                                                    gACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGt                 >  1:403082/1‑71 (MQ=255)
                                                    gACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGt                 <  1:3599026/71‑1 (MQ=255)
                                                    gACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGt                 >  1:904402/1‑71 (MQ=255)
                                                      cAACGCGCCAACCAGTACCGTCGGCATATTGTTCataata                                              >  1:2970419/1‑40 (MQ=255)
                                                      cAACGCGCCAACCAGTACCGTCGGCATATTGTTCataata                                              >  1:1952911/1‑40 (MQ=255)
                                                      cAACGCGCCAACCAGTACCGTCGGCATATTGTTCataata                                              >  1:1415244/1‑40 (MQ=255)
                                                       aaCGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGc               >  1:3508107/1‑70 (MQ=255)
                                                        gcgcaCCAATCAGCACCGTCGGCATATTGTTCATCACTGATGAGAGAAATGCCGTCAAGAAGCCGGTGCCg             <  1:1548561/70‑1 (MQ=255)
                                                        aCGCGCCAACCAGTACCGTCGGCATATTGTTCata                                                 >  1:925498/1‑35 (MQ=255)
                                                        aCGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagag                                        <  1:3312691/44‑1 (MQ=255)
                                                        aCGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAagag                                        <  1:958/44‑1 (MQ=255)
                                                         cgcgCCAACCAGTTCCGTCGGCATATTGTTCataata                                              >  1:61593/1‑37 (MQ=255)
                                                         cgcgCCAACCAGTACCGTCGGCATATTGTTCataata                                              >  1:797370/1‑37 (MQ=255)
                                                             ccAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAgg                         <  1:1364406/54‑1 (MQ=255)
                                                              cAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAg                          >  1:1403787/1‑52 (MQ=255)
                                                                aCCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAgg                                     <  1:1030686/39‑1 (MQ=255)
                                                                   aGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGc  >  1:3102491/1‑71 (MQ=255)
                                                                     |                                                                    
GGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATATTGTTCATAATAGAAGAGAGGAAGGCGGTGAGGAATCCGGTGCCGAGCGTCGCGGC  >  minE/2595506‑2595643

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: