Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,380 C→G 28.8% G406G (GGG→GGC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,3800CG28.8% 175.8 / 85.9 118G406G (GGG→GGCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (41/43);  new base G (18/16);  total (59/59)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 127 total aligned reads displayed.

ACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGC  >  minE/2595310‑2595449
                                                                      |                                                                     
aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATc                                                                       >  1:2094823/1‑71 (MQ=255)
aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATc                                                                       >  1:2982070/1‑71 (MQ=255)
 cAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATcc                                                                      <  1:543292/71‑1 (MQ=255)
 cAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATcc                                                                      >  1:1603561/1‑71 (MQ=255)
    aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                                                   <  1:1826622/71‑1 (MQ=255)
    aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                                                   <  1:4669639/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)
     gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                                                   >  1:5256164/1‑70 (MQ=255)
             gagaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                                                           >  1:1827920/1‑70 (MQ=255)
             gagaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                                                           >  1:825321/1‑70 (MQ=255)
             gagaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGt                                                          <  1:1696360/71‑1 (MQ=255)
              agaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTa                                                         <  1:1321587/71‑1 (MQ=255)
              agaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTa                                                         <  1:5410098/71‑1 (MQ=255)
                     gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATccc                                                  >  1:5591103/1‑71 (MQ=255)
                     gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATccc                                                  >  1:5127962/1‑71 (MQ=255)
                     acgcCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATccc                                                  >  1:3025161/2‑71 (MQ=255)
                      cgcCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATcccc                                                 <  1:3471919/71‑1 (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)
                              cagccagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                           >  1:3649646/1‑69 (MQ=255)
                               ggccagaGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAttg                                        >  1:3870623/2‑69 (MQ=255)
                                  cagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGCTAGTCa                                     <  1:5555/71‑1 (MQ=255)
                                   agcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa                                     <  1:1674147/70‑1 (MQ=255)
                                   agcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa                                     >  1:5646237/1‑70 (MQ=255)
                                    gcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCata                                   <  1:4893225/71‑1 (MQ=255)
                                     cGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCatat                                  <  1:4148128/71‑1 (MQ=255)
                                     cGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCatat                                  <  1:614030/71‑1 (MQ=255)
                                     cGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCatat                                  <  1:1256084/71‑1 (MQ=255)
                                      gTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCatat                                  <  1:2985680/70‑1 (MQ=255)
                                       tgacaAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttt                                <  1:1203956/69‑2 (MQ=255)
                                       tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                           <  1:4270512/60‑1 (MQ=255)
                                       tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                           <  1:1595552/60‑1 (MQ=255)
                                       tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc                                <  1:1791542/71‑1 (MQ=255)
                                        cacaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTc                                      >  1:5575767/1‑64 (MQ=255)
                                         acaAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttttt                              <  1:3852905/71‑1 (MQ=255)
                                         acaAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAtttt                               <  1:2152288/70‑3 (MQ=255)
                                         acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGaa                                                            >  1:2766781/1‑41 (MQ=255)
                                         acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGaa                                                            >  1:1231302/1‑41 (MQ=255)
                                         acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATcc                                                   >  1:24804/1‑50 (MQ=255)
                                         acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAg                                               >  1:2919375/1‑54 (MQ=255)
                                         acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                              >  1:4352398/1‑71 (MQ=255)
                                         acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                              >  1:4167808/1‑71 (MQ=255)
                                         acaAACAGCACAGGCAGGGTCATGATAATCCCTGAACGGAAGTAATATcc                                                   >  1:5395728/1‑50 (MQ=255)
                                          caAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAtttttc                             <  1:3068514/71‑1 (MQ=255)
                                           aaaCAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATttttct                            >  1:4119604/1‑71 (MQ=255)
                                             aCAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGc                          >  1:4643370/1‑71 (MQ=255)
                                              cagcagAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGc                          >  1:5555800/1‑70 (MQ=255)
                                               aGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTggg                         >  1:74608/1‑68 (MQ=255)
                                               aGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTggg                         >  1:4332141/1‑68 (MQ=255)
                                                  acaGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATtcttct                            >  1:34713/1‑64 (MQ=255)
                                                    aggcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGt                                       <  1:2549641/51‑1 (MQ=255)
                                                    aggcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTa                   >  1:1731795/1‑71 (MQ=255)
                                                     ggcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTAc                  <  1:1130225/71‑1 (MQ=255)
                                                     ggcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTAc                  <  1:4430395/71‑1 (MQ=255)
                                                      gcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACg                 <  1:5252685/71‑1 (MQ=255)
                                                       caggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGt                <  1:4120927/71‑1 (MQ=255)
                                                       caggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGt                <  1:1486229/71‑1 (MQ=255)
                                                        aggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCata                                   >  1:3604778/1‑51 (MQ=255)
                                                        aggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCata                                   >  1:2803758/1‑51 (MQ=255)
                                                        acagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCAc                  <  1:1809577/65‑1 (MQ=255)
                                                        acagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg               >  1:1389504/4‑71 (MQ=255)
                                                        acagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg               <  1:555069/68‑1 (MQ=255)
                                                        acagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg               <  1:5285263/68‑1 (MQ=255)
                                                         gggTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGc              <  1:1524422/71‑1 (MQ=255)
                                                         cagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCa                                     <  1:3979752/46‑1 (MQ=16)
                                                         cagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg               <  1:2483560/68‑1 (MQ=255)
                                                         cagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGc              <  1:725182/69‑1 (MQ=255)
                                                         cagTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGc              <  1:52761/69‑1 (MQ=255)
                                                          tgTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTAc                  <  1:2257996/65‑1 (MQ=255)
                                                          ggTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc                                <  1:3415756/52‑1 (MQ=255)
                                                            tCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGaa                       <  1:5615276/59‑1 (MQ=255)
                                                            tCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGc           >  1:4912920/1‑71 (MQ=255)
                                                            tCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttt                                <  1:4769993/50‑2 (MQ=255)
                                                            tCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAg                     >  1:2554648/1‑61 (MQ=255)
                                                            tCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAg                     >  1:774505/1‑61 (MQ=255)
                                                              aTGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTc                                      >  1:4134779/1‑42 (MQ=255)
                                                              aTGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTc                                      >  1:431553/1‑42 (MQ=255)
                                                              aTGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa         <  1:5321996/71‑1 (MQ=255)
                                                              aTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg               >  1:833614/1‑65 (MQ=255)
                                                              aTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCa         <  1:203005/71‑1 (MQ=255)
                                                               tGATAATCCCTGTACGGAAGTAATATCCCCAGCTg                                            <  1:3026791/35‑1 (MQ=255)
                                                               tGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCata                                   <  1:4453964/44‑1 (MQ=255)
                                                                 aTAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                              >  1:2099194/1‑47 (MQ=255)
                                                                 aTAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                              >  1:3303398/1‑47 (MQ=255)
                                                                  tAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGt                                       >  1:700295/1‑37 (MQ=255)
                                                                  tAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAg                     >  1:5424356/1‑55 (MQ=255)
                                                                  tAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAg                     >  1:4243826/1‑55 (MQ=255)
                                                                  tAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAg                     >  1:470386/1‑55 (MQ=255)
                                                                   aaTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTg                           >  1:2184508/1‑48 (MQ=19)
                                                                     tCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagc  >  1:4444805/1‑71 (MQ=255)
                                                                      cccTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta    >  1:5274802/1‑68 (MQ=255)
                                                                      cccTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta    >  1:2801405/1‑68 (MQ=255)
                                                                      cccTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta    >  1:4924942/1‑68 (MQ=255)
                                                                      |                                                                     
ACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGC  >  minE/2595310‑2595449

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: