SDgenerator is a simple program to simulate a FASTQ, the product of 'sequencing' a sample conformed of a wild type virus and a known population of DVGs.
wgsim
(reads simulator)
To get a local copy up and running follow these simple steps.
-
Clone the repo in the directory of your choice
git clone https://github.com/MJmaolu/SDgenerator.git
-
Go to the DVGfinder directory
cd SDgenerator
-
Give execution permission to all the scripts in the SDgenerator directory
chmod -R +x .
-
Create a new environment with conda with all the dependencies needed to run SDgenerator
conda env create -f sdgenerator_env.yaml
-
Activate SDgenerator environment
conda activate sdgenerator_env
python3 SDgenerator.py -t population_csv -N number_total_reads [-l length_reads. Default 100] -o output_basename [-d outer_distance. Default 300]
BP,RI,DVG_type,proportion
6194,8961,Deletion_forward,0.1
23000,120,Deletion_reverse,0.1
9478,9489,3cb/sb,0.1
10000,9800,5cb/sb,0.1
2200,2100,Insertion_forward,0.1
1000,1200,Insertion_reverse,0.1
The proportion of the wild type genome is calculated as the substraction to 1 of the sum of the DVG proportions
-
Deletion_forward
-
Deletion_reverse
-
3cb/sb
-
5cb/sb
-
Insertion_forward
-
Insertion_reverse
María José Olmo-Uceda - maolu@alumni.uv.es
PhD student
EvolSysVir Group, I2SysBio (CSIC-UV)