Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,761 G→T 17.9% P279P (CCC→CCA arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,7610GT17.9% 218.6 / 19.3 72P279P (CCC→CCAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base T (6/7);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGC  >  minE/2595693‑2595830
                                                                    |                                                                     
gTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAggg                                                                     >  1:3618205/1‑71 (MQ=255)
          cATTGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCAccag                                                             <  1:677764/69‑3 (MQ=255)
            ttGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCgcgc                                                           <  1:70289/69‑1 (MQ=255)
             tGCGCAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAg                                                         >  1:3425185/1‑70 (MQ=255)
                 cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                                    <  1:2898541/71‑1 (MQ=255)
                 cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                                    <  1:325150/71‑1 (MQ=255)
                 cAGGCCATAGACCACGATGTACATGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACtt                                                    <  1:3630922/71‑1 (MQ=255)
                 cAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                                                    <  1:2720411/71‑1 (MQ=255)
                 cAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACtt                                                    <  1:3579096/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: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)
                        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)
                         aaaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa                                            >  1:1050198/1‑71 (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:1150191/1‑71 (MQ=255)
                                          ccGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                           <  1:3724475/71‑1 (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)
                                                            aTCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                                        <  1:2015137/40‑1 (MQ=255)
                                                            aTCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa         <  1:649158/71‑1 (MQ=255)
                                                             tCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa         >  1:2263615/1‑70 (MQ=255)
                                                              cTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                                         >  1:424177/1‑37 (MQ=255)
                                                               tGCCAGGGGGCACCGCTCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat      >  1:1346485/1‑71 (MQ=255)
                                                               tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCAt                                  >  1:2794700/1‑43 (MQ=255)
                                                               tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa         >  1:339359/1‑68 (MQ=255)
                                                               tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat      >  1:3684755/1‑71 (MQ=255)
                                                               tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAatat      >  1:1679962/1‑71 (MQ=255)
                                                                gCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcga                   >  1:63436/1‑57 (MQ=255)
                                                                 ccAGTGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                    <  1:2647464/55‑1 (MQ=255)
                                                                 ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                    <  1:1266369/55‑1 (MQ=255)
                                                                 ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcg                    <  1:349721/55‑1 (MQ=255)
                                                                 ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat        <  1:1589338/67‑1 (MQ=255)
                                                                 ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAataata       <  1:4052833/68‑1 (MQ=255)
                                                                 ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCa    <  1:959771/71‑1 (MQ=255)
                                                                  cAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAg   <  1:936051/71‑1 (MQ=255)
                                                                    gggggCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAgc  <  1:70213/70‑1 (MQ=255)
                                                                    |                                                                     
GTTAATCCGGCATTGCGCAGGCCATAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGC  >  minE/2595693‑2595830

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: