Weighted metagenome mode about grodon HOT 1 CLOSED

Hi!

The read depth can be calculated using any read-mapping software to map the reads back to the genes output by prokka

For example, with the BWA aligner and using samtools to convert to bam output:

bwa index genes.ffn                                                                        
                          
bwa mem genes.ffn reads.fwd.fastq.gs reads.rev.fastq.gz | samtools view -bS - > mapped_reads.bam 

And then you can get a nice summary table listing the average coverage of each gene using bamcov:

bamcov -H mapped_reads.bam > genes_coverage.tbl

from grodon.