Question: Fastq screen: use bwa aln and not bwa mem
0
gravatar for Matteo
8 months ago by
Matteo0
Matteo0 wrote:

Hi All,

I'm using fastq screen to check some files. It runs fine, but I noticed that when using bwa as aligner it defaults to bwa mem, as mentioned in the docs for the --aligner parameter

... BWA with mem -a

As I have 42bp SE reads I'd rather use bwa aln, which isn't mentioned anywhere in the docs, so I was wondering if there's a way to hack this default and instruct fastq screen to use aln and not mem. I tried to play with the parameter

bwa "<"text">" : Specify extra parameters to be passed to BWA

but I don't seem to get an acceptable syntax. Does anyone know whether this
1) is possible at all
2) if yes, how to do it?

alignment • 307 views
ADD COMMENTlink modified 8 months ago by Pierre Lindenbaum129k • written 8 months ago by Matteo0
1
gravatar for Pierre Lindenbaum
8 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum129k wrote:

fastq screen is a perl script. How about editing the script and replacing

else {    #BWA
        $aligner_command = "$path_to_bwa mem -a $bwa_opts $illumina_flag -t $number_of_threads $library->[1] \'$file\' 2>$error_filename |";
    }

with a custom script that would wrap 'bwa aln' ?

ADD COMMENTlink written 8 months ago by Pierre Lindenbaum129k

That's not a bad idea, i'll give it a shot, thanks! Any clue if the .sai file produced by aln would be enough or if I need to wrap samse/sampe as well?

ADD REPLYlink written 8 months ago by Matteo0

Are you not seeing expected results? Since fastq screen is a quality check (and not serious analysis) you may be ok with the bwa default.

ADD REPLYlink written 8 months ago by genomax85k

I'm using the --tag --filter parameters to select reads based on alignment on multiple genomes, which I will next map for "serious analysis". It has to map all the files (without subsetting) to multiple genomes and mem is slower than aln for short reads, it's taking ages.

It is rather to optimize the poipeline for the next time I run it rather than for a serious issue with the result

ADD REPLYlink written 8 months ago by Matteo0

Are you looking to bin the reads subsequently? If you know the genomes you are interested in I suggest you use bbsplit.sh from BBMap suite to do this job once. No need for multiple passes.

ADD REPLYlink modified 8 months ago • written 8 months ago by genomax85k

you'll need same/sampe

ADD REPLYlink written 8 months ago by Pierre Lindenbaum129k
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: 729 users visited in the last hour