Hi, everyone! I assembled a fungal genome using abyss-pe (ABySS, version 1.9.0). I notice that there are some degenerate nucleotides in the resulted scaffolds, e.g. [RYKM]. I read the readme file carefully but cannot find which parameter is about this.
Does anyone know how to turn off the degenerate nucleotides? Thank you very much!
•
written
by
vtefnfqp •
2 answers
PathConsensus outputs the IUPAC ambiguity codes. Set either a=1 or p=1 to disable PathConsensus ambiguity codes. Note that this will also disable PopBubbles. You may prefer to simply replace ambiguity codes with an ACGT nucleotide or N after assembly, as Ben suggested.
See https://github.com/bcgsc/abyss/blob/master/bin/abyss-pe#L264-L267
•
written
by
Shaun Jackman
Unfortunately, I don't think there is an ABySS option to turn off degenerate nucleotides. (I just did a grep through the source code to make sure.)
If the ambiguity codes are causing problems in downstream tools, you will have to "flatten" the bases as a post-processing step after assembly. The unix 'tr' program might be helpful for that.
Thanks! I replace the degenerate nucleotides as "N", and then the replaced "Ns" can be replaced by "ATCG" in the gapfilling step.
•
written
by
vtefnfqp •
Log in to answer this question.
-
ABySS Process Killed
written by asarcillo •
-
PSMC -p parameter
written by Alba •
-
update genotypes in vcf
written by biosol
-
Multiple alignment: compare one sequence with many
written by mdihcroe •
-
Masking short sequences between gaps in genome assembly
written by pbigbig
-
The problem about tophat software -G parameter
written by Wayne Lee •
-
Trimmomatic : keepBothReads True vs False option
written by candida.vaz •
-
Reading single-end reads takes forever
written by RoseString •
-
Newbler: Merged Miseq Reads + Mate Pairs?
written by Adrian Pelin
-
Abyss Speed-Up Tricks
written by Darked89
