How to take reads from WGS fastq files in order to map it with Exome?
Hi
I have a fastq files from WGS but I want to subsample the reads and just runs it as WES, is there a way to do that?
Thanks Sara
wgs
• 2,011 views
•
link
written
by
jonessara770
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
download WGS from dbGAP authorized with prefetch tools
written by jonessara770Hi I would like to download some WGS from dbgap with prefetch.2.9.1. But I get the below message: 2018-08-27T23:17:32 prefetch.2.9.1 warn: Maximum file size download …
-
interpretation of intronic variant detected by WES
written by jonessara770Hi I detected a variant in the intron of a gene that is predicted by MutationFinder to affect splicing by activation of an intronic cryptic …
-
How to take in consideration UMI in WES pipeline?
written by jonessara770Hi I have received WES data treated with SureSelectXT HS Reagent that has UMI. I wonder how I should do the analysis for this? what …
-
What are the step for germline analysis with CNVKit?
written by jonessara770Hello, I would like to call CNV on germline samples with CNVKit. the documentation for CNV calling from germline is a bit confusing. could you …
-
best and easy Germline CNV caller with WES
written by jonessara770Hi all, I would like to check the presence of germline CNV in my samples. I wonder what are the best germline CNV callers using …
-
duplicate read names in TCGA raw bam file
written by jonessara770Hi, I downloaded some TCGA ExomeSeq bam file and would like to convert them to fatsq. Picard can not convert them because there are some …
-
How to systematically check if a bam file is truncated
written by jonessara770Hi, I have downloaded some specific regions from public database bam files. I would like to check if downloaded bam files are ok or truncated …
-
how to extract human from mature.fa?
written by jonessara770Hello, I downloaded mature.fa files from mirbase database. it contains all of species but i want only human. is there any script/ Rpackage ect to …
-
run our pipeline on TCGA data
written by jonessara770Hello all, I am wondering if it is possible to run our own pipeline (WES, WGS or RNA-Seq) on TCGA raw data without downloading TB …
-
Finding Common Reads Across Multiple Fastq Files
written by Abhi<p>Hi All </p> <p>We have some metagenome samples(multiple illumina lanes). What I would like to do is find out % of reads that are common …
Why not map as usual and find out which reads map to the exome?
Actually, because my pipeline is not optimized for WGS and it takes a lot of time to run the WGS so I want just subsample and map it to WES?
You don't know which reads come from exonic regions until you map them. As far as I know your only option would be to first map them all and then, if needed, extract the "exome" regions.