Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Time frame for BLASTN against NT jackfrost Bioinformatics 5 06-11-2015 12:08 PM
Matching a column from a data frame on the columns of another data frame zisis86 Bioinformatics 2 05-29-2014 05:32 AM
six frame translation Alagurajvelu Bioinformatics 1 08-12-2013 07:51 AM
When is Open reading frame=gene? ritzriya RNA Sequencing 4 10-06-2010 09:10 PM

Thread Tools
Old 03-09-2018, 05:51 AM   #1
Junior Member
Location: Utrecht

Join Date: Jan 2018
Posts: 1
Question Keep in-frame reads only

Hi all,

I have been working with Ribosome Sequencing datasets over the past few months. My ultimate gole is to determine the codon that's in the P-site of the Ribosome Protected Fragements. These fragments are roughly 29-30 nt in length. The P-site is located at position 12 (= nts 13, 14, 15).

We've gotten pretty far: standard adapter trimming, contaminant removal and alignment to the reference genome. Because for 'translating' the codon in the P-site to an amino acid it is very important that the read is in-frame (the first three nucleotides code for an amino acid), I would like to filter out any non-aligned or ambigious-frame reads from the alignment output (BAM).

Output preferably in fastq format, since this is compatible with the tool we use to map the P-site codons to amino acids.

Does anyone know of a tool or a way to do this?

Thanks a lot!

All the best,


Last edited by Flock; 03-09-2018 at 08:12 AM.
Flock is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 08:34 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO