OK, though so. Just wanted to be sure. Thanks
• 0 views
•
link
written
by
marongiu.luigi
Align to a reversed genome
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
• 1,464 views
•
link
written
by
marongiu.luigi
1 answer
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.
Log in to answer this question.
More posts like this
-
No hits for influenza virus in an RNA seq. data using FAST Q screen (from babraham institute)
written by Mohamed Samir •
-
Artificial reads - remove multiple mapped reads against reference genome, and only keep reads that …
written by Jalil Sharif •
-
Question about RNA-Seq data alignment
written by mohammedtoufiq91
-
Bowtie2 reference genome
written by Ada •
-
how to confirm that mapped read are accurate?
written by marongiu.luigi
-
BWA mapping to multiple reference genome
written by dsklsdoiwld •
-
How to interpret a read hit with IGV?
written by marongiu.luigi
-
CIGAR operations performed according to strand or to read?
written by Matteo Schiavinato
-
SAM Flags for illumina paired-end reads mapped on transcriptome
written by Matteo Schiavinato
-
WES reads mapping: RNASeq or DNASeq aligner?
written by user230613
Perhaps you could be a bit more specific here, by adding the commands you used?
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...