Convert abyss scaffolds to contigs
Entering edit mode
3.1 years ago
jm440 ▴ 10

Hi all,

I assembled a genome using paired end reads and abyss assembler with the following command:

abyss-pe k=128 in='R1.fq R2.fq'

I got a scaffold file, which I then filtered based on several criteria. I now want to get an updated contig file from this scaffold file. My strategy is to split the scaffolds into new lines based on the occurrence of 'N' similar to what is shown in this post: How to get contigs from scaffolds

I have two concerns though:

  1. I noticed in the scaffold file that there are also lower case Ns ("n") that appear (pretty rarely) in some scaffolds, see below for snippet from one of the scaffolds. What are these and is it safe to ignore these when splitting by "N"s?



2.Should I be splitting to contigs by a certain minimum # of "N"s or is 1 or more Ns the correct way to do this?

Thanks a bunch for any tips!

abyss contigs • 954 views
Entering edit mode
3.1 years ago
h.mon 34k

I am not very familiar with ABySS, but I am almost sure it should output a contigs file and a scaffolds file.

In any case, don't break the scaffolds at Ns, use the AGP file to correctly identify and break the scaffolds into contigs. See ABySS file formats and the NCBI AGP specification.

Entering edit mode

Yes ABySS does output a contig file along with a scaffold file, however I then filtered out scaffolds from the scaffold file using several criteria and so I want to produce an updated contig file based on this new scaffold file. Hopefully that makes sense..

I will check out AGP files - thank you!


Login before adding your answer.

Traffic: 2409 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6