[EMBOSS] Searching for repeats in fuzznuc

Phillip San Miguel pmiguel at purdue.edu
Thu Jun 21 16:46:39 UTC 2007


Peter Rice wrote:
> [...]
>   
>
> % fuzznuc -pattern @tg.pat
>
> Hope that helps,
>
> Peter
>   
Hi Peter,
    Along these lines, I've noticed that dreg wants the pattern file 
name to be all caps. Eg.

# cat test.pat
 >CCHC
TG[CT].{6}TG[CT].{12}CA[CT](?:...){4,5}TG[CT]
# dreg -rformat gff -stdout -auto -sequence HEX0238L06merge.fasta @test.pat
Error: Unable to open regular expression file 'TEST.PAT'
Error: Bad regular expression pattern:
   '@TEST.PAT'
Died: dreg terminated: Bad value for '-pattern' with -auto defined
# ln -s test.pat TEST.PAT
# dreg -rformat gff -stdout -auto -sequence HEX0238L06merge.fasta @test.pat
##gff-version 2.0
##date 2007-06-21
##Type DNA HEX0238L06Merged_Contigs
HEX0238L06Merged_Contigs        dreg    misc_feature    25288   25329   
0.000   +       .       Sequence "HEX0238L06Merged_Contigs.1" ; note 
"*pat CCHC"

Phillip




More information about the EMBOSS mailing list