This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to use 2 paired end .sam files in featureCounts?

Hello,

I split paired ends data from RNA high throughput sequencing data into two .fastq and then .sam files (using fastq-dump, then hisat2) and want to put that data into featureCounts, but the only sample code I could find in their documentation references only one file. Does anyone know how I can do this or how I should merge my two .sam files?

Thank you so much!

rna-seq sequencing alignment
• 3,453 views • link
updated by devbt15 • devbt15 • written by lauragouldgoetz • lauragouldgoetz

You will have to explain more carefully what you are doing - preferably, showing the commands used. Are the two fastq files from paired-reads? If yes, they should have been mapped together.

If you "put" several files into featureCounts, it will output a table of counts, with one column per bam file. Is that what you want?

• 0 views • link
written by h.mon h.mon

The original .sra file fileName.sra was from paired reads.

I first did: $ fastq-dump --split-files fileName.sra

Then (on the two generated files: fileName_1.fastq and fileName_2.fastq): $ hisat2 -x mm10/genome -U fileName_1.fastq -S fileName_1.sam $ hisat2 -x mm10/genome -U fileName_2.fastq -S fileName_2.sam

Now I have the two .sam files and don't know how to proceed from here to put them into featureCounts. Thanks!

• 0 views • link
written by lauragouldgoetz • lauragouldgoetz

Yes, the hisat2 command should use both the fastq files in case of a paired end read sequencing. The resultant .bam file has to be put into featureCounts to get read counts for the features...

• 0 views • link
written by devbt15 • devbt15

1 answer

I just realized that I should have used the hisat2 command:

hisat2 -x mm10/genome -1 fileName_1.fastq -2 fileName_2.fastq -S fileName.sam

Then the resultant file goes on to featureCounts.

Thank you for your help!

• 0 views • link
updated by Istvan Albert Istvan Albert • written by lauragouldgoetz • lauragouldgoetz

Log in to answer this question.

More posts like this