Entering edit mode
8.4 years ago
kcamnairb
▴
40
Hi, I have a paired end library that was split among two sequencing lanes. The left and right reads are in separate files so I have a total of 4 files for each library. How would I specify this library in abyss, or do I need to concatenate the reads from the two lanes? For example, would this be proper:
abyss-pe k=64 name=ecoli lib='pe200 pe500'
pe200='pe200_lane1_1.fa pe200_lane1_2.fa pe200_lane2_1.fa pe200_lane2_2.fa' pe500='pe500_lane1.fa pe500_lane1_2.fa pe500_lane2.fa pe500_lane2_2.fa'
I definitely think the easiest thing to do is to concatenate your reads - putting all your Left reads in one file, and all your Right reads in the other, for each library that you have.
However I have to say, I don't think I've ever heard of 1 library being split and sequenced between two lanes. I'd be careful with that data. If there are any disparities between the performance of the two lanes, your data might be a bit wacky. On the other hand, maybe it's not such a big deal. Just curious, is there any particular reason you sequenced the library in this fashion?
Thanks, I'll try concatenating the reads. I actually have two libraries, so I think pooling libraries across multiple lanes is supposed to reduce any lane bias.
When we have N samples to run across M lanes, it's normal to pool them across all to normalize any lane performance disparities across all samples. I haven't seen a lane completely fail , but the sequencing people claim it used to happen and this strategy is to make sure no samples are completely lost by a failed lane.