SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Samtools "is recognized as '*'" "truncated file" error axiom7 Bioinformatics 3 11-26-2014 02:53 AM
Inquery: not-stranded transcripts ("." in strand column) in cufflinks output gtf file sunfuhui Bioinformatics 0 11-15-2012 02:43 PM
pileup file symbols "<" and ">" geraldgsw Bioinformatics 3 12-30-2011 08:19 AM
The position file formats ".clocs" and "_pos.txt"? Ist there any difference? elgor Illumina/Solexa 0 06-27-2011 07:55 AM
Help: what a "*" indicates in the reads base column of pileup format? rainbow Illumina/Solexa 1 01-29-2010 03:59 AM

Reply
 
Thread Tools
Old 03-21-2014, 09:18 AM   #1
slp
Junior Member
 
Location: NY, USA

Join Date: Dec 2010
Posts: 9
Default Dealing with ">" and "<" signs in mapping column in RNA-Seq pileup file for RNA-edit

I have received a TopHat aligned bam file from sequencing core. This was for paired-end RNA-Seq and I wish to do the RNA-editing analysis. I am using SAMTools/bcf for calling the variants. But there seems to be a problem. When I look at the pileup file for the samples i see a lot of ">>>" and "<<<" symbols in the mapping column (apart from standard A ; C ; T ; G ; a ; c ; t ; g ; . ; , $ and ^S). I guess there are just the intervals between the reads that are mapped across exon-exon junctions.

e.g.
chr1 3995422 t 10 >><<<>>>>< CEFEJJJJJJ
So actually there are no reads mapped to this position but SAMTools still shows that this position has a coverage of 10X.

Now I assume that this will really be a problem in variants calling. Is there a way I can deal with this problem?

One way is of-course to generate the pileup file, remove all the "<>" signs and adjust depth accordingly and the use some tool like VarScan to call the variants. But I would like to use the bcftools itself.
slp is offline   Reply With Quote
Reply

Tags
rna-editing, rna-seq, tophat

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 07:45 PM.


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