I actually used STAR tool to align it to viral genomes.
• 0 views
•
link
written
by
ZuelTech •
RNA seq strandedness (reverse-strand)
Hi. I have an RNA seq data that is reverse-stranded. The library type is ISR. I want to count reads mapping to the genomes. I am new to bioinformatics. Can you check if my code is correct?
Plus strand reads (Read2 mapped to forward strand) ?
PLUS_COUNT=$(samtools view -f 0x80 -F 0x10 "$BAM" "$SEG" | wc -l)
Minus strand reads (Read2 mapped to reverse strand) ?
MINUS_COUNT=$(samtools view -f 0x80 -f 0x10 "$BAM" "$SEG" | wc -l)
• 4,457 views
•
link
written
by
ZuelTech •
2 answers
It depends on the alignment method you used. If you aligned against the transcriptome, you can do it that way you described. In case you aligned against the genome, you need to take into account on which DNA-strand the transcript is on. For that you can use RSeQC's infer-experiment tool.
• 0 views
•
link
written
by
michael.ante
You can use IGV to look at your alignment using your reference and its gene-annotation. You'll see how your reads orient themselves in regard to the genes.
Check also GenoMax' comment.
• 0 views
•
link
written
by
michael.ante
My purpose is to count reads mapping to the genomes. Was my code above in doing so correct?
• 0 views
•
link
written
by
ZuelTech •
Log in to answer this question.
More posts like this
-
Viral positive and negative strand with paired sequencing and bowtie
written by Lluís R.
-
paired-end reads, idxstats, and samtools and grep
written by Bogdan
-
Reads mapped on wrong strand
written by qstefano •
-
Splitting BAM file from paired end into forward and reverse strands?
written by Susmita Mandal
-
Orientation of reads in .fastq files from RNA-Seq
written by newbio17
-
Counting reads aligned to forward or reverse strand
written by caggtaagtat
-
BAM file header after splitting paired end reads
written by ilyco •
-
first and second mate reports by samtools 0x40 and 0x60 FLAGs
written by R •
-
Tutorial: How To Separate Illumina Based Strand Specific Rna-Seq Alignments By Strand
written by Istvan Albert
-
Samtools View: Only Forward Or Reverse Strand
written by Gregor Rot
Use a standard tool like
featureCountsas noted here --> RNA-seq analysis