New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | REL606 | = 242725 | 15 (0.380) | 4 (0.120) | 3/132 | 5.1 | 16.6% | intergenic (+343/‑88) | ECB_00212/yhhI | hypothetical protein/predicted transposase |
? | REL606 | = 719050 | 27 (0.780) | coding (66/762 nt) | ybfD | hypothetical protein | |||||
Rejected: Coverage evenness skew score above cutoff. | |||||||||||
Rejected: Frequency below/above cutoff threshold. |
ATGGTAGTGAGCGCAGGGGCAGGGAAAGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTCATTACCTGAAGCTGTTTTTTATTG > REL606/242650‑242768 | aTGGTAGTGAGCGCAGGGGCAGGGAAAGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGttt > 1:225634/1‑76 (MQ=255) tGGTAGTGAGCGCAGGGGCAGGGAAAGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTa > 1:890971/1‑76 (MQ=255) aGCGCAGGGGCAGGGAAAGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTga < 1:1099476/76‑1 (MQ=255) gcaggggcagggAAAGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGat < 1:822207/76‑1 (MQ=255) caggggcagggAAAGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGata < 2:713170/76‑1 (MQ=255) aggggcagggAAAGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGatat < 1:74411/76‑1 (MQ=255) gggcagggAAAGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTc < 2:962916/76‑1 (MQ=255) agggAAAGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTCa < 1:7464/76‑1 (MQ=255) ggAAAGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTCAtt < 1:757825/76‑1 (MQ=255) aGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTCATTACCt < 2:430596/76‑1 (MQ=255) aGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCCCATTACCt < 1:1209394/76‑1 (MQ=255) gATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTCATTACCTg < 1:689946/76‑1 (MQ=255) cACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTCATTACCTGAAg < 1:613900/76‑1 (MQ=255) cACGTAACGCTGAGTGTGTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTCATTACCTGAAg < 1:613926/76‑1 (MQ=14) cGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTCATTACCTGAAGCt > 1:1117866/1‑76 (MQ=255) ttttttGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTCATTACCTGAAGCTGTTTTTTATTg > 2:747992/1‑76 (MQ=31) | ATGGTAGTGAGCGCAGGGGCAGGGAAAGATCACGTAACGCTGATTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTCATTACCTGAAGCTGTTTTTTATTG > REL606/242650‑242768 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ 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. |