MethylCoder
Application data |
|
Created by | Pedersen BS |
---|---|
Biological application domain(s) | Genomics, Sequencing, DNA methylation, Epigenomics |
Principal bioinformatics method(s) | Read mapping, Bisulfite mapping |
Technology | Illumina |
Created at | UC Berkeley |
Maintained? | Yes |
Input format(s) | FASTQ, FASTA |
Output format(s) | binary, TXT, SAM |
Programming language(s) | Python, C |
Licence | BSD |
Operating system(s) | Linux, Linux 64, Mac OS X |
Summary: Pipeline for fast, simple processing of BiSulfite-treated reads into methylation data. Includes scripts for analysis and visualization. In addition to a binary output, the direct output of methylcoder is a text file that indicates per-nucleotide methylation context (CG/CHG/CHH) and methylation levels (both coverage and C-T conversions)
"Error: no local variable "counter" was set." is not a number.
Pipeline for fast, simple processing of BiSulfite-treated reads into methylation data. Includes scripts for analysis and visualization. In addition to a binary output, the direct output of methylcoder is a text file that indicates per-nucleotide methylation context (CG/CHG/CHH) and methylation levels (both coverage and C-T conversions). The direct output of methylcoder is a text file that looks like
#chrom mt bp c t 1 3 1354 0 1 1 3 1358 0 1 1 3 1393 0 1 1 3 1394 0 1 1 3 1402 0 1 1 3 1409 0 1
where columns are reference seqid methylation context (type) basepair location(bp) number of reads where a (c)ytosine was unconverted, number of reads where where a cytosine was converted to (t)hymine. Making methylation at every methylable basepair easily calculated as c / (c + t).
In addition, a SAM file with corrected positions and sequences is created for use with other bioinformatics packages.
Can use either bowtie or gsnap as the aligner. Example invocation:
methylcoder \ --bowtie /usr/local/src/bowtie-0.12.5/ \ --reference hu19.fa \ --extra-args "-m 1 -e 100" \ reads/r_1.fastq reads/r_2.fastq
Links
- MethylCoder [ edit link ]
References
To add a reference for MethylCoder, enter the PubMed ID in the field below and click 'Add'.
Search for "MethylCoder" in the SEQanswers forum / BioStar or:
Web Search | Wiki Sites | Scientific |
---|---|---|