mapping long-reads to a reference library
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?
• 27 views
Read more here: Source link