Hi,
I wonder if there is any option in, e.g., samtools to display regions (positions) with the base coverage >= given threshold. Something opposite to the mpileup option (when we have positions and want to see base coverage and quality). I have very small number of reads mapped to the genome and I want to see those tens or hundreds regions covered by, at least, 10-15 reads.
Thanks.
I wonder if there is any option in, e.g., samtools to display regions (positions) with the base coverage >= given threshold. Something opposite to the mpileup option (when we have positions and want to see base coverage and quality). I have very small number of reads mapped to the genome and I want to see those tens or hundreds regions covered by, at least, 10-15 reads.
Thanks.
Comment