1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
#include <seqan/bam_io.h>
using namespace seqan2;
int main()
{
CharString bamFileName = getAbsolutePath("demos/tutorial/sam_and_bam_io/example.sam");
BamFileIn bamFileIn(toCString(bamFileName));
BamHeader header;
readHeader(header, bamFileIn);
typedef FormattedFileContext<BamFileIn, void>::Type TBamContext;
TBamContext const & bamContext = context(bamFileIn);
for (unsigned i = 0; i < length(contigNames(bamContext)); ++i)
std::cout << contigNames(bamContext)[i] << '\t'
<< contigLengths(bamContext)[i] << '\n';
return 0;
}
|
