Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,306 C→A 20.7% intergenic (‑11/+2) arsC ← / ← arsB arsenate reductase/arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,3060CA20.7% 216.6 / 27.5 77intergenic (‑11/+2)arsC/arsBarsenate reductase/arsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base A (10/6);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

TCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGA  >  minE/2595237‑2595375
                                                                     |                                                                     
tcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa                                                                      >  1:3187360/1‑71 (MQ=255)
  cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAt                                                                     >  1:213189/1‑70 (MQ=255)
     cGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACa                                                                 <  1:2099097/71‑1 (MQ=255)
      gTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa                                                                 >  1:1597207/1‑70 (MQ=255)
       tATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa                                                                 >  1:2672352/1‑69 (MQ=255)
           aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACaa                                                                >  1:2202706/1‑66 (MQ=255)
           aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg                                                            >  1:1639357/1‑70 (MQ=255)
           aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg                                                            <  1:2030771/70‑1 (MQ=255)
            cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaa                                                          <  1:2384280/71‑1 (MQ=255)
            cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaa                                                          <  1:2880073/71‑1 (MQ=255)
            cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa                                                           <  1:2113878/70‑1 (MQ=255)
            cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa                                                           <  1:2634848/70‑1 (MQ=255)
             gcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa                                                         <  1:2336324/71‑1 (MQ=255)
              cgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa                                                           >  1:745275/1‑68 (MQ=255)
               gACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                                                       <  1:77634/71‑1 (MQ=255)
               gACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                                                       <  1:885400/71‑1 (MQ=255)
                aCGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagag                                                      <  1:3456575/71‑1 (MQ=255)
                aCGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagag                                                      <  1:3804033/71‑1 (MQ=255)
                   tGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                                                       <  1:427661/67‑1 (MQ=255)
                     ccGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagaga                                                     >  1:3346957/1‑67 (MQ=255)
                     ccGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagaga                                                     <  1:3826643/67‑1 (MQ=255)
                     ccGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACg                                                 <  1:2828871/71‑1 (MQ=255)
                      cGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAag                                                        >  1:1790879/1‑63 (MQ=255)
                      cGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAag                                                        >  1:2236511/1‑63 (MQ=255)
                        cAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACagagagc                                                   >  1:1754559/1‑65 (MQ=255)
                        cAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa                                                         >  1:1723584/1‑60 (MQ=255)
                        cAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa                                                         >  1:2392040/1‑60 (MQ=255)
                            ccGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGcc                                          <  1:1336311/71‑1 (MQ=255)
                             cGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAgcgc                                           >  1:3884315/1‑69 (MQ=255)
                                  tgtgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAag                                                        >  1:1906276/1‑51 (MQ=255)
                                  tgtgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAgagaga                                                   <  1:1168801/56‑1 (MQ=255)
                                    tgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa                                                                 <  1:3507236/40‑1 (MQ=255)
                                    tgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGcc                                          <  1:2358810/63‑1 (MQ=255)
                                    tgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcag                                   <  1:1479919/70‑1 (MQ=255)
                                     gATAAATGGTAATGTTGCTCATATCAGTATCTCATTACaa                                                                >  1:829721/1‑40 (MQ=255)
                                      aTAGATGGTAATGTTGCTCATATCAGTATCTAATTACaaa                                                               >  1:3230023/1‑40 (MQ=38)
                                      aTAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCa                                         >  1:3451772/1‑62 (MQ=255)
                                      aTAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagc                                  >  1:31842/1‑69 (MQ=255)
                                       tAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa                                                           <  1:135342/43‑1 (MQ=255)
                                        agaTGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGcggcca                                >  1:1307160/3‑69 (MQ=255)
                                         gaTGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGcggccagag                             >  1:417408/2‑69 (MQ=255)
                                         aaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCg                             >  1:1405204/1‑71 (MQ=255)
                                         aaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCg                             >  1:2351324/1‑71 (MQ=255)
                                         aaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCg                             <  1:3550584/71‑1 (MQ=255)
                                            ggTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTca                          >  1:4039555/1‑71 (MQ=255)
                                            ggTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTca                          >  1:4030391/1‑71 (MQ=255)
                                             gTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAgcgc                                     <  1:831615/59‑1 (MQ=255)
                                              tAATGTTGCTCATATCAGTATCTAATTACAAAGTGACaga                                                       >  1:3896011/1‑40 (MQ=38)
                                              tAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTgaca                        <  1:873995/71‑1 (MQ=255)
                                              tAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTgaca                        <  1:1124709/71‑1 (MQ=255)
                                               aaTGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAg                                              >  1:1842547/1‑48 (MQ=255)
                                                aTGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGt                                                <  1:607554/45‑1 (MQ=255)
                                                 tGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGCGACGTAGCGCCAGCGCAGCCAGCGTCACAAAc                     <  1:666367/71‑1 (MQ=255)
                                                 tGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAgcgc                                           <  1:786112/49‑1 (MQ=255)
                                                 tGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCCGCGCGGCCAGAGt                            <  1:4037393/64‑1 (MQ=255)
                                                  gTTGCTCATATCAGTATCTCATTACAAAGTGAAagaga                                                     >  1:1130018/1‑38 (MQ=255)
                                                  gTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACaaa                      >  1:1021749/1‑69 (MQ=255)
                                                   ttGCTCATATCAGTATCTCATTACAAAGTGAAagaga                                                     >  1:2393771/1‑37 (MQ=255)
                                                   ttGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAg                   <  1:2015386/71‑1 (MQ=255)
                                                   ttGCTCATATCAGTATCTAATTACAAAGTGACaga                                                       >  1:836820/1‑35 (MQ=37)
                                                   ttGCTCATATCAGTATCTAATTACAAAGTGACaga                                                       >  1:2224310/1‑35 (MQ=37)
                                                   ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCgc                                                >  1:2944192/1‑41 (MQ=37)
                                                   ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACaaa                      >  1:2576625/1‑68 (MQ=255)
                                                   ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACa                    <  1:3694442/70‑1 (MQ=255)
                                                   ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAg                   >  1:2906983/1‑71 (MQ=255)
                                                         atatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcagg             >  1:3357550/1‑71 (MQ=255)
                                                          tatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggc            <  1:3069705/71‑1 (MQ=255)
                                                          tatCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcagg             >  1:2013140/1‑70 (MQ=255)
                                                            tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag          >  1:2561007/1‑71 (MQ=255)
                                                            tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag          >  1:1528910/1‑71 (MQ=255)
                                                            tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag          >  1:1521584/1‑71 (MQ=255)
                                                            tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag          <  1:1319623/71‑1 (MQ=255)
                                                              aGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTc                           >  1:3974419/1‑52 (MQ=255)
                                                               gTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagccagc                              <  1:2570829/48‑1 (MQ=255)
                                                                tATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTc      <  1:3517368/71‑1 (MQ=255)
                                                                tATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTc      <  1:1207333/71‑1 (MQ=255)
                                                                 aTCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCa     >  1:1475760/1‑71 (MQ=255)
                                                                  tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGcaca                                    <  1:21487/39‑1 (MQ=255)
                                                                  tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagc                                  >  1:4059495/1‑41 (MQ=255)
                                                                  tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGca                                    <  1:3236470/39‑1 (MQ=255)
                                                                  tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGt                            <  1:289700/47‑1 (MQ=255)
                                                                  tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGt                            <  1:3392721/47‑1 (MQ=255)
                                                                  tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTcaca                        <  1:2182550/51‑1 (MQ=255)
                                                                  tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTcaca                        <  1:94865/51‑1 (MQ=255)
                                                                  tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTc                           >  1:961741/1‑48 (MQ=255)
                                                                  tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCAt    >  1:2220851/1‑71 (MQ=255)
                                                                  tctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCAt    >  1:1229331/1‑71 (MQ=255)
                                                                  tctaATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCAt    >  1:1143955/1‑71 (MQ=255)
                                                                   ctcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACaaa                      >  1:581078/1‑52 (MQ=255)
                                                                    tcATTACATAGTGAAAGAGATACGTAGCGCCAGCGCAGCCAGCGTcaca                        <  1:1696701/49‑1 (MQ=38)
                                                                    tcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTcaca                        <  1:1323027/49‑1 (MQ=255)
                                                                    tcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTcaca                        <  1:1005448/49‑1 (MQ=255)
                                                                    tcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTcaca                        <  1:1005006/49‑1 (MQ=255)
                                                                    tcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGa  >  1:1610239/1‑71 (MQ=255)
                                                                    tcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGa  >  1:87160/1‑71 (MQ=255)
                                                                     cATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAggg        >  1:3077170/1‑64 (MQ=255)
                                                                     cATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATg   >  1:3475973/1‑69 (MQ=255)
                                                                     cATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATg   >  1:3463881/1‑69 (MQ=255)
                                                                     |                                                                     
TCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGA  >  minE/2595237‑2595375

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: