Closed YCCHEN23 closed 2 years ago
It seems that there is some thing wrong in the sam2tsv steps, can you check the filesize of final/extract.sort.bam.tsv.gz
Hoping these messages help, thanks!
Can you send me your final/extract.sort.bam file to gaoyubang@qq.com for checking?
I convert the bam file using sam2tsv -r {1} {0}/extract.sort.bam|gzip -c >{0}/extract.sort.bam.tsv.gz Then i test the extract.sort.bam.tsv.gz file with follow python code:
import sys,os,re import gzip fl=sys.argv[1] for i in gzip.open(fl,"r"): i=i.decode("utf-8").rstrip()
if i.startswith("#"):
continue
ele=i.rstrip().split()
try:
if ele[3]=="." or ele[6]==".":
continue
except:
print(i.rstrip())
sys.exit()
But there is no wrong with the test code. Thus can you give me your extract_raw_and_feature_fast output for checking?
Okay, I sent the extract_raw_and_feature_fast output to gaoyubang@qq.com few minutes ago.
Thanks!
It seemed that I used bed12 format for "predict_sites", my fault, lol!
Maybe it could add some foolproof notification for bed format?
Thanks a lot for all your testing.
Hi,
I'm trying to predict the m6A site of Arabidopsis DRS data (published data, vir-1) with nanom6A, but I received the following message while performing "predict_sites".
What files should I check for? or any other information should I provide?
Thanks!