SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Consensus from mpileup for haploid sequences (forcing base calls - no ambiguities) ericarcher Bioinformatics 6 01-24-2014 09:52 AM
strand specific wig files frymor Bioinformatics 0 07-27-2011 11:47 PM
Similarity programs for SNP files leofixings Bioinformatics 1 05-03-2011 09:03 AM
Gene prediction programs that work with multi-fasta files Hobbe Bioinformatics 0 11-16-2010 02:31 AM
Getting specific Exon sequences from Transcriptome andreask RNA Sequencing 0 06-15-2010 05:27 AM

Reply
 
Thread Tools
Old 01-24-2012, 07:41 AM   #1
aperdomos
Junior Member
 
Location: Sydney

Join Date: Jan 2012
Posts: 8
Default How to subtract specific DNA sequences from *mpileup* files (programs and scripts)

Hello everyone,

I have around 100 mpileup files from different species. These files have some specific DNA regions that are of my interest, since they include some genes that I want to compare. According to the pileup format generated with Samtools, these files have these columns:
chr pos ref coverage

Therefore, using these pileup files I want to specifics parts of the sequences, but keeping the information about position, and coverage. For example: DNA sequence from Samples A, B, C, etc from 33566 to 38000; 56000 to 59000; 110000 to 112500 and that presented a coverage of more than 10.

Consequently, I want to know if there is some specific programs that help to perform this, or if it is necessary to design a specific script. Does anyone have any script already designed that can be useful to subtract this info?

Many thanks,

Varo
aperdomos is offline   Reply With Quote
Old 01-24-2012, 09:34 AM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Starting from pileups? There isn't going to be much software written that's designed to take pileups as input. You have lots more options if you had .bams. SAMtools or BEDTools can both take .bam files and .bed files with your regions, and filter the .bam to only contain reads from those region. BEDTools also has some programs that will calculate coverage for you. So go back to the .bams if you still have them. Otherwise, I think a simple script can do what you want.
swbarnes2 is offline   Reply With Quote
Reply

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 12:51 AM.


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