if ( !m_settings->HasInputBamFilename ) m_settings->InputFiles.push_back(Options::StandardIn());
// opens the BAM files without checking for indexes
-// BamMultiReader reader;
-// reader.Open(m_settings->InputFiles, false);
-
- BamReader reader;
- reader.Open(m_settings->InputFiles.at(0));
+ BamMultiReader reader;
+ reader.Open(m_settings->InputFiles, false, true);
// retrieve header & reference dictionary info
std::string mergedHeader = reader.GetHeaderText();
writer.Open(m_settings->OutputFilename, mergedHeader, references);
// store alignments to output file
-// BamAlignment bAlignment;
-// while (reader.GetNextAlignment(bAlignment)) {
-// writer.SaveAlignment(bAlignment);
-// }
-
BamAlignment bAlignment;
- while (reader.GetNextAlignment(bAlignment)) {
+ while (reader.GetNextAlignmentCore(bAlignment)) {
writer.SaveAlignment(bAlignment);
}
-
+
// clean & exit
reader.Close();
writer.Close();