I have mapped 6 paired illumina raw reads against a reference genome and filtered by mapping quality and chromosome number 6. Then sliced using predefined coordinates of a gene region and merged the resulting bam files. Now I have a bam file of newly mapped genome only within the coordinates. So I want to get this sequence in fasta format. How do I do that in galaxy or using any other software?