mapping long-reads to a reference library

mapping long-reads to a reference library

1

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

• 27 views

updated 1 hour ago by

54k

written 2 hours ago by

▴

30

Read more here: Source link