diff --git a/src/read_utils.cpp b/src/read_utils.cpp index c13a6ccaf9130945aca58608e60205c83da1499d..0639f4195a19fac9bb5a83f3dbfe4906220ce697 100644 --- a/src/read_utils.cpp +++ b/src/read_utils.cpp @@ -77,12 +77,12 @@ void getDNASeqstr(FqBaseBackend* fq_files_be [], f_id2=fid_stored &mask_fid; unsigned long offset1=j*UINT_MAX+sr.read_a1; - getFqRecord(fq_files_be,f_id1,offset1,dna_seqs[0].fq_record_buf); try { + getFqRecord(fq_files_be,f_id1,offset1,dna_seqs[0].fq_record_buf); processFqRecord(p_dna_seqs); } catch (int e) { cout<<"j="<<j<<" sr content: a1="<<sr.read_a1<<" a2= "<<sr.read_a2<<endl; - cout<<"pb reading at offset1="<<offset1<<endl; + cout<<"pb reading/seeking at offset1="<<offset1<<endl; cout<<"got: "<<endl; cout<<dna_seqs[0].fq_record_buf<<endl; } @@ -93,7 +93,7 @@ void getDNASeqstr(FqBaseBackend* fq_files_be [], getFqRecord(fq_files_be,f_id2,offset2,dna_seqs[1].fq_record_buf); } catch (int e) { cout<<"j="<<j<<" sr content: a1="<<sr.read_a1<<" a2= "<<sr.read_a2<<endl; - cout<<"pb reading at offset2="<<offset2<<endl; + cout<<"pb reading/seeking at offset2="<<offset2<<endl; cout<<"got: "<<endl; cout<<dna_seqs[1].fq_record_buf<<endl; }