Importing lowercase alignment file

[Maxy (Max Boeck)]

I had to use an online mafft program program because of memory issues. I am now trying to import the alignment file, but for some reason all the nucleotides were converted to lowercase. How do I import lowercase nucleotides? I could write a janky script to do it, but I am trying as much as possible to do this using qiime tools. My command was as follows:

qiime tools import
–input-path EXP_type1_rep-seqs-or-99/_out.18072802055553becwu7zMjlmSU9hBaKFHdlsflarge.pir
–output-path type1_rep-seqs-or-99_aligned-sequences.qza
–type ‘FeatureData[AlignedSequence]’

[Nicholas_Bokulich]

@Maxy

How do I import lowercase nucleotides?

lowercase nucleotides can cause issues with some tools downstream, so I'd recommend just converting to uppercase. I don't think a bash one-liner would be too janky, provided it does not cause issues with your seq IDs.

Good luck!

[Maxy (Max Boeck)]

Unfortunately it does interfere with naming, but if there isn’t an easy solution I’ll write some kind of masticated perl script.

[colinbrislawn (Colin J Brislawn)]

Try this command in bash:

awk '/^>/ {print($0)}; /^[^>]/ {print(toupper($0))}' file.fasta > file_upper.fasta

:qiime2: :

Colin