I am working on a project where I have bam file after alignment to reference genome. Now I want to filter out the mapped reads separately(forward and reverse read files) in fastq format. I have tried many commands but it shows the same reads in forward and reverse file. and I am unable to go ahead with the assembly step.
Please can anyone help me. Thanks