Question: Freebayes 'bamleftalign' function usage issues (Part of the 'Du Novo' Sequencing pipeline)
gravatar for jmath97
4 weeks ago by
jmath970 wrote:


I've been trying to use the 'Du Novo' sequencing pipeline to analyze a set of bacterial whole genome sequencing data. This is the paper outlining the method: Streamlined analysis of duplex sequencing data with Du Novo. At the 9th step, the pipeline suggests using a function from the package 'freebayes' called 'bamleftalign'. I've tried using this function as I've listed below.

bamleftalign duplex.filt.aln.sort.bam -f ref.fasta > output.bam

For some reason, it doesn't provide an output -- it leaves the output.bam file empty. Is there something wrong with my usage of the function?

Is there something else that I should be checking? I checked the function in my PATH and I've ran and passed the freebayes package test files.

Any and all input would be appreciated!

freebayes alignment du novo • 190 views
ADD COMMENTlink modified 4 weeks ago • written 4 weeks ago by jmath970
gravatar for jmath97
4 weeks ago by
jmath970 wrote:

For any other users that might have this question, I realized my usage was incorrect. The function requires piping the bam input file into the function as opposed to taking the bam file as a standard input.

cat duplex.filt.aln.sort.bam | bamleftalign -f ref.fasta > duplex.filt.aln.sort.leftalign.bam
ADD COMMENTlink written 4 weeks ago by jmath970
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: 1761 users visited in the last hour