New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | = 4166372 | 6 (0.170) | 6 (0.170) | 5/230 | 2.6 | 60.8% | intergenic (+87/‑287) | murI/rrsB | glutamate racemase/16S ribosomal RNA of rrnB operon |
? | NC_000913 | = 4166399 | 2 (0.060) | intergenic (+114/‑260) | murI/rrsB | glutamate racemase/16S ribosomal RNA of rrnB operon |
GCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAAC > NC_000913/4166281‑4166525 | gCTGATTTGGTAGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCgg > 7:93475/1‑147 (MQ=255) aTTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCagag > 3:132290/1‑148 (MQ=40) aGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGgtcag > 6:31390‑M2/1‑92 (MQ=255) ccacTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcaacattcaaccaaatcagcctaaaactgcca < 5:76145‑M2/139‑105 (MQ=255) ccacTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcaacattcaaccaaatcagcctaaaactgcca > 6:76145‑M2/1‑35 (MQ=255) ccacTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcaacattcaaccaaatc < 3:36397‑M2/124‑90 (MQ=255) ccacTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcaacattcaaccaaatc > 4:36397‑M2/1‑35 (MQ=255) acTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcaacattcaaccaaatcagcctaaaactgccagtttttcgagcg < 1:99753‑M2/149‑117 (MQ=255) aacaacGGCAAACACGCCGCCGTgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcaacattcaaccaaatcagcctaaaactgccagtttttcgagcgtttcgaagcc < 2:145737‑M2/149‑127 (MQ=255) | GCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAAC > NC_000913/4166281‑4166525 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 27 ≤ 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. |