Hello, I am trying to find a script that can help me check the depth of sequencing of the region. I want to see if the particular region was sequenced in all areas or is it different. Whole-exome sequencing FastQ file.
Thank you
•
written
by
nudrylab •
1 answer
Fastq files first have to be aligned to the genome, then you may apply mosdepth or something like samtools depth. Alignment software could be bwa mem or bowtie2 among other choices. Please google for popular alignment tools or read previous posts on that topic here. Since fastq files only contain the raw sequences without positional information you cannot quickly check any region from those files, and before you ask, no you cannot align against only the region you are interested in, it must be the entire genome that the WGS was based on to avoid false alignments.
Log in to answer this question.
-
Getting total sequenced megabases for each sample
written by zizigolu
-
VAAST (Variant Annotation, Analysis and Search Tool)
written by nudrylab •
-
paired-end fastq file
written by evelyn
-
Can I compute Tumour Mutational Burden from rnaseq fastq files?
written by Pin.Bioinf
-
extracting reads from bam file
written by KVC_bioinfo
-
Can I merge the exome interval lists from two different kits for identifying CNVs using XHMM?
written by bioinforesearchquestions
-
What are the typical IBD PIHAT values calculated by PLINK for duplicate samples?
written by anp375
-
MIN and MAX read depth for vcfutils.pl
written by nickmarinakis92 •
-
Whole exome length of exons from bam file
written by bioinformatics.cancer
-
Forum: TCGA: How to tell if the data is exome or WGS?
written by Les Ander
There is a program called
mosdepth(LINK) that will fit the bill. Install usingconda.