DRL/blobtools

Blobtools for symbiosis

Sandrine3 opened this issue · 2 comments

Sandrine3 commented

Hi,

I have been using blobtools to plot the diversity of symbiotic bacteria found in a plant. Everything seems to work great, the .bam file shows about 98% mapped reads depending on sample. However, when I get to plot the results, the coverage plot that is produced shows about 50% mapped reads. Why the difference?

Thank you very much for your support!
Sandrine

nostoc_blob_M2 blobDB json bestsum order p7 span 100 noscale blobplot cumulative read_cov bam0

Huangtianyu-caas commented

Hi,

I got the same problem that the Unmapped rate is so high, and I don't know why. Do you have any update?

My best,

Tianyu

DRL commented

Hi both,

this is how coverage parsing works: https://blobtools.readme.io/docs/coverage-parsing#1-bam-files

So i would assume that the difference is due to reads that map in supplementary/secondary alignments

cheers,

dom