More posts like this
-
Bacterial Seq Alignment
written by ravihansa82Hi Friends, I have illumina paired end seq reads for my bacterial species. Now I need to align it to a reference genome. What are …
-
Using paired-end date as single-end for mapping with bwa mem. Good, bad or ugly?
written by resug •Hi Biostars, I am trying to align my paired-end reads to my assembly with bwa mem (to be used for polishing with Racon afterwards). However …
-
How to calculate the mapping rate from BAM files generated from BWA
written by wangdp123Hi there, I used BWA MEM to map the paired-end reads onto a genome but it sees that BWA didn't report the mapping metrics such …
-
Extract mitochondrial DNA from paired-end RAD-seq data: bwa and stacks
written by CaffeSospeso •Hello, I'm using single digested and paired-end RAD-seq data. I would like to extract mitochondrial DNA data. I would expect to have some sequences derived …
-
Low mapping percentage after mapping RNA-seq reads to a closely related species
written by unawaz •Hi, I have Illumina 100 bp paired end RNA-Seq data from a non-model species. I mapped it to the closely related genome available and I …
-
How to reduce cross mapping alignments with paired-end reads ?
written by DavidHi, I´m trying to align paired-end reads (illumina) coming from a metagenomics experiment to a custom database of bacterial species. I´m using bwa-mem (althought could …
-
beginner question BWA program
written by rezahi i want to map my own reads (illumina hiseq) to reference using "bwa mem". i have 2 lane of paired end data (read1.fastq, read2.fastq, …
-
BWA option for FR sens only
written by PicasaHi, I am looking to `bwa mem` to map my reads. However, I would like to keep only the real Forward-reverse paired end and to …
-
Proper use of BWA-backtrack and BWA-MEM
written by cfarmeriHi, Biostars. I have a set of Illumina GA2 paired-end fastq files. These reads of the fastqs has 76b length. Before mapping, I trimmed barcode …
-
Better statergy for Gap Closing
written by second_exonHi Guys, Couple of months back, I sequenced (MiSeq) few BACs and assembled (Paired-end reads) using SOAPdenovo, but my assembly was fallen into many scaffolds. …
With the limited amount of information, this is difficult to answer (I admit, I don't know what I would need to give a good answer, anyway)
I would approach this in an iterative way, you can play with minimum seed length and bandwith, z-dropoff, then dive into match and mismatch scores and gap penalties. After all, much depends on how similar the genomes are. You will always align some xy percentage, but does that help you?
complementing Carambakaracho comment :
it's always advisable to start with the default settings and take it from there.
In principle, it's a bad idea to map to a different genome, but as suggested, start with the defaults and see how much can be used.