This is a test version of Biostars. For the public version, visit https://www.biostars.org.
frame frequency analysis per position per gene

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 of the codon) per gene. do you know how I can do that?

sequencing
• 920 views • link
written by alireza346 1

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • how treat FLAG 16 in ribosome footprint
    written by dwc06631 0
    hello. first, sorry for that I'm not good at this field. I aligned reads to reference, and now I'm finding each read's start location from …
  • single base resolution count at the genome level
    written by alireza346 1
    I am trying to get the single base resolution bedgraph file at the genome level . do you know how to do that? in fact, …
  • isolate the output of lignment in fastq format
    written by alireza346 1
    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 …
  • making the gene specific bam file
    written by alireza346 1
    I have RNAseq data and aligned the fastq files to the genome. so now I have bam file. I want to make the bam file …
  • extracting the longest transcript in a fasta file
    written by alireza346 1
    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 …
  • frequency of each nucleotide at each position in sequence
    written by Sara 28
    I have aligned RNAseq data and now I have bam files for control and treated conditions. I am trying to get the frequency of each …
  • Best way to normalize occurance of a nucleotide subsequence
    written by snamjoshi87 4
    If I have a set of nucleotide sequences of varying lengths and I want to find the number of times a certain subsequences occurs for …
  • CDS read counts per gene
    written by ashkan 16
    Hi Guys, I have aligned RNA-seq data and have bam files. I am goint to coundt the number of reads per gene but I need …
  • how to get read density for the first 1000 nt of each transcript
    written by Sara 28
    Hi, I have 2 RNA-seq files for two samples and aligned them to the reference genome so I have BAM files. I would like to …
  • Counting reads per codon
    written by fire_water 8
    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. …