New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NZ_CP009273 | 2837257 = | 27 (1.240) | 3 (0.150) | 3/152 | 1.6 | 12.5% | coding (313/768 nt) | hycG | formate hydrogenlyase subunit HycG |
? | NZ_CP009273 | 2837317 = | 18 (0.920) | coding (253/768 nt) | hycG | formate hydrogenlyase subunit HycG | |||||
Rejected: Frequency below/above cutoff threshold. |
GCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAG > NZ_CP009273/2837229‑2837395 | gCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTa < 1:388799/90‑1 (MQ=255) cACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTaa < 2:366687/90‑1 (MQ=255) cgtcatgcggatattttactgtttacgggcgcggtcacccgtgcaatgcgatcccctgcgctgcgtgcgtggcagtccgcgccggTAAAc < 1:395859‑M2/5‑1 (MQ=255) ggAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTaa < 1:274495/90‑1 (MQ=255) tACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAatat > 2:300526/1‑90 (MQ=255) aaaTTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCg > 1:176042/1‑90 (MQ=255) gtttaCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCCGTAAACAGTAAAATATCCGCATGAcg < 1:395859‑M1/85‑1 (MQ=255) gggTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGAcgc < 1:259355/90‑1 (MQ=255) gggTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGAcgc < 1:98979/90‑1 (MQ=255) ttaccGGCGCGGTCTGCCAAGCACTCCGCGACGGGGGTCGCGTTGCACGGGTGACCGCGCCGGTAAACAGGAAAATGTCCGCATGGcgcg > 1:48447/4‑90 (MQ=255) gcgatcccctgcgctgcgtgcgtggcagtccgcgccggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGaa > 2:94389‑M2/39‑90 (MQ=255) cagcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTACTGCCGaa > 1:337037/1‑90 (MQ=255) cagcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGaa > 2:601913/1‑90 (MQ=255) cgatcccctgcgctgcgtgcgtggcagtccgcgccggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAt < 1:94389‑M2/53‑1 (MQ=255) agcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAt < 1:430559/90‑1 (MQ=255) gcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATc < 1:358822/90‑1 (MQ=255) agGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTc < 1:337064/90‑1 (MQ=255) agGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTc < 2:405664/90‑1 (MQ=255) cGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggc > 2:264725/1‑90 (MQ=255) gTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGaa > 1:212122/1‑90 (MQ=255) gACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAg > 1:216587/1‑90 (MQ=255) gACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAg > 2:74799/1‑90 (MQ=255) | GCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAG > NZ_CP009273/2837229‑2837395 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ 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. |