Comments (4)
Hi,
I'm not sure about the reason. Could you provide me with the command you used for the two different references?
Thanks,
Chen
from nanosim.
Hi @cheny19 ,
This is the commands:
nanosim-h -n100 --max-len 15000 --min-len 25000 -p ecoli_R9_1D reference1.fa
nanosim-h -n100 --max-len 15000 --min-len 25000 -p ecoli_R9_1D reference2.fa
The references has same length.
(So NanoSim look for homopolymers length or the errors probabilities is the same for sequences of same length?)
Thanks,
from nanosim.
Sorry for the late reply. It seems you were asking for min_len to be larger than max_len based on the commands above.
from nanosim.
Dear @belaidH
NanoSim models the error rate, length distribution and other characteristics of the input read and then in the simulation stage, you give it a reference and NanoSim applies those profile and produces some synthetic reads.
Therefore, it means that when you use a predefined model, no matter if you change the reference file in the simulation stage, your simulated reads are going to have similar error rates and length distribution. That's because you use the same profiles for both of them!
I hope it explains your question. Besides, as Chen suggested, you are setting min length higher than max length which is not correct.
Finally, use the latest version, read the documentation and README file and let me know if I can be of more help. I am closing this issue now.
from nanosim.
Related Issues (20)
- Could not retrieve index file for alignments HOT 2
- Abnormal error distribution in reads HOT 3
- Simulating exact number of reads for each input -- transcriptome mode HOT 2
- Reads come mostly from one chromosome only HOT 1
- hangs at start of simulation HOT 5
- Can NanoSim be applied to simulate single cell Nanopore reads HOT 1
- Invalid fastq format HOT 4
- mis.hist missing when running characterization stage
- file has no sequences defined (mode='r') - is it SAM/BAM format? Consider opening with check_sq=False HOT 1
- Nanosim hangs in the middle HOT 18
- Infinite loop in function extract_reads in metagenome mode when length equals max length HOT 2
- Transcriptome mode error rate tsv explanation HOT 2
- Models for R10.3 or R10.4 flow cell
- Option to specify desired read coverage or sequencing depth HOT 2
- ValueError: Found array with 0 sample(s) (shape=(0, 1)) while a minimum of 1 is required. HOT 6
- Please specify the training reads and its reference genome! HOT 3
- Stuck at simulation stage HOT 4
- simulator.py genome FileNotFoundError: [Errno 2] No such file or directory: 'training_model_profile' HOT 1
- NanoSim for tuning Minimap2 parameters? HOT 2
- Models for newer versions of Guppy with sup basecalls 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 nanosim.