SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
SAM flag forward/reverse first/second nupurgupta Bioinformatics 3 01-20-2016 03:47 AM
rmdup can not move duplicates in forward and reverse strand for single-end reads ct586 Bioinformatics 4 03-11-2012 05:01 PM
Separate forward and reverse coverage in Artemis nivea Bioinformatics 4 12-12-2011 03:20 AM
how to define a forward or reverse read file poorphd Illumina/Solexa 3 11-22-2011 12:34 PM
forward and reverse sequance asankaf General 5 05-27-2009 07:48 AM

Reply
 
Thread Tools
Old 02-11-2013, 10:21 AM   #1
says_anova
Member
 
Location: California

Join Date: May 2011
Posts: 10
Default samtools mpileup forward vs. reverse strand?

Hi all,

Apologies if this has been addressed in an earlier post, I'm definitely a n00b. Just need a sanity check:

Working with RNASeq data, the mpileup output of reads assigns '.' to a reference match on the "forward strand" and a ',' for a reference match on the "reverse strand."

In this case, the "forward" and "reverse" terminology doesn't reflect the double strandedness of DNA, but the cDNA strands being sequenced (as complements to the mRNA from the sample), where the "forward" and "reverse" tells you something about primer alignment and sequencing.

Do I have that right??

Thanks in advance!
says_anova is offline   Reply With Quote
Old 02-12-2013, 09:14 AM   #2
says_anova
Member
 
Location: California

Join Date: May 2011
Posts: 10
Default

Can anyone out there confirm or deny?
says_anova is offline   Reply With Quote
Old 02-12-2013, 10:36 AM   #3
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Honestly, I don't get your question. And bumping it won't fix that, or make it more attractive to people who might answer it.

Unless you do a strand-specific RNA prep, you should have half forward reads, half reverse reads, at a given locus, even though your RNA is obviously stranded.
swbarnes2 is offline   Reply With Quote
Old 02-12-2013, 11:29 AM   #4
says_anova
Member
 
Location: California

Join Date: May 2011
Posts: 10
Default

swbarnes2 - Sorry if the question wasn't clear. i think I'm confused on some basic mol. bio.

When you say half forward and half reverse reads, what are they forward and reverse with respect to?

Thanks in advance.
says_anova is offline   Reply With Quote
Old 02-12-2013, 01:02 PM   #5
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Quote:
Originally Posted by says_anova View Post
swbarnes2 - Sorry if the question wasn't clear. i think I'm confused on some basic mol. bio.

When you say half forward and half reverse reads, what are they forward and reverse with respect to?

Thanks in advance.
In respect to your reference fasta. If your reference fasta says

Quote:
AGTAGTATAGGTA
and the fastq has a read like:

Quote:
AGTAGTATAGGTA
That's forward. If the fastq had this as the sequence:

Quote:
TACCTATACTACT
That aligns in the reverse direction.
swbarnes2 is offline   Reply With Quote
Reply

Tags
mpileup, rnaseq, samtools

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 01:47 PM.


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