This should have been posted as a separate question. Yes you can use the alignment files generated by either for DE. You will need to use featureCounts to get counts. STAR is able to generate gene counts with an option provided at run time, if you want to avoid having to run featureCounts.
RNA-seq read alignment
I tried mapping PE RNA-seq reads to the reference genome (chickpea) using hisat2 and tophat tools with default settings. I got 0.00% alignment rate. what are the main parameters that can be tweaked to get a good alignment rate.
PS: I am still a novice in this area.
2 answers
Could you post the commend line? Maybe some parameter is wrong.
• 0 views
•
link
written
by
Sparrow_kop
hisat2 -p 2 -q -c -5 15 -3 5 --phred33 -t -x ht2_Ca_v.1_index -1 P1_C1_CBBYRANXX_CGGCTATG-TATAGCCT_L003_R1_paired_1.fastq,P1_C1_CBBYRANXX_CGGCTATG-TATAGCCT_L004_R1_paired_1.fastq,P1_C1_CBBYRANXX_CGGCTATG-TATAGCCT_L005_R1_paired_1.fastq,P1_C1_CBBYRANXX_CGGCTATG-TATAGCCT_L006_R1_paired_1.fastq -2 P1_C1_CBBYRANXX_CGGCTATG-TATAGCCT_L003_R2_paired_2.fastq,P1_C1_CBBYRANXX_CGGCTATG-TATAGCCT_L004_R2_paired_2.fastq,P1_C1_CBBYRANXX_CGGCTATG-TATAGCCT_L005_R2_paired_2.fastq,P1_C1_CBBYRANXX_CGGCTATG-TATAGCCT_L006_R2_paired_2.fastq -S P1_C1_CBBYRANXX_CGGCTATG-TATAGCCT.sam --new-summary
You can remove "- c" option, and try again. From hisat2 help: "-c <m1>, <m2>, <r> are sequences themselves, not files", but in your command , -1 and -2 refer to sequence file, not the sequence themselves....
"- c" can used in the below demo:
hisat2 -c -x index -1 ATCTCGTCGTCGCT -2 ATGATGATCGTGTC -S out.sam
• 0 views
•
link
written
by
Sparrow_kop
Log in to answer this question.
More posts like this
-
Low mapping rate
written by yoshifumimiya •
-
miRNA-seq alignment bowtie2
written by imaparna27 •
-
HiSat2 vs BowTie2 for mapping miRNA-seq reads
written by Arindam Ghosh
-
Calculate Mapping Rate for an alignment to each sequence of a multi-fasta reference
written by kspata
-
RNA-seq biological replicates
written by amy16 •
-
RNA-seq quality assessment
written by amy16 •
-
Optimization of RNA-Seq data mapping with tophat2
written by Denis
-
Improving alignment rates for RNA-seq
written by pm2012
-
The comparison between HISAT2 and Tophat2
written by Peng Huang
-
DNA-seq PE150 BWA alignment
written by zhuozhu132 •
If the alignment was truly 0% then take a sample of reads (convert them to fasta) and blast at NCBI to verify that the data is really from chickpea and not something else.