@@ -49,25 +49,25 @@ one can calculate:
4949
5050 b. Make sure that the opam packages are up to date:
5151
52- $ opam update
52+ $ opam update
5353
5454 c. Make sure that you're on the relevant compiler:
5555
56- $ opam switch 4.05.0
57- $ eval `opam config env`
56+ $ opam switch 4.05.0
57+ $ eval `opam config env`
5858
5959 d. Get source:
6060
61- $ git clone https://github.com/hammerlab/prohlatype.git prohlatype
62- $ cd prohlatype
61+ $ git clone https://github.com/hammerlab/prohlatype.git prohlatype
62+ $ cd prohlatype
6363
6464 e. Install the dependent packages:
6565
66- $ make setup
66+ $ make setup
6767
6868 f. Build the programs (afterwards they'll be in ` _build/default/src/apps ` ):
6969
70- $ make
70+ $ make
7171
7272### Make sure that you have [ IMGT/HLA] ( https://github.com/ANHIG/IMGTHLA ) available:
7373
@@ -81,26 +81,26 @@ one can calculate:
8181 we normally do ** not** have sequence information will still align (in the
8282 next filtering step), albeit poorly:
8383
84- $ align2fasta path-to-imgthla/alignments -o imputed_hla_class_I.fasta
84+ $ align2fasta path-to-imgthla/alignments -o imputed_hla_class_I.fasta
8585
8686 This step needs to be performed only once, per each IMGT version.
8787 Run ` $align2fasta --help ` for further information.
8888
8989 2 . Filter your data against the reference, by first aligning. Ex:
9090
91- $ bwa mem imputed_hla_class_I.fasta ${SAMPLE}.fastq | \
92- samtools view -F 4 -bT imputed_hla_class_I.fasta -o ${SAMPLE}.bam
91+ $ bwa mem imputed_hla_class_I.fasta ${SAMPLE}.fastq | \
92+ samtools view -F 4 -bT imputed_hla_class_I.fasta -o ${SAMPLE}.bam
9393
9494 While fundamentally, the algorithms here are * alignment* based. They're
9595 too slow to run for all sequences. Sequences that do not originate from
9696 the HLA-region would just act as background noice.
9797
9898 3 . and then convert aligned reads back to FASTQ:
9999
100- $ samtools fastq ${SAMPLE}.bam > ${SAMPLE}_filtered.fastq
100+ $ samtools fastq ${SAMPLE}.bam > ${SAMPLE}_filtered.fastq
101101
102- 4 . Infer types:
102+ 4 . Infer types (see ` $ multi_par --help ` for further details):
103103
104- $ multi_par path-to-imgthla/aignments ${SAMPLE}_filtered.fast -o ${SAMPLE}_output.tsv
104+ $ multi_par path-to-imgthla/aignments ${SAMPLE}_filtered.fastq -o ${SAMPLE}_output.tsv
105105
106- See ` $ multi_par --help ` for further detail and optimizations.
106+
0 commit comments