I have 10 bacterial genomes FASTA sequences (few of them might have chromosome and plasmids together and rest only chromosome) . And I have another FASTA file with 20 different plasmids. My plan is to remove the plasmid part from those 10 bacterial genomes and keep only chromosomal part as FASTA form if any section of bacterial genome match 100% to any of my 20 plasmids list. What are the best ways to do it?
1 answer
I would take a look at the Biostrings R package for this. You can read your FASTA sequences in, then simply use the `maskMotif` method to mask your plasmid sequences and then you can use `injectHardMask` to replace your masked sequences with a letter of your choice. This way you can still identify where there were sequences, and the total length of your chromosome doesn't change.
Log in to answer this question.
-
Include plasmid sequences in RNA-seq differential expression analysis for bacteria?
written by MaxMin •
-
Bacterial plasmid analysis
written by nicole.kavanagh •
-
Should we remove plasmid sequences for Whole genome Alignment
written by Optimist
-
Plasmid and Phage identification package
written by 3335098459 •
-
Determine number and length of plasmids in bacterial genome (using something like Prokka)
written by suzuBell •
-
How to interprete PlasmidSeeker results when you have some plasmid clusters?
written by The Terminator •
-
How to interprete PlasmidSeeker results when you have some plasmid clusters?
written by The Terminator •
-
How to separate the plasmid protein files from the main chromosomal protein file into a GenBank bac…
written by melisechaves •
-
How to remove sequences from a scaffold/contig file?
written by mandigene94 •
-
Finding Plasmids
written by bioinfo
