New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4576863 =45 (1.620)4 (0.160) 4/246 2.9 8.8% intergenic (+8/+49) yjiV/mcrC pseudogene; conserved hypothetical protein/5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
?NC_000913 4576890 = 41 (1.600)intergenic (+35/+22) yjiV/mcrC pseudogene; conserved hypothetical protein/5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Rejected: Frequency below/above cutoff threshold.

ATAAAATTAATGGCTTCGATATTGGCTTGTGTACCGTCAATTTAGGTCAGGAATGGCCGTGTATACATCAAGAATTACTCGACATTTTCGATGAATATCTCAAATAAGCGGTAAGTTTTATAGAAAATACCGCTCCCGACACCACATCGGGAGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4577018‑4576863
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgggagcggTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGC  >  NC_000913/4576890‑4577034
                                                                                                                                                                                                                                                                                                             
ATAAAATTAATGGCTACGATATTGGCTTGTGTACCGTCAATTTAGGTCAGGAATGGCCGTGTATACATCAAGAATTACTCGACATTTTCGATGAATATCTCAAATAAGCGGTAAGTTTTATAGAAAATACCGCTCCCGACACCACATC                                                                                                                                                           <  2:101264/148‑1
  AAAATTAATGGCTTCGATATTGGCTTGTGTACCGTCAATTTAGGTCAGGAATGGCCGTGTATACATCAAGAATTACTCGACATTTTCGATGAATATCTCAAATAAGCGGTAAGTTTTATAGAAAATACCGCTCCCGACACCACATCGG                                                                                                                                                         >  2:23393/1‑148
                  ATATTGGCTTGTGTACCGTCAATTTAGGTCAGGAATGGCCGTGTATACATCAAGAATTACTCGACATTTTCGATGAATATCTCAAATAAGCGGTAAGTTTTATCGAAAATACCGCTCCCGATGTGGCGTCGGGAGCGGTATTTTCTATA                                                                                                                                        <  4:82266/149‑1
                                                                                                                          GAAAATACCGCTCCAGACACCACATCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTATTGATGTATA                                                                               >  6:102834/1‑102
                                                                                                                             AATACCGCTCCCGACACCACATCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAG                                                                                               <  1:114115/83‑1
                                                                                                                             AATACCGCTCCCGACACCACATCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAG                                                                                               >  2:114115/1‑83
                                                                                                                                                    GGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCA      <  6:126655/149‑1
                                                                                                                                                        GCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGC  <  4:88160/149‑1
                                                                                                                                                                                                                                                                                                             
ATAAAATTAATGGCTTCGATATTGGCTTGTGTACCGTCAATTTAGGTCAGGAATGGCCGTGTATACATCAAGAATTACTCGACATTTTCGATGAATATCTCAAATAAGCGGTAAGTTTTATAGAAAATACCGCTCCCGACACCACATCGGGAGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4577018‑4576863
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgggagcggTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGC  >  NC_000913/4576890‑4577034

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ 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.