Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,749 G→A 14.6% V283V (GTC→GTT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,7490GA14.6% 237.0 / 27.2 89V283V (GTC→GTTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (41/35);  new base A (6/7);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

AGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACA  >  minE/2595680‑2595818
                                                                     |                                                                     
aGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGa                                                                      >  1:2427970/1‑71 (MQ=255)
 gAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGa                                                                      >  1:1155118/1‑70 (MQ=255)
  aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGa                                                                      >  1:3267023/1‑69 (MQ=255)
  aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGAc                                                                     >  1:510221/1‑70 (MQ=255)
  aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                                    <  1:909375/71‑1 (MQ=255)
  aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                                    <  1:1492527/71‑1 (MQ=255)
   aaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                                                   <  1:2321938/71‑1 (MQ=255)
         ttCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                              <  1:2306123/70‑1 (MQ=255)
          tCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                              >  1:740359/1‑69 (MQ=255)
             gTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggg                                                         >  1:3618205/1‑71 (MQ=255)
              ttAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                              >  1:628440/1‑65 (MQ=255)
                aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAAt                                                                  >  1:572633/1‑59 (MQ=255)
                aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAAt                                                                  >  1:3752041/1‑59 (MQ=255)
                       cATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCAccag                                                 <  1:677764/69‑3 (MQ=255)
                         ttGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                                    <  1:2026522/48‑1 (MQ=255)
                         ttGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCt                                                                <  1:180973/52‑1 (MQ=255)
                         ttGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                              <  1:1668365/54‑1 (MQ=255)
                         ttGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCgcgc                                               <  1:70289/69‑1 (MQ=255)
                          tGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAg                                             >  1:3425185/1‑70 (MQ=255)
                           gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATg                                                                       >  1:1397900/1‑43 (MQ=255)
                           gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGa                                                                      >  1:3033549/1‑44 (MQ=255)
                           gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCa                                                            >  1:2406516/1‑54 (MQ=255)
                           gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCa                                                            >  1:2359650/1‑54 (MQ=255)
                              cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                        <  1:2898541/71‑1 (MQ=255)
                              cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                        <  1:325150/71‑1 (MQ=255)
                              cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                        <  1:3630922/71‑1 (MQ=255)
                              cAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                                        <  1:2720411/71‑1 (MQ=255)
                              cAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                                        <  1:3579096/71‑1 (MQ=255)
                               aGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACttt                                       <  1:2918027/71‑1 (MQ=255)
                               aGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                                        <  1:3526430/70‑1 (MQ=255)
                                ggCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCgcgc                                               <  1:3944700/62‑1 (MQ=255)
                                 gCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGa                                           >  1:37375/1‑65 (MQ=255)
                                   cATAAACCACCAGATACATGCCGAGCGAGAAGATGa                                                                      >  1:1662047/1‑36 (MQ=255)
                                   cATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                                   <  1:474058/71‑1 (MQ=255)
                                   cATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                                   >  1:161659/1‑71 (MQ=255)
                                   cATAAACAACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAgg                                                          >  1:358062/1‑48 (MQ=255)
                                    aTAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                  <  1:351432/71‑1 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                                                   <  1:611245/37‑1 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                                                   <  1:1701324/37‑1 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                                    <  1:2303278/36‑1 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                                    <  1:2568257/36‑1 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                                    <  1:3172504/36‑1 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCt                                                                >  1:2356609/1‑40 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCt                                                                >  1:4033837/1‑40 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggggg                                                       <  1:1820424/49‑1 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAgggg                                                        <  1:299629/48‑1 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAgg                                                          <  1:807049/46‑1 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                  <  1:1543030/70‑1 (MQ=255)
                                     tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                  >  1:4079911/1‑70 (MQ=255)
                                      aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTg                                                               >  1:2998659/1‑40 (MQ=255)
                                      aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTg                                                               >  1:3003348/1‑40 (MQ=255)
                                      aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                              >  1:1056650/1‑41 (MQ=255)
                                      aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                              >  1:307417/1‑41 (MQ=255)
                                      aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa                                >  1:1050198/1‑71 (MQ=255)
                                      aaaCCACCAGATACATGCCCAGCGAGAAGATGACAATCTGc                                                              >  1:1192800/1‑41 (MQ=255)
                                       aaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCa                                                     >  1:2956606/1‑49 (MQ=255)
                                        accaccAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTAtt                              <  1:3981654/71‑1 (MQ=255)
                                         ccaccaGATACATGCCGAGCGAGAAGATGACAATCTGCCAggg                                                         >  1:4042943/1‑43 (MQ=255)
                                         ccaccaGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                             <  1:3118844/71‑1 (MQ=255)
                                          caccaGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                            >  1:1221211/1‑71 (MQ=255)
                                             caGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                         <  1:873797/71‑1 (MQ=255)
                                             caGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                         >  1:3218807/1‑71 (MQ=255)
                                              aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa                                >  1:2292813/1‑63 (MQ=255)
                                               gATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                       >  1:479754/1‑71 (MQ=255)
                                                aTACATGCCGAGCGAGAAGATGACAATCTGCCAggggg                                                       <  1:19379/38‑1 (MQ=255)
                                                aTACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCa                                                     <  1:1026382/40‑1 (MQ=255)
                                                 tACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                            >  1:1446015/1‑64 (MQ=255)
                                                 tACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAAt                      <  1:3560186/70‑1 (MQ=255)
                                                 tACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATg                     <  1:1272004/71‑1 (MQ=255)
                                                  aCATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGa                    >  1:1143159/1‑71 (MQ=255)
                                                   cATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                         <  1:1546797/65‑1 (MQ=255)
                                                   cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAg                                             >  1:1187780/1‑45 (MQ=255)
                                                   cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAg                                             >  1:2411812/1‑45 (MQ=255)
                                                   cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                        >  1:3468886/1‑50 (MQ=255)
                                                   cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                        >  1:2056368/1‑50 (MQ=255)
                                                       ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg               >  1:1150191/1‑71 (MQ=255)
                                                       ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg               <  1:3724475/71‑1 (MQ=255)
                                                       ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg               >  1:3715574/1‑71 (MQ=255)
                                                       ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgc                <  1:1239249/70‑1 (MQ=255)
                                                        cgagcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgc                >  1:1707907/1‑69 (MQ=255)
                                                         gagcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGt                                 >  1:3120323/1‑51 (MQ=255)
                                                         gagcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGtt             >  1:2354371/1‑71 (MQ=255)
                                                          agcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcg                 <  1:542549/66‑1 (MQ=255)
                                                           gcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTa                                      >  1:575450/1‑44 (MQ=255)
                                                           gcgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTAtt                              >  1:4146121/1‑52 (MQ=255)
                                                            tgagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGttg                             <  1:1255131/51‑2 (MQ=255)
                                                            cgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                                   >  1:310340/1‑46 (MQ=255)
                                                             gagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCTTTTCGc         >  1:2720152/1‑71 (MQ=255)
                                                                  gaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                  <  1:2428047/41‑1 (MQ=255)
                                                                   aTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg        <  1:3595693/66‑1 (MQ=255)
                                                                     gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcgac      >  1:3410138/1‑66 (MQ=255)
                                                                     gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa  >  1:2972399/1‑70 (MQ=255)
                                                                     |                                                                     
AGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACA  >  minE/2595680‑2595818

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: