| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | 3643087 = | 22 (0.760) | 4 (0.150) | 4/236 | 3.2 | 14.9% | coding (46/753 nt) | rsmJ | 16S rRNA m(2)G1516 methyltransferase, SAM‑dependent |
| ? | NC_000913 | 3643107 = | 25 (0.930) | coding (26/753 nt) | rsmJ | 16S rRNA m(2)G1516 methyltransferase, SAM‑dependent | |||||
| Rejected: Coverage evenness skew score above cutoff. | |||||||||||
| Rejected: Frequency below/above cutoff threshold. | |||||||||||
CGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCA > NC_000913/3642973‑3643216 | cgccAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAg < 5:33975/149‑1 (MQ=255) gccAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCC‑‑AACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAg > 5:89145/1‑146 (MQ=255) aGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGAt < 4:119249/149‑1 (MQ=255) cGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCC‑‑AACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAg < 6:89145/149‑1 (MQ=255) gCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCTCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAAtt < 8:22552/125‑1 (MQ=255) gCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAAtt > 7:22552/1‑125 (MQ=255) cATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAAt > 3:115518/1‑84 (MQ=255) cATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAAt < 4:115518/84‑1 (MQ=255) gTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCa > 4:106336/1‑149 (MQ=255) tCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACgg > 7:61746/1‑148 (MQ=255) tGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCt < 8:2845/95‑1 (MQ=255) tGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCt > 7:2845/1‑95 (MQ=255) tCCAGCCGCCAGCGGGCCGCCAGGACGGATAAGGCACGGGCTGCGGTGCCTGTGTCAGCAATTAAGCAGATTTTCACTGAGTGGTCGGCCCTTAGTGCCGTAATGCTCCAGCATCGCATCCGGCTGCGGGTCAGGACCAGGGAAGCGtg > 8:85982/1‑148 (MQ=255) tCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGtt > 1:97156/1‑149 (MQ=255) cAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTACGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGttt < 6:119792/148‑1 (MQ=255) aGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGa < 8:61746/149‑1 (MQ=255) gCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCTGTTCACGACCACGGAAGCGTTTGaa < 7:85982/149‑1 (MQ=255) gCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGaa < 3:106336/149‑1 (MQ=255) gATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGttt > 7:16841/1‑123 (MQ=255) gATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGttt < 8:16841/123‑1 (MQ=255) cACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAgcatc < 6:114266/81‑1 (MQ=255) cACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAgcatc > 5:114266/1‑81 (MQ=255) aCCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAAt < 6:124239/47‑1 (MQ=255) aCCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAAt > 5:124239/1‑47 (MQ=255) ccGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCa < 4:26926/125‑1 (MQ=255) | CGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCA > NC_000913/3642973‑3643216 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ 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. |