I ran the four commands below.
python -m HTSeq.scripts.count accepted_hits.sam Mh.gtf > count.txt
samtools view accepted_hits.bam | python -m HTSeq.scripts.count - Mh.gtf > count.txt
samtools view accepted_hits.bam | python -m HTSeq.scripts.count -m intersection-strict - Mh.gtf > count.txt
samtools view accepted_hits.bam | python -m HTSeq.scripts.count -m intersection-nonempty - Mh.gtf > count.txt
However, all of the outputs show the same result:
no_feature 0
ambiguous 0
too_low_aQual 0
not_aligned 0
alignment_not_unique 409168
With this result, I don't think I can move to the next step, DESeq.
How could I make interpretation on this? Would you see what is wrong with my commands or result and what whould be fixed?
Thank you in advance.
python -m HTSeq.scripts.count accepted_hits.sam Mh.gtf > count.txt
samtools view accepted_hits.bam | python -m HTSeq.scripts.count - Mh.gtf > count.txt
samtools view accepted_hits.bam | python -m HTSeq.scripts.count -m intersection-strict - Mh.gtf > count.txt
samtools view accepted_hits.bam | python -m HTSeq.scripts.count -m intersection-nonempty - Mh.gtf > count.txt
However, all of the outputs show the same result:
no_feature 0
ambiguous 0
too_low_aQual 0
not_aligned 0
alignment_not_unique 409168
With this result, I don't think I can move to the next step, DESeq.
How could I make interpretation on this? Would you see what is wrong with my commands or result and what whould be fixed?
Thank you in advance.
Comment