Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,518 C→T 28.2% A360A (GCG→GCA arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5180CT28.2% 299.9 / 137.6 199A360A (GCG→GCAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (79/63);  new base T (42/14);  total (122/77)
Fisher's exact test for biased strand distribution p-value = 1.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 200 total aligned reads displayed.

AGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCG  >  minE/2595452‑2595582
                                                                  |                                                                
agcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGct                                                               >  1:117541/1‑70 (MQ=255)
agcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAACCGCAGCCAATCACATTGGCATAAACCATCGctt                                                              >  1:1749143/1‑71 (MQ=255)
 gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctc                                                             >  1:4111743/1‑71 (MQ=255)
 gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctc                                                             >  1:138197/1‑71 (MQ=255)
    aGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                                                          <  1:1136169/71‑1 (MQ=255)
      gCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGa                                                        >  1:44154/1‑71 (MQ=255)
       ctgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctct                                                            >  1:3569147/4‑66 (MQ=255)
       ctgccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctct                                                            >  1:445603/4‑66 (MQ=255)
          ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAAc                                                    >  1:3451052/1‑71 (MQ=255)
          ccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGacg                                                      <  1:530882/69‑3 (MQ=255)
            aaTCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                                                  <  1:2905130/71‑1 (MQ=255)
             atcggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                                                  <  1:97747/67‑1 (MQ=255)
              ttGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                                                  >  1:3177689/1‑69 (MQ=255)
              ttGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                                                  >  1:3748577/1‑69 (MQ=255)
                ggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgt                                              >  1:1923256/1‑70 (MQ=255)
                ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGcttctt                                                           >  1:2215367/1‑58 (MQ=255)
                ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                                                     >  1:2880212/1‑64 (MQ=255)
                ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                                                     >  1:623043/1‑64 (MQ=255)
                ggggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                                                  >  1:4738066/1‑67 (MQ=255)
                 gggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                                             >  1:3301054/1‑69 (MQ=255)
                 gggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                                             >  1:5121328/1‑69 (MQ=255)
                 gggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                              <  1:720008/70‑1 (MQ=255)
                  ggTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGAt                                                       <  1:1616517/60‑1 (MQ=255)
                   gTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTg                                                         <  1:4444327/57‑1 (MQ=255)
                    tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                          >  1:4639318/1‑71 (MQ=255)
                     aaTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg                                         >  1:1316500/1‑71 (MQ=255)
                       ttttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGta                                       >  1:870299/1‑70 (MQ=255)
                        tttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATa                                                      <  1:1283530/55‑1 (MQ=255)
                        tttCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtac                                      >  1:4252284/1‑69 (MQ=255)
                             gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                                                          >  1:2205122/1‑46 (MQ=255)
                             gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                                                          <  1:1582934/46‑1 (MQ=255)
                              tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt                                >  1:5619693/1‑71 (MQ=255)
                               cctaaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTc                               <  1:2033291/68‑1 (MQ=255)
                               cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                          >  1:3994363/1‑60 (MQ=255)
                               cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                          >  1:5393428/1‑60 (MQ=255)
                               cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                          >  1:50638/1‑60 (MQ=255)
                               cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt                                >  1:3396011/1‑70 (MQ=255)
                                ccAAATCGCAGCCAATCACATTGGCATAAACCATCGcttct                                                            <  1:3713674/41‑1 (MQ=255)
                                ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt                                                          <  1:5252944/43‑1 (MQ=255)
                                ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                               <  1:3367170/70‑1 (MQ=255)
                                ccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                               <  1:3949922/70‑1 (MQ=255)
                                 cAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                             <  1:3422265/71‑1 (MQ=255)
                                  aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                            >  1:2927406/1‑71 (MQ=255)
                                   aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                           <  1:334155/71‑1 (MQ=255)
                                   aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                                             >  1:269945/1‑51 (MQ=255)
                                   aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAtc                           <  1:911476/71‑2 (MQ=255)
                                   aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAtc                           <  1:5437133/71‑2 (MQ=255)
                                     tcacAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                         >  1:2392320/4‑71 (MQ=255)
                                     tCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                         >  1:4288879/1‑71 (MQ=255)
                                     tCGCAGCCAATCACATTGGCATAAACAATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                         >  1:167266/1‑71 (MQ=255)
                                      cacAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                                               >  1:5465245/3‑48 (MQ=255)
                                       gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                        <  1:3255432/70‑1 (MQ=255)
                                       gCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGAc                        <  1:4587681/70‑1 (MQ=255)
                                       acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGAc                        <  1:1526969/69‑1 (MQ=255)
                                       acAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACa                       <  1:214923/70‑1 (MQ=255)
                                        cAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                            <  1:84031/65‑1 (MQ=255)
                                         aGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                        <  1:4159363/52‑1 (MQ=255)
                                              aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                         >  1:595489/1‑62 (MQ=255)
                                              aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                         >  1:2777726/1‑62 (MQ=255)
                                              aTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                         >  1:2176540/1‑62 (MQ=255)
                                              aTCACATTGGCATAAACCATCGCTTCTTTGATAACGcc                                                 <  1:553610/38‑1 (MQ=255)
                                              aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                <  1:2079645/71‑1 (MQ=255)
                                               tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                         >  1:2594187/1‑61 (MQ=255)
                                               tCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGa                         >  1:4010419/1‑61 (MQ=255)
                                               tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGa                                              >  1:2980187/1‑40 (MQ=255)
                                                cacaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa              >  1:5013667/1‑71 (MQ=255)
                                                 acaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgca                  >  1:3832121/1‑65 (MQ=255)
                                                 acaTTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgca                  >  1:223661/1‑65 (MQ=255)
                                                 acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa              >  1:3466013/1‑70 (MQ=255)
                                                 acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa              <  1:466127/70‑1 (MQ=255)
                                                 acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc             <  1:1037860/71‑1 (MQ=255)
                                                 acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc             >  1:2308503/1‑71 (MQ=255)
                                                  caTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAc             <  1:1966156/70‑1 (MQ=255)
                                                   aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                         >  1:823567/1‑57 (MQ=255)
                                                   aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa           <  1:4849085/71‑1 (MQ=255)
                                                    ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtactc                                    >  1:2718459/1‑44 (MQ=255)
                                                    ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCaa                             >  1:79368/1‑52 (MQ=255)
                                                    ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                            >  1:2988091/1‑53 (MQ=255)
                                                    ttGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                            >  1:3608310/1‑53 (MQ=255)
                                                    ttGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAg          <  1:3909527/71‑1 (MQ=255)
                                                     tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                             >  1:2678707/1‑51 (MQ=255)
                                                     tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                         >  1:2541286/1‑55 (MQ=255)
                                                     tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                         >  1:1689067/1‑55 (MQ=255)
                                                     tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCaa              >  1:3055542/1‑66 (MQ=255)
                                                      ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                             >  1:3308327/1‑50 (MQ=255)
                                                      ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                             >  1:3142306/1‑50 (MQ=255)
                                                      ggCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGt         <  1:1822723/70‑1 (MQ=255)
                                                       gCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTa        <  1:2396974/70‑1 (MQ=255)
                                                          tAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                        <  1:1312794/35‑1 (MQ=255)
                                                          tAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                        <  1:2900094/35‑1 (MQ=255)
                                                          tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                           <  1:135212/48‑1 (MQ=255)
                                                           aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAc       >  1:1249439/1‑67 (MQ=255)
                                                           aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTAc       >  1:1088319/1‑67 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                      >  1:2762815/1‑51 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg  >  1:190643/1‑71 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg  >  1:4762948/1‑71 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCg  >  1:4005558/1‑71 (MQ=255)
                                                            aaCCATCGCTTCTTTGATAACGCCAGATGCCGAGCTGCCATCAATGGACAACGCGCCAAcc            >  1:4243777/1‑61 (MQ=255)
                                                            aaCCATAGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgcg                  >  1:3215243/1‑55 (MQ=255)
                                                               cATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                               <  1:769114/39‑1 (MQ=255)
                                                                  |                                                                
AGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCG  >  minE/2595452‑2595582

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: