Skip to content

length of the reads #54

@parsboy66

Description

@parsboy66

Thanks for the nice job!
I was trying to simulate 100K cDNA reads which each read consist of:
BC + UMI + Polya(t)+ 1500 (random seq)
at first, i couldn't ran it all over 100K in one fasta file and at the end I spliced reads in 10 batches.
so each read consists of 1570 ~ 1590 (variable length cos of the poly a/t) i set the -l or mean read length to 800
but by using this option, I got the error that says the input genome is null, while by just not using -l, it ran correctly.

i am wondering to know while the mean read length in my sample is around 800 while default in deepsimulatre is 8000 , does it affect the output??? and also is there anyways to run 100k in one fasta file?? (error was like list was too long ) .

i ran the instruction like this:
./deep_simulator.sh -B 2 -l 800 -i ../nano_fasta_small_9.fasta -n -1 -o ../output_9

Thanks

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions