| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913_3_vch_pgi | = 4171803 | 22 (1.000) | 4 (0.190) +TTGATGT |
4/276 | 2.5 | 16.7% | intergenic (+47/‑254) | rrfB/murB | 5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding |
| ? | NC_000913_3_vch_pgi | = 4213155 | 20 (0.910) | noncoding (116/120 nt) | rrfE | 5S ribosomal RNA of rrnE operon | |||||
| Rejected: Frequency below/above cutoff threshold. | |||||||||||
GCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAA > NC_000913_3_vch_pgi/4213013‑4213290 | gCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGACGGTAGAGTGGGGTCTCCCCATGCGAGCGAAGGGAACTACCAGccatca < 2:223455‑M2/149‑7 (MQ=255) gctgTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAg < 1:119359/146‑1 (MQ=17) tCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaaaa > 1:310256/1‑149 (MQ=35) tCTGATAAAACAGAATTTGCCTGGCGGCAGTAAAGTTTTGAGCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGCTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaaaa > 1:384615/1‑149 (MQ=11) ggaTAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaa > 2:256867‑M2/2‑131 (MQ=255) aCAGGATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGCCTCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTc < 2:343191/149‑1 (MQ=255) ggactttGCCTGGCGGCAGCAGTGCGGTGGTCCCACCTGACCCCATACCGAATTCAGAAGTGAAACGCCGTATCGCCTATGGTAGTGTGTGGTCTCCCCATGCGAGTGTAGGGAACTGCCAGTCATGAAATTAGAAAACCCCCGGTCCa < 2:242626/145‑1 (MQ=14) gcgcGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggctcagtctttcgactgagcctttcgttt > 1:324995‑M2/1‑101 (MQ=255) cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTACCGCCGATGGTAGTGAGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAAAAAACCCAGGTCACAAAGGCCGGGGTTTTTTGCATATCAatcaat > 2:185107/1‑145 (MQ=17) cACCTGACCCCATGCCGAACTCAGAAGTGAAACGACGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTGGAAAGACCCCGGTCCATAAGGCCGGGGTTTTTTGCGTATCAATTTTTTg > 2:308419/1‑149 (MQ=18) ccTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCa < 1:198775/149‑1 (MQ=255) ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCat < 1:230182/128‑1 (MQ=255) ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCat > 2:230182/1‑128 (MQ=255) ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATAAATTATTTGAATGATGaa > 1:261255/1‑149 (MQ=255) agcagAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGGGTGTGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAggg < 2:110425/145‑1 (MQ=255) tGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGgcatcagataaaacgaaaggcccagtctttcgactgagcctttcgt < 1:143174‑M2/108‑47 (MQ=255) tGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGgcatcagataaaacgaaaggcccagtctttcgactgagcctttcgt > 2:143174‑M2/1‑62 (MQ=255) aaCGCCGTAGGGTGCGTGGTAGTCGGGGGTCTCCCCAGGCGAGAGTAGGGAACTGCCAGGAATCAAATTAGAAAATCCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTg < 2:382433/149‑1 (MQ=37) tGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTAAGTTCTGTATATCAAAATtcgtc > 1:386952/1‑149 (MQ=255) ggggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCaaa > 1:47640/1‑149 (MQ=255) ggggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCaaa < 1:342339/149‑1 (MQ=255) cATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTACAAAAGGCCGGGGTTTTTTGCATATCAATTATTTGCCTGATAAAGGGACTCTCATGTGAGTTCTTTATATCCAAATGCATCGCAGTCAAATTACtttggg > 1:43240/1‑144 (MQ=37) gagTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGCATCTCATGTCAGTTCTGTATATCCAACTCCGACCTATTCAAATTACTGGTGCCGATcgt > 2:33190/1‑147 (MQ=255) gTAGGGAACTGCGAGGCATCGAATTAGAAAACCCCCGGTCCATAAGGCCGGTGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGa < 2:285811/149‑1 (MQ=255) tAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATACGCCAGGGGTTTTTTGCATATCAGTTATTTGCATGATGAAGGGAATCTCAGGTCAGGTCTGTATATCCAAATTCGTAGTATTGTAATTCCTATTAGCGATCTTgta > 2:57063/1‑147 (MQ=255) | GCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAA > NC_000913_3_vch_pgi/4213013‑4213290 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ 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. |