coverage.pyΒΆ

usage: coverage.py

Construct contig coverage table given an input `assembly` and provided files.
Provided files may include one from the list below:
1. `fwd_reads` and/or `rev_reads` and/or `se_reads`
2. `sam` - alignment of `assembly` and `reads` in SAM format
3. `bam` - alignment of `assembly` and `reads` in BAM format
4. `bed` - alignment of `assembly` and `reads` in BED format

optional arguments:
  -h, --help            show this help message and exit
  -f ASSEMBLY, --assembly ASSEMBLY
                        </path/to/metagenome.fasta>
  -1 [FWD_READS [FWD_READS ...]], --fwd-reads [FWD_READS [FWD_READS ...]]
                        </path/to/forwards-reads.fastq>
  -2 [REV_READS [REV_READS ...]], --rev-reads [REV_READS [REV_READS ...]]
                        </path/to/reverse-reads.fastq>
  -U [SE_READS [SE_READS ...]], --se-reads [SE_READS [SE_READS ...]]
                        </path/to/single-end-reads.fastq>
  --sam SAM             </path/to/alignments.sam>
  --bam BAM             </path/to/alignments.bam>
  --lengths LENGTHS     Path to tab-delimited lengths table with columns of
                        contig & length.
  --bed BED             </path/to/alignments.bed>
  --cpus CPUS           Num processors to use. (default: 2)
  --from-spades         Extract k-mer coverages from contig IDs. (Input
                        assembly is output from SPAdes)
  --out OUT             Path to write a table of coverages