Hi, I have a fasta file representing a reference genome and I would like to modify it to mask variable sites when I map variable sites. I'm interested in doing this because I have bisulfite reads from several related species, but BSmap and Bismark don't offer an option to mask variable sites while mapping.
The initial genome is in a fasta file. The sites I would like masked in a vcf file.
Thank you!
1 answer
If you can convert your VCF to BED format (see Converting a VCF with SNPs and indels to BED format) you can use the pyfaidx faidx command to mask your FASTA file with a special character, or as lowercase letters:
vcf2bed < variable_sites.vcf | faidx genome.fasta --bed - -m
Note that the -m and -M options will modify your FASTA file in-place, so you probably want to make a copy first.
•
written
by
Matt Shirley
Log in to answer this question.
-
Denovo assembly in Bisulphite sequencing
written by Smilesky •
-
STAR bisulfite mapping
written by noodle
-
News: Comprehensive benchmarking of software for mapping whole genome bisulfite data
written by David Langenberger
-
About Mapping Efficiency of Bisulfite aligners
written by chaudharyc61 •
-
Reference based assembly
written by deepti1rao •
-
more than 28 million methylation sites?
written by ccagg
-
BSMap to Bismark bam file formats
written by Jautis
-
Genotype calls for a list of variable sites
written by Jautis
-
Getting reads that end at a digest site from a sam/bam file
written by Jautis
-
Publicly Available Pipeline For Analysis Of Whole Genome Bisulfite Sequencing Data
written by skm770 •
