bwa mem to feature Counts
1
0
Entering edit mode
ajajoo • 0
@ajajoo-21737
Last seen 3.1 years ago
United States

bwa mem -t $Ncpu -k 14 -T 1 -L 3,3 -O 6,6 -E 3,3 ../humangenomes/hg38.fa $filename.R1.fastq.gz $filename.R2.fastq.gz | samtools view -Sbh > $filename.bam

Can I use above bam file for feature counts, how will it handle chimeric alignment and what is the quality for linear match with the above parameters. Is there any away I can filter bam file from above command to pass it to feature count. This is what I use right now.

featureCounts -p -T 5 -a hg38GENCODE.gtf -o $filename.txt $filename.bam &

I am running data in bulk, that requires bwa mem with above parameters so do not want to use resources for e.g. STAR alignment etc for passing it to featureCount.

I am very new to this field, so please direct me to good literature for all these things. I am still understanding bwa mem algo and how these parameters will affect the results.

bwa mem featureCounts • 1.5k views
ADD COMMENT
2
Entering edit mode
swbarnes2 ★ 1.4k
@swbarnes2-14086
Last seen 2 days ago
San Diego

I'm pretty sure bwa-mem is not splice aware, so it's not suitable for RNASeq. You really should use STAR.

ADD COMMENT
2
Entering edit mode

Or use Bioconductor software Rsubread::align and Rsubread::featureCounts. This is the Bioconductor forum after all!

ADD REPLY

Login before adding your answer.

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