3,045
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to searchno edit summary
SamValidationErrors samInvalidErrors;
SamValidationErrors samInvalidErrors;
// Set returnStatus to success. It will be changed
// to the failure reason if any of the writes fail.
SamStatus::Status returnStatus = SamStatus::SUCCESS;
// Keep reading records until ReadRecord returns false.
while(samIn.ReadRecord(samHeader, samRecord))
{
// Successfully read a record from the file, so write it.
if(!samOut.WriteRecord(samHeader, samRecord))
{
// Failed to write a record.
fprintf(stderr, "%s\n", samOut.GetStatusMessage());
returnStatus = samOut.GetStatus();
}
}
if(samIn.GetStatus() != SamStatus::NO_MORE_RECS)
{
// Failed to read a record.
fprintf(stderr, "%s\n", samIn.GetStatusMessage());
// Set the return status to the reason why
// the read failed.
returnStatus = samIn.GetStatus();
}
std::cout << std::endl << "Number of records read = " <<
samIn.GetCurrentRecordCount() << std::endl;
std::cout << "Number of records written = " <<
samOut.GetCurrentRecordCount() << std::endl;
// Since the reads were successful, return the status based
// on the status of the writes. If any failed, return
// their failure status.
return(returnStatus);
}
}
</source>
</source>
