Hi all,
I'm studying some documents about countOverlaps() (I never work with this tool), and come across with the below sentence: "The first run of the countOverlaps()marks the individual reads with running numbers. The reads with the flag one map to a single gene only and these are mapped to the genes they map to during the second run of the function countOverlaps()." I got confused with the bold sentence. Could you please kindly to explain a bit what's its concept and meaning?
Many thanks in advance
2 answers
The first run of the countOverlaps()marks the individual reads with running numbers.
First run marks how many genes the reads overlap ...
The reads with the flag one map to a single gene only
... then takes the reads overlapping only one gene (excluding multi-mappers or overlapping genes) ...
these are mapped to the genes they map to during the second run of the function countOverlaps().
... then the second run associates these reads with the corresponding genes.
If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted.
Actually, these sentences refer specifically to the following two commands:
hits <- countOverlaps(h1b, txdb)
ol <- countOverlaps(txdb, h1b[hits==1])
The first command counts overlaps with txdb elements (e.g. genes) per each read in h1b.
In the second command, only h1b items characterized by one hit: h1b[hits==1] are applied to countOverlaps. The output of the second command is a named integer vector reporting, per each gene, the number of reads which overlap that gene.
Log in to answer this question.
-
why hisat2 mapping sam file has many sam flags
written by szp770 •
-
Please help to understand the concept of the sentence about cufflinks
written by seta
-
inner and outer genomic coordinates concept
written by seta
-
Please help me with the fate of multimapping reads in the RNA-seq context
written by seta
-
PE read flag 69 and 137 meaning
written by epigene
-
Question about non-polymorphic site and flanking sequence.
written by mangfu100
-
I didn't understand a term about non-template DNA sequence
written by mangfu100
-
Question about Somatic and Germline with homozygous and heterozygous
written by mangfu100
-
Best way to map paired-end, uniquely mapped reads with Bowtie2?
written by Ian
-
The Reference Genome Of The Reading In Bam File
written by shl198
Can you link to that documentation, so we know which package you're referring to?