Average Coverage after Assembly (Spades)
0
0
Entering edit mode
5 months ago
davidmaimoun ▴ 50

Hi, I would like to know to average coverage of my assemblies (gotten by spades). I ran bbmap to mapping my reads (the same ones I used to make the assembly) to the assembly. But I am not really sure that is the way and the answers on the different forums didn't completely help.

Command

bbmap.sh in=./trimmed/NM167_S34_R1.fastq in2=./trimmed/NM167_S34_R2.fastq ref=./assembly/NM167_S34.fasta covstats=covstats.txt

I got output like this:

#ID Avg_fold    Length  Ref_GC  Covered_percent Covered_bases   Plus_reads  Minus_reads Read_GC Median_fold Std_Dev

NODE_1_length_79421_cov_85.482448   261.1709    79421   0.5185  100.0000    79421   60378   60230   0.5259  253 111.34

....

I will really appreciate your help!

Thank you

coverage average assembly • 231 views
ADD COMMENT
0
Entering edit mode

I see now I got the average, but it's printed on the terminal and doesn't appear in the covstats.txt file

Average coverage: 164.398 Standard deviation: 155.626 Percent scaffolds with any coverage: 84.91 Percent of reference bases covered: 99.56

And I don't know how to set it in a file I tried with the redirect output ">" but it is not working

Any Idea ?

Thank you

ADD REPLY

Login before adding your answer.

Traffic: 779 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6