InSilicoSeq icon indicating copy to clipboard operation
InSilicoSeq copied to clipboard

coverage file seems not working properly in setting --coverage

Open kaysahu opened this issue 9 months ago • 1 comments

I want to apply coverage file in the command to override the number of read settings. I have tested with two genomes to generate hiseq styled read files. In one coverage.csv file, I put both coverage at 100, and in another coverage.csv file, I double the coverage at 200... However, in both cases, there gives me the same size for R1 and R2. Increasing the coverage doesn't help to increase the data amount. When I look at the error file, it clearly says that the coverage file disabled the number of reads, but later it still shows the number of reads is by default 1 million.

Kindly please advice regarding this issue, and if there is an efficient way to use the coverage file, without using 'reads' as benchmark.

Image

kaysahu avatar Apr 09 '25 22:04 kaysahu

I also witnessed the similar problem, i have setted my coverage file as follows:

Image

however,after generating the sequence , i count the read number of different genome the file used, and the results showed below:

Image

and it seems that the generator only caught the top two line of coverage in my coverage file.

Attached file coverm_utf8.txt is the --coverage_file i use; and counts_summary.txt is the file for read number of different genome the file used.

coverm_utf8.txt

counts_summary.txt

thank you!

heyinghui22 avatar May 28 '25 06:05 heyinghui22