Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,347 C→A 23.4% T417T (ACG→ACT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,3470CA23.4% 218.6 / 43.1 98T417T (ACG→ACTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (35/40);  new base A (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 = 9.97e-01

Only 100 of 119 total aligned reads displayed.

AAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAT  >  minE/2595277‑2595417
                                                                      |                                                                      
aaaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagccagc                                                                        >  1:3958727/1‑71 (MQ=255)
aaaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagccagc                                                                        >  1:3462899/1‑71 (MQ=255)
aaaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagccagc                                                                        >  1:3469632/1‑71 (MQ=255)
  aTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGt                                                                      >  1:3468783/1‑71 (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)
           ttGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGGGCAGCCAGCGTCACAAACAg                                                             >  1:5265201/1‑71 (MQ=255)
           ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACa                                                              <  1:4281743/70‑1 (MQ=255)
           ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACa                                                              <  1:5067645/70‑1 (MQ=255)
           ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACa                                                              <  1:1086304/70‑1 (MQ=255)
              cTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGcac                                                          <  1:844464/71‑1 (MQ=255)
              cTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGca                                                           <  1:1061871/70‑1 (MQ=255)
                 atatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACa                                                              <  1:2459377/64‑1 (MQ=255)
                 atatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACa                                                              <  1:2809995/64‑1 (MQ=255)
                 atatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACa                                                              <  1:5688686/64‑1 (MQ=255)
                 atatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACa                                                              <  1:541780/64‑1 (MQ=255)
                 atatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcagg                                                       >  1:4622462/1‑71 (MQ=255)
                   atCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACa                                                              >  1:2065172/1‑62 (MQ=255)
                    tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTc                                                                     >  1:2403905/1‑54 (MQ=255)
                    tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag                                                    >  1:4892038/1‑71 (MQ=255)
                    tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggca                                                     <  1:1859326/70‑1 (MQ=255)
                    tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggca                                                     <  1:1831062/70‑1 (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)
                          tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag                                                    <  1:4290259/65‑1 (MQ=255)
                          tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag                                                    <  1:3355540/65‑1 (MQ=255)
                          tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCAt                                              >  1:65953/1‑71 (MQ=255)
                          tctaATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCAccggcaca                                                   >  1:3592595/1‑64 (MQ=255)
                          tctaATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCAt                                              >  1:1160607/1‑71 (MQ=255)
                          tctaATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCAt                                              >  1:1518678/1‑71 (MQ=255)
                             aaTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGa                                            >  1:1664135/2‑70 (MQ=255)
                              aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGac                                           >  1:1737157/1‑69 (MQ=255)
                              aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGac                                           >  1:4377325/1‑69 (MQ=255)
                              aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGac                                           <  1:3521973/70‑2 (MQ=255)
                              aTTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATg                                             <  1:2284194/68‑1 (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)
                                 aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATc                                       >  1:2982070/1‑71 (MQ=255)
                                  cAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAAtgc                                      <  1:3516961/71‑3 (MQ=255)
                                  cAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATcc                                      <  1:543292/71‑1 (MQ=255)
                                  cAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATcc                                      >  1:1603561/1‑71 (MQ=255)
                                   aaaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGcac                                                          <  1:4110322/50‑1 (MQ=255)
                                   aaaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGcac                                                          <  1:1261978/50‑1 (MQ=255)
                                    aaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcagg                                                   >  1:3923040/1‑56 (MQ=255)
                                     aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                   <  1:3202795/71‑1 (MQ=255)
                                     aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                   <  1:1826622/71‑1 (MQ=255)
                                     aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                   <  1:4786190/71‑1 (MQ=255)
                                      gTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATgccgg                                   >  1:5242204/1‑68 (MQ=255)
                                      gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                   >  1:1153276/1‑70 (MQ=255)
                                      gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGt                                  <  1:969119/71‑1 (MQ=255)
                                              gagaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                           >  1:1827920/1‑70 (MQ=255)
                                              gagaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGt                          <  1:1696360/71‑1 (MQ=255)
                                               agaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTa                         <  1:5410098/71‑1 (MQ=255)
                                               agaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTa                         <  1:1321587/71‑1 (MQ=255)
                                                     cacgcCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATcc                   >  1:5377867/3‑71 (MQ=255)
                                                     aGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAAt                                        >  1:5171048/1‑50 (MQ=255)
                                                     aGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAAt                                        >  1:993475/1‑50 (MQ=255)
                                                     aGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAAt                                        >  1:1784388/1‑50 (MQ=255)
                                                      gcgcCCGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCAt                                              >  1:1979192/1‑43 (MQ=255)
                                                      gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCAt                                              >  1:83857/1‑43 (MQ=255)
                                                      gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGa                                            >  1:3131452/1‑45 (MQ=255)
                                                      gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATccc                  >  1:5127962/1‑71 (MQ=255)
                                                      gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATccc                  >  1:5591103/1‑71 (MQ=255)
                                                      acgcCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATccc                  >  1:2590656/2‑71 (MQ=255)
                                                      acgcCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATccc                  >  1:3025161/2‑71 (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:29933/71‑1 (MQ=255)
                                                             cgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa           <  1:513441/71‑1 (MQ=255)
                                                             cgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa           <  1:1243161/71‑1 (MQ=255)
                                                               cagccagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCt                                    <  1:2041747/44‑1 (MQ=255)
                                                               cagccagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa           >  1:3649646/1‑69 (MQ=255)
                                                                ggccagaGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAttg        >  1:3870623/2‑69 (MQ=255)
                                                                 gccagaGTGACAAACAGCACCGGCACAGTCATGac                                           <  1:584133/35‑2 (MQ=12)
                                                                 gccagaGTGACAAACAGCACCGGCACAGTCATGac                                           <  1:2503828/35‑2 (MQ=12)
                                                                 gccagaGTGACAAACAGCACCGGCACAGTCATGac                                           <  1:257711/35‑2 (MQ=12)
                                                                 gccagaGTGACAAACAGCACCGGCACAGTCATGac                                           <  1:2147441/35‑2 (MQ=12)
                                                                   cagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGCTAGTCa     <  1:5555/71‑1 (MQ=255)
                                                                   cagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa     <  1:214073/71‑1 (MQ=255)
                                                                    agcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa     <  1:396040/70‑1 (MQ=255)
                                                                    agcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa     >  1:5646237/1‑70 (MQ=255)
                                                                    agcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa     <  1:1674147/70‑1 (MQ=255)
                                                                     gcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCata   <  1:4893225/71‑1 (MQ=255)
                                                                      cGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCatat  <  1:614030/71‑1 (MQ=255)
                                                                      cGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCatat  <  1:4148128/71‑1 (MQ=255)
                                                                      |                                                                      
AAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAT  >  minE/2595277‑2595417

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: