I have a bed file of some custom regions I have interest in, e.g.
chr1 1000 1015 my_region 0 + ...
I believe these regions may be translated, but are not necessarily annotated genes. So I would like to map some reads to these regions to confirm my suspicions.
The dataset GSE109313 has 1 Trillion reads from various human tissues.
Clicking on the run selector I can see the 37 SRR samples.
Using
prefetch
I downloaded these files.
Often reads are mapped to a GTF file, e.g. the human GTF file can be found on ensembl.
So now I have bed file or regions I would like to map to, the human GTF file, and the ~400Gb of SRR files.
Where do I go from here?
•
written
by
Mag.ds •
1 answer
As far as I know reads are never mapped to a GTF (or GFF) file, they are always mapped to a fasta file (or an indexed version of a fasta file), so that's a file with the actual sequence in it. The GTF/GFF do come into play when you want to link your mappings to annotated features (eg. genes, promoters, lncRNA, ... ) .
The best you can do is to map all your reads to the human genome.
As an alternative, you could subselect the human genome to only contain the regions you might be interested in. Based on you bed file you could do this with eg getFasta from the bedtools package . DO keep in mind that this has a huge potential in biasing your analysis.
no worries, actually there is no need to convert your bed to GFF or GTF, follow genomax suggestion on how to do this without converting your bed file.
Log in to answer this question.
-
Struggling to download using SRA toolkit from an Accession list
written by SethJ •
-
Same GEO Accession, different SRR number, how to download this RNA-seq paired-end data?
written by ev97 •
-
retrieving paired end sequencing data with fasts-dump
written by v.berriosfarias •
-
ChIP-seq Peak calling from BED or WIG file
written by gspirito •
-
Getting Paired end SRA
written by ovariohisterectomia •
-
Mapping multiple SRR ID from same experiment
written by singlecell •
-
Mapping short 26-32nt RPF reads to human genome
written by Adrian Pelin
-
NCBI GEO/GSM/SRA conventions: Is raw data a part of the SRA or a supplementary file?
written by abe •
-
Downloaing SRA using Aspera
written by always_learning
-
Tutorial: How to download raw sequence data from GEO/SRA
written by Obi Griffith
Map to the whole genome and then look at regions of your interest. Don't do it any other way, especially if original data is from entire genome.
so for clarity, use bowtie on file with hg38 and then what do you suggest to for looking at regions of my interest
You can use
mosdepth(LINK) orbedtools coverage(LINK) to look at your regions of interest.