single base resolution count at the genome level
I am trying to get the single base resolution bedgraph file at the genome level . do you know how to do that? in fact, I have DNAseq data and trying to get the read count per base per gene.
sequencing
• 572 views
•
link
written
by
alireza346 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Single-base bins for finding enrichment at promoters
written by sobethon •Using R, how would I be able to create single-base bins from a GenomicRanges object (transcripts/promoters) of the hg19 genome? I have a BAM file …
-
read counts per chromosome arm
written by SaraI have `DNAseq` data and trying to get the read `counts per chromosome arm`. but I do not know the correct way. shall I align …
-
isolating genomic sequence based on RNAseq data
written by alireza346 •using `RNAseq` data I have the coordinates of `5p end` of all expressed transcripts and I want to get the the sequence of that site …
-
isolate the output of lignment in fastq format
written by alireza346 •I have `RNAseq` data and trying to align the `fastq` files to human `transcriptome` (not `genome`). I want to isolate the reads that map to …
-
frame frequency analysis per position per gene
written by alireza346 •I have ribosome profiling data and I have aligned the files. I want to get the frame frequency for each codon (basically for each position …
-
coverage of the reads over coding sequence only for one gene
written by alireza346 •I have `RNAseq` data and trying to get the coverage of the reads over coding sequence but only for one gene for example `GAPDH`. and …
-
Using Samtools depth to get number of forward and reverse reads per base from a .bam file
written by jspainhour3 •I am trying to build a tab delimited table of count and read data from a deduplicated realigned bam file. I am looking at mitochondrial …
-
extracting the longest transcript in a fasta file
written by alireza346 •I have a fasta file like this example: >ENST00000511961.1|ENSG00000013561.13|OTTHUMG00000129660.5|OTTHUMT00000370661.3|RNF14-003|RNF14|278 MSSEDREAQEDELLALASIYDGDEFRKAESVQGGETRIYLDLPQNFKIFVSGNSNECLQNSGFEYTICFLPPLVLNFELPPDYPSSSPPSFTLSGKWLSPTQLSALCKHLDNLWEEHRGSVVLFAWMQFLKEETLAYLNIVSPFELKIGSQKKVQRRTAQASPNTELDFGGAAGSDVDQEEIVDERAVQDVESLSNLIQEILDFDQAQQIKCFNSKLFLCSICFCEKLGSECMYFLECRHVYCKACLKDYFEIQIRDGQVQCLNCPEPKCPSVATPGQ >ENST00000506822.1|ENSG00000013561.13|OTTHUMG00000129660.5|OTTHUMT00000370662.1|RNF14-004|RNF14|132 MSSEDREAQEDELLALASIYDGDEFRKAESVQGGETRIYLDLPQNFKIFVSGNSNECLQNSGFEYTICFLPPLVLNFELPPDYPSSSPPSFTLSGKWLSPTQLSALCKHLDNLWEEHRGSVVLFAWMQFLKE >ENST00000513019.1|ENSG00000013561.13|OTTHUMG00000129660.5|OTTHUMT00000370663.1|HAS-0|HAS|99 MSSEDREAQEDELLALASIYDGDEFRKAESVQGGETRIYLDLPQNFKIFVSGNSNECLQNSGFEYTICFLPPLVLNFELPPDYPSSSPPSFTLSGKWLS >ENST00000356143.1|ENSG00000013561.13|OTTHUMG00000129660.5|-|HAS-202|HAS|474 MSSEDREAQEDELLALASIYDGDEFRKAESVQGGETRIYLDLPQNFKIFVSGNSNECLQNSGFEYTICFLPPLVLNFELPPDYPSSSPPSFTLSGKWLSPTQLSALCKHLDNLWEEHRGSVVLFAWMQFLKEETLAYLNIVSPFELKIGSQKKVQRRTAQASPNTELDFGGAAGSDVDQEEIVDERAVQDVESLSNLIQEILDFDQAQQIKCFNSKLFLCSICFCEKLGSECMYFLECRHVYCKACLKDYFEIQIRDGQVQCLNCPEPKCPSVATPGQVKELVEAELFARYDRLLLQSSLDLMADVVYCPRPCCQLPVMQEPGCTMGICSSCNFAFCTLCRLTYHGVSPCKVTAEKLMDLRNEYLQADEANKRLLDQRYGKRVIQKAL I m trying to take the longest transcript of …
-
How to merge bedGraph records which are next to each other and have the same score?
written by James AshmoreIs anyone aware of software to merge bedGraph records if the score is the same? I have a bedGraph file calculated at single base pair …
-
coding sequnce read coun
written by ashkanHi Guys, I asked a question yesterday but did not get my answer and I really need it. I have RNA-seq data and its bam …
Hello alireza346!
We believe that this post does not fit the main topic of this site.
Please use google and the search function. This has been asked many times before. Look for genomecov (bedtools) or bamCoverage (deeptools). For the latter, check the -bs parameter.
For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.
If you disagree please tell us why in a reply below, we'll be happy to talk about it.
Cheers!
PER_BASE_COVERAGEparameter