Questions tagged [coverage]

When referring to the coverage of a base in a reference genome, the number of reads covering that base in an alignment. When referring to the coverage of a genome in a sequencing experiment, the average number of times that a base was sequenced. = total number of sequenced bases / genome size. For example, a 100x coverage means that each base was sequenced on average 100 times. For the human genome, that would be ~ 300 Gigabases.

6 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
2
votes
1answer
109 views

Coverage required

I was came across a problem during an exercise in a book and I don't really know how to solve it. I feel like something's missing. "coverage, c = $NL/G$ (N=number of reads, L=read length, G=genome ...
1
vote
0answers
38 views

Analyzing a blast result of NGS data

I am analyzing BLAST results from total RNA NGS, pair ends. This BLAST results was performed against the ssRNA virus database download from NCBI. I filtered those read data whose paired ends ...
1
vote
0answers
18 views

How to choose reference to calculate coverage of WGS of organism that has high intra-species variation?

Calculation of coverage involves mapping reads with a reference and depending on the choice of reference we can get very different result. If a species has a genome that has very low size variation ...
1
vote
0answers
69 views

bedtools coverage - Report the depth at each position in each A feature

I am trying to use bedtools coverage to compute the sequencing depth at every positions of a chromosome but it didn't work as I expected. It reported 0 coverage at every positions. This is how I did ...
1
vote
0answers
100 views

CollectHsMetrics base coverage output has overlapping targets

I am trying to get base coverage information, and am using the --PER_BASE_COVERAGE output from Picard tools to get a text file of base level coverage. However, some ...
1
vote
0answers
167 views

bedtools single-nucleotide coverage in BED-specified regions for multiple BAMs

I have a BED6 (BED + name, score, strand info) file that defines some regions of interest. I also have a set of BAMs corresponding to different samples. I would like to obtain output similar to <...