Question: Pileup Haploid Base Calling, Creating Consensus Sequence
gravatar for Lee Katz
8.0 years ago by
Lee Katz3.0k
Atlanta, GA
Lee Katz3.0k wrote:

Hi all, I couldn't find if this has been asked here. I realize that samtools assumes that you are working with a diploid and will output heterozygous snps. How could I modify my script to make it not output heterozygous snps? In other words, how do I trick it into thinking that it is looking at a haploid genome? My snippet of code is in Perl and uses a command() subroutine which runs a string and prints the command to the screen.

command("samtools mpileup -uf $$settings{assembly} $bam > $out1");
command("bcftools view -cg - < $out1 > $out2");
command(" vcf2fq -d $minDepth -D $maxDepth < $out2 > $fastqOutNonstandard");
command(" varFilter -d $minDepth -D $maxDepth < $out2 > $variantsFile");
samtools pileup • 2.3k views
ADD COMMENTlink written 8.0 years ago by Lee Katz3.0k

Is there perhaps something like GATK which could do it instead?

ADD REPLYlink written 8.0 years ago by Lee Katz3.0k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 948 users visited in the last hour