Mercurial > repos > simon-gladman > snippy

annotate test-data/test/snps.log @ 22:577b1b17a0ff draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added snippy-core tool
author simon-gladman
date Thu, 09 Jun 2016 09:44:50 -0400
parents c6ecf655fde1
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
1
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
2 ### samtools faidx reference/ref.fa
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
3
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
4
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
5 ### bwa index reference/ref.fa
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
6
2
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
7 [bwa_index] Pack FASTA... 0.04 sec
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
8 [bwa_index] Construct BWT for the packed sequence...
2
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
9 [bwa_index] 2.34 seconds elapse.
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
10 [bwa_index] Update BWT... 0.04 sec
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
11 [bwa_index] Pack forward-only FASTA... 0.03 sec
2
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
12 [bwa_index] Construct SA from BWT and Occ... 0.61 sec
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
13 [main] Version: 0.7.12-r1039
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
14 [main] CMD: bwa index reference/ref.fa
2
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
15 [main] Real time: 3.079 sec; CPU: 3.060 sec
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
16
2
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
17 ### (bwa mem -v 2 -M -R '@RG ID:snps SM:snps' -t 1 reference/ref.fa /Users/Simon/code/galaxy-wrappers/snippy/test-data/reads_1.fq /Users/Simon/code/galaxy-wrappers/snippy/test-data/reads_2.fq | samtools view -@ 1 -q 60 -F 3844 -S -b -u -T ref.fa - | samtools sort -@ 1 - snps)
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
18
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
19 [samfaipath] fail to read file ref.fa.
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
20 [M::mem_pestat] skip orientation FF as there are not enough pairs
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
21 [M::mem_pestat] analyzing insert size distribution for orientation FR...
2
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
22 [M::mem_pestat] (25, 50, 75) percentile: (429, 488, 536)
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
23 [M::mem_pestat] low and high boundaries for computing mean and std.dev: (215, 750)
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
24 [M::mem_pestat] mean and std.dev: (479.53, 81.93)
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
25 [M::mem_pestat] low and high boundaries for proper pairs: (108, 857)
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
26 [M::mem_pestat] skip orientation RF as there are not enough pairs
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
27 [M::mem_pestat] skip orientation RR as there are not enough pairs
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
28 [main] Version: 0.7.12-r1039
2
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
29 [main] CMD: bwa mem -v 2 -M -R @RG ID:snps SM:snps -t 1 reference/ref.fa /Users/Simon/code/galaxy-wrappers/snippy/test-data/reads_1.fq /Users/Simon/code/galaxy-wrappers/snippy/test-data/reads_2.fq
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
30 [main] Real time: 0.268 sec; CPU: 0.262 sec
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
31
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
32 ### samtools index snps.bam
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
33
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
34
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
35 ### samtools depth -q 20 snps.bam | bgzip > snps.depth.gz
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
36
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
37
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
38 ### tabix -s 1 -b 2 -e 2 snps.depth.gz
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
39
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
40
2
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
41 ### fasta_generate_regions.py reference/ref.fa.fai 1397543 > reference/ref.txt
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
42
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
43
2
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
44 ### freebayes-parallel reference/ref.txt 1 -p 1 -q 20 -m 60 -F 0.9 -C 10 -V -f reference/ref.fa snps.bam > snps.raw.vcf
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
45
2
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
46 WARNING(freebayes): Could not find any mapped reads in target region gi|15829254|ref|NC_002695.1|:1397543..2795086
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
47
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
48 ### vcffilter -f 'DP > 9' -f 'QUAL > 10' snps.raw.vcf > snps.vcf
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
49
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
50
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
51 ### bgzip -c snps.vcf > snps.vcf.gz
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
52
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
53
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
54 ### tabix -p vcf snps.vcf.gz
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
55
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
56
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
57 ### /usr/local/Cellar/snippy/2.6/bin/snippy-vcf_to_tab --gff reference/ref.gff --ref reference/ref.fa --vcf snps.vcf > snps.tab
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
58
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
59 Loading reference: reference/ref.fa
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
60 Loaded 1 sequences.
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
61 Loading features: reference/ref.gff
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
62 Parsing variants: snps.vcf
2
c6ecf655fde1 planemo upload commit 09dd886ffd8ab6348cb002c9383d7bf342c465a9-dirty
simon-gladman
parents: 0
diff changeset
63 Converted 0 SNPs to TAB format.
0
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
64
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
65 ### vcf-consensus snps.vcf.gz < reference/ref.fa > snps.consensus.fa
da9e55d788b8 planemo upload commit 4c7d4bfdccecf4acaf3ba612166182e41e5f398f-dirty
simon-gladman
parents:
diff changeset
66