Skip to content

Script that takes samtools pileup output and write consensus fasta file as output.

License

Notifications You must be signed in to change notification settings

MeHelmy/pileup_consensus

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

17 Commits
 
 
 
 
 
 
 
 

Repository files navigation

pileup consensus

Script takes samtools pileup output and write consensus fasta file as output.

usage: pileup_consensus.py [-h] [-v] [-n BASE] [-r READS] [input] [output]

Produce consensus fatsa from samtools pileup output

positional arguments:
  input                 pileup file 
  output                Output file if no file result will be directed to
                        stander output

optional arguments:
  -h, --help            show this help message and exit
  -v, --version         show program's version number and exit
  -n BASE, --base BASE  Number of bases in line. (default: 80)
  -r READS, --reads READS
                        Minimum number of reads to accept nucleotide as a
                        variant. (default: 2)

pileup_consensus.py version 0.01. use command -h for info.

About

Script that takes samtools pileup output and write consensus fasta file as output.

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages