diff --git a/src/rock.cpp b/src/rock.cpp index 2461fb65d306432a98140a23868cd3d4a4759183..fbb53023a5ad46314d58000ba5320c61164534d3 100644 --- a/src/rock.cpp +++ b/src/rock.cpp @@ -128,6 +128,7 @@ void lowFilterCMS(FqBaseBackend* map_id_backend[],int nb_be,int k, CountMinSketc unsigned int j=it_offs->first; for (it_struct=it_offs->second.begin();it_struct!=it_offs->second.end();it_struct++) { // read dna string corresponding to fastq record + if (it_struct->fileid) { DnaSeqStr a_seqs[2]; init_DnaSeqStr(&a_seqs[0]); init_DnaSeqStr(&a_seqs[1]); @@ -137,6 +138,7 @@ void lowFilterCMS(FqBaseBackend* map_id_backend[],int nb_be,int k, CountMinSketc ret=pcms->isBeneathMinKappa(nbrKmerDecompo); if (ret) it_struct->fileid=0; nbrKmerDecompo.clear(); + } } } }