bedtools multicov problem
bedtools multicov problems:
I am counting number of mapped reads by using bedtools multicov:
bedtools multicov -s -D -bams AA.sorted.bam -bed yyy.bed > zzz_low_s.txt
Read counts for all bed file coordinates are zero.
When I use uppercase -S , there are reads in all coordinates.
bedtools multicov -S -D -bams AA.sorted.bam -bed yyy.bed > zzz_up_s.txt
When I do not use any stand option, there are reads in all coordinates, which are same number when I use uppercase -S
These coordinates have reads from both strands (I checked in IGV). Do you have any idea why lowercase s is not giving any count?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
what does the bed file in bedtools multicov command refer to ?
written by rania.hamdy1 •
-
coordinate of BAM and BED
written by yuabrahamliu •
-
bedtools multicov displays only the last BED feature
written by ThePresident
-
Bedtools multicover how to submit a bam list
written by benoit.tessoulin •
-
bedtools error: could not find indexes
written by Katie D •
-
Alignment with STAR produces unique bam files, but identical output from bedtools multicov
written by dec986
-
bedtools multicov not running
written by Bnf83
-
Identify genes in a tophat aligned bam/sam file
written by gchaves •
-
bedtools multicov and fraction overlap flag
written by fitziano •
-
Using Gnu Parallel For Bedtools
written by GouthamAtla
Can you post a subset of the BED file. I would guess that either the BED has "*" or "." in the strand column or no strand column at all. Therefore, neither "+", nor "-" would be specified and your result.