Question: Generate consensus from BAM file
0
gravatar for blur
4.8 years ago by
blur110
European Union
blur110 wrote:

Hi,

I wish to get the consensus sequence of a certain section of my bowtie alignment.

I have the bam file, I thought perhaps I need to do faidx - is this possible in Galaxy?

Thanks,

 

ADD COMMENTlink modified 4.7 years ago by thackl2.7k • written 4.8 years ago by blur110
4
gravatar for Devon Ryan
4.8 years ago by
Devon Ryan92k
Freiburg, Germany
Devon Ryan92k wrote:

I have no idea if someone has made a galaxy pipeline to do this. The general steps are as follows:

  1. Sort and index the BAM file (this can be done with samtools)
  2. Use samtools mpileup like in this post: Questions Regarding Consensus Sequence Calling With Samtools / Bcftools / Vcfutils.Pl Note that I think you can skip the bcftools part step and just give samtools mpileup the -v option.
ADD COMMENTlink written 4.8 years ago by Devon Ryan92k
1
gravatar for thackl
4.7 years ago by
thackl2.7k
MIT
thackl2.7k wrote:

You can use the bam2cns module of proovread (https://github.com/BioInf-Wuerzburg/proovread). It's a consensus based PacBio-Illumina correction software. An advantage might be that the module let's you set a coverage cutoff and it only considers the best alignments (by score) for each location up this cutoff. This is paricularly useful, if you are dealing with repetitive regions, where suboptimal alignments introduce a lot of noise.

ADD COMMENTlink written 4.7 years ago by thackl2.7k
Please log in to add an answer.

Help
Access

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