Open thedam opened 10 months ago
Hey,
I run: seqToProfile -B 300 -b $bam -t $bed -v $vcf -r $ref > test.profile
seqToProfile -B 300 -b $bam -t $bed -v $vcf -r $ref > test.profile
then in the output test.profile I see:
test.profile
#model created at Fri Nov 17 09:56:43 2023 #reads: /path/to/my/test.bam bases: ACTG readLength: 54 binCount: 54 kmer: 3
Why readLength resluts in 54? Reads in my bam files are 101nt long.
After running simuReads my.config I get two fastq files with reads that ar 54nt long... I want them to be 101nt long
simuReads my.config
How to change the read length in the output fastq files?
Hey,
I run:
seqToProfile -B 300 -b $bam -t $bed -v $vcf -r $ref > test.profile
then in the output
test.profile
I see:Why readLength resluts in 54? Reads in my bam files are 101nt long.
After running
simuReads my.config
I get two fastq files with reads that ar 54nt long... I want them to be 101nt longHow to change the read length in the output fastq files?