Chop off sequences before aligment
Hello!
I have a lot of sequences that have been obtained using different sets of primers. I want to align them, but first I have to remove that are not similar and won't align. Is there any way to do it automatically and not by hand? Thank you!!!
1 answer
To remove primer sequences, try cutadapt.
Log in to answer this question.
More posts like this
-
Missing primers in multi amplicon sequencing
written by madalton •
-
Meta-analysis using 16S rRNA NGS data
written by bala •
-
Realign previous alignments
written by QLFblaireau •
-
Is blacklisted regions filtering in ChIP-seq still needed? (let's talk about the state of the art)
written by msimmer92
-
What are the available automated microscopy options in 2018?
written by ceruleanivy •
-
What can I use to identify primer sequences in my fasta files?
written by c.older •
-
Finding adapters in sequences
written by genomics Newbie •
-
Trimming out primer sequences in the middle of reads
written by s.kyungyong64 •
-
Command Line Degenerate Primer Design tools
written by adp7 •
-
Automated Primer-Blast like funcionnality
written by Eric Normandeau
It is not clear what you want to do here.
As @b.nota says are you worried about parts of the reads that are not expected to align with the reference or are you referring to contaminants that may be full length reads?
In either case an aligner should be able to tackle those reads and handle them appropriately as long as you provide a correct reference.