From e990c5b852c62cd1237919ccc75cd3ea049b915f Mon Sep 17 00:00:00 2001 From: Derek Barnett Date: Wed, 24 Jul 2013 15:56:40 -0400 Subject: [PATCH] Added model description to stats file header in 'bamtools resolve' tool. (issue #76) --- src/toolkit/bamtools_resolve.cpp | 20 +++++++++++++++++++- 1 file changed, 19 insertions(+), 1 deletion(-) diff --git a/src/toolkit/bamtools_resolve.cpp b/src/toolkit/bamtools_resolve.cpp index cb42f5b..9e5fb84 100644 --- a/src/toolkit/bamtools_resolve.cpp +++ b/src/toolkit/bamtools_resolve.cpp @@ -2,7 +2,7 @@ // bamtools_resolve.cpp (c) 2011 // Marth Lab, Department of Biology, Boston College // --------------------------------------------------------------------------- -// Last modified: 14 October 2011 +// Last modified: 24 July 2013 (DB) // --------------------------------------------------------------------------- // Resolves paired-end reads (marking the IsProperPair flag as needed). // *************************************************************************** @@ -73,6 +73,20 @@ static const string OPTION_FORCEMARKREADGROUPS = "ForceMarkReadGroups"; static const string RG_FIELD_DESCRIPTION = "# "; +static const string MODEL_DESCRIPTION = + "# ------------- Model Types Description ---------------\n" + "#\n" + "# ID Position Orientation \n" + "# 1 mate1 < mate2 mate1:forward, mate2:forward \n" + "# 2 mate1 < mate2 mate1:forward, mate2:reverse \n" + "# 3 mate1 < mate2 mate1:reverse, mate2:forward \n" + "# 4 mate1 < mate2 mate1:reverse, mate2:reverse \n" + "# 5 mate2 < mate1 mate2:forward, mate1:forward \n" + "# 6 mate2 < mate1 mate2:forward, mate1:reverse \n" + "# 7 mate2 < mate1 mate2:reverse, mate1:forward \n" + "# 8 mate2 < mate1 mate2:reverse, mate1:reverse \n" + "# -----------------------------------------------------\n"; + // -------------------------------------------------------------------------- // unique readname file constants // -------------------------------------------------------------------------- @@ -731,9 +745,13 @@ void ResolveTool::StatsFileWriter::WriteHeader(void) { << BAMTOOLS_VERSION_BUILD; // # bamtools resolve (vX.Y.Z) + // # + // # MODEL DESCRIPTION - see above for actual text // \n m_stream << COMMENT_CHAR << " bamtools resolve (" << versionStream.str() << ")" << endl + << COMMENT_CHAR << endl + << MODEL_DESCRIPTION << endl; } -- 2.39.2