I have data from an experiment which is looking into differential expression of small RNA. It seems that a fraction of the reads are mapped outside of the annotated regions of the mouse genome (as listed in the Ensemble GTF file). To be honest, I am not expecting to find new genes in such a popular model organism but I am still curious where these reads map. I would like to produce statistics for each chromosome listing the number of reads that are mapped outside of the annotated regions. What's the easiest way to produce such statistics? I am using bowtie1, samtools and htseq-count in my current pipeline but so far I have always limited my analysis to the annotated regions so I am not sure how to best approach this - any advice is most welcome!
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
If you have read alignments in bam format on the one hand and a gene annotation in gtf/gff/bed format on the other hand it should be simple to compare both with bedtools.
For instance something like "bedtools intersect -v -abam mapped_reads.bam -b exons.bed" should give you the reads that do not intersect any annotated exon (you can use entire gene regions if you prefer to discard intronic reads). Pipe this to "samtools view -c -" to get the number of reads.
More precisely:
Code:bam1=mapped_reads.bam bam2=intergenic_reads.bam annotation=annotated_genes.bed bedtools intersect -v -abam $bam1 -b $annotation > $bam2 samtools index $bam2
Code:for chr in `cat list_of_chromosomes.txt` ; do N=`samtools view -c $bam1 $chr` ; n=` samtools view -c $bam2 $chr` ; echo $chr $n $N ; done
Code:samtools idxstats $bam2
Latest Articles
Collapse
-
by seqadmin
Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.
Somatic Genomics
“We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...-
Channel: Articles
05-24-2024, 01:16 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 06-03-2024, 06:55 AM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
06-03-2024, 06:55 AM
|
||
Started by seqadmin, 05-30-2024, 03:16 PM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
05-30-2024, 03:16 PM
|
||
Comprehensive Sequencing of Great Ape Sex Chromosomes Yields Insights into Evolution and Genetic Variability
by seqadmin
Started by seqadmin, 05-29-2024, 01:32 PM
|
0 responses
29 views
0 likes
|
Last Post
by seqadmin
05-29-2024, 01:32 PM
|
||
Started by seqadmin, 05-24-2024, 07:15 AM
|
0 responses
216 views
0 likes
|
Last Post
by seqadmin
05-24-2024, 07:15 AM
|
Comment