How can I get polyA position and coverage from polyA(+)RNASeq data?
1
2
Entering edit mode
4 months ago
Ri ▴ 30

Hi everyone!

As polyA(+) RNAseq library produce reads with polyA tails, is there any way I can extract the position and coverage informations of reads that contain polyA from mapped sam/bam file?

Thank you!

NGS polyA RNASeq • 242 views
ADD COMMENT
1
Entering edit mode
4 months ago

Just grep the FASTQ for 'AAAAAAA$' and pull the read names. There is no way that the reads aligned without soft-clipping anyway.

ADD COMMENT

Login before adding your answer.

Traffic: 2550 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6