Read alignment anomaly

0

Entering edit mode

3.5 years ago

daewowo ▴ 80

I have aligned metagenomic reads to a genome, and find that at the 5' ends of the reads often there is no alignment to the genome - this is indicated by the multicoloured regions in the plot blow. Solid blue and solid pink are perfectly aligned sections of reads (fwd and rev).

Any ideas why parts of reads would not align in such a consistent manner?

aligned reads plotted in igv

igv reads rna-seq • 904 views

ADD COMMENT • link 3.5 years ago by daewowo ▴ 80

0

Entering edit mode

Are these soft-trimmed by the aligner? Did you do a quick quality check? (fastqc). Metagenomic sample and one reference genome alignment? Could also be some sample preparation bias? Just ideas

ADD REPLY • link 3.5 years ago by Gregor Rot ▴ 540

1

Entering edit mode

Could be that adapters not clipped, need to check workflow

ADD REPLY • link 3.5 years ago by daewowo ▴ 80

0

Entering edit mode

Same result after trimming with fastp. These are metagenomic sample and one ref genome. fastqc didnt turn up any obvious anomalies to me. Suspect there is some odd sample prep. issue

ADD REPLY • link 3.5 years ago by daewowo ▴ 80

Login before adding your answer.