Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,405 C→G 9.3% S398T (AGC→ACC) ‡ arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,4050CG 9.3% 253.1 / 13.3 86S398T (AGC→ACC) ‡arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (36/42);  new base G (2/6);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAAT  >  minE/2595336‑2595475
                                                                     |                                                                      
aGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCt                                                                       <  1:3509221/71‑1 (MQ=255)
aGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATACCCAGCt                                                                       <  1:2398680/71‑1 (MQ=255)
   gcggccagaGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAt                                                                    >  1:304701/4‑71 (MQ=255)
   gcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                                                     >  1:1051765/1‑70 (MQ=255)
   gcAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGa                                                                     >  1:345183/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:3696913/1‑70 (MQ=255)
                                   cATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGc                                     >  1:2458076/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:2359636/1‑69 (MQ=255)
                                      gATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCa                                   >  1:2284123/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:3287712/54‑1 (MQ=255)
                                                gTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGc                                        <  1:438968/54‑1 (MQ=255)
                                                gTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgt                               >  1:436095/1‑63 (MQ=255)
                                                 tACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgtagc                       <  1:1214801/70‑1 (MQ=255)
                                                 tACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgtagc                       >  1:3381718/1‑70 (MQ=255)
                                                 tACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgtagc                       <  1:528318/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:2274053/71‑1 (MQ=255)
                                                         tAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTAcc              <  1:3748884/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:478390/1‑53 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgta                              >  1:1923009/1‑53 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTa                <  1:1208785/67‑1 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCaa            >  1:148406/1‑71 (MQ=255)
                                                           atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCaa            >  1:3269888/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)
                                                                     cTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAgtagcgta                         <  1:1840950/48‑1 (MQ=255)
                                                                     cTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAAt  <  1:1007563/71‑1 (MQ=255)
                                                                     |                                                                      
AGCGCAGCCAGCGTCACAAACAGCACAGGCAGGGTCATGATAATCCCTGTACGGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAAT  >  minE/2595336‑2595475

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: