isolating sequence from a specific site of the genes
I am trying to get the sequence at both sides of the coding sequence stop site for all of the genes for the alignment. 40 nt before stop site (in the CDS) and 50 nt at the downstream of CDS stop site(in the 3p UTR) for all of the genes. do you know how I can do that correctly
• 1,283 views
•
link
written
by
alireza346 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
extracting gencode 3utr and 5utr as part of R pipeline
written by Jalil Sharif •
-
How to locate and extract data from a df in R
written by margo •
-
SnpEff START/STOP codon not found error
written by thomasbersez
-
isolating genomic sequence based on RNAseq data
written by alireza346 •
-
coverage of the reads over coding sequence only for one gene
written by alireza346 •
-
counting the reads that map to a special window
written by Christleemans •
-
Transcriptomic features of peaks at mature mRNA level
written by Ezhil La •
-
alignment to a special part
written by alirezamomeni707 •
-
How to extract whole gene sequence with 5' and 3' UTR from genome
written by Bioinfonext
-
Conservation Of Coding And Non Coding Exons
written by Puthier
It's unclear which data format you have.
I have RNAseq and want to align that to this part of transcriptome.
That is important information you should have mentioned in your question. Why do you think this is a good idea?
Anyway, to get this done you would first need a GTF/GFF of your organism and isolate the stop sites. Then you would make a bed file with -40nt and + 50nt intervals for each stop site. Then you could use
bedtools getfastato get the nucleotide sequence of those intervals from the reference genome.