Extract Negative chromosome(gene) location from VCF files
1
0
Entering edit mode
8.0 years ago
MAPK ★ 2.1k

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"
VCF extract vcf chromosome location gene location • 1.9k views
ADD COMMENT
0
Entering edit mode
8.0 years ago

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 COMMENT

Login before adding your answer.

Traffic: 1960 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6