Repository navigation
question about “--sams=bam.list” #157
Description
Activity
This is my command.thanks
Hi, the bamlist format looks right to me. Can you check if those are the the correct read paths, there are no whitespaces/tabs in the bamlist, and that you have permissions to read the data by using i.e. samtools
samtools view /data/..../2.bam/I1.bam Chr1:1-1200000 | less -SBest,
HannesThanks for your reply.
I tried the following command as you suggested, these results do not seem to be problematic:
The newline style needs to be UNIX (
\n) but it looks you have the DOS style (\r\n) from yourcat -Aoutput. You can convert to UNIX style withsed $'s/\r$//' tst.bam.list > unix.tst.bam.list cat -A unix.tst.bam.list # Should not print the ^M in output
and use
unix.tst.bam.listin graphtyperDamn, I didn't catch the formatting problem, thanks for your careful answer. The program is running.
One more small question, how many samples does the program support at a time? I have 5,000 samples
I ran graphtyper with my samples and got this error.alright, well that error message doesn't explain much. Only that some std::map had an access problem, somewhere. Is there really nothing else helpful in your logs?
And it would help me a lot if you could demonstrate your problem in a way for me to reproduce it
graphtyper genotype_sv ${reference_fasta} 5k.svimmer.vcf.gz --sams=5k.bam.list --region=Chr${1}
--threads=52 --vverbose --log=5k.Chr${1}.log --output=5k.graphtyper.SV
5k.Chr12.log
0_graphtyper_4481436.log
Hi, this is my command line and all the logs I can get. I hope this can help.
hi,
I am using the genotype_sv function of V2.7.x, my bamlist is shown in the picture, but it will report an error, is my bamlist format wrong? When I use "--sam" to give a separate bam file program can run.
bamlist:
err.log:
looking forward your reply
best wishes
ding