RNA-seq analysis


 
Thread Tools Search this Thread
Top Forums UNIX for Dummies Questions & Answers RNA-seq analysis
Prev   Next
# 1  
Old 02-11-2013
RNA-seq analysis

I am processing RNA-seq data files that have been aligned using RUM. One of the output files is a *.sam that includes:

Unique alignments
Non-unique alignments
original read files

I want to extract only the unique alignments by pulling out alignments that have "IH:i:1" (indicates this read had only one alignment).

I have tried.... grep "IH:i:1" file.sam > filtered.sam

but this will also return "IH:i:11 IH:i:12" etc etc. I have also tried "IH:i:1 " which returns nothing (I believe it's tab delimited)

Any suggestions?
 
Login or Register to Ask a Question

Previous Thread | Next Thread

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

Easy seq Question

Hi! I'm trying to do this: 1 - 2 - 3 - 4 - 5 - I'm using seq for this: seq 1 20 > filename.txt How do I get the "-"? I've tried -f per man but can't get anything to work. Also, is there an easier or better way than using sequence? Thanks! (6 Replies)
Discussion started by: TonyBe
6 Replies

2. Shell Programming and Scripting

awk command to find seq lines

I have a below file FILE.cfg JAN_01 VAR1=4 VAR2=SUM VAR3=PRIVATE JAN_10 VAR1=44 VAR2=GUN VAR3=NATUR JAN_20 VAR1=3 VAR2=TQN VAR3=COMMA code: (JAN_10 is argument passed from script) (6 Replies)
Discussion started by: Roozo
6 Replies

3. Shell Programming and Scripting

How to show first 0 using seq or +1 count?

Greetings, Using linux based OS and KSH. I m trying to make a simple script to parse some logs to show a count per hour on a specific alarm starting from midnight to the current hour. So I format my "HOUR" variable to show the current time and so I can use it in the following bit of code.... (6 Replies)
Discussion started by: Sekullos
6 Replies

4. Shell Programming and Scripting

sorting a fixed width seq file

I have a file like this... 2183842512010-11-25 15379043 453130325 2386225062010-11-30 4946518 495952336 2386225062010-11-30 4946518 495952345 2386225062010-11-25 262066688 -516224026 2679350512010-11-25 262066688 -516224124 3196089062010-11-25 262066688 203238229... (5 Replies)
Discussion started by: issaq84mohd
5 Replies

5. UNIX for Dummies Questions & Answers

Help with seq (print a series of dates)

Assuming one does not have such luxuries as bash, zsh, jot, rs, perl, etc. what is the most elegant way to print out a formatted date series like this: 01-01-2010 01-02-2010 01-03-2010 ... 02-01-2010 02-02-2010 ... Can I accomplish this with just basic shell builtins and seq, or... (3 Replies)
Discussion started by: uiop44
3 Replies

6. Shell Programming and Scripting

Using seq (Or alternative)

I usually just browse the forum/google for answers, however I've been stuck on a problem for a number of hours now and I've decided to join up and actually ask I've searched the forum ad naseum in an attempt to find answer to my query, however so far I have been unsuccessful. I'm no expert... (3 Replies)
Discussion started by: gtc
3 Replies

7. Shell Programming and Scripting

Using Seq As A Variable With Padded Digits

Hi all. Im trying to use a sequence in a while loop like this below. I need it for navigating a year, month, day folder structure where a user can input the start date and have it go to the desired end date. The script will grab a certain file on each day then move onto the next. Ive got all that... (3 Replies)
Discussion started by: Grizzly
3 Replies

8. Shell Programming and Scripting

declaring variable with for $(seq)

Hi guys. i have the following script: 1 #!/bin/bash 2 linkcount=$(grep "/portal" tickets | wc -l) 3 grep "/portal" tickets > links 4 for i in $(seq 1 $linkcount); do 5 echo "BLYAT" 6 let link$i=$(sed -n "$i"p links) 7 echo $ 8 done the problem is, that "let" can`t... (1 Reply)
Discussion started by: neverhood
1 Replies

9. Shell Programming and Scripting

Move and rename files in seq. with padded digits

Greetings, I am new to scripting, but find if I can see the code working for a given problem, then I can eventually figure it out. (9 Replies)
Discussion started by: rocinante
9 Replies

10. Shell Programming and Scripting

script to loop and check jumping seq.

Hi, Normally, I will manually to use "ll" command to list the following file from \FILE\CACHE\ directory and check the jump seq. Can I write a script to loop or/and check jump seq file (if jumped seq and show "missing seq no" message for me) -rw-rw----+ 1 user develop 14012 Sep 4... (1 Reply)
Discussion started by: happyv
1 Replies
Login or Register to Ask a Question
Bio::SeqFeature::SiRNA::Oligo(3pm)			User Contributed Perl Documentation			Bio::SeqFeature::SiRNA::Oligo(3pm)

NAME
Bio::SeqFeature::SiRNA::Oligo - Perl object for small inhibitory RNAs. SYNOPSIS
use Bio::SeqFeature::SiRNA::Oligo; my $oligo = Bio::SeqFeature::SiRNA::Oligo-> new( -seq => 'AUGCCGAUUGCAAGUCAGATT', -start => 10, -end => 31, -strand => -1, -primary => 'SiRNA::Oligo', -source_tag => 'Bio::Tools::SiRNA', -tag => { note => 'A note' }, ); # normally two complementary Oligos are combined in an SiRNA::Pair # object $pair->antisense($oligo); DESCRIPTION
Object methods for single SiRNA oligos - inherits Bio::SeqFeature::Generic. Does not include methods for designing SiRNAs - see Bio::Tools::SiRNA for that. SEE ALSO
Bio::Tools::SiRNA, Bio::SeqFeature::SiRNA::Pair. FEEDBACK
Mailing Lists User feedback is an integral part of the evolution of this and other Bioperl modules. Send your comments and suggestions preferably to the Bioperl mailing list. Your participation is much appreciated. bioperl-l@bioperl.org - General discussion http://bioperl.org/wiki/Mailing_lists - About the mailing lists Support Please direct usage questions or support issues to the mailing list: bioperl-l@bioperl.org rather than to the module maintainer directly. Many experienced and reponsive experts will be able look at the problem and quickly address it. Please include a thorough description of the problem with code and data examples if at all possible. Reporting Bugs Report bugs to the Bioperl bug tracking system to help us keep track of the bugs and their resolution. Bug reports can be submitted via the web: https://redmine.open-bio.org/projects/bioperl/ AUTHOR
Donald Jackson (donald.jackson@bms.com) APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _ new Title : new Usage : my $sirna_oligo = Bio::SeqFeature::SiRNA::Oligo->new(); Function : Create a new SiRNA::Oligo object Returns : Bio::Tools::SiRNA object Args : -seq sequence of the RNAi oligo. Should be in RNA alphabet except for the final TT overhang. -start start position -end end position -strand strand -primary primary tag - defaults to 'SiRNA::Oligo' -source source tag -score score value -tag a reference to a tag/value hash -seq_id the display name of the sequence -annotation the AnnotationCollectionI object -location the LocationI object Currently passing arguments in gff_string or gff1_string is not supported. SiRNA::Oligo objects are typically created by a design algorithm such as Bio::Tools::SiRNA seq Title : Seq Usage : my $oligo_sequence = $sirna_oligo->seq(); Purpose : Get/set the sequence of the RNAi oligo Returns : Sequence for the RNAi oligo Args : Sequence of the RNAi oligo (optional) Note : Overloads Bio::SeqFeature::Generic seq method - the oligo and parent sequences are different. Note that all but the last 2 nucleotides are RNA (per Tuschl and colleagues). SiRNA::Pair objects are typically created by a design algorithm such as Bio::Tools::SiRNA. perl v5.14.2 2012-03-02 Bio::SeqFeature::SiRNA::Oligo(3pm)