EFieldDelimiter delim = eTab,
83 boolparse_local_ids =
false);
107 void SetScores(
intscore,
doublebit_score,
doubleevalue);
114 voidSetEndpoints(
intq_start,
intq_end,
ints_start,
ints_end);
121 voidSetCounts(
intnum_ident,
intlength,
intgaps,
intgap_opens,
122 intpositives =0,
intquery_frame = 1,
123 intsubject_frame = 1);
126 voidSetBTOP(
stringbtop_string);
132 intSetFields(
constobjects::CSeq_align& sal,
133objects::CScope& scope,
136x_SetFieldDelimiter(eCustom,customDelim);
139 virtual void Print(
void);
150 voidPrintHeader(
const string& program,
151 constobjects::CBioseq& bioseq,
154 unsigned intiteration =
156 constobjects::CSeq_align_set* align_set=0,
163 voidPrintNumProcessed(
intnum_queries);
166list<string> GetAllFieldNames(
void);
181 voidSetNoFetch(
boolnofetch);
195 voidPrintFieldNames(
bool is_csv=
false);
198 voidPrintFieldSpecs(
void);
210 voidx_AddDefaultFieldsToShow(
void);
213 voidx_SetFieldsToShow(
const string&
format);
215 voidx_ResetFields(
void);
218 voidx_SetFieldDelimiter(EFieldDelimiter delim,
stringcustomDelim =
"");
223 voidx_PrintQuerySeqId(
void)
const;
225 voidx_PrintQueryGi(
void);
227 voidx_PrintQueryAccession(
void);
229 voidx_PrintQueryAccessionVersion(
void);
231 voidx_PrintQueryAndDbNames(
const string& program,
232 constobjects::CBioseq& bioseq,
235 unsigned intiteration,
238 voidx_PrintSubjectSeqId(
void);
240 voidx_PrintSubjectAllSeqIds(
void);
242 voidx_PrintSubjectGi(
void);
244 voidx_PrintSubjectAllGis(
void);
246 voidx_PrintSubjectAccession(
void);
248 voidx_PrintSubjectAccessionVersion(
void);
250 voidx_PrintSubjectAllAccessions(
void);
252 voidx_PrintQuerySeq(
void);
254 voidx_PrintSubjectSeq(
void);
256 voidx_PrintQueryStart(
void);
258 voidx_PrintQueryEnd(
void);
260 voidx_PrintSubjectStart(
void);
262 voidx_PrintSubjectEnd(
void);
264 voidx_PrintEvalue(
void);
266 voidx_PrintBitScore(
void);
268 voidx_PrintScore(
void);
270 voidx_PrintAlignmentLength(
void);
272 voidx_PrintPercentIdentical(
void);
274 voidx_PrintNumIdentical(
void);
276 voidx_PrintMismatches(
void);
278 voidx_PrintNumPositives(
void);
280 voidx_PrintGapOpenings(
void);
282 voidx_PrintGaps(
void);
284 voidx_PrintPercentPositives();
286 voidx_PrintFrames();
287 voidx_PrintQueryFrame();
288 voidx_PrintSubjectFrame();
291 voidx_PrintQueryLength();
293 voidx_PrintSubjectLength();
295 voidx_PrintSubjectTaxIds();
296 voidx_PrintSubjectSciNames();
297 voidx_PrintSubjectCommonNames();
298 voidx_PrintSubjectBlastNames();
299 voidx_PrintSubjectSuperKingdoms();
300 voidx_PrintSubjectTaxId();
301 voidx_PrintSubjectSciName();
302 voidx_PrintSubjectCommonName();
303 voidx_PrintSubjectBlastName();
304 voidx_PrintSubjectSuperKingdom();
305 voidx_PrintSubjectTitle();
306 voidx_PrintSubjectAllTitles();
307 voidx_PrintSubjectStrand();
308 voidx_PrintSeqalignCoverage();
309 voidx_PrintSubjectCoverage();
310 voidx_PrintUniqSubjectCoverage();
314 voidx_SetQueryCovSubject(
constobjects::CSeq_align & align);
315 voidx_SetQueryCovUniqSubject(
constobjects::CSeq_align & align);
316 voidx_SetQueryCovSeqalign(
const CSeq_align& align,
intquery_len);
382 return(m_CustomDelim ==
","|| m_FieldDelimiter ==
",");
446 doubleperc_positives =
532 SIgDomain(
const string&
n,
ints,
inte,
intss,
intse):
533name(
n), start(s), end(e),
534s_start(ss), s_end(se), length(0),
535num_match(0), num_mismatch(0), num_gap(0) {};
549 void Set(
const string id,
ints,
inte) {
550 if(
id.substr(0,4) ==
"lcl|") {
551sid =
id.substr(4,
id.
size());
583 const string& program,
584 constobjects::CBioseq& bioseq,
586 const string& domain_sys,
588 unsigned intiteration =
590 constobjects::CSeq_align_set* align_set=0,
595 intSetMasterFields(
constobjects::CSeq_align& align,
596objects::CScope& scope,
597 const string& chain_type,
598 const string& master_chain_type_to_show,
602 int SetFields(
constobjects::CSeq_align& align,
603objects::CScope& scope,
604 const string& chain_type,
605 const string& master_chain_type_to_show,
609 virtual void Print(
void);
614 voidSetAirrFormatData(
CScope& scope,
620 voidPrintAirrRearrangement(
CScope& scope,
622 const string& program_version,
625 const string& domain_sys,
627 unsigned intiteration,
631 boolprint_airr_format_header,
644m_IsMinusStrand =
minus;
654 ints_start=-1,
ints_end=-1) {
655 if(start <0 || end <= start)
return;
657x_ComputeIgDomain(*domain);
658m_IgDomains.push_back(domain);
663 void SetVGene(
const string&
id,
ints,
inte) {
664m_VGene.Set(
id, s,e);
668 void SetDGene(
const string&
id,
ints,
inte) {
669m_DGene.Set(
id, s,e);
673 void SetJGene(
const string&
id,
ints,
inte) {
674m_JGene.Set(
id, s,e);
678 void SetCGene(
const string&
id,
ints,
inte) {
679m_CGene.Set(
id, s,e);
695 string& master_chain_to_show,
698 string& productive)
const{
703master_chain_to_show = m_MasterChainTypeToShow;
704cdr3_nuc = m_Cdr3Seq;
705cdr3_aa = m_Cdr3SeqTrans;
706productive = m_OtherInfo[4];
717 voidx_ResetIgFields();
718 voidx_PrintIgGenes(
boolisHtml=
false,
const string& header=
"# ")
const;
719 voidx_ComputeIgDomain(SIgDomain &domain);
720 voidx_PrintIgDomain(
constSIgDomain &domain)
const;
721 voidx_PrintIgDomainHtml(
constSIgDomain &domain)
const;
722 voidx_PrintPartialQuery(
intstart,
intend,
boolisHtml=
false)
const;
742 static const intnum_otherinfo = 5;
743 stringm_OtherInfo[num_otherinfo];
User-defined methods of the data storage class.
void SetScores(objects::CSeq_align &seq_align, objects::CScope &scope, const string &matrix_name="BLOSUM62")
Class containing information needed for tabular formatting of BLAST results.
void SetQueryGeneticCode(int q_gc)
ESeqIdType
In what form should the sequence identifiers be shown?
@ eAccession
Show only best accession.
@ eAccVersion
Show only best accession.version.
list< CRef< objects::CSeq_id > > m_QueryId
List of query ids for this HSP.
string m_BitScore
Bit score of this HSP, in appropriate format.
void SetParseSubjectDefline(bool val)
Should subject deflien be parsed for id or not?
void x_PrintSubjectSeq(void)
Print aligned part of subject sequence.
void x_PrintQueryLength()
Print the query sequence length.
void x_PrintNumIdentical(void)
Print number of identical matches.
map< string, ETabularField > m_FieldMap
Map of field enum values to field names.
void x_PrintSubjectEnd(void)
Print subject end.
pair< string, int > m_QueryCovUniqSubject
void SetQueryId(const objects::CBioseq_Handle &bh)
Set query id from a Bioseq handle.
void x_PrintGapOpenings(void)
Print number of gap openings.
set< string > m_SubjectSuperKingdoms
int m_AlignLength
Alignment length of this HSP.
set< string > m_SubjectBlastNames
bool m_ParseSubjectDefline
Parse subejct defline?
string m_SubjectBlastName
EFieldDelimiter
What delimiter to use between fields in each row of the tabular output.
void SetCustomDelim(string customDelim)
int m_SubjectStart
Starting offset in subject.
void x_PrintPercentIdentical(void)
Print percent of identical matches.
void x_PrintScore(void)
Print raw score.
void x_PrintQueryStart(void)
Print query start.
void x_PrintBitScore(void)
Print bit score.
TSeqPos m_SubjectLength
Length of subject sequence.
void x_PrintGaps(void)
Print total number of gaps.
void x_PrintNumPositives(void)
Print number of positive matches.
int m_SubjectEnd
Ending offset in subject.
string m_SubjectSuperKingdom
string m_SubjectSeq
Aligned part of the subject sequence.
void SetQueryRange(TSeqRange &q_range)
Set query range.
void SetSubjectId(list< CRef< objects::CSeq_id > > &id)
Set subject id from a objects::CSeq_id.
CNcbiOstream & m_Ostream
Stream to write output to.
CRef< CBlast_def_line_set > m_SubjectDefline
void SetParseLocalIds(bool val)
Should local IDs be parsed or not?
void x_PrintSubjectStart(void)
Print subject start.
int m_QueryEnd
Ending offset in query.
int m_SubjectFrame
subject frame
int m_NumGapOpens
Number of gap openings in this HSP.
vector< string > m_SubjectCommonNames
void SetSubjectId(const objects::CBioseq_Handle &bh)
Set subject id from a Bioseq handle.
pair< string, int > m_QueryCovSubject
void x_PrintSubjectLength()
Print the subject sequence length.
list< ETabularField > m_FieldsToShow
Which fields to show?
virtual void Print(void)
Print one line of tabular output.
int SetFields(const objects::CSeq_align &sal, objects::CScope &scope, CNcbiMatrix< int > *matrix=0)
Set all member fields, given a Seq-align.
bool m_ParseLocalIds
Should the query deflines be parsed for local IDs?
void x_PrintEvalue(void)
Print e-value.
void x_PrintMismatches(void)
Print number of mismatches.
void SetNoFetch(bool nofetch)
Avoid fetching sequence (if possible) If the sequence is needed (e.g., will be formatted,...
bool m_NoFetch
program as a string
void SetDbGeneticCode(int db_gc)
void x_PrintPercentPositives()
Print percent positives.
vector< list< CRef< objects::CSeq_id > > > m_SubjectIds
All subject sequence ids for this HSP.
string m_FieldDelimiter
Delimiter character for fields to print.
bool GetNoFetch()
Avoid fetch of sequence if true returned.
const list< CRef< CSeq_id > > & GetQueryId() const
Get query seqid list.
set< TTaxId > m_SubjectTaxIds
Blast-traceback-operations.
bool x_IsFieldRequested(ETabularField field)
int m_NumPositives
Number of positives in this HSP.
bool x_IsCsv() const
Convenience method to determine if this object is dealing with CSV output formats.
int m_QueryStart
Starting offset in query.
void SetQueryId(list< CRef< objects::CSeq_id > > &id)
Set query id from a objects::CSeq_id.
void x_PrintQuerySeq(void)
Print aligned part of query sequence.
void x_PrintSubjectFrame()
void x_PrintFrames()
Print frames.
void x_PrintAlignmentLength(void)
Print alignment length.
int m_Score
Raw score of this HSP.
void PrintHeader(const string &program, const objects::CBioseq &bioseq, const string &dbname, const string &rid=kEmptyStr, unsigned int iteration=numeric_limits< unsigned int >::max(), const objects::CSeq_align_set *align_set=0, CConstRef< objects::CBioseq > subj_bioseq=CConstRef< objects::CBioseq >(), bool is_csv=false)
Print the tabular output header.
vector< string > m_SubjectSciNames
int m_QueryFrame
query frame
list< CRef< objects::CSeq_id > > m_SubjectId
string m_SubjectCommonName
void x_PrintQueryEnd(void)
Print query end.
int m_NumGaps
Total number of gaps in this HSP.
int m_NumIdent
Number of identities in this HSP.
string m_QuerySeq
Aligned part of the query sequence.
string m_Evalue
E-value of this HSP, in appropriate format.
TSeqPos m_QueryLength
Length of query sequence.
Class containing information needed for tabular formatting of BLAST results.
void SetMinusStrand(bool minus=true)
Set strand information.
vector< SIgDomain * > m_IgDomains
~CIgBlastTabularInfo()
Destructor.
map< string, string > m_AirrData
CIgBlastTabularInfo(CNcbiOstream &ostr, const string &format=kDfltArgTabularOutputFmt, EFieldDelimiter delim=eTab)
What delimiter to use between fields in each row of the tabular output.
CRef< CSeq_align > m_TopAlign_V
void SetSeqType(bool isNucl)
Set sequence type.
void SetJGene(const string &id, int s, int e)
Set gene info.
const vector< SIgDomain * > & GetIgDomains() const
Get Ig domain.
void SetDGene(const string &id, int s, int e)
Set gene info.
void AddIgDomain(const string &name, int start, int end, int s_start=-1, int s_end=-1)
Set domain info.
void SetFrame(const string &frame="N/A")
Set out-of-frame information.
void SetCGene(const string &id, int s, int e)
Set gene info.
CRef< CSeq_align > m_TopAlign_J
CRef< CSeq_align > m_TopAlign_D
CRef< CSeq_align > m_TopAlign_C
void SetVGene(const string &id, int s, int e)
Set gene info.
void GetIgInfo(string &v, string &d, string &j, string &c, string &master_chain_to_show, string &cdr3_nuc, string &cdr3_aa, string &productive) const
Getter.
string m_AirrCdr3SeqTrans
string m_MasterChainTypeToShow
void Print(const CCompactSAMApplication::AlignInfo &ai)
ETabularField
Enumeration for all fields that are supported in the tabular output.
const char * kDfltArgTabularOutputFmt
Default value for tabular and comma-separated value output formats.
unsigned int TSeqPos
Type for sequence locations and lengths.
SStrictId_Tax::TId TTaxId
Taxon id type.
#define END_NCBI_SCOPE
End previously defined NCBI scope.
#define END_SCOPE(ns)
End the previously defined scope.
#define BEGIN_NCBI_SCOPE
Define ncbi namespace.
#define BEGIN_SCOPE(ns)
Define a new scope.
IO_PREFIX::ostream CNcbiOstream
Portable alias for ostream.
static string DoubleToString(double value, int precision=-1, TNumToStringFlags flags=0)
Convert double to string.
#define NCBI_ALIGN_FORMAT_EXPORT
char * dbname(DBPROCESS *dbproc)
Get name of current database.
Declares CIgBlast, the C++ API for the IG-BLAST engine.
The blob sat and sat key Both must be positive integers</td > n< td > Non empty string The interpretation of the blob id depends on a processor Cassandra n processor expects the following format
const struct ncbi::grid::netcache::search::fields::SIZE size
NCBI C++ stream class wrappers for triggering between "new" and "old" C++ stream libraries.
Defines BLAST database access classes.
struct containing annotated domain information
SIgDomain(const string &n, int s, int e, int ss, int se)
struct containing annotated gene information
void Set(const string id, int s, int e)
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4