Question: How to map the reads from a Bam file to another reference genome????
0
gravatar for Cindy
23 months ago by
Cindy20
Cindy20 wrote:

Dear all,

I have got BAM files of mitochondrial DNA sequencing, which were aligned to mtDNA reference genome. Now I want to look at whether those mapped reads also have some overlap with nuclear genome. What should I do? I have extracted mapped reads from bam files, should I just use bwa to align to nuclear reference genome?

Thanks

Cindy

sequencing • 1.8k views
ADD COMMENTlink modified 23 months ago by WouterDeCoster38k • written 23 months ago by Cindy20
1

Hi Cindy, If you have extract a fastq yes you can simply realign your reads against your new reference using bwa :)

Best

ADD REPLYlink written 23 months ago by Titus810
3
gravatar for WouterDeCoster
23 months ago by
Belgium
WouterDeCoster38k wrote:

You can convert the bam file back to fastq (samtools fastq), and use that for mapping just like usual mapping.

ADD COMMENTlink written 23 months ago by WouterDeCoster38k

Thanks. I have used samtools to remap.

ADD REPLYlink written 23 months ago by Cindy20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2029 users visited in the last hour