Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,728 A→G 11.9% Y290Y (TAT→TAC arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,7280AG11.9% 220.0 / 17.5 84Y290Y (TAT→TACarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (33/41);  new base G (7/3);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGC  >  minE/2595659‑2595796
                                                                     |                                                                    
cGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGa                                                                      >  1:3093729/1‑70 (MQ=255)
cGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGa                                                                      >  1:3222186/1‑70 (MQ=255)
  ccAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATa                                                                    >  1:2179942/1‑70 (MQ=255)
         gTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATa                                                                    >  1:271931/1‑63 (MQ=255)
         gTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcga                                                            <  1:381778/71‑1 (MQ=255)
         gTTGAGTACGCCAGAAAGATATTCCGTTAATCAGGCATTGCGCAGGCCATAAACCACCAGATACATGcc                                                              >  1:337736/1‑69 (MQ=255)
           tGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATg                                                                >  1:3548705/1‑65 (MQ=255)
             aGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcg                                                         <  1:1017314/70‑1 (MQ=255)
              gcaCGCCAGACAGGTACTCCGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGcc                                                              <  1:2402821/62‑1 (MQ=255)
              gTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcgag                                                       <  1:2219218/71‑1 (MQ=255)
                 cGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGAGagaag                                                    <  1:1168868/71‑1 (MQ=255)
                   ccAGACAGGTACTCCGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGagaaaa                                                   >  1:1892086/1‑68 (MQ=255)
                   ccAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                  <  1:3838769/71‑1 (MQ=255)
                   ccAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                  <  1:2421712/71‑1 (MQ=255)
                    cAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAACCCACCAGATACATGCcg                                                             <  1:3641748/59‑1 (MQ=255)
                    cAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcg                                                             <  1:636218/59‑1 (MQ=255)
                     aGACAGGTACTCCGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAAtaa                                                >  1:346431/1‑69 (MQ=255)
                       aaaGATATTCCGTTAATCCGTCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGAc                                               >  1:2643353/1‑70 (MQ=255)
                       aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCgaga                                                      >  1:3468182/1‑63 (MQ=255)
                       aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCgaga                                                      >  1:3053704/1‑63 (MQ=255)
                       aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCgaga                                                      >  1:182216/1‑63 (MQ=255)
                       aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaa                                                     <  1:1519251/64‑1 (MQ=255)
                       aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGAc                                               >  1:1498331/1‑70 (MQ=255)
                       aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGAc                                               >  1:2225303/1‑70 (MQ=255)
                        aaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                             <  1:582541/71‑1 (MQ=255)
                        aaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                             >  1:449962/1‑71 (MQ=255)
                        aaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                             <  1:406610/71‑1 (MQ=255)
                        aaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                             <  1:34104/71‑1 (MQ=255)
                        aaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                             <  1:1708659/71‑1 (MQ=255)
                         aGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaag                                                    >  1:1738566/1‑63 (MQ=255)
                         aGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAAt                                            <  1:1266934/71‑1 (MQ=255)
                         aGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAAt                                            <  1:3062096/71‑1 (MQ=255)
                         aGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAAt                                            <  1:3340653/71‑1 (MQ=255)
                         aGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGAGAGAAGATGACAAt                                            >  1:1714121/1‑71 (MQ=255)
                            tatTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCt                                          <  1:3226815/70‑1 (MQ=255)
                            tactCCGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTg                                         >  1:2525585/4‑71 (MQ=255)
                            tactCCGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTg                                         >  1:1357401/4‑71 (MQ=255)
                              ttCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGcc                                       <  1:2541932/71‑1 (MQ=255)
                              ttCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                        <  1:1969698/70‑1 (MQ=255)
                                 cGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCatg                                    >  1:3096342/1‑69 (MQ=255)
                                 cGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCa                                      <  1:551761/69‑1 (MQ=255)
                                  gTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggg                                   <  1:3734046/71‑1 (MQ=255)
                                  gTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAg                                     <  1:405584/69‑1 (MQ=255)
                                   ttAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                          <  1:1228813/47‑1 (MQ=255)
                                   ttAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                          <  1:603853/47‑1 (MQ=255)
                                    tAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggggg                                 >  1:2311898/1‑71 (MQ=255)
                                    tAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggggg                                 >  1:472871/1‑71 (MQ=255)
                                    tAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAgggg                                  >  1:1300510/1‑70 (MQ=255)
                                    tAATACGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAgggg                                  >  1:1199977/1‑70 (MQ=255)
                                     aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACa                                                                  <  1:2162589/37‑1 (MQ=255)
                                     aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACa                                                                  <  1:3586714/37‑1 (MQ=255)
                                     aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACa                                                                  <  1:3484224/37‑1 (MQ=255)
                                     aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggggg                                 >  1:1630935/1‑70 (MQ=255)
                                     aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggggg                                 <  1:500776/70‑1 (MQ=255)
                                       tCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCAc                              <  1:3048356/71‑1 (MQ=255)
                                         cGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCAcc                             <  1:2819862/70‑1 (MQ=255)
                                          ggCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACcgc                           >  1:2607462/1‑71 (MQ=255)
                                          ggCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACcgc                           <  1:551492/71‑1 (MQ=255)
                                          ggCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACcgc                           >  1:548508/1‑71 (MQ=255)
                                            cATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                        <  1:1164393/56‑1 (MQ=255)
                                             aTTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                  <  1:2340247/45‑1 (MQ=255)
                                             aTTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                  <  1:1163289/45‑1 (MQ=255)
                                                gcgcAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGAt                                            >  1:353962/1‑48 (MQ=255)
                                                gcgcAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGAt                                            >  1:2624037/1‑48 (MQ=255)
                                                gcgcAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAgg                      <  1:3652944/70‑1 (MQ=255)
                                                  gcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAgg                                    <  1:1644523/54‑1 (MQ=255)
                                                  gcAGGCCATAAACCACCAGATAAATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACt                   <  1:2089511/71‑1 (MQ=255)
                                                   cAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                  >  1:1379005/1‑71 (MQ=255)
                                                    aGGCCATAAACCACCAGATACATGCGGAGCGagaaga                                                   <  1:673405/37‑1 (MQ=255)
                                                    aGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACttt                 <  1:3576747/71‑1 (MQ=255)
                                                    aGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTCCCAGGGGGCACCGCGCAGGACttt                 <  1:3056390/71‑1 (MQ=255)
                                                      gCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                             <  1:2817311/41‑1 (MQ=255)
                                                        cATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc             <  1:3393576/71‑1 (MQ=255)
                                                          tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg            >  1:2769848/1‑70 (MQ=255)
                                                          tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg            >  1:3753967/1‑70 (MQ=255)
                                                          tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg            >  1:3751767/1‑70 (MQ=255)
                                                          tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg            >  1:3745906/1‑70 (MQ=255)
                                                           aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa          >  1:3490036/1‑71 (MQ=255)
                                                           aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa          >  1:1123326/1‑71 (MQ=255)
                                                            aaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAg                                     >  1:2332685/1‑43 (MQ=255)
                                                             accaccAGATACATGCCGAGCGAGAAGATGACAATCTGcc                                       >  1:1237171/1‑40 (MQ=255)
                                                              ccaccaGATACATGCCGAGCGAGAAGATGACAATCTGCCa                                      >  1:86078/1‑40 (MQ=255)
                                                              ccaccaGATACATGCCGAGCGAGAAGATGACAATCTGCCa                                      >  1:1073704/1‑40 (MQ=255)
                                                                accaGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa      >  1:3168596/1‑70 (MQ=255)
                                                                  caGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg   >  1:44901/1‑71 (MQ=255)
                                                                   aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGc  <  1:2151999/71‑1 (MQ=255)
                                                                     |                                                                    
CGCCAGCACGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGC  >  minE/2595659‑2595796

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: