mapping long-reads to a reference library
1
0
Entering edit mode
11 months ago
mthm ▴ 40

Hi,

I have long, pacbio, reads and I have a reference library of only repeats, I want to map the long reads on the repeats library using bwa mem, is this command correct?

bwa index mmm.pacbio.fastq.gz
bwa mem mmm.pacbio.fastq.gz  repeat-library.fasta | samtools sort -@16 -  -o mm.vs.lib.sorted.bam

Do I need reference genome for mapping?

bwa TE repeat pacbio • 448 views
ADD COMMENT
2
Entering edit mode
11 months ago
ATpoint 65k

You need to index the reference genome, not the fastq files. I think for long reads people prefer minimap2 though these days.

what is the purpose of indexing a genome

Please see the manuals of bwa-mem or minimap2, it explains what to do.

https://lh3.github.io/minimap2/minimap2.html

http://bio-bwa.sourceforge.net/bwa.shtml

ADD COMMENT

Login before adding your answer.

Traffic: 638 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6