New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4459872104 (1.140)5 (0.060) 5/242 NT 4.8% intergenic (+17/+28) cybC/nrdG pseudogene, truncated cytochrome b562; cytochrome b(562)/anaerobic ribonucleotide reductase activating protein
?NC_000913 = 4459886 103 (1.220)intergenic (+31/+14) cybC/nrdG pseudogene, truncated cytochrome b562; cytochrome b(562)/anaerobic ribonucleotide reductase activating protein
Rejected: Frequency below/above cutoff threshold.

CGCCACGGTTTCGACATTCTGGTCGGTCAGATTGACGACGCGCTGAAGCTGGCAAATGAAGGTAAAGTAAAAGAAGCGCAGGCTGCTGCAGAGCAACTGAAAACGACCCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTTCCCTGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4459718‑4459872
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttccctgttACCTGAGTGCAGGCAACAGGGAAATGAGGAATTAACGATACTTCTTGTGATATGAATTGCGGGTCGTTTTCAGTTGCTCTGCAGCAGCCTGCGCTTCTTTTACTTTACCTTCATTTGCCAGCTTCAGCGCGTCGTCAATCTGAC  <  NC_000913/4459886‑4459743
                                                                                                                                                                                                                                                                                                           
CGCCACGGTTTCGACATTCTGGTCGGTCAGATTGACGACGCGCTGAAGCTGGCAAATGAAGGTAAAGTAAAAGAAGCGCAGGCTGCTGCAGAGCAACTGAAAACTACACGCAATTCATATCACAAGAAGTATCGTTAATTCATCATTTC                                                                                                                                                        >  6:191370/1‑149
 GCCACGGTTTCGACATTCTGGTCGGTCAGATTGACGACGCGCTGAAGCTGGCAAATGAAGGTAAAGTAAAAGAAGCGCAGGCTGCTGCAGAGCAACTGAAAACGACCCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTTCC                                                                                                                                                       <  1:201210/149‑1
 GCCACGGTTTCGACATTCTGGTCGGTCAGATTGACGACGCGCTGAAGCTGGCAAATGAAGGTAAAGTAAAAGAAGCGCAGGCTGCTGCAGAGCAACTGAAAACGACCCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTTCC                                                                                                                                                       <  5:75080/149‑1
  CCACGGTTTCGACATTCTGGTCGGTCAGATTGACGACGCGCTGAAGCTGGCAAATGAAGGTAAAGTAAAAGAAGCGCAGGCTGCTGCAGAGCAACTGAAAACGACCCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTTCCC                                                                                                                                                      <  8:150843/149‑1
   CACGGTTTCGACATTCTGGTCGGTCAGATTGACGACGCGCTGAAGCTGGCAAATGAAGGTAAAGTAAAAGAAGCGCAGGCTGCTGCAGAGCAACTGAAAACGACCCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTTCCCT                                                                                                                                                     >  4:299838/1‑149
   CACGGTTTCGACATTCTGGTCGGTCAGATTGACGACGCGCTGAAGCTGGCAAATGAAGGTAAAGTAAAAGAAGCGCAGGCTGCTGCAGAGCAACTGAAAACGACCCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTTCCCT                                                                                                                                                     >  6:336306/1‑149
                                                     AAATGAAGGTAAAGTAAAAGAAGCGCAGGCTGCTGCAGAGCAACTGAAAACGACCCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTTCCCTGTTACCTGAGTGCAGGC                                                                                                                                    >  7:64986/1‑116
                                                     AAATGAAGGTAAAGTAAAAGAAGCGCAGGCTGCTGCAGAGCAACTGAAAACGACCCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTTCCCTGTTACCTGAGTGCAGGC                                                                                                                                    <  8:64986/116‑1
                                                                                                    AAACGACCCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTTCCCTGTTGCCTGCACTCAGGTAACAGGGAAATGAGGAATTAACGATACTTCTTGTGATATGAATTGCGGGTCGTTTTCAGTTGCTCTGCAGCAGCCTGCGC                                                    >  6:73556/1‑149
                                                                                                        GACCCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTTCCCTGTTGCCTGCACTCAGGTAACAGGGAAATGAGGAATTAACGATACTTCTTGTGATATGAATTGCGGGTCGTTTTCAGTTGCTCTGCAGCAGCCTGCGCTTCT                                                <  5:73556/149‑1
                                                                                                           CCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTACCCTGTTGCCTGCACTCAGGTAACAGGGAAATGAGGAATTAACGATACTTCTTGTGATATGAATTGCGGGTCGTTTTCAGTTGCTCTGCAGCAGCCTGCGCTTCTTT                                              >  2:171453/1‑148
                                                                                                                                                     CCTGTTACCTGAGTGCAGGCAACAGGGAAATGAGGAATTAACGATACTTCTTGTGATATGAATTGCGGGTCGTTTTCAGTTGCTCTGCAGCAGCCTGCGCTTCTTTTACTTTACCTTCATTTGCCAGCTTCAGCGCGTCGTCAATCTGA   <  4:292254/149‑1
                                                                                                                                                      CTGTTACCTGAGTGCAGGCAACAGGGAAATGAGGAATTAACGATACTTCTTGTGATATGAATTGCGGGTCGTTTTCAGTTGCTCTGCAGCAGCCTGCGCTTCTTTTACTTTACCTTCATTTGCCAGCTTCAGCGCGTCGTCAATCTGAC  <  3:368125/149‑1
                                                                                                                                                                                                                                                                                                           
CGCCACGGTTTCGACATTCTGGTCGGTCAGATTGACGACGCGCTGAAGCTGGCAAATGAAGGTAAAGTAAAAGAAGCGCAGGCTGCTGCAGAGCAACTGAAAACGACCCGCAATTCATATCACAAGAAGTATCGTTAATTCCTCATTTCCCTGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4459718‑4459872
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttccctgttACCTGAGTGCAGGCAACAGGGAAATGAGGAATTAACGATACTTCTTGTGATATGAATTGCGGGTCGTTTTCAGTTGCTCTGCAGCAGCCTGCGCTTCTTTTACTTTACCTTCATTTGCCAGCTTCAGCGCGTCGTCAATCTGAC  <  NC_000913/4459886‑4459743

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