Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,596,109 A→G 9.6% V163V (GTT→GTC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,596,1090AG 9.6% 199.7 / 10.2 73V163V (GTT→GTCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (28/38);  new base G (3/4);  total (31/42)
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.74e-01

CGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCA  >  minE/2596040‑2596176
                                                                      |                                                                   
cgcgGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                                    >  1:3327534/1‑71 (MQ=255)
cgcgGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGaaa                                                                    >  1:762590/1‑71 (MQ=255)
 gcggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                                   <  1:3281810/71‑1 (MQ=255)
 gcggcgATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                                   <  1:3357981/71‑1 (MQ=255)
      gATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGatat                                                              <  1:1127373/71‑1 (MQ=255)
      gATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGatat                                                              <  1:212430/71‑1 (MQ=255)
       atatCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATAtt                                                             <  1:162433/71‑1 (MQ=255)
        tatCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAAc                                                                   <  1:2776925/64‑1 (MQ=255)
        tatCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATAtt                                                             >  1:1732809/1‑70 (MQ=255)
           ccACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                                          >  1:3179704/1‑70 (MQ=255)
             cccGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                                          <  1:1863816/67‑1 (MQ=255)
              ccGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAgg                                                      <  1:606113/71‑1 (MQ=255)
               cGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACg                                                                  >  1:2931477/1‑58 (MQ=255)
               cGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGata                                                               <  1:465757/61‑1 (MQ=255)
               cGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAg                                                       <  1:3276202/69‑1 (MQ=255)
                 gCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTg                                                   <  1:2632470/71‑1 (MQ=255)
                   aCCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                                 <  1:2427126/71‑1 (MQ=255)
                     cATAACGGAGGCGTACTGCGCAAAGCTCAGACCGAAGAAGTCCGCCGAGACGATATTCACCAGGTTAgaaa                                               >  1:1756316/1‑69 (MQ=255)
                      aTCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                              <  1:855573/71‑1 (MQ=255)
                       tCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGa                                                                 >  1:1560873/1‑51 (MQ=255)
                       tCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                              >  1:228468/1‑70 (MQ=255)
                         aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                                 >  1:3771641/1‑65 (MQ=255)
                         aCCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGga                                                 >  1:2201873/1‑65 (MQ=255)
                          ccGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATa                                          <  1:1929019/71‑1 (MQ=255)
                           cGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATaa                                         <  1:453881/71‑1 (MQ=255)
                            gAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAg                                        >  1:1513322/1‑71 (MQ=255)
                                  tATTCGCTAAAGCCGAGGCCAAAGAAATCAGCGGAAACGat                                                                <  1:644957/41‑1 (MQ=255)
                                  tATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                                          <  1:2234769/47‑1 (MQ=255)
                                   aTTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAgg                                 >  1:1371377/1‑71 (MQ=255)
                                    ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGatat                                                              >  1:1864844/1‑41 (MQ=255)
                                    ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGa                                                                 <  1:1952124/38‑1 (MQ=255)
                                    ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATAtt                                                             <  1:1329411/42‑1 (MQ=255)
                                    ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCAc                                                          <  1:3294984/45‑1 (MQ=255)
                                    ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGt                                                     <  1:2267484/50‑1 (MQ=255)
                                    ttCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGt                                                     <  1:496861/50‑1 (MQ=255)
                                      cgcgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCa                                                           >  1:633443/1‑42 (MQ=255)
                                      cgcgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCa                                                        >  1:3599048/1‑45 (MQ=255)
                                      cgcgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAgg                                                      <  1:389335/47‑1 (MQ=255)
                                      cgcgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAgg                                                      <  1:3072550/47‑1 (MQ=255)
                                      cgcgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTg                              >  1:274491/1‑71 (MQ=255)
                                       gcgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTgg                             >  1:1117304/1‑71 (MQ=255)
                                       gcgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTgg                             <  1:1896383/71‑1 (MQ=255)
                                        cgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                              <  1:2733741/53‑1 (MQ=255)
                                        cgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                              <  1:2492156/53‑1 (MQ=255)
                                        cgAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGAc                                              <  1:2396547/53‑1 (MQ=255)
                                         gAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATAtt                                                             <  1:876438/37‑1 (MQ=255)
                                         gAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCg                           <  1:2789780/71‑1 (MQ=255)
                                           aaaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGc                                       >  1:744467/1‑58 (MQ=255)
                                            aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                          >  1:2877939/1‑70 (MQ=255)
                                            aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                          >  1:2715622/1‑70 (MQ=255)
                                            aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCgg                          >  1:635676/1‑70 (MQ=255)
                                            aaGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGt                         <  1:203151/71‑1 (MQ=255)
                                              gCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTAt                       <  1:2529413/71‑1 (MQ=255)
                                               ccGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATc                      >  1:3677093/1‑71 (MQ=255)
                                                cGAGGCCAAAGAAATCAGCGTAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCg                     <  1:3030136/71‑1 (MQ=255)
                                                  aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTgg                                                  >  1:3900979/1‑40 (MQ=255)
                                                  aGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGc                   <  1:1731717/71‑1 (MQ=255)
                                                    gCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCaa                 >  1:2084086/1‑71 (MQ=255)
                                                     ccAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAAt                <  1:2061612/71‑1 (MQ=255)
                                                     ccAAAGAAATAAGCGGAAACGATATTCACCAGGTTGGAGAc                                              <  1:1967548/41‑1 (MQ=255)
                                                      cgaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTa                        <  1:2746529/60‑1 (MQ=255)
                                                       gaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACaa                                            >  1:239345/2‑41 (MQ=12)
                                                       gaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAAtg                                          <  1:881907/42‑2 (MQ=16)
                                                       gaagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAAtaa              <  1:1599426/70‑3 (MQ=255)
                                                       aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATaa                                         <  1:2609176/44‑1 (MQ=255)
                                                       aaaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATaa                                         <  1:2669200/44‑1 (MQ=255)
                                                        aagaagTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATCTGCAATaaa             <  1:2122001/71‑1 (MQ=255)
                                                        aaGAAATCAGCGGAAACGATATTCACCAGGTTGgaga                                               >  1:1524877/1‑37 (MQ=255)
                                                        aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa             >  1:333578/1‑71 (MQ=255)
                                                        aaGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGaa             >  1:3449102/1‑71 (MQ=255)
                                                         aGAAGTCCGCCGAGACGATATTCACCAGGTTAGAAACAATGAGCGGCAGGCTGGCCGTATct                     >  1:675806/1‑61 (MQ=255)
                                                         aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCg                           >  1:3522283/1‑56 (MQ=255)
                                                         aGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAAt            >  1:3510425/1‑71 (MQ=255)
                                                                  cGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGgccgc    >  1:70191/1‑70 (MQ=255)
                                                                   ggAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGTTATCGGCAATGAATCCGGCCGCCa  >  1:3568100/1‑71 (MQ=255)
                                                                      |                                                                   
CGCGGCGATATCCACCGGCACCATCACCGAGGCGTATTCGCGAAAGCCGAGGCCAAAGAAATCAGCGGAAACGATATTCACCAGGTTGGAGACAATAAGCGGCAGGCTGGCGGTATCGGCAATGAATCCGGCCGCCA  >  minE/2596040‑2596176

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: