BAM file download for just a few genes
1
0
Entering edit mode
7 weeks ago
abdi.zhr26 • 0

Hi all,

I want to download the BAM file just for few genes and because Don’t need the whole data of patients, I can download it by my interest gene position. I tried to find my gene position in the UCSC and RefSeq, but they showed me different positions -very small difference- (I set my Reference genome on hg19 in both sites). Now I don’t know which one should I choose for download?

Could I randomly download 200bp upstream and downstream of my gene position? Does anyone have the same experience?

RefSeq BAM position RNA-Seq UCSC • 304 views
ADD COMMENT
0
Entering edit mode

That is most likely going to mess up your analysis. The reads in an RNA-seq experiment are not independent of each other. It is the composition of the entire library that determines the distribution of reads across all genes. Generally you always want to download the entire file, analyze it as a whole, and then later filter the results for genes of interest. I strongly suggest to stick to that advise (and standard workflows) unless you have a good idea what you're doing.

ADD REPLY
0
Entering edit mode

Thank you for your kind suggestion.

ADD REPLY
2
Entering edit mode
7 weeks ago

samtools view is capable of pointing to a remote file and download only reads of interest. I used in the past to retrieve mitochondrial reads from 1000g bams instead of downloading the whole bunch of bam files available and extracting them locally, and it worked flawlessly:

samtools view -b http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/phase1/data/HG00096/alignment/HG00096.mapped.ILLUMINA.bwa.GBR.low_coverage.20101123.bam MT > HG00096.mt.bam
ADD COMMENT

Login before adding your answer.

Traffic: 2090 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