--legalizing=1 .
--maskShortPortion=0 .
--noLeadingN=1 .
--removeShortSeq=200 .
Problems with the fasta file; no symbol > found !
I'm sure that the entries in my genome file start with ">", like this:
>1
TTGCACTGGTCTGACTGAATAGTTTCATTTCACCATGAATGACTTTGTGGAGcaacaaaa
gaaatattttgCTGAACATGAacgtttttttttgttattccgAGCTGATTGTCCCGTTGA
TATGCATGCTCTTGTTCAACTCATTTTCactgttattattattattataaaaaaactttt
Do you know why the > symbols in my file are not recognized?
Dear Shengfeng,
I got expected results when I ran scripts A3 and B1 with the genome bbv18wm. However, I got the same error when I ran A3 and B1 with my own genome:
I'm sure that the entries in my genome file start with ">", like this: >1 TTGCACTGGTCTGACTGAATAGTTTCATTTCACCATGAATGACTTTGTGGAGcaacaaaa gaaatattttgCTGAACATGAacgtttttttttgttattccgAGCTGATTGTCCCGTTGA TATGCATGCTCTTGTTCAACTCATTTTCactgttattattattattataaaaaaactttt
Do you know why the > symbols in my file are not recognized?
Best regards, Haoran