Been a while since we had a code golf challenge (shortest script wins, but expect to receive upvotes for cleverness)!
Simplifying slightly, bisulfite conversion means that Cs get converted to Ts, but some bases may be protected by methylgroups and thus remain unconverted. (This happens on the opposite strand too, so some complimentary Gs could get converted to As)
So, given an input sequence read (assume a [ACTGN] string of no more than 150 bp), output all reads that might result from partial or complete bisulfite conversion of that sequence.
EDIT: To be clear, C>T happens on one read, G>A happens on the other. Reads will not exist with both conversions! See this screenshot:
Small example: Input:
AATGTGAA AATGCGAA AACGTGAA AACACAAA AACACGAA AACGCAAA AACGCGAA
(the original string could be in there as well, since all bases could be protected!)