I am trying to count single-end reads.
Based on amel_OGSv3.2.gene.gff
(Group1.1 amel_OGSv3.1 gene 507599 515039 1 - . ID=GB42155),
I written the following command:
python -m HTSeq.scripts.count -s no -t gene -i ID CK_2.sam /media/wenfu/LaCie/my_rnaseq_dat/amel_OGSv3.2.gene.gff > CK_2.gff
Is there any one who can check if the command was written correctly, because I don't know whether “-s no -t gene -i ID” are right.
Thanks a lot!!
Richard
Based on amel_OGSv3.2.gene.gff
(Group1.1 amel_OGSv3.1 gene 507599 515039 1 - . ID=GB42155),
I written the following command:
python -m HTSeq.scripts.count -s no -t gene -i ID CK_2.sam /media/wenfu/LaCie/my_rnaseq_dat/amel_OGSv3.2.gene.gff > CK_2.gff
Is there any one who can check if the command was written correctly, because I don't know whether “-s no -t gene -i ID” are right.
Thanks a lot!!
Richard
Comment