bowtie mate 2 error

Hello everyone,

I'm working on ChIP-seq files and I'm getting this type of error

mem: invalid option -- '2'

here is the code :

bwa mem hg38.fa -1 R1_001.fastq.gz -2 R2_001.fastq.gz | samtools sort > file.bam

I'm working on paired end ChIP-seq experiment so I have to use both files, I'm new to it so I wonder also if I have to unzip first ,

thank you for you support

chip-seq

• 1,163 views • link

updated 6.1 years ago by h.mon

• written 6.1 years ago by Morris_Chair

The current version of bwa mem does not have any -1 or -2 option.

bwa mem [options] <idxbase> <in1.fq> [in2.fq]

• 0 views • link

written 6.1 years ago by Arup Ghosh

1 answer

You are mixing up bowtie2 and bwa mem: you are calling bwa mem, but -1 and -2 are bowtie / bowtie2 arguments.

• 0 views • link

written 6.1 years ago by h.mon

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 keeps failing to locate the index files
written 1.4 years ago by pairedttest •

Okay, so I am at my wits end here. I just want to run bwa-mem on some fastq files but bwa-mem does not seem to …
It is weird there are no reads less than 100 bp in qualimap.
written 4.1 years ago by Dan •

For ATAC-seq, I trimmed the reads using ``` trim_galore --illumina --paired --fastqc \ fastq/${FILE}/*R1_001.fastq.gz fastq/${FILE}/*R2_001.fastq.gz \ -o trim/${FILE} ``` When I run Bowtie2 using ``` …
skip orientation FF/RF/RR by bwa mem
written 5.9 years ago by 14301010010 •

Hi everyone I extracted DNA fragments from human somatic tissue, sequenced them, and then try to align data to the ref genome (hg38) using bwa …
Skip orientation as there are not enough pairs by bwa mem
written 7.2 years ago by seta

Hi all, I extracted the specific region from a BAM file generated by aligning fastq (Illumina/whole genome sequencing) on hg38 and converted to corresponding fastq …
remove duplicates with markdup
written 7.5 years ago by dimitrischat

Hello all. I am trying to do a paired-end analysis. Ran the 2 fastq files with hisat2: hisat2 -p 10 -x '..index_hg19/indexed' -1 R1_001.fastq.gz -2 …
Is there something wrong with my Indel Calling protocol?
written 8.6 years ago by lishen0709 •

hello everyone, I'm using BWA+SAMtools+Varscan to call indels to check my sgRNA cutting efficiency. The protocol I used is shown below: **1> set index:** bwa …
BWA mem error
written 9.0 years ago by gingergeiger22 •

Hi, I am trying to run BWA using the following command: module load bwa/0.7.15 time bwa mem -aCMP -t 12 -R '@RG\tID:MiSeq\tSM:VA-002\tPL:illumina\tLB:XT40 PU:Batch22' referencebasename R1_001.fastq.gz …
Loop tophat to align multiple RNA-seq files
written 10.2 years ago by mmrcksn

Hello, I'm a beginner! trying to use a loop to align multiple paired-end RNA-seq samples with tophat. Here's what I have that isn't working (basically …
BWA batch script does not execute
written 10.3 years ago by madkitty

I'm running an alignment with BWA on multiple exomes files on our server. For some reason it does not output anything, and doesn't even seem …