Hi All,
I have a bam file and I calculated read depth using samtools depth and I now want to filter the bam file to have only the contigs that have a depth between a certain value. I was able to filter the samtools output to create a txt file that only contains the depth region im looking for, but Im struggling to then filter the bam file with this information. The depth file has three columns, the contig name, the POS, and the depth
I attempted to remove the depth column and use grep to search the bam file, but my output was empty. I also tried to find answers with awk, but quickly got confused as I am still a fairly novice coder. Any ideas on how to proceed would be greatly appreciated.
•
written
by
Dak •
1 answer
awk -F '\t' '(int($3)> 123 ) {printf("%s\t%d\t%s\n",$1,int($2)-1,$2);}' depth.output | sort -t $'\t' -k1,1 -k2,2n | bedtools merge > select.bed
samtools view -M -L select.bed in.bam
Log in to answer this question.
-
Bam File Average Coverage Depth Changes After Fixmate
written by jdru •
-
samtools depth -r
written by lz1903 •
-
Samtools depth to get the read depth
written by bandanaschapagain •
-
samtools index error
written by evelyn
-
Errors with samtools depth on merged bam file
written by rbronste
-
no output after running samtools depth
written by Ana
-
compute depth of coverage for certain regions of the genome using samtools depth
written by Ana
-
How to filter a multisample VCF based on individual sample depth?
written by jtwalker •
-
bam-readcount tool: output zero coverage bases
written by user230613
-
Samtools Depth Can Only Have Stdout?
written by Chris •
