samtools depth -a in1.bam > depth_in1_both.tsv

samtools view -b -F 0x10 in1.bam | samtools depth -a > depth_in1_fwd.tsv
samtools view -b -f 0x10 in1.bam | samtools depth -a > depth_in1_rev.tsv

bedtools intersect -a in1.bam -b regions.bed |
  samtools view -b -F 0x10 in1.bam |
  samtools depth -a > depth_in1_regions_fwd.tsv

paste <(samtools depth -aa in1.bam)
      <(samtools view -b -F 0x10 in1.bam | samtools depth -aa)
      <(samtools view -b -f 0x10 in1.bam | samtools depth -aa) |
  cut -f 1,2,3,6,9 > depth_in1_bfr.tsv

paste depth_in1_both.tsv depth_in1_fwd.tsv depth_in1_rev.tsv |
  cut -f 1,2,3,6,9 > depth_in1_bfr.tsv