Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,404 G→A 9.4% S398S (AGC→AGT) ‡ arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,4040GA 9.4% 259.7 / 13.3 85S398S (AGC→AGT) ‡arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (35/42);  new base A (2/6);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

AGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTA  >  minE/2595336‑2595473
                                                                    |                                                                     
aGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCt                                                                     <  1:3509221/71‑1 (MQ=255)
aGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATACCCAGCt                                                                     <  1:2398680/71‑1 (MQ=255)
   gcggccagaGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAt                                                                  >  1:304701/4‑71 (MQ=255)
   gcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                                                   >  1:345183/1‑70 (MQ=255)
   gcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                                                   >  1:1051765/1‑70 (MQ=255)
       ccagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTg                                                                    <  1:3610481/65‑1 (MQ=255)
        cagcGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCa                                                             >  1:542016/1‑71 (MQ=255)
             tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc                                                        >  1:2031137/1‑71 (MQ=255)
             tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc                                                        <  1:561201/71‑1 (MQ=255)
             tCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttc                                                        <  1:1822834/71‑1 (MQ=255)
               acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTg                                                                    <  1:484389/57‑1 (MQ=255)
               acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTg                                                                    <  1:1752835/57‑1 (MQ=255)
               acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                                                      >  1:2165363/1‑71 (MQ=255)
               acaAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttctt                                                      >  1:2506286/1‑71 (MQ=255)
                caAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttcttc                                                     <  1:1889340/71‑1 (MQ=255)
                caAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttcttc                                                     >  1:2525046/1‑71 (MQ=255)
                caAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttcttc                                                     <  1:2508316/71‑1 (MQ=255)
                      gCACAGGCAGGGTCATGTTAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGaa                                               <  1:3516601/71‑1 (MQ=255)
                        acaGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAg                                             <  1:1201508/71‑1 (MQ=255)
                           ggcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTc                                                              <  1:3288447/51‑1 (MQ=255)
                           ggcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTc                                                              <  1:848417/51‑1 (MQ=255)
                            gcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACg                                         >  1:330065/1‑71 (MQ=255)
                            gcaggGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACg                                         >  1:51357/1‑71 (MQ=255)
                               gggTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATAttcttc                                                     <  1:974816/56‑1 (MQ=255)
                                ggTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTg                                       >  1:76264/1‑69 (MQ=255)
                                  tCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGc                                   <  1:770669/71‑1 (MQ=255)
                                   cATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGc                                   >  1:2458076/1‑70 (MQ=255)
                                   cATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGc                                   >  1:3696913/1‑70 (MQ=255)
                                   cATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttttt                                                      <  1:315477/51‑1 (MQ=255)
                                    aTGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGc                                   >  1:2529828/1‑69 (MQ=255)
                                    aTGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa                                 <  1:3659215/71‑1 (MQ=255)
                                    aTGATAATCACTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTa                                           >  1:4099395/1‑61 (MQ=255)
                                    aTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCa                                 <  1:3072840/71‑1 (MQ=255)
                                    aTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCa                                 >  1:1819255/1‑71 (MQ=255)
                                     tGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCag                                <  1:2091029/71‑1 (MQ=255)
                                      gATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCat                                                            >  1:1226182/1‑42 (MQ=255)
                                      gATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa                                 >  1:2284123/1‑69 (MQ=255)
                                      gATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa                                 >  1:2359636/1‑69 (MQ=255)
                                       aTAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCagag                              >  1:1817348/1‑71 (MQ=255)
                                       aTAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa                                 >  1:2737528/1‑68 (MQ=255)
                                        tAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgt                             <  1:528847/71‑1 (MQ=255)
                                          aTCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta                            <  1:2374533/70‑1 (MQ=255)
                                           tCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagc                          <  1:3396641/71‑1 (MQ=255)
                                             ccTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcg                         >  1:2036715/1‑70 (MQ=255)
                                              cTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCaga                               <  1:3315332/63‑1 (MQ=255)
                                                gTACGGAAGTAATATCCCGAGCTGATAGTCATATTCTTCTGCGAAAg                                             <  1:3914880/47‑1 (MQ=255)
                                                gTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGc                                      <  1:438968/54‑1 (MQ=255)
                                                gTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGc                                      <  1:3287712/54‑1 (MQ=255)
                                                gTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgt                             >  1:436095/1‑63 (MQ=255)
                                                 tACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgtagc                     <  1:528318/70‑1 (MQ=255)
                                                 tACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgtagc                     >  1:3381718/1‑70 (MQ=255)
                                                 tACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgtagc                     <  1:1214801/70‑1 (MQ=255)
                                                 tACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGcc                    <  1:1973345/71‑1 (MQ=255)
                                                   ctgAAGTAATATCCCCTGCTGATAGTCATATTCTTCTGCGAAAGt                                            <  1:3977680/43‑1 (MQ=255)
                                                   cGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGt                                            <  1:2361278/45‑1 (MQ=255)
                                                   cGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAg                  >  1:3162509/1‑71 (MQ=255)
                                                   cGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCa                   <  1:783470/70‑1 (MQ=255)
                                                    ggAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAg                                             <  1:2964930/43‑1 (MQ=255)
                                                     gAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGc                >  1:3318851/1‑71 (MQ=255)
                                                     gAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCa                   <  1:778668/68‑1 (MQ=255)
                                                      aaGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGc                >  1:213149/1‑70 (MQ=255)
                                                         tAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTAcc            <  1:3748884/71‑1 (MQ=255)
                                                         tAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTAcc            <  1:2274053/71‑1 (MQ=255)
                                                         tAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAgc                                            <  1:2155694/39‑2 (MQ=12)
                                                          aaTATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTa              <  1:1499285/68‑1 (MQ=255)
                                                          aaTATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTa              <  1:2452325/68‑1 (MQ=255)
                                                          aaTATCCCCAGCTGATAGTCATATTCTTCTGCGAAACTACg                                         >  1:442012/1‑41 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta                            >  1:1923009/1‑53 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta                            >  1:478390/1‑53 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTa              <  1:1208785/67‑1 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCaa          >  1:3269888/1‑71 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCaa          >  1:148406/1‑71 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCa           <  1:473574/70‑1 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCa           <  1:3014649/70‑1 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAAAGTAGCGTAGCCAGGCTACCaa          >  1:844001/1‑71 (MQ=255)
                                                              tCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTa                                           <  1:3823627/35‑1 (MQ=255)
                                                              tCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta                            <  1:140861/50‑1 (MQ=255)
                                                              tCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGcc                    >  1:548470/1‑58 (MQ=255)
                                                               ccccAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgtagc                     <  1:3824497/56‑1 (MQ=255)
                                                               ccccAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTgg      >  1:2661469/1‑71 (MQ=255)
                                                                cccAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTggg     >  1:3862579/1‑71 (MQ=255)
                                                                 ccAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTgggg    >  1:1200205/1‑71 (MQ=255)
                                                                  cAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgta                       <  1:698790/51‑1 (MQ=255)
                                                                  cAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgta                       <  1:1175047/51‑1 (MQ=255)
                                                                   aGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTa  <  1:3205980/71‑1 (MQ=255)
                                                                    gCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACaaa          >  1:1958845/1‑62 (MQ=255)
                                                                    |                                                                     
AGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTA  >  minE/2595336‑2595473

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: