Biostar Beta. Not for public use.
Question: Align to a reversed genome
1
Entering edit mode

Dear all,

to test the alignment procedure, I inserted a genome in reverse sense into an artificial set of fastq files and aligned against the reference genome in normal orientation. The resulting alignment does not show any mapping.

Is it normal that the aligner could not map the reads because they were in the reverse orientation? I thought that orientation could not affect the alignment procedure...

Thank you

ADD COMMENTlink 12 months ago marongiu.luigi • 380
Entering edit mode
0

I inserted a genome in reverse sense into an artificial set of fastq files and aligned against the reference genome in normal orientation.

Perhaps you could be a bit more specific here, by adding the commands you used?

ADD REPLYlink 12 months ago
WouterDeCoster
39k
Entering edit mode
0

Well, nothing much to add: I reversed the reference fasta file, then used ART to generate the fastq files, then used BWA to align the fatsq to the reference file but there is no alignment. Do you really need all the commands? This is more of a theoretical question...

ADD REPLYlink 12 months ago
marongiu.luigi
• 380
4
Entering edit mode

Reference genome orientation of course matters. What is interchangable would be the reference genome in its standard representation and its reverse complement but not its reverse. That would completely ignore that DNA has a 5'->3' polarity that cannot simply be changed at will.

ADD COMMENTlink 12 months ago ATpoint 17k
Entering edit mode
0

OK, though so. Just wanted to be sure. Thanks

ADD REPLYlink 12 months ago
marongiu.luigi
• 380

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0