How can I get polyA position and coverage from polyA(+)RNASeq data?
1 answer
Just grep the FASTQ for 'AAAAAAA$' and pull the read names. There is no way that the reads aligned without soft-clipping anyway.
• 0 views
•
link
written
by
benformatics
Log in to answer this question.
More posts like this
-
Matching read IDs to mapped transcripts/genes
written by joshcylee •
-
Extract stderr from mapped BAM files of bowtie2
written by Farzaneh •
-
Extarct specific bases from BAM/SAM files
written by Ankit
-
Extract primer-containing reads from RNA-Seq
written by Ruixuan •
-
Get genes and coverage/depth per gene from sam/bam
written by prishly •
-
tool for a pattern of repeats
written by 2nelly
-
Extract the mapped contig sequences from SAM/BAM file
written by BDK_compbio
-
Extracting reads with homopolymer errors
written by aparnakrish89 •
-
Nucletotide distribution, at each position, in a .sam/.bam file ?
written by remi.maglione •
-
Simulation Of Rnaseq Read Mapping With Reads Having Different Polya Tail Lengths.
written by Gaima Kushwaha •
