Odd reads from RAD seq
2
0
Entering edit mode
6.8 years ago

Hello all, I´m trying to align RAD seq reads to de novo in Bowtie2, but I get the next problem:

Warning: Could not open read file "Spa.vir_CAR147_1.fastq" for reading; skipping...
Warning: Could not open read file "Spa.vir_CAR153_1.fastq" for reading; skipping...
Warning: Could not open read file "Spa.vir_CAR154_1.fastq" for reading; skipping...
Error: No input read files were valid
(ERR): bowtie2-align exited with value 1

I saw the reads with the command head, and the reads looks like this

@SN747:516:C8M5TACXX:1:2210:2138:56104 1:N:0:1
TGCAGGGATCTTTCTCTCCCAGCAGACCCTCTGTTATTTTTACAAAATATACTCTCCAGAGACTTTGATGTGTAGGATTTCAGTTATTGTG
+
HHHJJJJJJJIJJJJJJJJJJIJIJJJJJJJJJJIJJJJJJJJIJJJJJJJJJJJJGHHHHFFFFDFEEEFEDDCDCCDDDDDFEDEDDEC
@SN747:516:C8M5TACXX:1:2210:2166:56214 1:N:0:1
TGCAGGTTCAGCGTGGAGTCTGGGTACAGCTCCTCTCCCTGTGAGCCGGCCCACCTCACTGCTGCTGTACCGGACTGACTCCTCAGACTCC
+
HHHJJJGJJJJJJHIIJJHHIIJIDHIJJIJHIIJJJJJIJGIIJJJJEHFDDDDDDDDDDCDDDACDDDDDDDDDDDDDDDDDDDDDDDD

I´m very new in NGS. Please somebody could explain me, what is exactly the issue with the reads? Thank you so much

alignment • 934 views
ADD COMMENT
0
Entering edit mode
6.8 years ago
GenoMax 141k

The reads are fine/normal.

It appears that bowtie2 is unable to read your sequence files. This may simply be a problem of providing correct file paths (full or relative). Since you did not include the original command line used in this case try adding correct paths to that line.

ADD COMMENT
0
Entering edit mode
6.8 years ago

Thank you! I will try it, again

ADD COMMENT

Login before adding your answer.

Traffic: 1990 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6