Mercurial > repos > dfornika > snippy

annotate test-data/out/snps.log @ 1:e7e5da92d0ab draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
author dfornika
date Fri, 18 Jan 2019 19:35:46 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
1 ### cd /home/dfornika/Code/tools-iuc/tools/snippy/test-data
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
2
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
3 ### /home/dfornika/miniconda3/envs/snippy-4.3.6/bin/snippy --outdir out --ref wildtype.fna --R1 mutant_R1.fastq --R2 mutant_R2.fastq
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
4
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
5 ### samtools faidx reference/ref.fa
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
6
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
7
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
8 ### bwa index reference/ref.fa
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
9
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
10 [bwa_index] Pack FASTA... 0.00 sec
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
11 [bwa_index] Construct BWT for the packed sequence...
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
12 [bwa_index] 0.03 seconds elapse.
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
13 [bwa_index] Update BWT... 0.00 sec
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
14 [bwa_index] Pack forward-only FASTA... 0.00 sec
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
15 [bwa_index] Construct SA from BWT and Occ... 0.02 sec
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
16 [main] Version: 0.7.17-r1188
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
17 [main] CMD: bwa index reference/ref.fa
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
18 [main] Real time: 0.060 sec; CPU: 0.053 sec
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
19
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
20 ### mkdir -p reference/genomes && cp -f reference/ref.fa reference/genomes/ref.fa
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
21
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
22
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
23 ### ln -sf reference/ref.fa .
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
24
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
25
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
26 ### ln -sf reference/ref.fa.fai .
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
27
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
28
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
29 ### mkdir -p reference/ref && gzip -c reference/ref.gff > reference/ref/genes.gff.gz
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
30
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
31
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
32 ### bwa mem -Y -M -R '@RG\tID:out\tSM:out' -t 8 reference/ref.fa /home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R1.fastq /home/dfornika/Code/tools-iuc/tools/snippy/test-data/mutant_R2.fastq | samclip --max 10 --ref reference/ref.fa.fai | samtools sort -n -l 0 -T /tmp/snippy.19210. --threads 8 -m 1000M | samtools fixmate -m - - | samtools sort -l 0 -T /tmp/snippy.19210. --threads 8 -m 1000M | samtools markdup -T /tmp/snippy.19210. -r -s - - > snps.bam
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
33
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
34 READ 24871 WRITTEN 24858
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
35 EXCLUDED 0 EXAMINED 24871
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
36 PAIRED 24782 SINGLE 89
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
37 DULPICATE PAIR 10 DUPLICATE SINGLE 3
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
38 DUPLICATE TOTAL 13
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
39
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
40 ### samtools index snps.bam
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
41
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
42
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
43 ### fasta_generate_regions.py reference/ref.fa.fai 13379 > reference/ref.txt
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
44
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
45
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
46 ### freebayes-parallel reference/ref.txt 8 -p 2 -P 0 -C 10 --min-repeat-entropy 1.5 --strict-vcf -q 13 -m 60 --min-coverage 10 -F 0.05 -f reference/ref.fa snps.bam > snps.raw.vcf
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
47
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
48
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
49 ### bcftools view --include 'FMT/GT="1/1" && QUAL>=100 && FMT/DP>=10 && (FMT/AO)/(FMT/DP)>=0' snps.raw.vcf | vt normalize -r reference/ref.fa - | bcftools annotate --remove '^INFO/TYPE,^INFO/DP,^INFO/RO,^INFO/AO,^INFO/AB,^FORMAT/GT,^FORMAT/DP,^FORMAT/RO,^FORMAT/AO,^FORMAT/QR,^FORMAT/QA,^FORMAT/GL' > snps.filt.vcf
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
50
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
51
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
52 ### cp snps.filt.vcf snps.vcf
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
53
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
54
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
55 ### /home/dfornika/miniconda3/envs/snippy-4.3.6/bin/snippy-vcf_to_tab --gff reference/ref.gff --ref reference/ref.fa --vcf snps.vcf > snps.tab
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
56
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
57 Loading reference: reference/ref.fa
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
58 Loaded 1 sequences.
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
59 Loading features: reference/ref.gff
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
60 Parsing variants: snps.vcf
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
61 Converted 16 SNPs to TAB format.
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
62
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
63 ### /home/dfornika/miniconda3/envs/snippy-4.3.6/bin/snippy-vcf_extract_subs snps.filt.vcf > snps.subs.vcf
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
64
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
65
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
66 ### bcftools convert -Oz -o snps.vcf.gz snps.vcf
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
67
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
68
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
69 ### bcftools index -f snps.vcf.gz
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
70
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
71
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
72 ### bcftools consensus -f reference/ref.fa -o snps.consensus.fa snps.vcf.gz
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
73
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
74 Note: the --sample option not given, applying all records regardless of the genotype
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
75
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
76 ### bcftools convert -Oz -o snps.subs.vcf.gz snps.subs.vcf
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
77
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
78
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
79 ### bcftools index -f snps.subs.vcf.gz
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
80
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
81
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
82 ### bcftools consensus -f reference/ref.fa -o snps.consensus.subs.fa snps.subs.vcf.gz
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
83
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
84 Note: the --sample option not given, applying all records regardless of the genotype
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
85
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
86 ### rm -f snps.subs.vcf.gz snps.subs.vcf.gz.csi snps.subs.vcf.gz.tbi
e7e5da92d0ab planemo upload for repository https://github.com/tseemann/snippy commit d9e4e18619d5b9378f4fcbe518856fc1fa202d73-dirty
dfornika
parents:
diff changeset
87