libsequence
1.9.5
valid_dna.cc
#include <
Sequence/Fasta.hpp
>
#include <
Sequence/SeqRegexes.hpp
>
#include <
Sequence/SeqAlphabets.hpp
>
#include <fstream>
#include <iostream>
#include <algorithm>
int
main(
int
argc,
char
**argv)
{
std::ifstream in(argv[1]);
Sequence::Fasta
seq;
while
(in >> seq)
{
std::cout <<
Sequence::validSeq
(seq.
begin
(),seq.
end
())
<<
'\t'
<<
Sequence::validSeq
(seq.
begin
(),seq.
end
(),
Sequence::full_dna_alphabet
)
<<
'\t'
<< (std::find_if(seq.
begin
(),seq.
end
(),
Sequence::ambiguousNucleotide
())
!= seq.
end
())
<<
'\n'
;
}
}
Generated on Thu Sep 13 2018 15:07:05 for libsequence by
1.8.13