This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Abyss results

I am trying to use abyss, but even if I get the results and any error, the final file just have a few lines.

sample41-contigs.fa

>86 92 510 0+,13+
TGAGTCTCGTGTACCCATTCATTCATTTATTCATTCAGGTTAAGTTCCTATTTGAGGGCAGACCTTTGAGTCCAAGGAACTTGCCTTTTCCG
>87 100 510 1+,19+
TGACATGCTAGAAATGATATGAGTTGCTTCAACCTAGTTAAAGTGAAAATTTTAGTAGGAAATGAGGAAGGTGGCCTAGCCCATGACAGCATCAGGAACT
>88 94 510 16+,31+
TAGATGCTGGAATAAAACATAAAAACAAAACAAAGACCTTGGTGTAATGTACTATACAGGGAAGGGAATCAGATATTAAATACATATGTAGTAA
>89 97 510 20+,35+
CTGATTCTACATTATGATGAGTTGTATAATTATTTCATTATGAATTACAATGTAATAATAATAGAAATAAACTGCACAACAAATGTAATGCACTTGA
>90 100 510 33+,60+
AATAAAGACACATACTTCATGTAGTCAGAAGAGTGGTCTGAAGCAAAATATTCAAAGCCTATTGAAGAAATCATTTTAGAATTTTTCTACAAGTTTTTTC
>91 100 510 37+,59+
TTCTGTATTGCATTTCCGAACACTGGTGAGGAGTGCTTGGACCTGTGCCCTTTGCTGTCCAGGCTGCAATTTTCTCGTTGAGTATCACGTGAATATTTAT
>92 100 510 61+,85+
TAGGTGGGAGAATGCTTGAACCCAGGAGTTCAAGTCAAGCCAGGCAACACAGTGAAACCCTGTCTCTAAAACATAAATAAATAAAATAAAAATGGGACTT
>93 99 510 62+,84+
GTTTAAAATGTATATTTTACATACTTATTCTTCATTAGTGTTATGTAATATTAATGTCCTAAAGTCAGTTAACTGTAGCAACAACTTATGAGTTTTTCA

The samples looks like this:

@SRR062641.1 HWI-EAS110_103327062:5:1:1091:7885/1
GAAAGAAAGAAAGTCAACTGTATGCTTAAAAATCCAAGTTGTGGGTGGGAAGCTGATTGAATTTTTTACTACGGTTCATAAAAAAACACAAGACTCACAT
+
@7:6+)=0577''3=?->>;:A>A############################################################################
@SRR062641.2 HWI-EAS110_103327062:5:1:1092:19656/1
AGACTTTATCAAATTATAACTGGGATTAGATAGACAGTCCAGAATAATTTTAATAATTCATGTACATGCTTCATGTATTTTCCTTGCTTACACTTGCCTA
+
=@>@6?>A;AA?5?DDB5DD=DAD5-=?6@::A-ADDCD?C55,CB:AD?BD:5DD?BD=:5C@5C##################################
@SRR062641.3 HWI-EAS110_103327062:5:1:1092:17219/1
TTGAAGGATACAAAGTATTGATCCTGGGTGTGTCTGTGTGGGTGTTGCCAAAGGAGACTAACAGTGGAGTCAGGGGGCGGGGAAAGGCAGACCACCTGCT

and the 2 file

@SRR062641.1 HWI-EAS110_103327062:5:1:1091:7885/2
TCAATGAATGCACTTTAATTTTTTAAAATNNNNNNANATNATGATNTANTTTANNNNTNCAANAGAAGTGTTGGGTTTTTTTATGAACCGTTTTATCAAA
+
5DD=-B#######################!!!!!!#!##!#####!##!####!!!!#!###!#####################################
@SRR062641.2 HWI-EAS110_103327062:5:1:1092:19656/2
TCAATAACTTTATGGATAGTGTATCATCTNNNNNNCNAGNGAAATNTTNATACNNNNGNAACNACTGAAAATCTAGTTCATTGTTAGGCAAGATTCACCA
+
5DB-D-BDDD=D=D:--<9:-??6@3738!!!!!!6!15!;0?A9!;5!,782!!!!;!16*!25?25+=?<+@<CC=-:A?><A)4;->A:55:+=466
@SRR062641.3 HWI-EAS110_103327062:5:1:1092:17219/2
TGGGAGGCAAGGCCAGACTCACGTTCTTANNNNNNANTANCGGCTNTCNTCCGNNNNGNCGGNCAGCTGGGTAGGGTGTGTCCTCCCCCCTCAGGGGTGG

I run abyss like this:

abyss-pe k=88 name=sample41 B=50 in='SRR062641_1.fastq.gz SRR062641_2.fastq.gz'

The final report of abyss is this:

Best scaffold N50 is 100 at n=15 s=100. 



n       n:100   L50     min     N75     N50     N25     E-size  max     sum     name 

8       4       4       100     100     100     100     100     100     400     n=15 s=100 

abyss-stack-size 65536 PathConsensus  --dot -k88  -p0.9  -s sample41-7.fa -g sample41-7.dot -o sample41-7.path sample41-6.fa sample41-6.dot sample41-6.path 

Running with max stack size of 65536 KB: PathConsensus --dot -k88 -p0.9 -s sample41-7.fa -g sample41-7.dot -o sample41-7.path sample41-6.fa sample41-6.dot sample41-6.path 

Ambiguous paths: 0 

Merged:          0 

No paths:        0 

Too many paths:  0 

Too complex:     0 

Dissimilar:      0 

cat sample41-6.fa sample41-7.fa \ 

        | MergeContigs   -k88 -o sample41-8.fa - sample41-7.dot sample41-7.path 

ln -sf sample41-8.fa sample41-scaffolds.fa 

PathOverlap --overlap   -k88 --dot sample41-7.dot sample41-7.path >sample41-8.dot 

ln -sf sample41-8.dot sample41-scaffolds.dot 

abyss-fac   sample41-unitigs.fa sample41-contigs.fa sample41-scaffolds.fa |tee sample41-stats.tab 

n       n:500   L50     min     N75     N50     N25     E-size  max     sum     name 

8       0       0       0       0       0       0       0       0       0       sample41-unitigs.fa 

8       0       0       0       0       0       0       0       0       0       sample41-contigs.fa 

8       0       0       0       0       0       0       0       0       0       sample41-scaffolds.fa 

ln -sf sample41-stats.tab sample41-stats 

tr '\t' , <sample41-stats.tab >sample41-stats.csv 

abyss-tabtomd sample41-stats.tab >sample41-stats.md 

I dont understand the problem

assembly abyss

Highlight text you want to show as code first before clicking the 101010 button in the edit bar. Blockquotes are used for quoting someone else or a specific set of lines. I have edited the post this time.

1 answer

I suppose the mistake is in "B=50". This is a Bloom filter memory size, usually it should be larger (see https://github.com/bcgsc/abyss).

I dont think so because I am using the advise of abyss documentation: H. sapiens (~3.1Gbp): B=50G

Note the letter "G" after "50". It stands for gigabytes. The command in your post has "B=50", without "G". This means, ABySS used 50 bytes.

Totally right, I am gonna try with the G and see what happens

Log in to answer this question.