-
Notifications
You must be signed in to change notification settings - Fork 0
/
README
33 lines (23 loc) · 718 Bytes
/
README
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
Usage:
falco
input:
- bam file
- manifest file
- ref + fai
- locifilt
- manifest
output:
- *.vcf file
- *.qc.ann.txt -> catalog of each assayed bp
- *.qc2.ann.txt -> catalog of transitions
- *.qc.targets.txt -> amplicon centered information, eg depth, start/end positions.
example:
OUTBASE="HCT116-TEST"
## variant calling
falco --bam data/HCT-116/HCT116.bam --output $OUTBASE --ref genome.fa
file1=$OUTBASE".qc.ann.qual.txt"
file2=$OUTBASE".qc2.ann.txt"
file3=$OUTBASE".qc.targets.txt"
file4=$OUTBASE".res.filtered.tsv"
## variant filtering and reporting
falco-filter-report --vcf snpeff.vcf --output $OUTBASE --qc_ann_qual_txt $file1 --qc2_ann_txt $file2 --qc_targets_txt $file3 --res_filtered_tsv $file4