Interpret coverage from mapping back de-novo assembly
0
0
Entering edit mode
3.9 years ago
nisrinalulu ▴ 10

Dear all,

I have assembly my reads with MEGAHIT to get final contigs. Then I mapping back my contigs with raw reads to get coverage rate with

bbmap/bbwrap.sh ref=megahit/IA_result/final.contigs.fa in=raw_data/IA_1.fq in2=raw_data/IA_2.fq out=IA.sam.gz kfilter=22 subfilter=15 maxindel=80

and

bbmap/pileup.sh in=IA.sam.gz out=IA_mapping.txt hist=IA_histogram.txt

I got result like this:

Average coverage:                       4.25
Percent scaffolds with any coverage:    99.91
Percent of reference bases covered:     98.83

How to interpret the result? Is the result better enough?

Thank you, Nisrina

assembly metagenomic next-gen • 603 views
ADD COMMENT

Login before adding your answer.

Traffic: 2680 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