[Seqan-dev] reading fasta with non-DNA5 characters


Hi,

I have a couple of genome seqeunces that contain characters other than ACTGN (i.e. Y, M,...)...

Is there a way to read those sequences in as well and automatically convert those non conforming letters to N?

Thanks,

Bernd

PS:

I am using:

RecordReader<String<char, MMap<> >, DoublePass<Mapped> > refReader(seqMMapString);
        int read2out = read2(seqIds, faSeqs, refReader, Fasta());

for reading in the data and get an INVALID_FORMAT error...