one command line for getting consensus sequences from bam file

samtools mpileup -uf ref.fa aln.bam | bcftools view -cg – | vcfutils.pl vcf2fq > cns.fq

Leave a Reply

  

  

  

You can use these HTML tags

<a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>