I am trying to create a smaller BAM file from the original BAM file. The smaller BAM file should have info on only about 10 loci.
Can I not just append the
mpileup() output of the region of interest to the header (
samtools view -H <bam>)?
Any pointers on how I can accomplish this would be appreciated!