This is a test version of Biostars. For the public version, visit https://www.biostars.org.
I am not able to let R read my file

I proceed to import this alignment file to save it into a variable in R

x <- read.alignment("humanSeqAligned.fasta", "fasta")

and I get this error:

Error in read.alignment("humanSeqAligned.fasta",  : 
  File humanSeqAligned.fasta is not readable
Execution halted

Can anyone tell me what I am doing wrong?

r fasta
• 4,616 views • link
updated by Arup Ghosh Arup Ghosh • written by trejomarco6 • trejomarco6

It looks like you might have two different quotation mark types in your command (“ versus "). Try manually typing the command into your terminal instead of copy-pasting it from somewhere if you did so.

• 0 views • link
written by rpolicastro rpolicastro

Python code above is not writing an aligned file. Why are you trying to read it as an alignment?

• 0 views • link
written by GenoMax GenoMax

Sorry i did an alignment on my terminal already.

• 0 views • link
written by trejomarco6 • trejomarco6

Then please remove the python code in original question. It is not related to issue you are experiencing.

• 0 views • link
written by GenoMax GenoMax

Okay fixed the question and fix what rpolicastro suggested to do. Now I have updated the question with what R is telling me.

• 0 views • link
written by trejomarco6 • trejomarco6

Is the file in the same working directory as R? You can check the R working directory with getwd(), and list the files with list.files().

• 0 views • link
written by rpolicastro rpolicastro

Are you able to view the file with less/cat? Is it in the current working directory?

• 0 views • link
written by GenoMax GenoMax

Yes I able to run less -S humanSeqAlign.Fasta in the working directory and im able to see all of the ids and base pairs in this file. Which makes me confused on why is not being read through R.

• 0 views • link
written by trejomarco6 • trejomarco6

Is your file in one of the following alignment formats?

Read a file in mase, clustal, phylip, fasta or msf format

Can you show us head -5 yourfile?

• 0 views • link
written by GenoMax GenoMax

Is in a fasta format. Is a long file but when I run head -5.

>NC_001643.1
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------GTTTATGTAGCTTACCCCCTCAAAGCAATACACTGAAAATGTTTCGACGGGTTTACATCACCCCATAAACAAACAGGTTTGGTCCTAGCCTTTCTATTAGCTCTTAGTAAGATTACACATGCAAGCATCCCCG-CCCCGTGAG-TCACCCTCTAAATCGCCATGATCAAAAGGAACAAGTATCAAGCACGCAGCAATGCAGCTCAAAACGCTTAGCCTAGCCACACCCCCACGGGAGACAGCAGTGATAAACCTTTAGCAATAAACGAAAGTTTAACTAAGCCATACTAACCTCAGGGTTGGTCAATTTCGTGCTAGCCACCGCGGTCATACGATTAACCCAAGTCAATAGAAACCGGCGTAAAGAGTGTTTTAGATCACC-----CCCCCATAAAGCTAAAATTCACCTGAGTTGTAAAAAACTCCAGCTGATACAAAATAAACTACGAAAGTGGCTTTAACACATCTGAATACACAATAGCTAAGACCCAAACTGGGATTAGATACCCCACTATGCTTAGCCCTAAACTTCAACAGTTAAATTAACAAAACTGCTCGCCAGAACACTACGAGCCACAGCTTAAAACTCAAAGGACCTGGCGGTGCTTCATATCCCTCTAGAGGAGCCTGTTCTGTAATCGATAAACCCCGATCAACCTCACCGCCTCTTGCTCAGCCTATATACCGCCATCTTCAGCAAACCCTGATGAAGGTTACAAAGTAAGCACAAGTACCCACGTAAAGACGTTAGGTCAAGGTGTAGCCTATGAGGTGGCAAGAAATGGGCTACATTTTCTACCCCAGAAAATTACGATAACCCTTATGAAACCTAAGGGTCAAAGGTGGATTTAGCAGTAAACTAAGAGTAGAGTGCTTAGTTGAACAGGGCCCTGAAGCGCGTACACACCGCCCGTCACCCTCCTCAAGTATACTTCAAAGGATACTTAACTTAAACCCCCTACGTATTTATATAGAGGAGATAAGTCGTAACATGGTAAGTGTACTGGAAAGTGCACTTGGACGAACCAGAGTGTAGCTTAACATAAAGCACCCAACTTACACTTAGGAGATTTCAACTCAACTTGACCACTCTGAGCCAAACCTAGCCCCAAACCCCCTCCACCCTACTACCAAACAACCTTAACCAAACCATTTACCCAAATAAAGTATAGGCGATAGAAATTGTAAACCGGCGCAATAGACATAGTACCGCAAGGGAAAGATGAAAAATTATACCCAAGCATAATACAGCAAGGACTAACCCCTGTACCTTTTGCATAATGAATTAACTAGAAATAACTTTGCAAAGAGAACCAAAGCTAAGACCCCCGAAACCAGACGAGCTACCTAAGAACAGCTAAAAGAGCACACCCGTCTATGTAGCAAAATAGTGGGAAGATTTATAGGTAGAGGCGACAAACCTACCGAGCCTGGTGATAGCTGGTTGTCCAAGATAGAATCTTAGTTCAACTTTAAATTTACCTACAGAACCCTCTAAATCCCCTTGTAAACTTAACTGTTAGTCCAAAGAGGAACAGCTCTTTAGACACTAGGAAAAAACCTTGTAAAGAGAGTAAAAAATTTAACACCCATAGTAGGCCTAAAAGCAGCCACCAATTAAGAAAGCGTTCAAGCTCAACACCCACAACCT-TAAAGATCCCAAACATACAACCGAACTCCTTACACCCAATTGGACCAATCTATTACCCCATAGAAGAACTAATGTTAGTATAAGTAACATGAAAACATTCTCCTCCGCATAAGCCTACATCAGACCAAAATATTAAACTGACAATTAACAGCCTAATATCTACAATCAACCAACAAGCCATTATTACCCCCGCTGTTAACCCAACACAGGCATGCCCACAAGGAAAGGTTAAAAAAAGTAAAAGGAACTCGGCAAATCTTACCCCGCCTGTTTACCAAAAACATCACCTCTAGCATTACCAGTATTAGAGGCACCGCCTGCCCGGTGACATATGTTTAACGGCCGCGGTACCCTAACCGTGCAAAGGTAGCATAATCACTTGTTCCTTAAATAGGGACTTGTATGAATGGCTCCACGAGGGTTTAGCTGTCTCTTACTTTCAACCAGTGAAATTGACCTACCCGTGAAGAGGCGGGCATAACATAACAAGACGAGAAGACCCTATGGAGCTTTAATTCATTAATGCAAACAATACTTAACAAACCTACAGGTCCTAAACTATTAAACCTGCATTAAAAATTTCGGTTGGGGCGACCTCGGAGCACAACCCAACCTCCGAGCAATACATGCTAAGACCTCACCAGTCAAAGCGAATTACTACATCCAATTGATCCAATGACTTGACCAACGGAACAAGTTACCCTAGGGATAACAGCGCAATCCTATTCCAGAGTCCATATCAACAATAGGGTTTACGACCTCGATGTTGGATCAGGACATCCCGATGGTGCAGCCGCTATTAAAGGTTCGTTTGTTCAACGATTAAAGTCCTACGTGATCTGAGTTCAGACCGGAGTAATCCAGGTCGGTTTCTATCTGT-TCTAAATTTCTCCCTGTACGAAAGGACAAGAGAAATGAGGCCTACTTCACAAAGCGCCTTCCCCAATAAATGATATTATCTCAATTTAGCGCCATGCCAACACCCACTCAAGAACAGAGTTT-GTTAAGATGGCAGAGCCCGGTAATTGCATAAAACTTAAAACTTTACAATCAGAGGTTCAATTCCTCTTCTTGACAACACACCCATGACCAACCTCCTACTCCTCATTGTACCCATCCTAATCGCAATAGCATTCCTAATGCTAACCGAACGAAAAATTCTAGGCTACATACAACTACGCAAAGGTCCCAACATTGTAGGTCCTTACGGGCTATTACAGCCCTTCGCTGACGCCATAAAACTCTTCACTAAAGAACCCTTAAAACCCTCCACTTCAACCATTACCCTCTACATCACCGCCCCAACCCTAGCCCTCACCATTGCCCTCTTACTATGAACCCCCCTCCCCATACCCAACCCCCTAGTCAATCTTAACTTAGGCCTCCTATTTATTCTAGCCACCTCCAGCCTAGCCGTTTACTCAATCCTCTGATCAGGGTGAGCATCAAACTCGAACTACGCCTTAATCGGTGCACTACGAGCAGTAGCCCAAACAATCTC
• 0 views • link
updated by GenoMax GenoMax • written by trejomarco6 • trejomarco6

aligned fasta format is not the same as plain fasta format. Just making sure. example of an aligned fasta file below

>APE_PIG 
MRVLWVALVVTLLAGCRTEDEPGPPPE.VHVWWEESKWQGSQP.WEQALGRFWDYLRWVQ
SLSDQVQEELLSTKVTQELTELIEESMKEVKAYREELEAQLGPVTQETQARLSKELQAAQ
ARVGADMEDVRNRLVLYRSEVHNMLGQTTEELRSRLASHLRNVRKRLVRDTEDLQKRLAV
YQAGLREGAERSVSALRERLGPLVEQGRLRAATLSTRAGQPLRERAEAWGQKLRGRLEEM
GSRTRDRLDEMRDELEEVRTKVEEQGSQLRLQAEAFQARLKGWFEPLVEDMRRQWAGLVE
RMQSAVSISSSTSA..PSDNQ
>APE_MOUSE 
MKALWAVLLVTLLTGCLA........EGEPEVTDQLEWQSNQP.WEQALNRFWDYLRWVQ
TLSDQVQEELQSSQVTQELTALMEDTMTEVKAYKKELEEQLGPVAEETRARLGKEVQAAQ
ARLGADMEDLRNRLGQYRNEVHTMLGQSTEEIRARLSTHLRKMRKRLMRDAEDLQKRLAV
YKAGAREGAERGVSAIRERLGPLVEQGRQRTANLGAGAAQPLRDRAQAFGDRIRGRLEEV
GNQARDRLEEVREHMEEVRSKMEEQTQQIRLQAEIFQARLKGWFEPIVEDMHRQWANLME
KIQASVATNPII.TPVAQENQ
• 0 views • link
written by GenoMax GenoMax

Is a long file so I can't post the full output.

• 0 views • link
written by trejomarco6 • trejomarco6

As long as all sequences are same length and leading and trailing positions are indicated by tildes (~) and gap positions are indicated by periods (.) you should have the correct format.

• 0 views • link
written by GenoMax GenoMax

Yeah I checked to see if they were all the same length and if there's any suspicious characters, but I couldn't find any.

• 0 views • link
written by trejomarco6 • trejomarco6

I accidentally spammed this post which caused it to be removed and for you to be suspended. I approved the post and reinstated you. Sorry!

• 0 views • link
written by rpolicastro rpolicastro

As mentioned by rpolicastro , make sure the file is present in the R work directory.

• 0 views • link
written by Arup Ghosh Arup Ghosh

0 answers

No answers yet.

Log in to answer this question.

More posts like this