weird behaviour on bedtools
Hi All,
I want to extract the counts on defined regions (my.bed). But i'm not able to understand the results. The command says extract the counts form the bam files using bed file also the region of the bed region has to have >=80%.
multiBamCov -f 0.80 -bams 1.bam 2.bam -bed my.bed > output.bed
sample result.
14 21853003 21853010 1 2 ENST00000390435 5 5
According to the multiBamCov i have 5 reads the regions i defined but the IGV shows different count.
Why is that?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
problem with gene abundance calculation
written by v.berriosfarias •
-
quantify counts from sorted bam files using a bed file containing regions of interest
written by bioyas •
-
Tutorial: Calculation of TSS Scores for ATAC-seq Data
written by Sael •
-
Getting different results from bedtools and samtools
written by wstla27 •
-
How to compute target aligned reads/total aligned reads (read enrichment) with Python?
written by vw •
-
cwlexec and cwltool not interpreting baseCommand the same way
written by drkennetz
-
bedtools multicov displays only the last BED feature
written by ThePresident
-
ngs.plot.r heatmap for single-ntd bed entries
written by rbpdee •
-
Convert paired end bam file into bed of fragments
written by lkmklsmn
-
Extract Reads From A Bam File That Fall Within A Given Region
written by abi •
Maybe IGV not displaying all reads? (downsample display).