Comments (6)
Hi @natir
I'm affraid that, in this particular case, yacrd is not correctly identifiyng the chimeras that I have already identified manually
I've been trying several values for -c and I get many "not covered" reads and not all the chimeras that it should identify.
Thanks anyway for all your useful help
Regards
from yacrd.
Thanks for using yacrd
Important point:
- yacrd has been tested only on DNA direct sequencing data
- you used parameters optimized for read scrubbing not chimera detection
Can you give me a little more information on the sequencing protocol, PCR amplifications, read selection, nanopore sequencing cell version?
from yacrd.
Hello again,
Thanks for your quick response.
The data we work with have been obtained from an amplification of the complete 16RNA gene from a known mixture of bacteria,( Gut Zymobiomics standard, (about 15 different species)).
The amplicon obtained has been sequenced with the R10 4 1 ONT flowcell and then classified with EMU (https://www.nature.com/articles/s41592-022-01520-4).
In a detailed analysis we have seen that out of about 20000 reads 5% are chimerical, being able to identify the origin of the chimerical reads. We have seen that Latorre (https://github.com/adlape95/Spaghetti) says that in a similar assay he is able to identify and eliminate chimeras using yacrd, but we have not been able to reproduce the results with our data, as we did not detect any chimera.
We have also already tried this instructions, as you suggested, without results.
minimap2 -x ava-ont reads.fq reads.fq > overlap.paf
yacrd -i overlap.paf -o reads.yacrd
Thanks in advance for your help.
from yacrd.
I tried to contact the authors of Spaghetti to tell them that they didn't use yacrd with a good parameter for chimera detection in my opinion. It's interesting to scrub reads, but it's not the same operation as chimera detection. But I haven't been able to get an answer.
Since you're multiplying reads by PCR, I'm assuming you've got a very high coverage rate, but you may also have duplicating chimeras.
Yacrd analyzes the coverage of each read to determine whether it is chimeric or not. You can try increasing the value of the -c
parameter, which corresponds to the minimum coverage value at which a read is considered to have a poor quality region. If a read has a poor-quality region in the middle, it's considered chimeric.
I recommend you play with this parameter on a subset of reads, including reads you're sure are chimeric, to find the best -c
values for your dataset.
This should enable you to discover reads that have become chimeric after or during PCR multiplication. But for reads that were "generated" before PCR multiplication and therefore multiplied by it, yacrd won't be able to help you.
from yacrd.
Hi @natir ,
First of all, thanks a lot for all your help and efforts with yacrd.
I've been making some trials modifying the -c parameter. I've tried values going from 50 to 10000 but, although now yacrd finds some chimeras, I'm not getting the amount of chimeric detections needed.
I'll try some more values to check if I can get all the chimeras from the reads.
Thanks again for your quick responses and for your help.
from yacrd.
Be careful with very large values of -c
, there may be many false positives.
I'm interested in feedback to possibly add recommendations to the Readme.
I close the issue, but do not hesitate to repost a message with your conclusion.
from yacrd.
Related Issues (20)
- Kernel panic at 'called `Option::unwrap()` HOT 2
- Reduce memory usage HOT 1
- escape characters `/` and `\` in ondisk filename HOT 1
- thread 'main' panicked at 'slice index starts at 8092 but ends at 7507' HOT 12
- Error: --input provided more than once.... HOT 2
- yacrd parameters optimization HOT 4
- Filter from .yacrd report HOT 8
- Faster yacrd HOT 5
- Mismatched types HOT 2
- interpreting results HOT 5
- What pipeline should I use? HOT 7
- raw fastq file sequences ID did not fully appear in the *report.yacrd HOT 5
- how about using -X when running minimap2? HOT 1
- Compilation error due to incompatible clap version? HOT 3
- Why the overlap? HOT 2
- Error while running yacrd HOT 6
- Order of Chimera detection and scrubbing HOT 2
- yacrd for detecting chimeras in amplicon sequences HOT 3
- many metagenomic samples has redundant reads but not errors HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from yacrd.