# To remove reads with poor mapping quality: samtools view -b -q30 SRRXXXXXXX_sort.bam > SRRXXXXXXX_sort_q30.bam