Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,737 G→C 35.9% 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,7370GC35.9% 118.3 / 128.0 131L287L (CTC→CTGarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (44/40);  new base C (23/24);  total (67/64)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

Only 100 of 139 total aligned reads displayed.

CGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGT  >  minE/2595667‑2595806
                                                                      |                                                                     
cGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcg                                                                       <  1:2684351/71‑1 (MQ=255)
   tGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                                    >  1:1370406/1‑71 (MQ=255)
        aCGCCAGACAGGTACTCCGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGagaa                                                               <  1:3268683/71‑1 (MQ=255)
         cGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagcgag                                                                 >  1:3286343/1‑68 (MQ=255)
          gCCAGACAGGTACTCCGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGagaaaa                                                             <  1:951590/71‑3 (MQ=255)
           ccAGACAGGTACTCCGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAAt                                                            <  1:328670/71‑1 (MQ=255)
           ccAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                            <  1:5096957/71‑1 (MQ=255)
            cAGACAGGTACTCCGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAAt                                                            <  1:3390875/70‑1 (MQ=255)
            cAGACAGGTACTCCGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAAt                                                            >  1:4179180/1‑70 (MQ=255)
               aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaa                                                               >  1:3888093/1‑64 (MQ=255)
               aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaa                                                               >  1:1902425/1‑64 (MQ=255)
               aaaGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaa                                                               >  1:5593260/1‑64 (MQ=255)
                 aGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAAt                                                      >  1:5290268/1‑71 (MQ=255)
                   atatTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaag                                                              <  1:3366957/61‑1 (MQ=255)
                      ttCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCcgagc                                                                    <  1:343597/52‑1 (MQ=255)
                      ttCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGcc                                                 <  1:5425102/71‑1 (MQ=255)
                      ttCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGcc                                                 <  1:2489466/71‑1 (MQ=255)
                      ttCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                  <  1:639752/70‑1 (MQ=255)
                       tCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                  >  1:4215380/1‑69 (MQ=255)
                        ccGTGAACCCTGCATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGc                                                  >  1:1502745/1‑68 (MQ=255)
                         cGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaaga                                                             >  1:3009797/1‑56 (MQ=255)
                         cGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCa                                                <  1:2437817/69‑1 (MQ=255)
                             aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggggg                                           <  1:5190734/70‑1 (MQ=255)
                             aaTCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggggg                                           <  1:1520409/70‑1 (MQ=255)
                               tCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGcc                                                 <  1:4721886/62‑1 (MQ=255)
                               tCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGc                                          <  1:395201/69‑1 (MQ=255)
                                 cGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCt                                                    <  1:4706119/57‑1 (MQ=255)
                                   gCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATc                                                     >  1:5226937/1‑54 (MQ=255)
                                    cATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATc                                                     >  1:5413957/1‑53 (MQ=255)
                                    cATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATc                                                     >  1:4682640/1‑53 (MQ=255)
                                    cATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATc                                                     >  1:5480005/1‑53 (MQ=255)
                                    cATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACcgcg                                    <  1:4448805/70‑1 (MQ=255)
                                    cATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggg                                             >  1:1810563/1‑61 (MQ=255)
                                    cATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACcgcg                                    >  1:3923888/1‑70 (MQ=255)
                                    cATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACcgcg                                    >  1:5508891/1‑70 (MQ=255)
                                      ttGCGCAGGCCATAGACCACGATGTACATGCCCAGTgaga                                                                <  1:628327/40‑1 (MQ=255)
                                      ttGCGCAGGCCATAGACCACGATGTACATGCCCAGTgaga                                                                <  1:4334479/40‑1 (MQ=255)
                                      ttGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCa                                  >  1:1753101/1‑70 (MQ=255)
                                      ttGCGCAGGCCATAAACCACCAGATACATGCCGAGCGagaag                                                              <  1:5193117/42‑1 (MQ=255)
                                      ttGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAg                                 <  1:4337317/71‑1 (MQ=255)
                                       tGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAg                                 >  1:4959400/1‑70 (MQ=255)
                                       tGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAg                                 >  1:2682276/1‑70 (MQ=255)
                                       tGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                        <  1:4538676/47‑1 (MQ=255)
                                       tGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                        <  1:4774993/47‑1 (MQ=255)
                                       tGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                        <  1:1414141/47‑1 (MQ=255)
                                        gcgcAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGAt                                                      >  1:4017813/1‑48 (MQ=255)
                                        gcgcAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATc                                                     >  1:5028236/1‑49 (MQ=255)
                                        gcgcAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCAtggc                                            >  1:5300690/1‑57 (MQ=255)
                                        gcgcAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCAtggc                                            >  1:4981360/1‑57 (MQ=255)
                                        gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGAt                                                            >  1:4180699/1‑42 (MQ=255)
                                        gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                                       >  1:1478264/1‑47 (MQ=255)
                                        gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACaa                                                       >  1:5657266/1‑47 (MQ=255)
                                        gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATc                                                     >  1:3537021/1‑49 (MQ=255)
                                        gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATc                                                     >  1:838677/1‑49 (MQ=255)
                                        gcgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGa                               >  1:3144657/1‑71 (MQ=255)
                                         cgcAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGAc                              <  1:3725958/71‑1 (MQ=255)
                                           cAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                            >  1:117966/1‑71 (MQ=255)
                                            aGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                            <  1:5676939/70‑1 (MQ=255)
                                             ggCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                            <  1:2847733/69‑1 (MQ=255)
                                              gCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAc                         >  1:3936944/1‑71 (MQ=255)
                                               ccATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtttt                          >  1:3117455/1‑68 (MQ=255)
                                               ccATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtttt                          >  1:176605/1‑68 (MQ=255)
                                               ccATAAACCACCAGATACATGCCGAGCGAGAAGATGa                                                          >  1:1987143/1‑37 (MQ=255)
                                               ccATAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                        <  1:5257019/39‑1 (MQ=255)
                                                cATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAgg                                              >  1:1163787/1‑48 (MQ=255)
                                                 aTAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGAc                              <  1:4301593/63‑1 (MQ=255)
                                                 aTAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGAc                              <  1:786349/63‑1 (MQ=255)
                                                 aTAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGAc                              <  1:192852/63‑1 (MQ=255)
                                                 aTAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGAc                              <  1:3387899/63‑1 (MQ=255)
                                                 aTAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                  <  1:880141/43‑1 (MQ=255)
                                                 aTAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                  <  1:3435103/43‑1 (MQ=255)
                                                 aTAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGc                                                  <  1:3889169/43‑1 (MQ=255)
                                                  tAAACCACCAGATATATGCCGAGCGAGAAGATGACAATCTGCCAggggg                                           <  1:4049902/49‑1 (MQ=255)
                                                  tAAACCACCAGATACATGCCGAGCGAGAAGATGACa                                                        <  1:1484990/36‑1 (MQ=255)
                                                  tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                       <  1:299542/69‑1 (MQ=255)
                                                   aaaCCACCAGATACATGCCGAGCGAGAAGATGACAAt                                                      >  1:2089697/1‑37 (MQ=255)
                                                    aaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCa                                                >  1:5545859/1‑42 (MQ=255)
                                                    aaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCAcc                                       >  1:1532783/1‑51 (MQ=255)
                                                      ccaccaGATACATGCCGAGCGAGAAGATGACAATCTGc                                                  >  1:1349049/1‑38 (MQ=255)
                                                      ccaccaGATACATGCCGAGCGAGAAGATGACAATCTGc                                                  >  1:1025090/1‑38 (MQ=255)
                                                      ccaccaGATACATGCCGAGCGAGAAGATGACAATCTGc                                                  >  1:457586/1‑38 (MQ=255)
                                                      ccaccaGATACATGCCGAGCGAGAAGATGACAATCTGc                                                  >  1:945133/1‑38 (MQ=255)
                                                       caccaGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                <  1:257140/71‑1 (MQ=255)
                                                         ccaGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATc              <  1:2879642/71‑1 (MQ=255)
                                                          caGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg             <  1:3722434/71‑1 (MQ=255)
                                                           aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAg                                 >  1:1036405/1‑50 (MQ=255)
                                                           aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAg                                 >  1:234510/1‑50 (MQ=255)
                                                           aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg             >  1:4186020/1‑70 (MQ=255)
                                                             aTACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGc            <  1:5510319/69‑1 (MQ=255)
                                                             aTACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCAt          <  1:1668657/71‑1 (MQ=255)
                                                               aCATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGa        >  1:5107652/1‑71 (MQ=255)
                                                                cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGtg                    >  1:5689818/1‑57 (MQ=255)
                                                                cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGtg                    >  1:4378702/1‑57 (MQ=255)
                                                                cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAc       >  1:180481/1‑71 (MQ=255)
                                                                cATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAc       >  1:1060550/1‑71 (MQ=255)
                                                                  tGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                 <  1:2975736/59‑1 (MQ=255)
                                                                  tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACttt                           <  1:4840439/49‑1 (MQ=255)
                                                                  tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACttt                           <  1:4282590/49‑1 (MQ=255)
                                                                   gCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGa                               >  1:4394685/1‑44 (MQ=255)
                                                                    ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg   <  1:2361138/71‑1 (MQ=255)
                                                                      |                                                                     
CGTTGAGTACGCCAGAAAGATATTCCGTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGT  >  minE/2595667‑2595806

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: