I have generated multiple read mapping files (.sam) generated with bowtie2 "--no-unal" parameter. As the number of files are huge (7 reference genomes mapped against >6000 reads) i want a loop to count the % (or total number) mapped reads for each file.
HTSeq could be an option but don't know how to loop it for multiple files.
Any help would be much appreciated.
HTSeq could be an option but don't know how to loop it for multiple files.
Any help would be much appreciated.
Comment