Entering edit mode
2.5 years ago
gubrins
▴
290
Heys,
I'm working with samtools view to check a specific region from my bam file and I would like to convert that region into fasta piping it to angsd. However, I'm not being able to do it. How should I indicate I want the output as the input? I tried several things but they didn't work:
samtools view -b my.bam chrom | angsd -i - -doFasta 3 -out my.fasta
Any help would be welcomed!