Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Sim-it not producing any long reads #20

Open
Oieswarya opened this issue Nov 13, 2023 · 1 comment
Open

Sim-it not producing any long reads #20

Oieswarya opened this issue Nov 13, 2023 · 1 comment

Comments

@Oieswarya
Copy link

Oieswarya commented Nov 13, 2023

I am trying to simulate a set of HiFi reads 10x coverage of Pea Aphid genome. I am getting this in the error script:

rm: cannot remove `Sim-it-master/Output/PeaA/Long_reads_to_merge_PeaA_NW_021758228.1_Acyrthosiphon_pisum_isolate_AL4f_unplaced_genomic_scaffold,_pea_aphid_22Mar2018_4r6ur_Scaffold_10000': No such file or directory

rm: cannot remove `Sim-it-master/Output/PeaA/Long_reads_to_merge_PeaA_NW_021758229.1_Acyrthosiphon_pisum_isolate_AL4f_unplaced_genomic_scaffold,_pea_aphid_22Mar2018_4r6ur_Scaffold_10001': No such file or directory

rm: cannot remove `Sim-it-master/Output/PeaA/Long_reads_to_merge_PeaA_NW_021758230.1_Acyrthosiphon_pisum_isolate_AL4f_unplaced_genomic_scaffold,_pea_aphid_22Mar2018_4r6ur_Scaffold_10002': No such file or directory

Can you please help me.

This is my config.txt:

Project:

Project name = PeaA
Reference sequence = PeaAphid_Genome.fa
Replace ambiguous nts(N) =
Max threads = 16

Structural variation:

VCF input =
Foreign sequences =

Deletions = 0
Length (bp) = 30-150000

Insertions = 0
Length (bp) = 30-100000

Tandem duplications = 0
Length (bp) = 50-10000
Copies = 1-20

Inversions = 0
Length (bp) = 150-1300000

Complex substitutions = 0
Length (bp) = 30-100000

Inverted duplications = 0
Length (bp) = 150-350000

Heterozygosity = 60%

Long Read simulation:

Sequencing depth = 10
Median length = 10000
Length range = 10000-15000
Accuracy = 99.9%
Error profile = error_profile_PB_Sequel_CCS_hifi.txt

This is my job script to run:
perl Sim-it1.3.4.pl -c t_config.txt -o Sim-it-master/Output/PeaA/

@ndierckx
Copy link
Owner

Sorry for the late response
It is probably too late now, but I uploaded a new version that should have resolved your problem.
The read ids of your assembly file were too long to be used as a file name

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants