Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,740 C→T 13.8% S286S (TCG→TCA arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,7400CT13.8% 251.7 / 27.0 94S286S (TCG→TCAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (39/41);  new base T (5/8);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

Only 100 of 112 total aligned reads displayed.

TGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTT  >  minE/2595670‑2595807
                                                                      |                                                                   
tGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                                     >  1:1370070/1‑71 (MQ=255)
tGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                                     >  1:1962964/1‑71 (MQ=255)
 gAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                                     >  1:3869406/1‑70 (MQ=255)
 gAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                                     >  1:795533/1‑70 (MQ=255)
 gAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                                     <  1:2370377/70‑1 (MQ=255)
   gTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcgag                                                                  <  1:3264086/71‑1 (MQ=255)
    tACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAAGCACCAGATACATGCcgagcgag                                                                  >  1:2466920/1‑70 (MQ=255)
    tACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcgag                                                                  >  1:1309749/1‑70 (MQ=255)
      cGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcgag                                                                  >  1:3842745/1‑68 (MQ=255)
      cGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcgag                                                                  >  1:3686725/1‑68 (MQ=255)
       gCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaaga                                                              <  1:3723524/71‑1 (MQ=255)
        ccAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                             >  1:3782293/1‑71 (MQ=255)
        ccAGAAAGATAATCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                             <  1:3180948/71‑1 (MQ=255)
         cAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaaga                                                              <  1:2127303/69‑1 (MQ=255)
          aGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGa                                                           >  1:2427970/1‑71 (MQ=255)
           gAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGa                                                           >  1:1155118/1‑70 (MQ=255)
            aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGAc                                                          >  1:510221/1‑70 (MQ=255)
            aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                         <  1:1492527/71‑1 (MQ=255)
             aaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                                        <  1:2321938/71‑1 (MQ=255)
                atatTCCGTTAATCCGGCATTGCTCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                             <  1:4053749/63‑1 (MQ=255)
                   ttCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                   <  1:2306123/70‑1 (MQ=255)
                    tCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaag                                                               <  1:856557/57‑1 (MQ=255)
                    tCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                   >  1:740359/1‑69 (MQ=255)
                       gTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggg                                              >  1:3618205/1‑71 (MQ=255)
                          aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAAt                                                       >  1:3752041/1‑59 (MQ=255)
                          aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAAt                                                       >  1:572633/1‑59 (MQ=255)
                                 cATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCAccag                                      <  1:677764/69‑3 (MQ=255)
                                   ttGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAAta                                                            <  1:1532692/45‑2 (MQ=255)
                                   ttGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAAta                                                            <  1:377836/45‑2 (MQ=255)
                                   ttGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaa                                                                <  1:1272313/41‑1 (MQ=255)
                                   ttGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                             <  1:2983377/44‑1 (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)
                                     gcgcAGGCCATAAACCACCAGATACATGCCGAGCGagaaga                                                              >  1:3223350/1‑41 (MQ=255)
                                     gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                             >  1:3388879/1‑42 (MQ=255)
                                     gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATg                                                            >  1:1397900/1‑43 (MQ=255)
                                     gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGa                                                           >  1:3033549/1‑44 (MQ=255)
                                     gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCa                                                 >  1:2359650/1‑54 (MQ=255)
                                     gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCa                                                 >  1:2406516/1‑54 (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)
                                           gCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                             <  1:824554/36‑1 (MQ=255)
                                           gCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGa                                >  1:37375/1‑65 (MQ=255)
                                             cATAAACCACCAGATACATGCCGAGCGAGAAGATGa                                                           >  1:1662047/1‑36 (MQ=255)
                                             cATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                        >  1:161659/1‑71 (MQ=255)
                                             cATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                        <  1:474058/71‑1 (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:3172504/36‑1 (MQ=255)
                                               tAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                         <  1:2303278/36‑1 (MQ=255)
                                               tAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                         <  1:2568257/36‑1 (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)
                                                aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTg                                                    >  1:2998659/1‑40 (MQ=255)
                                                aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTg                                                    >  1:3003348/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: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)
                                                          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:3715574/1‑71 (MQ=255)
                                                                 ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg    >  1:1150191/1‑71 (MQ=255)
                                                                 ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg    <  1:3724475/71‑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)
                                                                      cgagAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                        >  1:310340/1‑46 (MQ=255)
                                                                      |                                                                   
TGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTT  >  minE/2595670‑2595807

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: