snippy icon indicating copy to clipboard operation
snippy copied to clipboard

snps.tab file not created

Open jennifercalawa opened this issue 1 year ago • 0 comments

Hello, I am trying to run snippy-multi on a MacBook Pro M2 and am running into an error. Here is the snps log output:

cd /Users/jennifercalawa/Desktop/FASTAforRoary/FNAonly

/Users/jennifercalawa/anaconda3/envs/snippyenv/bin/snippy --outdir CM1.1 --ctgs /Users/jennifercalawa/Desktop/FASTAforRoary/FNAonly/ncbi_dataset/CM1.1/CM1.1.fna --ref Vfischeri_ES114.all.gbk --cpus 6

samtools faidx reference/ref.fa

bwa index reference/ref.fa

[bwa_index] Pack FASTA... 0.04 sec [bwa_index] Construct BWT for the packed sequence... [bwa_index] 0.93 seconds elapse. [bwa_index] Update BWT... 0.03 sec [bwa_index] Pack forward-only FASTA... 0.03 sec [bwa_index] Construct SA from BWT and Occ... 0.30 sec [main] Version: 0.7.17-r1188 [main] CMD: bwa index reference/ref.fa [main] Real time: 1.329 sec; CPU: 1.330 sec

mkdir reference/genomes && cp -f reference/ref.fa reference/genomes/ref.fa

mkdir reference/ref && bgzip -c reference/ref.gff > reference/ref/genes.gff.gz

snpEff build -c reference/snpeff.config -dataDir . -gff3 ref

WARNING: All frames are zero! This seems rather odd, please check that 'frame' information in your 'genes' file is accurate.

(bwa mem -v 2 -M -R '@RG\tID:snps\tSM:snps' -t 6 reference/ref.fa fake_reads.fq | samtools view -u -T reference/ref.fa -F 3844 -q 60 | samtools sort --reference reference/ref.fa > snps.bam)

[main] Version: 0.7.17-r1188 [main] CMD: bwa mem -v 2 -M -R @RG\tID:snps\tSM:snps -t 6 reference/ref.fa fake_reads.fq [main] Real time: 16.675 sec; CPU: 92.473 sec

samtools index snps.bam

samtools depth -aa -q 20 snps.bam | bgzip > snps.depth.gz

tabix -s 1 -b 2 -e 2 snps.depth.gz

fasta_generate_regions.py reference/ref.fa.fai 181047 > reference/ref.txt

freebayes-parallel reference/ref.txt 6 -p 1 -q 20 -m 60 --min-coverage 10 -V -f reference/ref.fa snps.bam > snps.raw.vcf

env: python: No such file or directory

And the terminal output says there is no .tab file (what I think is meant to be snps.tab) in my output folder. From what I can tell, all the other output files are being made, but the snps VCF is empty. Any help is greatly appreciated!

jennifercalawa avatar May 16 '23 17:05 jennifercalawa