Counting reads per codon
Hello. So I have some ribosome profiling sequencing data which I've clipped, trimmed, mapped with Bowtie, and converted to BAM files for viewing in IGV. The data looks good in IGV but now I want to know the number of reads per codon. Does a tool exist for doing that?
Thanks!
• 1,913 views
•
link
updated
by
Jorge Amigo
•
written
by
fire_water
•
written
by
fire_water
1 answer
none that I'm aware of, but codon information can be considered as any other annotation, so if you have all the codons you're interested in described in a bed file you can obtain such information using bedtools intersect for instance:
bedtools intersect -c -a codons.bed -b reads.bam
• 0 views
•
link
written
by
Jorge Amigo
Log in to answer this question.
More posts like this
-
STAR alignment speed
written by Gilles •
-
miRDeep2 - identifying miRNAs from deep sequencing data
written by Hasan •
-
CHIP-Seq data analysis
written by prs •
-
STAR aligner - multimaps vs. strandedness
written by cjag •
-
Reads that map to two different genomes, but with different map coordinates
written by Susmita Mandal
-
Using reads aligned on the Negative strand of the genome to select soft clipped nucleotides
written by fusion.slope
-
Mapping Ribosome Profiling Data
written by fire_water
-
low number mapped reads in ribosome profiling
written by laura.fancello •
-
Mapping Start Codons With Ribosome Profiling
written by beauripert90 •
-
Tajima'S D Test-Ngs Data
written by jackuser1979
