Multiple entries in a some columns on FeatureCounts output
1
0
Entering edit mode
13 days ago
myoui3122010 ▴ 20

Performed alignment with Bowtie2 using

bowtie2 -p 4 -x ~/shared/NGS/P4/genome/bowtie2_ind/      bowtie2_ind -1 /mnt/shared/NGS/P4/control/SRR30712601_1.fastq.gz -2 /mnt/shared/NGS/P4/control/SRR30712601_2.fastq.gz -S aligned_SRR30712601.sam

Sam to Bam using

samtools view -bS results/aligned_SRR30712601.sam -o results/aligned_SRR30712601.sorted.bam

samtools sort results/aligned_SRR30712601.sorted.bam -o results/aligned_aligned_SRR30712601.sorted.bam

samtools index results/aligned_aligned_SRR30712601.sorted.bam 

Then featureCounts was run using

 featureCounts -a ~/shared/NGS/P4/genome/Homo_sapiens.GRCh38.108.gtf \ 

-o aligned_SRR30712601_counts.txt \

-T 4 \

-p \

aligned_aligned_SRR30712601.sorted.bam 

This .txt file has multiple entries in the same rows in chr, start, end columns.

Is this normal? Can I just calculate TPM, FPKM or RPKM ?
enter image description here

FeatureCounts • 180 views
ADD COMMENT
0
Entering edit mode
13 days ago
biofalconch ★ 1.3k

Yes, this is normal as the output contains the chromosome, start and end positions for all exons. But you can just take the first and seventh (last) column which contain the Gene ID and its respective counts :) You can also take the 6th column (Length) to calculate FPKM as well.

ADD COMMENT

Login before adding your answer.

Traffic: 2257 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