I need to mask phages, mobile elements, and transposons before calling SNP's in my bacterial genome. How should I go about doing this? I have found the coordinates to mask by annotating my reference genome with RAST. I was planning on trying to mask my genome before mapping but I have found from this link: Which Aligners Recognize Soft-Masked Repeats In Reference Sequences? that this will lead to wrongly mapped sequences and spurious SNPs/indels.
I'm using the CFSAN SNP calling pipeline to call SNP's. Can I mask with samtools? Does anyone have a recommendations on how to mask phages, mobile elements, and transposons before calling SNP's?
•
written
by
taffners •
1 answer
You should be able to use bbmask.sh from BBMap suite for doing the masking. Take a look at bbmap/docs/guides/BBMaskGuide.txt once you download the software for instructions.
Log in to answer this question.
-
Code Review: Nextflow pipeline for phage-bacteria assembly from Nanopore data
written by eyaggelia146 •
-
Estimation of mobile elements insertions frequencies in heterogeneous samples
written by semion.leyn •
-
Mobile element discovery with MELT brings different results when chromosomes are analysed separate
written by optimistsso4co3 •
-
delete several ranges from a DNA sequence
written by oinkost •
-
Structural variants call from circular genome
written by rthapa •
-
Determine number and length of plasmids in bacterial genome (using something like Prokka)
written by suzuBell •
-
Comparing the WG sequenced isolate against the other same genus isolates on the public database
written by bioinforesearchquestions
-
Bacterial mobile elements database
written by hed.robin •
-
Masking repeats in your alignment
written by biogirl
-
Dealing with bacterial transposable elements during assembly
written by bioinfo
