Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,458 G→A 12.5% S380S (AGC→AGT arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,4580GA12.5% 288.8 / 24.8 112S380S (AGC→AGTarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (52/46);  new base A (8/6);  total (60/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 118 total aligned reads displayed.

GGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGA  >  minE/2595388‑2595528
                                                                      |                                                                      
ggAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAgg                                                                        >  1:2068162/1‑71 (MQ=255)
  aaGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGc                                                                       >  1:1061447/1‑70 (MQ=255)
  aaGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGc                                                                       >  1:4607899/1‑70 (MQ=255)
   aGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACtg                                                                     >  1:2813062/1‑70 (MQ=255)
   aGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACt                                                                      <  1:1140045/70‑1 (MQ=255)
   aGTAATATCCCCAAGTGATTGTCATATTTTTCAGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACt                                                                      <  1:878483/70‑1 (MQ=255)
    gTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACtgc                                                                    <  1:1763288/71‑3 (MQ=255)
     tAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACtgc                                                                    >  1:3296041/1‑68 (MQ=255)
       atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAgg                                                                        >  1:292915/1‑64 (MQ=255)
       atatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCa                                                                  <  1:5022952/70‑1 (MQ=255)
        tatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCATAGTAGCGTAGCCAGGCTACCaa                                                                 <  1:4347204/70‑1 (MQ=255)
        tatCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAAt                                                                <  1:2615888/71‑1 (MQ=255)
         atCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAAtt                                                               >  1:2975719/1‑71 (MQ=255)
          tCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTg                                                              >  1:3338200/1‑71 (MQ=255)
           ccccAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTgg                                                             >  1:4586071/1‑71 (MQ=255)
            cccAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTggg                                                            >  1:4068552/1‑71 (MQ=255)
            cccAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTggg                                                            >  1:1531737/1‑71 (MQ=255)
                gCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTaa                                                        <  1:1272417/71‑1 (MQ=255)
                gCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTaa                                                        <  1:2839072/71‑1 (MQ=255)
                  tGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCTAGGCTACCa                                                                  <  1:2456311/59‑1 (MQ=255)
                    aTAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAAtt                                                      <  1:4232981/69‑1 (MQ=255)
                      aGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAAt                                                                >  1:4653318/1‑57 (MQ=255)
                      aGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCg                                                  >  1:3891123/1‑71 (MQ=255)
                      aGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCg                                                  >  1:3583917/1‑71 (MQ=255)
                       gTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCg                                                  <  1:5434796/70‑1 (MQ=255)
                        tCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCggg                                                <  1:1234240/71‑2 (MQ=255)
                          atatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTc                                                   >  1:61438/1‑66 (MQ=255)
                          atatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTc                                                   >  1:3968983/1‑66 (MQ=255)
                          atatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTc                                                   >  1:866379/1‑66 (MQ=255)
                          atatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTc                                                   >  1:2489482/1‑66 (MQ=255)
                          atatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCg                                                  >  1:1747697/1‑67 (MQ=255)
                          atatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTcc                                              >  1:3251462/1‑71 (MQ=255)
                          atatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTcc                                              >  1:3141656/1‑71 (MQ=255)
                          atatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTcc                                              >  1:4548297/1‑71 (MQ=255)
                           tatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTcc                                              >  1:4974575/1‑70 (MQ=255)
                           tatTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTcc                                              >  1:3336725/1‑70 (MQ=255)
                             ttcttcTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCaa                                           <  1:5141489/71‑1 (MQ=255)
                              tcttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCa                                                                  <  1:1290473/47‑1 (MQ=255)
                               cttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCaaa                                          <  1:4285372/70‑1 (MQ=255)
                                ttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGc                                                                       <  1:1569405/40‑1 (MQ=255)
                                ttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGc                                                                       <  1:1002196/40‑1 (MQ=255)
                                ttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGc                                                                       >  1:1011515/1‑40 (MQ=255)
                                ttctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATc                                        >  1:696216/1‑71 (MQ=255)
                                 tctGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCg                                       >  1:855357/1‑71 (MQ=255)
                                    gCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAg                                    >  1:2687148/1‑71 (MQ=255)
                                    gCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAg                                    >  1:3120286/1‑71 (MQ=255)
                                     cGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCt                                                                      >  1:4534146/1‑36 (MQ=255)
                                     cGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGc                                      <  1:3526813/68‑1 (MQ=255)
                                     cGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGc                                      <  1:2964023/68‑1 (MQ=255)
                                      gAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGc                                   >  1:3121949/1‑70 (MQ=255)
                                       aaaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCa                                 <  1:1807966/71‑1 (MQ=255)
                                        aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTa                                                         >  1:2402409/1‑46 (MQ=255)
                                        aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCgg                                                 >  1:1273099/1‑54 (MQ=255)
                                        aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCg                                                  >  1:3799103/1‑53 (MQ=255)
                                        aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCaa                                <  1:3251813/71‑1 (MQ=255)
                                        aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCaa                                <  1:4062936/71‑1 (MQ=255)
                                        aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCaa                                <  1:1582911/71‑1 (MQ=255)
                                        aaGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACAAATTGGGGTAATTTTc                                                   >  1:1483323/1‑52 (MQ=255)
                                         aGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAAt                               >  1:3615097/1‑71 (MQ=255)
                                         aGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAAt                               >  1:4120768/1‑71 (MQ=255)
                                          gTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCaa                                                                 >  1:1796600/1‑36 (MQ=255)
                                          gTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCaa                                                                 >  1:4497734/1‑36 (MQ=255)
                                           tACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCaa                                                                 >  1:3714587/1‑35 (MQ=255)
                                           tACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAAt                                                       <  1:5153649/45‑1 (MQ=255)
                                           tACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATca                             >  1:3863892/1‑71 (MQ=255)
                                           tACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATca                             <  1:2751941/71‑1 (MQ=255)
                                           tACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATca                             <  1:2745415/71‑1 (MQ=255)
                                            acatGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATca                             >  1:1514912/4‑70 (MQ=255)
                                              gTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACAt                          <  1:4778729/71‑1 (MQ=255)
                                              gTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACAt                          <  1:3790476/71‑1 (MQ=255)
                                              atGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACAt                          <  1:3814260/70‑1 (MQ=255)
                                               tGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACAt                          >  1:665170/1‑70 (MQ=255)
                                                gCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTg                        >  1:528571/1‑71 (MQ=255)
                                                  aGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTgg                       >  1:2543930/1‑70 (MQ=255)
                                                     cAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                   <  1:1972875/71‑1 (MQ=255)
                                                      agagTAGCGTAGCCAGGCTACCAATTGGGGTAAttt                                                     >  1:3833663/1‑36 (MQ=255)
                                                      agagTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATaa                  <  1:4839785/71‑1 (MQ=255)
                                                      agagTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                   >  1:2716251/1‑70 (MQ=255)
                                                        agTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAc                <  1:3884190/71‑1 (MQ=255)
                                                          tagcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCa                                            <  1:698790/41‑1 (MQ=255)
                                                          tagcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCa                                            <  1:892936/41‑1 (MQ=255)
                                                          tagcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCa                                            <  1:507659/41‑1 (MQ=255)
                                                               ttgcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATcac                            <  1:5693837/50‑1 (MQ=255)
                                                               ttgcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCAt             >  1:308917/3‑67 (MQ=255)
                                                               ttgcCAGACTGCCAACCGGGGTGATTTTCGGGCCTAAATCACAGCCAATcac                            <  1:5169098/50‑1 (MQ=255)
                                                               tagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATcac                            <  1:3148925/52‑1 (MQ=255)
                                                               tagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATcac                            <  1:3109399/52‑1 (MQ=255)
                                                               tagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCAt             <  1:5261264/67‑1 (MQ=255)
                                                               tagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCAt             <  1:2479251/67‑1 (MQ=255)
                                                                agcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGct         >  1:117541/1‑70 (MQ=255)
                                                                agcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAACCGCAGCCAATCACATTGGCATAAACCATCGctt        >  1:1749143/1‑71 (MQ=255)
                                                                 gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTg                        >  1:830136/1‑54 (MQ=255)
                                                                 gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctc       >  1:4111743/1‑71 (MQ=255)
                                                                 gcCAGACTGCCAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcctc       >  1:138197/1‑71 (MQ=255)
                                                                  ccAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAAcc               <  1:1278617/62‑1 (MQ=255)
                                                                   cAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATaaa                 >  1:725078/1‑59 (MQ=255)
                                                                    aGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCttt    <  1:1136169/71‑1 (MQ=255)
                                                                      gCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                   <  1:5621723/54‑1 (MQ=255)
                                                                      gCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATa                   <  1:5274453/54‑1 (MQ=255)
                                                                      gCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGa  >  1:44154/1‑71 (MQ=255)
                                                                      |                                                                      
GGAAGTAATATCCCCAGCTGATAGTCATATTCTTCTGCGAAAGTACGTGCAGCCAGAGTAGCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGA  >  minE/2595388‑2595528

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: