How to extract mapping features of multiple samples in one output?
1
0
Entering edit mode
3.6 years ago
star ▴ 290

I have a multiple samples and I run mapping with Bowtie2 for them and now I have .bam files. I have tried to run Multiqc to see how is mapping rate, duplication rate and other features but It did not work with bowtie output I think. I would like to know is there any graphical program like Multiqc for observing aligning features?

bowtie2 mapping multiqc ChIP-Seq RNA-Seq • 3.1k views
0
Entering edit mode

but It did not work with bowtie output I think

what happened ?

0
Entering edit mode

If I remember well, in bowtie2 informations about the alignment (alignment rate, not aligned rate...) are in the standard output (terminal).

If you want to check reads quality, GC content... You will have to run a tool like fastqc on your data and use multiqc to compile them

4
Entering edit mode
3.6 years ago
h.mon 34k

MultiQC doesn't extract information from the bam files, rather, it extracts from log files. Bowtie2 logs are output to stderr, you have to redirect them to a file and pass this file to MultiQC.

2
Entering edit mode

Correct! If you forget, you can always run samtools stats on your aligned BAM file, which will give you a file to feed to MultiQC for a % alignment plot.