New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | USA300TCH1516_ALE | = 800928 | 33 (1.260) | 4 (0.160) | 4/266 | 3.0 | 14.7% | intergenic (+62/‑289) | hisC_1/USA300HOU_RS03915 | Histidinol‑phosphate aminotransferase/Putative 5'(3')‑deoxyribonucleotidase |
? | USA300TCH1516_ALE | 801948 = | 15 (0.600) | intergenic (+189/+435) | USA300HOU_RS03915/USA300HOU_RS03920 | Putative 5'(3')‑deoxyribonucleotidase/Putative lipid kinase | |||||
Rejected: Frequency below/above cutoff threshold. |
GTGACTGAGATGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGT > USA300TCH1516_ALE/801818‑802077 | gtgACTGAGATGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCAGGGCCCCAACATAGAGAAATTGGAAAGCGAATTTCTACAGACCAAGCAAGTTGGAGTAGGCCACAACATAAAGAAATa > 2:100815/1‑141 (MQ=37) acaattggttttaaagaacaaaatgataaaatgttagaagttttatcaaactttaaatacgaatagtaagtggggagtgggacagaaatgatactttcgcaaaatttatttcgtcgtcccaccccaacATAAAGAAATACt < 2:164382‑M2/13‑1 (MQ=255) ttttATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCGAATTGCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGaaa < 1:254874/141‑1 (MQ=255) aTGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAg > 2:146819/1‑141 (MQ=255) aGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCt < 1:62084/141‑1 (MQ=255) tGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTCGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCACTCata > 2:138203/1‑141 (MQ=255) gggCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAg > 1:80333/1‑141 (MQ=255) gCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATACATTGATCTTTTTAAGTa > 1:151415/1‑141 (MQ=255) agtggggagtgggacagaaatgatactttcgcaaaatttatttcgtcgtcccaccccaacATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACAt > 1:200136‑M2/61‑141 (MQ=255) aCATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACAtt < 1:157273/141‑1 (MQ=255) gggagtgggacagaaatgatactttcgcaaaatttatttcgtcgtcccaccccaacATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGc < 1:59501‑M2/85‑1 (MQ=255) ggAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAAt < 1:58058/141‑1 (MQ=255) aTTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATCCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGaa < 1:291399/141‑1 (MQ=255) aCAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTa > 2:288569/1‑141 (MQ=255) ttcgtcgtgccacccaaacCTAGATAAGTACTTTTTGTGTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAatcaatca < 2:33622‑M2/122‑1 (MQ=255) gTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCAt > 2:108891/1‑141 (MQ=255) ggTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGCTCAAGAATCGCTACATCATCTATCATTAGTc > 1:68609/1‑141 (MQ=255) gTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCg < 2:179659/141‑1 (MQ=255) tGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGt < 1:275067/141‑1 (MQ=255) | GTGACTGAGATGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGT > USA300TCH1516_ALE/801818‑802077 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
Reads not counted as support for junction |
read_name Not counted due to insufficient overlap past the breakpoint. |
read_name Not counted due to not crossing MOB target site duplication. |