Tweet
I know that the read order of xxx_F3.csfasta and xxx_F3_QV.qual are the same, but I can't find any information that describes how the read order of xxx_F3.csfasta and xxx_R3.csfasta are related. I'd appreciate any pointers to documentation that describe the order relationship.
-
-
In practice, they will be sorted based on read name, with the numbers indicating panel, x-position, y-position. If a read has no mate, then it will only be present in one of the files. See solid2fastq programs/scripts like the ones in BFAST or MAQ that use the above properties.Originally posted by smarkel View Post -
The files I'm using are from AB's site (http://solidsoftwaretools.com/gf/project/ecoli2x50/) and aren't sorted. I looked at Rosalind_20080729_2_Chris5_F3.csfasta.zip and Rosalind_20080729_2_Chris5_R3.csfasta.zip. My understanding is that the sorting happens when the reads are mapped. Maybe the posted files aren't representative.Comment
-
I wasn't too descriptive in my explanation.Originally posted by smarkel View Post
Take a look at these reads from real data:
I claim they are sorted based on read name. They have the form:Code:>1_6_55_F3 T01100000000201002010120012300200011000100.01101131 >1_6_64_F3 T01203010110102003000000101100111000010100.01100131 >1_6_69_F3 T01031320200103032011110221111112110020111.11111131 >1_6_97_F3
>%d_%d_%d_F3
where %d stands for some integer. It is sorted by the right-most integer, then middle integer, then left-most integer. The equivalent read (the mate) in the R3 file will be
>%d_%d_%d_R3
The "Rosalind" file follows the same pattern:
Code:>469_26_42_F3 T12113310031232112221003120021221223320222122212122 >469_26_379_F3 T31202223003310000130302323312223212011000010033200 >469_26_540_F3 T11012313031030123033113130100223110001231232303210 >469_26_560_F3
Comment
-
Thank you. Yes, you're right. Both files are individually sorted. I should have worded my original question differently. I expected xxx_F3.csfasta and xxx_R3.csfasta to have the same number of entries, with the nth entry in the F3 file being the mate of the nth entry in the R3 file. What does it mean for a paired read's mate to not exist?Comment
-
It means that given some quality threshold, a proper the mate at that location could not be identified. I use "unpaired' reads without any problem. After >10 slides those extra reads really add up!Originally posted by smarkel View PostComment
-
Thank you for the explanation.Comment
-
Genomics studies in neuroscience face a special challenge due to the brain’s complexity and scarcity of samples. Mapping changes in cell type and state using conventional next-generation sequencing methods remains challenging. Advances in technologies like single-cell sequencing, spatial transcriptomics, and long-read sequencing have opened the door to deeper studies of the brain and diseases like Alzheimer’s, amyotrophic lateral sclerosis (ALS), and schizophrenia.
...07-09-2026, 11:10 AM -
Cancer survival rates have significantly increased in the last few decades in the United States, reaching a combined 70% 5-year survival rate by 2021. Behind this number, there are years of research to find new therapies, drug targets, and early detection methods. But there is one core challenge that keeps slowing down these advances, and it’s about drug resistance.
There is no single reason why many patients don’t respond to treatment as expected. Cancer is...07-08-2026, 05:17 AM -
Love this - good data definitely starts from good input, and poor input can only give relatively poor data. I particularly like the mention of Nanodrop/absorbance based methods for quantification. It's such a toss up if you'll get an accurate reading or what amounts to a randomly generated number, and a lot of library/sequencing related issues can be traced back to poor quant.07-01-2026, 11:43 AM
Comment