adVNTR
adVNTR copied to clipboard
sqlite3.OperationalError: disk I/O error
Hi Mehrdad,
I'm trying to run advntr on some PacBio reads mapped to hg38 and I'm using the hg38 models (hg38_selected_VNTRs_Illumina.db) however I get a sqlite3.OperationalError
when running with the following options:
advntr genotype --alignment_file HG00733.sorted.bam --working_directory $PWD --pacbio -r hg38.no_alt.fa --models vntr_data/hg38_selected_VNTRs_Illumina.db --outfile advntr.vcf -
-outfmt vcf
Using Theano backend.
Traceback (most recent call last):
File "/software/anaconda3/4.9.2/lssc0-linux/envs/advntr-1.4.1/bin/advntr", line 11, in <module>
sys.exit(main())
File "/software/anaconda3/4.9.2/lssc0-linux/envs/advntr-1.4.1/lib/python3.6/site-packages/advntr/__main__.py", line 134, in main
genotype(args, genotype_parser)
File "/software/anaconda3/4.9.2/lssc0-linux/envs/advntr-1.4.1/lib/python3.6/site-packages/advntr/advntr_commands.py", line 104, in genotype
reference_vntrs = load_unique_vntrs_data()
File "/software/anaconda3/4.9.2/lssc0-linux/envs/advntr-1.4.1/lib/python3.6/site-packages/advntr/models.py", line 140, in load_unique_vntrs_data
right_flanking, repeats, scaled_score FROM vntrs''')
sqlite3.OperationalError: disk I/O error
All the files referenced in the command exist. I tried redownloading the VNTR models in case they were corrupted. I get the same error if I use the hg19 ones. Any suggestions? Thanks.