Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,350 C→G 7.4% V416V (GTG→GTC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,3500CG 7.4% 290.6 / 9.7 95V416V (GTG→GTCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (45/43);  new base G (4/3);  total (49/46)
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.93e-01

Only 100 of 110 total aligned reads displayed.

GGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTC  >  minE/2595281‑2595419
                                                                     |                                                                     
ggTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTca                                                                      >  1:4030391/1‑71 (MQ=255)
ggTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTca                                                                      >  1:4039555/1‑71 (MQ=255)
  tAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTgaca                                                                    <  1:873995/71‑1 (MQ=255)
  tAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTgaca                                                                    <  1:1124709/71‑1 (MQ=255)
     tGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGCGACGTAGCGCCAGCGCAGCCAGCGTCACAAAc                                                                 <  1:666367/71‑1 (MQ=255)
      gTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACaaa                                                                  >  1:1021749/1‑69 (MQ=255)
       ttGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAg                                                               <  1:2015386/71‑1 (MQ=255)
       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:1521584/1‑71 (MQ=255)
                tCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag                                                      >  1:1528910/1‑71 (MQ=255)
                  aGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTc                                                                       >  1:3974419/1‑52 (MQ=255)
                    tATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTc                                                  <  1:1207333/71‑1 (MQ=255)
                    tATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTc                                                  <  1:3517368/71‑1 (MQ=255)
                     aTCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCa                                                 >  1:1475760/1‑71 (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:1005448/49‑1 (MQ=255)
                        tcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTcaca                                                                    <  1:1005006/49‑1 (MQ=255)
                        tcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTcaca                                                                    <  1:1323027/49‑1 (MQ=255)
                        tcATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGa                                              >  1:87160/1‑71 (MQ=255)
                         cATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAggg                                                    >  1:3077170/1‑64 (MQ=255)
                         cATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATg                                               >  1:3463881/1‑69 (MQ=255)
                         cATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATg                                               >  1:3475973/1‑69 (MQ=255)
                          aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGac                                             <  1:2803491/70‑2 (MQ=255)
                          aTTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCAt                                                <  1:427955/67‑1 (MQ=255)
                          aTTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGAt                                             >  1:1450432/1‑70 (MQ=255)
                          aTTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATa                                            <  1:14396/71‑1 (MQ=255)
                            tACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTcaca                                                                    <  1:3973420/45‑1 (MQ=255)
                            tACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcag                                                      <  1:4028717/59‑1 (MQ=255)
                            tACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAAt                                          >  1:1778553/1‑71 (MQ=255)
                             aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcaggcagg                                                     <  1:1522056/59‑1 (MQ=255)
                             aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGa                                              >  1:2419528/1‑66 (MQ=255)
                             aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATc                                         >  1:2465341/1‑71 (MQ=255)
                             aCAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATc                                         >  1:1923899/1‑71 (MQ=255)
                               aaaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATccc                                       <  1:1333661/71‑1 (MQ=255)
                               aaaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATccc                                       <  1:3610673/71‑1 (MQ=255)
                                aaGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTcaca                                                                    <  1:3864270/41‑1 (MQ=255)
                                 aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCa                                                 >  1:741238/1‑59 (MQ=255)
                                 aGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                     <  1:3023071/71‑1 (MQ=255)
                                  gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                     >  1:1961298/1‑70 (MQ=255)
                                  gTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTg                                     >  1:2408914/1‑70 (MQ=255)
                                    gAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTgatgat                                             <  1:1408265/60‑1 (MQ=255)
                                     aaaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCAcag                                                          >  1:915525/1‑46 (MQ=255)
                                     aaaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATg                                               <  1:3317553/57‑1 (MQ=255)
                                     aaaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGa                                              >  1:1269738/1‑58 (MQ=255)
                                     aaaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACg                                 >  1:2924354/1‑71 (MQ=255)
                                      aaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                <  1:1728026/71‑1 (MQ=255)
                                      aaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                <  1:2684986/71‑1 (MQ=255)
                                      aaGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                <  1:3705634/71‑1 (MQ=255)
                                       agagagACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                >  1:2619540/1‑70 (MQ=255)
                                       agagagACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                >  1:1833044/1‑70 (MQ=255)
                                       agagagACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGa                               <  1:3008348/71‑1 (MQ=255)
                                          gagaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCa                                                 >  1:3803614/1‑50 (MQ=255)
                                           agaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                             >  1:816757/1‑69 (MQ=255)
                                           agaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                             >  1:821700/1‑69 (MQ=255)
                                            gaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                          <  1:3628161/71‑1 (MQ=255)
                                            gaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                          >  1:3399759/1‑71 (MQ=255)
                                            gaCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                          <  1:778476/71‑1 (MQ=255)
                                             aCGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                          >  1:537398/1‑70 (MQ=255)
                                                tAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAata                        >  1:2637145/1‑69 (MQ=255)
                                                 aGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATaa                                           >  1:3272811/1‑49 (MQ=255)
                                                 aGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATcc                                        >  1:2510901/1‑52 (MQ=255)
                                                  gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATccc                    >  1:2700483/1‑71 (MQ=255)
                                                  gcgcCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATcc                     <  1:75801/70‑1 (MQ=255)
                                                   cgcCTGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGt                                                   <  1:2854584/39‑1 (MQ=255)
                                                     ccAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGa                                              <  1:2345536/42‑1 (MQ=255)
                                                     ccAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGa                                              <  1:3652998/42‑1 (MQ=255)
                                                     ccAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                          <  1:631034/62‑1 (MQ=255)
                                                       aGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCt               <  1:3509221/71‑1 (MQ=255)
                                                       aGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATACCCAGCt               <  1:2398680/71‑1 (MQ=255)
                                                         cgcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                          <  1:2804827/58‑1 (MQ=255)
                                                          gcggccagaGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAt            >  1:304701/4‑71 (MQ=255)
                                                          gcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa             >  1:1051765/1‑70 (MQ=255)
                                                             gccagaGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATcccc                   <  1:849417/61‑1 (MQ=255)
                                                              ccagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTg              <  1:3610481/65‑1 (MQ=255)
                                                               cagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa       >  1:542016/1‑71 (MQ=255)
                                                                    tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                <  1:394732/41‑1 (MQ=255)
                                                                    tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACgg                                <  1:1036500/41‑1 (MQ=255)
                                                                    tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGaa                              >  1:3783589/1‑43 (MQ=255)
                                                                    tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                             <  1:84994/44‑1 (MQ=255)
                                                                    tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAg                             <  1:2691768/44‑1 (MQ=255)
                                                                    tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGt                            <  1:1568847/45‑1 (MQ=255)
                                                                    tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGt                            <  1:3844983/45‑1 (MQ=255)
                                                                    tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTaa                          <  1:3121005/47‑1 (MQ=255)
                                                                    tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc  <  1:561201/71‑1 (MQ=255)
                                                                    tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc  <  1:1822834/71‑1 (MQ=255)
                                                                    tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc  >  1:2031137/1‑71 (MQ=255)
                                                                     cacaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGt                            <  1:849909/44‑1 (MQ=255)
                                                                     |                                                                     
GGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTC  >  minE/2595281‑2595419

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: