I have download the hiC data from NCBI and nd instead of mapping them on reference genome, I mapped the data on list of genes separately. it means I mapped the R1 and R2 separately.
Now I want quantify how many reads are their hows R1 part are mapped on gene1 (suppose) and R2 part on gene 2.
for example
if total number of reads are 100 (50 R1 and 50 R2). total number of genes are 5 (suppose). out of 100 reads, 30 reads (R1/R2) are mapped on gene1. then find out counter part of reads (R2/R1) which mapped apart from gene1..
positive suggestions are always welcome
2 answers
Try running BLAT. You should be able to compare the table output.
I suggest that with current state of Hi-C technology, it is better to use the normalized data, not the raw one. For the resolution considerations, it is also better to bin Hi-C data by 1Mb/500kb bins and for a pair of genes compare the Hi-C value for pairs of bins in which those genes are positioned.
You can download normalized Hi-C contact maps here:
Log in to answer this question.
-
Determining total nucleotides for paired end metagenomic sequences
written by cssulliv •
-
Tophat problem - R1 and R2 reads mapped to different strandeds of viral genome
written by WZOUXM •
-
Realign Aligned.sortedByCoord.out.bam to Aligned.toTranscriptome.out.bam for RSEM input
written by mkvalluru •
-
Weird fragments when aligning PE reads to duplicated genes.
written by bsksln •
-
How to properly combine two bam files of a paired-end data
written by izzy.yichao.cai
-
How to obtain insert sequence from paired-end mapping result (SAM/BAM)
written by gundalav
-
filtering paired end mapped reads form SAM/BAM file
written by Raghav
-
how to extract those reads which mapped at 1kb of specific location on chromosomes
written by Raghav
-
samtools + map
written by elfeusilva •
-
Filtering Bed Files By Using Bedops
written by Raghav
