Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,491 G→A 25.0% C369C (TGC→TGT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,4910GA25.0% 229.3 / 78.8 128C369C (TGC→TGTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (54/42);  new base A (21/11);  total (75/53)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

Only 100 of 134 total aligned reads displayed.

CTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGA  >  minE/2595422‑2595559
                                                                     |                                                                    
ctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCg                                                                      >  1:2872330/1‑70 (MQ=255)
ctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCg                                                                      <  1:2438638/70‑1 (MQ=255)
 tGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCg                                                                      >  1:1985086/1‑69 (MQ=255)
    gAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGc                                                                  >  1:3683788/1‑70 (MQ=255)
     aaaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGc                                                                  >  1:1781988/1‑69 (MQ=255)
     aaaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGc                                                                  >  1:744154/1‑69 (MQ=255)
      aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCaa                                                               >  1:1406051/1‑71 (MQ=255)
          aCGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATca                                                            >  1:1281937/1‑70 (MQ=255)
           catGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATcaca                                                          >  1:751980/3‑71 (MQ=255)
           cGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATcaca                                                          >  1:1926630/1‑71 (MQ=255)
           cGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATcaca                                                          <  1:272622/71‑1 (MQ=255)
             tGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACAtt                                                        <  1:2537102/71‑1 (MQ=255)
              gCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTg                                                       >  1:3045573/1‑71 (MQ=255)
              gCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTg                                                       >  1:3239388/1‑71 (MQ=255)
               cAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGTGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTg                                                       <  1:1087227/70‑1 (MQ=255)
               cAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTg                                                       <  1:831074/70‑1 (MQ=255)
               cAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTg                                                       <  1:1043247/70‑1 (MQ=255)
                aGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTg                                                       <  1:2735754/69‑1 (MQ=255)
                  ccAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCa                                                    >  1:2097381/1‑70 (MQ=255)
                   cAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                                  <  1:3680750/71‑1 (MQ=255)
                    agagTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                                                  >  1:927278/1‑70 (MQ=255)
                     gagTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATcaca                                                          <  1:2372881/61‑1 (MQ=255)
                     gagTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACAtt                                                        >  1:533963/1‑63 (MQ=255)
                     gagTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACAtt                                                        >  1:875123/1‑63 (MQ=255)
                      agTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAc                                               >  1:65735/1‑71 (MQ=255)
                       gtagcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAcc                                              <  1:2796306/71‑1 (MQ=255)
                        tagcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCa                                             <  1:1573024/71‑1 (MQ=255)
                         agcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCa                                             <  1:1378962/70‑1 (MQ=255)
                          gcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATc                                           >  1:2131489/1‑71 (MQ=255)
                           cgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATc                                                             <  1:300475/52‑1 (MQ=255)
                               gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCAtt                                           <  1:3134047/66‑2 (MQ=255)
                                 cAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctctt                                    <  1:3613558/71‑1 (MQ=255)
                                  aGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCg                                                                      <  1:2198133/36‑1 (MQ=255)
                                     cTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGa                                 >  1:3749195/1‑70 (MQ=255)
                                      tACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGAt                                >  1:1270227/1‑70 (MQ=255)
                                        ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                              <  1:496167/70‑1 (MQ=255)
                                        ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                              <  1:3592426/70‑1 (MQ=255)
                                        ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAAc                             >  1:2744289/1‑71 (MQ=255)
                                        ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAAc                             <  1:686166/71‑1 (MQ=255)
                                        ccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGacga                              <  1:1759982/70‑4 (MQ=255)
                                         cAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACg                            >  1:3597233/1‑71 (MQ=255)
                                          aaTCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                           <  1:2056505/71‑1 (MQ=255)
                                          aaTCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                           >  1:1142962/1‑71 (MQ=255)
                                           aTTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGcc                          <  1:448498/71‑1 (MQ=255)
                                           aTTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                           >  1:1826369/1‑70 (MQ=255)
                                            ttGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCa                         <  1:966097/71‑1 (MQ=255)
                                            tcggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                           >  1:3146126/3‑69 (MQ=255)
                                             cggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGc                                                     >  1:3469152/2‑42 (MQ=25)
                                             cggggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGc                                                     >  1:3645109/2‑42 (MQ=25)
                                               gggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                      <  1:1039489/71‑3 (MQ=255)
                                                 gTAATTTTCTGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGc                    >  1:3629754/1‑71 (MQ=255)
                                                 gTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATc                                           <  1:45937/48‑1 (MQ=255)
                                                  tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAcc                                              <  1:3684269/44‑1 (MQ=255)
                                                  tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                   <  1:1716718/71‑1 (MQ=255)
                                                   gaTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCg                  >  1:2310222/2‑71 (MQ=255)
                                                   aaTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCa                                                    >  1:920003/1‑37 (MQ=255)
                                                   aaTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg                  >  1:3548375/1‑71 (MQ=255)
                                                    aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATc                                           >  1:3104766/1‑45 (MQ=255)
                                                    aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATc                                           >  1:1823567/1‑45 (MQ=255)
                                                    aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAAGGCCAGATGCCGt                 >  1:3657593/1‑71 (MQ=255)
                                                    aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg                  <  1:3632583/70‑1 (MQ=255)
                                                    aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                 >  1:3045419/1‑71 (MQ=255)
                                                    aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                 >  1:2282745/1‑71 (MQ=255)
                                                    aTTTTCGGTCCCAAATCGCACCCAATCACATTGGCATAAACCATc                                           >  1:3530414/1‑45 (MQ=255)
                                                    aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCttt                                   >  1:1575010/1‑53 (MQ=255)
                                                    aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCttt                                   >  1:3659092/1‑53 (MQ=255)
                                                    aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                        <  1:2776004/64‑1 (MQ=255)
                                                    aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                        <  1:373726/64‑1 (MQ=255)
                                                     ttttCGGGCCTAAATCACAGCCAATCACATTGGCATaaatc                                              <  1:636744/41‑3 (MQ=255)
                                                      tttCGGTCCCAAATCGCAGCCAATCACATTGGCATaa                                                 >  1:152663/1‑37 (MQ=255)
                                                       ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct              >  1:1318215/1‑71 (MQ=255)
                                                       ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct              >  1:1834918/1‑71 (MQ=255)
                                                        tCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctg             >  1:1079007/1‑71 (MQ=255)
                                                        tCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                        >  1:1243003/1‑60 (MQ=255)
                                                        tCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtactc             >  1:1335928/1‑70 (MQ=255)
                                                         cGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtactc             <  1:1014757/70‑2 (MQ=255)
                                                          ggTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctgc            <  1:2603902/70‑1 (MQ=255)
                                                          ggTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctgc            >  1:1098596/1‑70 (MQ=255)
                                                           gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCa          <  1:3379826/71‑1 (MQ=255)
                                                           gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCa          <  1:3081292/71‑1 (MQ=255)
                                                            tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTg                                  <  1:803192/46‑1 (MQ=255)
                                                            tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTg                                  <  1:2384325/46‑1 (MQ=255)
                                                            tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCat                        <  1:3070565/56‑2 (MQ=255)
                                                            tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtg                <  1:3699881/64‑1 (MQ=255)
                                                            tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtg                <  1:1221579/64‑1 (MQ=255)
                                                             cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctg             >  1:3733789/1‑66 (MQ=255)
                                                             cccAAATCGAAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt         >  1:556/1‑70 (MQ=255)
                                                               cAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGAt                      >  1:1126735/1‑55 (MQ=255)
                                                                aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt     >  1:2797251/1‑71 (MQ=255)
                                                                 aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc               >  1:96736/1‑60 (MQ=255)
                                                                 aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc               >  1:3524820/1‑60 (MQ=255)
                                                                 aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg    >  1:1132875/1‑71 (MQ=255)
                                                                 aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg    >  1:1673561/1‑71 (MQ=255)
                                                                 aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg    >  1:899769/1‑71 (MQ=255)
                                                                 aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg    <  1:2865012/71‑1 (MQ=255)
                                                                 aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg    >  1:2405114/1‑71 (MQ=255)
                                                                 aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg    >  1:2105044/1‑71 (MQ=255)
                                                                 aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAtc    >  1:453421/1‑70 (MQ=255)
                                                                  aTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCa          >  1:3412466/1‑64 (MQ=255)
                                                                    cGCAGCCAATCACATTGGCATAAACCATCGGTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa  >  1:810901/1‑70 (MQ=255)
                                                                     |                                                                    
CTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGA  >  minE/2595422‑2595559

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: