Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Robin Erich Muench
metaSNV
Commits
2e0a25ae
Commit
2e0a25ae
authored
Jan 21, 2017
by
Robin Erich Muench
Browse files
compilation fix
parent
4ae32f43
Changes
1
Hide whitespace changes
Inline
Side-by-side
src/snpCaller/call_vC.cpp
View file @
2e0a25ae
...
...
@@ -53,7 +53,7 @@ static int print_usage()
fprintf
(
stderr
,
"Contact: Paul Costea <costea@embl.de>,
\n\t
Robin Muench <rmuench@embl.de>
\n\n
"
);
fprintf
(
stderr
,
"Usage: snpCall [options] <stdin.mpileup>
\n
"
);
fprintf
(
stderr
,
"Options:
\n
"
);
fprintf
(
stderr
,
" -b, TODO: list of input BAM filenames, one per line (for header)
\n
"
);
//
fprintf(stderr, " -b, TODO: list of input BAM filenames, one per line (for header) \n ");
fprintf
(
stderr
,
" -f, faidx indexed reference metagenome
\n
"
);
fprintf
(
stderr
,
" -g, gene annotation file [NULL].
\n
"
);
fprintf
(
stderr
,
" -i, individual SNPs output file [NULL].
\n\n
"
);
...
...
@@ -603,7 +603,7 @@ std::map<std::string, char> CodonMap(codons, codons + sizeof codons / sizeof cod
sToWrite
=
&
s
;
}
else
{
//May be a individual variant
for
(
int
i
=
1
;
i
<=
nrSamples
;
++
i
)
{
int
s
=
getSum
(
map
,
snp
,
i
);
int
s
=
getSum
(
map
,
check
,
i
);
if
(
s
>=
CALLING_THRESHOLD
)
{
//We observe the variant "CALLING_THRESHOLD" times in one sample.
writeThis
=
true
;
sToWrite
=
&
indiv
;
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment