Are there any tools that can report all areas of very low coverage or areas of high coverage in NGS data?
-
-
I've written some custom code to programatically do this, but it only works on a few assembly formats (don't have time or funding to add more formats yet). But several assembler programs have ways to output the coverage levels of the contigs. It would then be up to you to parse that output and find any locations that are below or above your thresholds. -
A commonly used method is to map all your reads to contigs using bowtie2/bwa mem/bwa and calculate the coverage with samtools or bedtools.Comment
-
Can you give us some more information? What type of data do you have? Have you mapped reads to a reference yet?Originally posted by willMD View Post
If you do have a bam file of aligned reads to a reference sequence then you will be able to generate a plot of the per-base coverage of the reference sequence using bedtools genomeCoveragebed.Comment
-
Human whole genome data aligned to hg19 in bam file. Thanks. I'll look at bedtools.Comment
-
CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).
Despite this, “CRISPR helped turn genome editing from a specialized technique into...07-31-2026, 11:01 AM
Comment