Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
how to eliminate reads that map to multiple places in genome ammoon Bioinformatics 7 08-23-2017 05:43 AM
Differential Expression: is it better to map reads to genome or transcriptome? dalesan Bioinformatics 26 04-22-2014 03:02 AM
Tophat RNASeq mapping - right reads map and left reads map 50% less airad22 Bioinformatics 2 08-14-2013 07:26 AM
How to map short reads to a distant genome? ynwh Illumina/Solexa 5 08-03-2011 05:56 AM
How to map 454 reads/contigs to a mitochondrial genome? fruktimport Bioinformatics 2 03-28-2011 10:35 AM

Thread Tools
Old 01-09-2019, 09:54 AM   #1
Location: Oklahoma

Join Date: Jan 2012
Posts: 35
Default Deciding how many reads map to a plasmid or to a genome

Hi all,
I have illumina genome reads for an E. coli a collaborator is studying. The prep used had a plasmid as well. I have two questions:

1) For the reads that match both (e.g. lacI gene), how can I tell which came from the plasmid and which came from the genome?
2) Should I normalize by dividing on the total reads per library or only the mapped reads in a library?

Thanks in advance...

Last edited by fznajar; 01-09-2019 at 12:44 PM.
fznajar is offline   Reply With Quote
Old 01-09-2019, 02:09 PM   #2
Senior Member
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,101

#1 = Unless there is/are SNP's which are captured by that particular read you would not be able to tell where the read came from. I assume the sequence of the gene is otherwise identical?
GenoMax is offline   Reply With Quote
Old 01-09-2019, 02:10 PM   #3
Location: Oklahoma

Join Date: Jan 2012
Posts: 35

You are correct GenoMax. They are identical. Appreciate it.
fznajar is offline   Reply With Quote
Old 01-29-2019, 11:39 AM   #4
Senior Member
Location: Québec, Canada

Join Date: Jul 2008
Posts: 260

If you'are using "bwa sampe" for aligning, the documentation says:

Generate alignments in the SAM format given paired-end reads. Repetitive read pairs will be placed randomly.
I don't know if it is the case with "bwa mem" or with "bwa aln".
seb567 is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 12:38 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO