Question: Extract Negative chromosome(gene) location from VCF files
gravatar for MAPK
4.2 years ago by
MAPK1.5k wrote:

I have about 3000 genes list and their chromosome locations(start and end). Some of the locations have negative positions as you can see below in the table. I think this mean these genes are position referenced from negative strand. I need to extract these genes from VCF file, but there aren't any negative location in VCF file. Can someone please explain me how this issue could be addressed when extracting these gene locations from VCF files.

   gene_id symbol chromosome start_location end_location
1  "1"     "A1BG" "19"       " -58858171"   " -58864865"
2  "10"    "NAT2" "8"        "  18248754"   "  18258723"
3  "100"   "ADA"  "20"       " -43248162"   " -43280376"
4  "1000"  "CDH2" "18"       " -25530926"   " -25616549"
ADD COMMENTlink modified 4.2 years ago by Devon Ryan95k • written 4.2 years ago by MAPK1.5k
gravatar for Devon Ryan
4.2 years ago by
Devon Ryan95k
Freiburg, Germany
Devon Ryan95k wrote:

The negative positions are typically used to indicate the - strand. Yes this is incredibly unintuitive. So " -58858171" " -58864865" should be 58858171-58864865 on the - strand.

ADD COMMENTlink written 4.2 years ago by Devon Ryan95k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1212 users visited in the last hour