New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | = 2333629 | 15 (0.490) | 4 (0.150) | 4/112 | 2.1 | 18.7% | pseudogene (774/4605 nt) | yfaS | pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein |
? | NC_000913 | = 2334464 | 22 (0.850) | coding (496/624 nt) | yfaT | DUF1175 family protein, putative host defense protein | |||||
Rejected: Frequency below/above cutoff threshold. |
CTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTAC > NC_000913/2333562‑2333692 | cttcACCCTTTTTTTTACCCGCTGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCTCG < 2:552555‑M1/68‑1 (MQ=255) aCCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACTGTAt < 1:389268/68‑1 (MQ=255) tttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaac < 2:341454/68‑1 (MQ=255) ttttACCAGCGGTCAACACCAGAAGCTCTTAGCCTGACACTTTGCATAGCGCAACGGTAACGGaaaca > 2:861952/1‑68 (MQ=255) aCCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAaacaaaca > 1:354812/1‑68 (MQ=255) aCCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAaacaaaca > 2:1106960/1‑68 (MQ=255) gTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACcgt < 1:186741/68‑1 (MQ=255) tCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACcgtc > 1:687110/1‑68 (MQ=255) cacCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGaaaactgccaacggcatgcgcgcagtca > 1:549548‑M1/1‑40 (MQ=255) acCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGc < 2:490153/68‑1 (MQ=255) ccAGAGGCTCTTTGCCTGCCACTTTGCTACGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGcc > 2:308020/1‑68 (MQ=255) ccAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGaaaactgacaacggaatgcgcgca < 2:637355‑M1/62‑25 (MQ=255) ccAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGaaaactgacaacggaatgcgcgca > 1:637356‑M1/1‑38 (MQ=255) aGCTCTTTGCCTGACACTTTGCTAAGCGCCACGaaaactgacaacggaatgcgcgcagtcagtctgcg > 2:392541‑M1/1‑33 (MQ=255) ctctTTGCCTGACGCTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATAcc > 2:584857/1‑68 (MQ=255) acTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGAAACCCACCAACCATCGc > 2:234117/1‑68 (MQ=255) tGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTcg > 1:883108/1‑68 (MQ=255) tGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTcg > 1:485940/1‑68 (MQ=255) cttaGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTcgac < 2:533645/65‑1 (MQ=255) aaGCGCCACGaaaactgacaacggaatgcgcgcagtcagtctgcaacaacttatgacatggaaggaca < 1:658663‑M1/68‑59 (MQ=255) aGCGCCACGaaaactgacaacggaatgcgcgcagtcagtctgcaacaacttatgacatggaaggacac > 1:439592‑M1/1‑9 (MQ=255) gcgcCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAgg < 2:935120/68‑1 (MQ=255) ccACGaaaactgacaacggaatgcgcgcagtcagtctgcaacaacttatgacatggaaggacacccga < 1:1153544‑M1/68‑64 (MQ=255) ccACGaaaactgacaacggaatgcgcgcagtcagtctgcaacaacttatgacatggaaggacacccga < 1:1152030‑M1/68‑64 (MQ=255) | CTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTAC > NC_000913/2333562‑2333692 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 37 ≤ 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. |