Fantastic thank you this worked. However, I am trying to use humann3 for peptide identification. However, I think my fasta file is parsed incorrectly:
1
LLLARKLKLNLLLLLVLAAK
2
DEARRRMRRWWNARHNHPRRRLWVSRMPR
3
HDMCKWGEWWYQYLPTAYYDCMMAK
4
VKLLKKKKKGAKLLVSVCK
Do you have a required fasta format that is required? I would greatly appreciate knowing what this is!
Best,
Carla
#CPU threads: 1
Scoring parameters: (Matrix=BLOSUM62 Lambda=0.267 K=0.041 Penalties=11/1)
Temporary directory: /Users/curanga/Desktop/humann3/C1peps/C1novor_humann_temp/tmp4bngxi3p
Opening the database… [0s]
Percentage range of top alignment score to report hits: 1
Reference = /Users/curanga/opt/anaconda3/envs/metaphlan3.0/lib/python3.7/site-packages/humann/data/uniref_DEMO/uniref90_demo_prots_v201901.dmnd
Sequences = 44
Letters = 18842
Block size = 2000000000
Opening the input file… [0.001s]
Opening the output file… [0s]
Loading query sequences… [0s]
Error: Error reading input stream at line 2: Invalid character (L) in sequence