New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | 3079260 = | 35 (1.350) | 3 (0.140) | 3/112 | 2.1 | 10.4% | coding (71/444 nt) | cmtB | putative mannitol‑specific enzyme IIA component of PTS |
? | NC_000913 | 3080360 = | 22 (1.000) | coding (1276/1992 nt) | tktA | transketolase 1, thiamine triphosphate‑binding | |||||
Rejected: Frequency below/above cutoff threshold. |
TAATTTATGCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAGGAAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ < NC_000913/3079336‑3079260 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcaggaagcCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGT > NC_000913/3080360‑3080425 TAATTTATGCGGCTTAGTGATTATTTTCCAGAATCATCAGTCTCAGTTATACATTCAGCAAAAGATTG > 2:465674/1‑68 AATTTATGCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGG < 2:352968/68‑1 TTATGCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAG > 1:729900/1‑68 TTATGCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAG > 2:430360/1‑68 GCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCTGGAAG > 1:218784/1‑68 CGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAGGAAGC < 2:60518/68‑1 CGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAGGAAGC > 2:94663/1‑68 CAATCTCAGTTATACATTCAGCAAAAGATTGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAAT > 1:81549/1‑68 TATACATTCAGCAAAAGATTGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATA < 2:81549/68‑1 GCAAAAGATTGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGC < 1:226822/68‑1 GCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACC > 1:53966/1‑64 GCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACC < 2:53966/64‑1 AGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTG > 1:429342/1‑68 AAGCCACCGTGCAGGGAGTTACCGTTAGCAATCGCGGTCATACCGAACTCTCGAACTCCGTTGTGGAT < 1:258585/68‑1 GCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGT < 2:962479/68‑1 TAATTTATGCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAGGAAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ < NC_000913/3079336‑3079260 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcaggaagcCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGT > NC_000913/3080360‑3080425 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ 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. |