Pull FASTQ reads and map them with a gene
Hello,
I want to pull fastq reads those who are mapping with a gene. After that, also need to find out the coverage and depth of the reads those who are mapped/aligned with the gene.
How to carry on the work ?
Please share your research ideas/bioinformatics technqiues for the above problem.
thanks, in-advance.
• 183 views
•
link
written
by
bioinf_data-analyst9449 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Haplotype Phased Assembly Contigs to Chromosome Annotations
written by turcoa1 •
-
cover depth estimation with samtools depth and mpileup
written by Ирина •
-
How to extend assembled PACBIO contigs for scaffolding
written by K •
-
Reads mapped on wrong strand
written by qstefano •
-
How bedtools genomecov handles reads without both ends mapped in paired sequencing
written by rimelof •
-
very confused with GC bias
written by 9521ljh •
-
Should I include ambiguously mapped reads in the realignment around InDels?
written by stu111538
-
How to find nu. of variants present in a gene for my Exome seq. sample
written by bioinf_data-analyst9449 •
-
What is the definition of "read depth" vs "coverage"? (again...)
written by ariel.balter
-
All the reads aligned using STAR have low mapping quality (0-3). What is happening?
written by kirannbishwa01
samtools view input.bam "Chr10:100-200" > output.bamThis will extract the reads from bam file that are in the region 100:200 in chr10 as example (make sure to know how your chromosomes are named). But this is the reads only (not fastq) from the reads name you can write script or maybe there is already some tools to do that extract there corresponding fastq reads
a lot of examples for this issue already exist here
The answer for the coverage section is here A: Depth of coverage for specific target from bam file