Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,737 G→C 14.8% L287L (CTC→CTG arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,7370GC14.8% 236.8 / 26.7 94L287L (CTC→CTGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (39/41);  new base C (6/8);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 110 total aligned reads displayed.

GTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTT  >  minE/2595668‑2595807
                                                                     |                                                                      
gTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcga                                                                       <  1:3101453/71‑1 (MQ=255)
 ttGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgag                                                                      <  1:637127/71‑1 (MQ=255)
  tGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                                     >  1:1962964/1‑71 (MQ=255)
  tGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                                     >  1:1370070/1‑71 (MQ=255)
   gAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                                     >  1:795533/1‑70 (MQ=255)
   gAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                                     <  1:2370377/70‑1 (MQ=255)
     gTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcgag                                                                  <  1:3264086/71‑1 (MQ=255)
      tACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcgag                                                                  >  1:1309749/1‑70 (MQ=255)
        cGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcgag                                                                  >  1:3096479/1‑68 (MQ=255)
        cGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcgag                                                                  >  1:3686725/1‑68 (MQ=255)
        cGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcgag                                                                  >  1:3842745/1‑68 (MQ=255)
         gCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaaga                                                              <  1:3723524/71‑1 (MQ=255)
         gCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaaga                                                              <  1:3607875/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)
                  atatTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGAGagaaga                                                              >  1:3770257/1‑62 (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)
                          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)
                                     ttGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAAta                                                            <  1:377836/45‑2 (MQ=255)
                                     ttGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAAta                                                            <  1:1532692/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:2406516/1‑54 (MQ=255)
                                       gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCa                                                 >  1:2359650/1‑54 (MQ=255)
                                          cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                             <  1:325150/71‑1 (MQ=255)
                                          cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                             <  1:2898541/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)
                                               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: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:4079911/1‑70 (MQ=255)
                                                 tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                       <  1:1543030/70‑1 (MQ=255)
                                                  aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTg                                                    >  1:3003348/1‑40 (MQ=255)
                                                  aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTg                                                    >  1:2998659/1‑40 (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)
                                                           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)
                                                              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:3468886/1‑50 (MQ=255)
                                                               cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                             >  1:2056368/1‑50 (MQ=255)
                                                                   ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg    <  1:3724475/71‑1 (MQ=255)
                                                                   ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg    >  1:1150191/1‑71 (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)
                                                                     |                                                                      
GTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTT  >  minE/2595668‑2595807

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: