How to remove all recombinant viral sequences from given multi fasta file
Hello,
I tried RDP (Recombination Detection Program) to remove recombinant sequences from a multifasta file. after excluding recombinant sequences when i re-ran RDP on sequences kept after first run, I again get some recombinant sequences.
Is there any way to remove all recombinant sequences in one run or I will have to re-run RDP until it doens't detect recombinant sequences anymore?
Thank You.
• 1,693 views
•
link
updated
by
Biostar
•
written
by
rajeshkumar_vinod •
•
written
by
rajeshkumar_vinod •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How should recombinant haplotypes from vg haplotypes personalized GBZ be exported as FASTA?
written by Etenesh •
-
SlingShot and Subsetting data
written by carolofharvest •
-
Trying to use hmmer to search genes over metagenome assembled genomes
written by v.berriosfarias •
-
Index building failed
written by Gaye •
-
How to remove multiple adapter sequences by Cutadapt
written by okiedokie1208 •
-
Can I use FLEXBAR to detect and remove barcodes out of my short reads without demultiplexing ?
written by giulia.trauzzi •
-
help with vcftools filtering on individuals
written by schisler •
-
Run transmembrane prediction with one multifasta file and multiple MSA sequences.
written by enkh.tug •
-
How to reduce protein sequence dataset size for machine learning?
written by rajeshkumar_vinod •
-
Best Program For Recombination Detection In Msa
written by GR
Have you tried something very simple such as cutadapt (http://cutadapt.readthedocs.io/en/stable/guide.html)?
In its most basic syntax, you can trim any recognisable sequence from your sequence files:
cutadapt -a AACCGGTT -o output.fastq input.fastqIn the options, I think that you can specify the program to remove any reads that contain the identified sequence.
If you have a known list of recombinant sequences, I think that cutadapt also allows you to specify a file with a list of sequences to remove.
Kevin