How can I possibly extract the reads from a BAM file and put them into a FASTQ file for simulation (maq simutrain, then maq simulate)?
Should I just extract col. 1, 10 and 11 from a BAM file and put them in a text file along with a '+'? That is the output by read simulators.
But does it not happen that DNA sequences and base quality sequences are reversed and/or transformed depending on the direction and strand of reads relative to reference genomes? I would have to fix that too in that case.
Should I just extract col. 1, 10 and 11 from a BAM file and put them in a text file along with a '+'? That is the output by read simulators.
But does it not happen that DNA sequences and base quality sequences are reversed and/or transformed depending on the direction and strand of reads relative to reference genomes? I would have to fix that too in that case.
Comment