This is a test version of Biostars. For the public version, visit https://www.biostars.org.
featurecounts -p Paried ends which count total to use as read1 &read 2 give different totals.

Hi

First post so please be kind, its my first year.

I have run featurecounts -p on my files and the output has a count total for read 1 and read 2. My simple question is which total count should i use? read 1 or read 2 since they are not the same.

I tried finding the answer on the forumn but no luck.

Thanks for your help

R

rna-seq alignment r gene sequencing
• 1,403 views • link
written by ray.liv.wong • ray.liv.wong

featurecounts -p counts fragments if the reads are paired, I have never seen it reporting separated counts for R1 and R2... What is the exact command you use ? And can you show the summary results ?

• 0 views • link
written by Carlo Yague Carlo Yague

To generate Bam file i used hisat2

hisat2 -x "./refs/hg38_hisat2_index/hg38_tran/genome_tran" -1 "${file1}" -2 "${file2}" -S "./PRJNA316673_sam/${acc}.sam"
samtools sort -l 102 -o "./PRJNA316673_bam/${acc}.bam" "./PRJNA316673_sam/${acc}.sam"

Then for quantification

ref="./refs/hg38_ensembl_gtf/ensembl_gtf/hg38.gtf.gz"
list=(./PRJNA316673_bam/*)

featureCounts -p -T 16 -t exon -g gene_id -a $ref --tmpDir "./temp" -o "./PRJNA316673_3_counts.txt" "${list[@]}"

output is like this 

Geneid          Chr                   Strand            Length  read1   read2 
ENSG00000223972 1;1;1;1;1;1;1;1;1     +;+;+;+;+;+;+;+;+ 1735    0         0 
ENSG00000227232 1;1;1;1;1;1;1;1;1;1;1   ;-;-;-;-;-;-;-;-;-;-;-  1351    42        49    
ENSG00000268903         1                                -                           755             50        33
• 0 views • link
updated by Carlo Yague Carlo Yague • written by ray.liv.wong • ray.liv.wong

what is the output of ls PRJNA316673_bam/*?

• 0 views • link
written by Carlo Yague Carlo Yague

Can you post the full command you used? Also, can you post the command used to generate your alignment (i.e bam) file?

• 0 views • link
written by Rogerio Ribeiro • Rogerio Ribeiro

0 answers

No answers yet.

Log in to answer this question.

More posts like this