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 08:52 AM
Generate consensus sequence with Samtools mpileup? Heisman Bioinformatics 4 10-28-2013 04:57 PM
mpileup consensus behaviour rwness Bioinformatics 5 07-24-2012 07:12 AM
Differences in Consensus Sequences bhan Genomic Resequencing 0 08-01-2011 07:36 AM
consensus from mpileup prm36 Bioinformatics 3 07-14-2011 08:11 AM

Thread Tools
Old 05-05-2012, 12:15 PM   #1
Location: london

Join Date: Mar 2012
Posts: 12
Default My unreliable consensus sequences with mpileup?

Dear All,

I've used (what I think are) default settings to align paired end Illumina data using bwa then samtools and finishing with mpileup.

Before obtaining my sequences I PCR'd all our strains in a specific area of the genome to confirm their family according to the presence or absence of certain SNP's.

The problem is that the consensus fasta that is derived from mpileup always has the family of the reference strain and not the query strains that I've lined up (despite having good coverage in the region that defines the family).

Somehow the reference sequence is being output rather than the query sequence at this part of the genome.

Can anybody help?

Best wishes and thanks
lg36 is offline   Reply With Quote
Old 05-06-2012, 07:24 AM   #2
Senior Member
Location: St. Louis

Join Date: Dec 2010
Posts: 535

Why don't you just call SNPs and compare the SNP calls?
Heisman is offline   Reply With Quote
Old 05-06-2012, 04:53 PM   #3
Location: london

Join Date: Mar 2012
Posts: 12

Thanks for your reply.

I was wanting to somehow calibrate the SNP filters based on the correctness of the consensus sequence. I figured (perhaps incorrectly) that a reasonable test of the consensus sequence (and hence the SNP calls) would be to check the consensus sequence against confirmed oligonucleotides in the lab. Does this make any sense.

Can it be that the SNP's called are totally different from the consensus sequence.
lg36 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 05:38 PM.

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