BWA MEM -K option

In some alignment scripts, I find that they use the command bwa mem -K 10000000. However, I cannot find any description of the option -K (captial K, not -k). Does anyone know what -K really means?

Thanks!

alignment bwa

• 1,360 views • link

updated 3.1 years ago by DBScan •

• written 3.1 years ago by Yihang •

2 answers

"fixed chunck size" https://github.com/lh3/bwa/blob/139f68fc4c3747813783a488aef2adc86626b01b/fastmap.c#L193C22-L193C38

used in https://github.com/lh3/bwa/blob/master/bwa.c#L79

used in https://github.com/lh3/bwa/blob/master/bwa.c#L104

related to the number of bytes read (fastq) for each iteration

• 0 views • link

written 3.1 years ago by Pierre Lindenbaum

If I am not mistaken, this flag is usually used for reproducibility reasons.

• 0 views • link

written 3.1 years ago by DBScan •

More posts like this

BWA-MEM with an array of files
written 10 months ago by garcesj •

Hi, I'm trying to run BWA-MEM with an array of files rather than manually inputting each single one... but it seems BWA is not getting …
bwa mem skip orientation when aligning a small subset of reads
written 3.1 years ago by Yihang •

For pair-end read files, `read1.fastq`, `read2.fastq`, I use the following command `cat read1.fastq | head -n 4000 > read1_first1000.fastq` (and `cat read2.fastq | head -n …
BWA MEM -P option
written 4.2 years ago by chisqga •

Hello everyone I am trying to align my sequencing reads (short reads, paired-end) to reference using BWA MEM. I am wondering what is `-P` option …
Inconsistency between alignment rates of BWA-mem and bowtie2
written 5.3 years ago by jesvinjoy35 •

I have been trying to align a few bacterial species with the human genome to check the alignment rate (using BWA-mem and Bowtie2). However, both …
how to understand the -k parameter in bwa mem method?
written 6.2 years ago by dongguiying2017 •

I use the bwa-mem to align the reads to reference genome. The official website says the parameter -k INT represents Minimum seed length. Matches shorter …
bwa mem -K 100000000
written 7.9 years ago by gwotto •

Hi, in Regier et al.: "Functional equivalence of genome sequencing analysis pipelines enables harmonized variant calling across human genetics projects", bioRxiv, doi 10.1101/269316 The use …
SAM to BAM conversion error
written 9.3 years ago by jitendoshi2006 •

I am trying to convert SAM file to BAM file after genome alignment using bwa-mem on a cluster. I perform genome alignment on a cluster …
bwa mem algorithm filter secondary alignment
written 9.9 years ago by xudongyang •

Hello,everyone I want to use the bwa mem algorithm to analysis my PE 100 seqdata (BWA Version: 0.7.12-r1039). I used "-M" option to mark the …
bwa aln for paired-end sequencing data in one .fastq file
written 10.2 years ago by ciemanek

Hi! I have data from paired-end Illumina RNA-Seq in one .fastq file. I want to make comparison of alignment quality for bwa aln and mem …
BWA mem how to know if a read is mapped uniquely?
written 10.7 years ago by steven.lin.pur •

<p>I have read many posts regarding this question. The XT, XA tag doesn't seem to be used anymore by BWA mem, so I can't use …