I have 11 mpileups generated by novoalign from whole genome sequencing of Daphnia, which has a genome size of ~200MB. Each is the result of alignment of about 130M paired-end illumina reads. I find that all of the mpileups have about 110M lines, implying that we have info for about 55% of the reference. I want to quantify the degree of overlap in coverage among the mpileups. Is there an easy way to do this with samtools or something similar? Also, is there a way to construct an mpileup file that has a line for every site in the reference, rather than just those with read coverage?
Announcement
Collapse
No announcement yet.
X
Latest Articles
Collapse
-
by seqadmin
The recent pandemic caused worldwide health, economic, and social disruptions with its reverberations still felt today. A key takeaway from this event is the need for accurate and accessible tools for detecting and tracking infectious diseases. Timely identification is essential for early intervention, managing outbreaks, and preventing their spread. This article reviews several valuable tools employed in the detection and surveillance of infectious diseases.
...-
Channel: Articles
11-27-2023, 01:15 PM -
-
by seqadmin
Microbiome research has led to the discovery of important connections to human and environmental health. Sequencing has become a core investigational tool in microbiome research, a subject that we covered during a recent webinar. Our expert speakers shared a number of advancements including improved experimental workflows, research involving transmission dynamics, and invaluable analysis resources. This article recaps their informative presentations, offering insights...-
Channel: Articles
11-09-2023, 07:02 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 08:26 AM
|
0 responses
11 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:26 AM
|
||
Started by seqadmin, Yesterday, 08:12 AM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:12 AM
|
||
Started by seqadmin, 11-27-2023, 08:12 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
11-27-2023, 08:12 AM
|
||
Started by seqadmin, 11-22-2023, 09:29 AM
|
1 response
60 views
0 likes
|
Last Post
![]()
by VilliamPast
11-25-2023, 02:47 AM
|
Comment