Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,341 A→G 23.5% A419A (GCT→GCC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,3410AG23.5% 167.2 / 52.7 98A419A (GCT→GCCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (35/40);  new base G (15/8);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

Only 100 of 121 total aligned reads displayed.

GTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGA  >  minE/2595272‑2595408
                                                                     |                                                                   
gtgATAAATGCTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGca                                                                     <  1:4087390/70‑1 (MQ=255)
  gATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGcggc                                                                   >  1:4396535/1‑70 (MQ=255)
    tAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagccag                                                                <  1:2410854/71‑1 (MQ=255)
    tAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagcca                                                                 <  1:574231/70‑1 (MQ=255)
     aaaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagccagc                                                               >  1:3958727/1‑71 (MQ=255)
     aaaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagccagc                                                               >  1:3469632/1‑71 (MQ=255)
     aaaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagccagc                                                               >  1:3462899/1‑71 (MQ=255)
      gaTGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGcggccaga                                                               <  1:3367669/69‑2 (MQ=255)
       aTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGt                                                             <  1:5387995/71‑1 (MQ=255)
        tGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTc                                                            <  1:5228961/71‑1 (MQ=255)
         ggTAATGTTGCTCATTTCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGtga                                                           >  1:4247100/1‑69 (MQ=255)
         ggTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGtga                                                           >  1:2836887/1‑69 (MQ=255)
         ggTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGtga                                                           >  1:4797135/1‑69 (MQ=255)
            aaTGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTcaca                                                         >  1:1251148/1‑70 (MQ=255)
               gTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagc                                                                   >  1:2083160/1‑57 (MQ=255)
               gTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagc                                                                   >  1:2162843/1‑57 (MQ=255)
               gTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagc                                                                   >  1:4469462/1‑57 (MQ=255)
                ttGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGGGCAGCCAGCGTCACAAACAg                                                    >  1:5265201/1‑71 (MQ=255)
                ttGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAg                                                    >  1:287879/1‑71 (MQ=255)
                ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACa                                                     <  1:3899196/70‑1 (MQ=255)
                ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACa                                                     <  1:4281743/70‑1 (MQ=255)
                ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACa                                                     <  1:1086304/70‑1 (MQ=255)
                ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACa                                                     <  1:5067645/70‑1 (MQ=255)
                   cTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGcac                                                 <  1:844464/71‑1 (MQ=255)
                   cTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGcac                                                 <  1:2572734/71‑1 (MQ=255)
                   cTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGca                                                  <  1:1061871/70‑1 (MQ=255)
                      atatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACa                                                     <  1:5688686/64‑1 (MQ=255)
                      atatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACa                                                     <  1:2809995/64‑1 (MQ=255)
                      atatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACa                                                     <  1:2459377/64‑1 (MQ=255)
                      atatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcagg                                              >  1:4622462/1‑71 (MQ=255)
                       tatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGcggcca                                                                 <  1:4603532/51‑1 (MQ=255)
                       tatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGcggcca                                                                 <  1:4323806/51‑1 (MQ=255)
                        atCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGcggcca                                                                 >  1:2786411/1‑50 (MQ=255)
                         tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTc                                                            >  1:2403905/1‑54 (MQ=255)
                         tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag                                           >  1:4892038/1‑71 (MQ=255)
                         tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggca                                            <  1:1831062/70‑1 (MQ=255)
                         tCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGcggccaga                                                               >  1:124132/1‑50 (MQ=255)
                           aGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAggg                                         <  1:977500/71‑1 (MQ=255)
                             tATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGca                                                  <  1:3053832/60‑1 (MQ=255)
                               tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTcaca                                                         <  1:2242025/51‑1 (MQ=255)
                               tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACaaa                                                       >  1:4709551/1‑53 (MQ=255)
                               tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAg                                                    >  1:705235/1‑56 (MQ=255)
                               tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAg                                                    >  1:3134278/1‑56 (MQ=255)
                               tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag                                           <  1:3355540/65‑1 (MQ=255)
                               tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag                                           <  1:4290259/65‑1 (MQ=255)
                               tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCAt                                     >  1:65953/1‑71 (MQ=255)
                               tctaATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCAccggcaca                                          >  1:3592595/1‑64 (MQ=255)
                               tctaATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCAt                                     >  1:1518678/1‑71 (MQ=255)
                               tctaATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCAt                                     >  1:1160607/1‑71 (MQ=255)
                                  aaTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGa                                   >  1:1664135/2‑70 (MQ=255)
                                   aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGac                                  <  1:3521973/70‑2 (MQ=255)
                                   aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGac                                  >  1:1737157/1‑69 (MQ=255)
                                   aTTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGAt                                  <  1:4851837/70‑1 (MQ=255)
                                   aTTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATa                                 <  1:5540100/71‑1 (MQ=255)
                                    ttACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACa                                                     <  1:4860878/50‑1 (MQ=255)
                                    ttACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATaa                                <  1:308626/71‑1 (MQ=255)
                                     tACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAAt                               >  1:5577178/1‑71 (MQ=255)
                                      aCAAAGTGAAAGAGAGACGTAGCGCCAGCTCAGCCAGCGTCACAAAc                                                      >  1:1909327/1‑47 (MQ=255)
                                      aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGa                                   >  1:4675637/1‑66 (MQ=255)
                                      aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATc                              >  1:2094823/1‑71 (MQ=255)
                                       cAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAAtgc                             <  1:3516961/71‑3 (MQ=255)
                                       cAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATcc                             >  1:1603561/1‑71 (MQ=255)
                                       cAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATcc                             <  1:543292/71‑1 (MQ=255)
                                        aaaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGcac                                                 <  1:1261978/50‑1 (MQ=255)
                                        aaaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGcac                                                 <  1:4110322/50‑1 (MQ=255)
                                         aaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcagg                                          >  1:3923040/1‑56 (MQ=255)
                                         aaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcagg                                          >  1:2154353/1‑56 (MQ=255)
                                          aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                          <  1:3202795/71‑1 (MQ=255)
                                          aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                          <  1:4786190/71‑1 (MQ=255)
                                          aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                          <  1:1826622/71‑1 (MQ=255)
                                           gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                          >  1:1153276/1‑70 (MQ=255)
                                           gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                          >  1:5256164/1‑70 (MQ=255)
                                           gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGt                         <  1:969119/71‑1 (MQ=255)
                                                agagagACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATa                                 >  1:2000060/1‑58 (MQ=255)
                                                   gagaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                  >  1:1827920/1‑70 (MQ=255)
                                                   gagaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                  >  1:825321/1‑70 (MQ=255)
                                                   gagaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGt                 <  1:163546/71‑1 (MQ=255)
                                                   gagaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGt                 <  1:1696360/71‑1 (MQ=255)
                                                   gagaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGt                 <  1:5453997/71‑1 (MQ=255)
                                                    agaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTa                <  1:1321587/71‑1 (MQ=255)
                                                    agaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTa                <  1:5410098/71‑1 (MQ=255)
                                                         tAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggc                                             <  1:1791825/37‑1 (MQ=255)
                                                          cacgcCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATcc          >  1:5377867/3‑71 (MQ=255)
                                                          aGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAAt                               >  1:1784388/1‑50 (MQ=255)
                                                          aGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAAt                               >  1:993475/1‑50 (MQ=255)
                                                          aGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAAt                               >  1:5171048/1‑50 (MQ=255)
                                                           gcgcCCGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCAt                                     >  1:1979192/1‑43 (MQ=255)
                                                           gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATccc         >  1:5591103/1‑71 (MQ=255)
                                                           gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATccc         >  1:5127962/1‑71 (MQ=255)
                                                           acgcCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATccc         >  1:2590656/2‑71 (MQ=255)
                                                           acgcCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATccc         >  1:3025161/2‑71 (MQ=255)
                                                            cgcCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATcccc        <  1:3471919/71‑1 (MQ=255)
                                                              ccAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAg      >  1:4473558/1‑71 (MQ=255)
                                                                aGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCt    <  1:130726/71‑1 (MQ=255)
                                                                 gcgcGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATcccc        <  1:4704743/66‑1 (MQ=255)
                                                                 gcgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCt    >  1:2839713/1‑70 (MQ=255)
                                                                  cgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa  <  1:513441/71‑1 (MQ=255)
                                                                  cgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa  <  1:29933/71‑1 (MQ=255)
                                                                    cagccagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCt                           <  1:2041747/44‑1 (MQ=255)
                                                                    cagccagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa  >  1:5070701/1‑69 (MQ=255)
                                                                     |                                                                   
GTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGA  >  minE/2595272‑2595408

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: