Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,752 A→G 16.7% I282I (ATT→ATC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,7520AG16.7% 183.8 / 29.9 78I282I (ATT→ATCarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (39/26);  new base G (6/7);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

AAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAA  >  minE/2595682‑2595820
                                                                      |                                                                    
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:3752041/1‑59 (MQ=255)
              aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAAt                                                                    >  1:572633/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)
                         gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCa                                                              >  1:2406516/1‑54 (MQ=255)
                         gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCa                                                              >  1:2359650/1‑54 (MQ=255)
                            cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                          <  1:2898541/71‑1 (MQ=255)
                            cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                          <  1:3630922/71‑1 (MQ=255)
                            cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                          <  1:325150/71‑1 (MQ=255)
                            cAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                                          <  1:3579096/71‑1 (MQ=255)
                            cAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                                          <  1:2720411/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)
                                 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:2568257/36‑1 (MQ=255)
                                   tAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                                      <  1:3172504/36‑1 (MQ=255)
                                   tAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                                      <  1:2303278/36‑1 (MQ=255)
                                   tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCt                                                                  >  1:4033837/1‑40 (MQ=255)
                                   tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCt                                                                  >  1:2356609/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:3003348/1‑40 (MQ=255)
                                    aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTg                                                                 >  1:2998659/1‑40 (MQ=255)
                                    aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                                >  1:307417/1‑41 (MQ=255)
                                    aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                                >  1:1056650/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:3218807/1‑71 (MQ=255)
                                           caGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                           <  1:873797/71‑1 (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:2411812/1‑45 (MQ=255)
                                                 cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAg                                               >  1:1187780/1‑45 (MQ=255)
                                                 cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                          >  1:2056368/1‑50 (MQ=255)
                                                 cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                          >  1:3468886/1‑50 (MQ=255)
                                                     ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                 >  1:3715574/1‑71 (MQ=255)
                                                     ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                 <  1:3724475/71‑1 (MQ=255)
                                                     ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                 >  1:1150191/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)
                                                                    aCAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa  >  1:3746971/1‑71 (MQ=255)
                                                                     cAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaaa  >  1:3127947/1‑70 (MQ=255)
                                                                      |                                                                    
AAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAA  >  minE/2595682‑2595820

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: