Hi everyone,
This is probably a stupid/easy question but I'm feeling a bit lost and any help would be much appreciated.
I extracted a particular genomic region of interest (centered around a gene) from a chromosome-level genomic assembly. I have a reference genome (not as complete) from a closely related species. I would just like to map the data that I have from this closely related species to that particular genomic region of interest.
Should I just use bwa-mem with the sequence of the region of interest as the reference or will this cause undue problems?
Thanks in advance!
1 answer
e or will this cause undue problems?
this will cause problems see: Exome Sequencing: Masking The Non-Genic Sequences ?
Log in to answer this question.
-
Using other individuals and related species to improve a de novo genome assembly
written by Ed 991 •
-
Mapping difference with BWA
written by rj.rezwan •
-
Mapping a small number of contigs to a reference subsequence
written by alexandru.bologa.marian •
-
Why do all of my PCR reads not map to the reference sequence?
written by Saran •
-
Skip orientation as there are not enough pairs by bwa mem
written by seta
-
BWA alignment and very low properly paired
written by seta
-
How to set up a reference file for BWA so that highly variable sites between conserved regions are …
written by DNAngel
-
Extract mitochondrial DNA from paired-end RAD-seq data: bwa and stacks
written by CaffeSospeso •
-
What to use as Close-species for miRDeep2
written by GBC_Zonatos •
-
HLA typing WGS data
written by emblake •
If you are using the selected region as a query against a different genome then it should not cause any issues. Do you have sequence data or just a sequence (what length) that you have extracted. If you have latter using
bwm memwould likely not be appropriate.So, you're saying that using the region of interest as a query in BLAST (for instance) should be okay? For the region of interest, I have a FASTA file with 10kb either side of the gene.
Yes it should be fine to use BLAST. If it will work well or not will depend on how closely related the second genome is and the quality of sequence available.
You could also try to use a long read aligner like
minimap2is the genomes are reasonably close.