This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Guessing The File Format And Information

Hi everyone,

I have no information on this large sequence file(1.2Gb) nor its method of sequencing which i received as .txt. Is it advisable to proceed further with analysis. How do you go from here.

For starts, can anyone identify the format/background a little about the following few lines: How you did it or can do it is firstly great to know.

Thanks.

>S111:32:A03SG:1:1:12484:2206
NTCATTTAGATATCTGGCTTACACGCATAGCATTCTCAAGACATGGTACGCTTAATAAGTGATATNAATNTTTCAATTAAANCNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
>S111:32:A03SG:1:1:13579:2206
NTGTCGAGTCGATGTCTGATGGACGAGCATGAGTGACGCTCTGTTGTTTGTGCAGATTTGGCTGCNTTTNCGTTTTGNGTTNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
>S111:32:A03SG:1:1:12338:2206
NGGTCCAGAAATGTACTTTTGAGGGTTGTTTCAAGACCAACATGACTTTCAAACATATTCTGGAANATTNCTGGGTCATTCNGNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNCNNNNNNNNNNNNN
next-gen sequencing
• 3,335 views • link
updated by Biostar Biostar • written by Sequer • Sequer

3 answers

Just a guess, but this looks like an Illumina GAIIx single-ended run converted from fastq to fasta based on its 150 bp read length and 7-colon seperated identifier string: http://en.wikipedia.org/wiki/FASTQ_format#Illumina_sequence_identifiers

• 0 views • link
updated by Ram Ram • written by Casey Bergman Casey Bergman

The format is fasta (although technically, the lines with sequence should be wrapped; typically they contain 60 characters per line and no more than 80).

As raygozak said, the header line seems to contain some information in the ID. It looks like each sequence might be part of a longer sequence (since they all begin S111). However, since the content of header lines is completely arbitrary and determined by whoever created the file, it's impossible to say anything definitive.

I suggest, if possible, that you contact whoever generated the file for clarification.

• 1 views • link
written by Neilfws Neilfws

would the trailing consecutive NNNN's mean anything useful?

• 0 views • link
written by Sequer • Sequer

Yes, it means they're useless data. You want to pursue the originator of these files for the FASTQ files from which this has been generated. The FASTQ files will have the base qualities associated with each base, and will be much more amenable to analysis with NGS tools.

• 0 views • link
written by User 59 User 59

From what i can see this is a raw sequence file, with the identifier in one row and sequence in another, you can convert this to FASTA file by adding the > symbol on the identifier. The id line seems to contain annotation information about the sequence, separated by colon.

I'd recommend looking up the strings in the position of "A03SG", it looks like some sort of feature (gene, cds,mrna, etc.) id or name, and then compare the feature information with the last two numbers as they look as location information.

• 0 views • link
written by Raygozak Raygozak

the > was clipped by the biostar editor i'm not sure why

• 0 views • link
written by Sequer • Sequer

Because ">" at the start of a line is marked up in BioStar to be a block quote. I've edited your sequence so as lines are indented by 4 spaces.

• 0 views • link
written by Neilfws Neilfws

now that i look at it, it looks like the result of a whole genome shotgun sequencing.

• 0 views • link
written by Raygozak Raygozak

Log in to answer this question.

More posts like this