You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Hi,
I'm using SURVIVOR (v1.0.7) for the first time and trying out some test data.
I've generated a params file using SURVIVOR simSV test_params.param
and then simulated the SVs onto a reference genome. SURVIVOR simSV ref.fasta test_params.param 0.1 0 simulated
Upon viewing the generated simulated.vcf file, I see that the VCF file is using zero-based coordinates when VCF uses 1-based coordinates.
Here are the first 16 characters of my reference sequence:
>SJAP_Chr1
CACCAAAAACCCTAAG
Here are the first 16 characters of my simulated sequence:
>SJAP_Chr1
AACAAAAAAGCCTAAA
and here are the lines in the VCF file that relate to the variants (the VCF header states 'source=Sniffles').
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT
SJAP_Chr1 0 SNP0SURVIVOR C A PRECISE;SVMETHOD=SURVIVOR_sim;SVLEN=1 GT:GL:GQ:FT:RC:DR:DV:RR:RV 1/1
SJAP_Chr1 3 SNP1SURVIVOR C A PRECISE;SVMETHOD=SURVIVOR_sim;SVLEN=1 GT:GL:GQ:FT:RC:DR:DV:RR:RV 1/1
SJAP_Chr1 9 SNP2SURVIVOR C G PRECISE;SVMETHOD=SURVIVOR_sim;SVLEN=1 GT:GL:GQ:FT:RC:DR:DV:RR:RV 1/1
SJAP_Chr1 15 SNP3SURVIVOR G A PRECISE;SVMETHOD=SURVIVOR_sim;SVLEN=1 GT:GL:GQ:FT:RC:DR:DV:RR:RV
As you can see, the first variant is at position zero and then the subsequent co-ordinates are also zero-based. The actual 'POS' values should be 1, 4, 10, and 16.
Based on my initial use of SURVIVOR, I have two other comments:
The 'simulated.bed' file is not in BED format, so perhaps it would be better to give it a different file extension
I can't find any command line help for the simSV method. Typing SURVIVOR simSV --help gives the stdout of Parameter file generated and produces a file called --help
The text was updated successfully, but these errors were encountered:
Hi,
I'm using SURVIVOR (v1.0.7) for the first time and trying out some test data.
I've generated a params file using
SURVIVOR simSV test_params.param
and then simulated the SVs onto a reference genome.
SURVIVOR simSV ref.fasta test_params.param 0.1 0 simulated
Upon viewing the generated
simulated.vcf
file, I see that the VCF file is using zero-based coordinates when VCF uses 1-based coordinates.Here are the first 16 characters of my reference sequence:
Here are the first 16 characters of my simulated sequence:
and here are the lines in the VCF file that relate to the variants (the VCF header states 'source=Sniffles').
As you can see, the first variant is at position zero and then the subsequent co-ordinates are also zero-based. The actual 'POS' values should be 1, 4, 10, and 16.
Based on my initial use of SURVIVOR, I have two other comments:
SURVIVOR simSV --help
gives the stdout ofParameter file generated
and produces a file called--help
The text was updated successfully, but these errors were encountered: