I have a RNA-seq data and want to get the proportion of reads that map to a window of 100pb between 5'UTR and CDS (50 bp at 5'UTR and 50 bp at the CDS). since I am using Tophat to align and bowtie to build indexes, I guess I have to get the sequence of that part to build indexes. do you guys know how to find the sequence of that part?
•
written
by
alirezamomeni707 •
1 answer
I would advise against creating a custom genome index that only contains the intervals of interest. The alignments are most likely to be negatively affected.
Instead, you should align your reads against the full genome then intersect these resulting alignment with a custom interval file that contains only the regions of interests.
Tools like bedtools or bedops will allow you to do both, create upstream regions relative to the intervals of interest as well as to intersect the alignment files.
Log in to answer this question.
-
Mapping sequencing reads to 5-UTR, 3-UTR and CDS using Ensembl GTF files
written by Beginner •
-
counting the reads that map to a special window
written by Christleemans •
-
isolating sequence from a specific site of the genes
written by alireza346 •
-
counting the reads that mapped to a part of transcriptome
written by Sara
-
Transcriptomic features of peaks at mature mRNA level
written by Ezhil La •
-
problem with tRNAseq alignment
written by alirezamomeni707 •
-
alignment of sequencing data to the tRNA transcriptome
written by alirezamomeni707 •
-
Aligning reads to genome/transcriptome
written by EVR
-
ERCC and hg19 Bowtie Indexes
written by bgraphit •
-
Tophat Error IOError: [Errno 13] Permission denied
written by punjabdhaputar •
You have BAM files and a GTF or something like that already, yes?
yes I have bam file. I can also download GTF FILE.