New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4627662 =94 (0.940)6 (0.060) 4/170 NT 6.3% coding (348/1233 nt) nadR trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
?NC_000913 4627688 = 89 (0.930)coding (374/1233 nt) nadR trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Rejected: Frequency below cutoff threshold.

ATACGGCTCCATGCCCTCTTCGTTGAAAGCATGAATGCGAATATTTTTCTGATATTTAAAAGTTTGCAATAACCAACGCAGACGATCCGGCACGGTCGGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4627764‑4627662
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgcGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGA  >  NC_000913/4627688‑4627786
                                                                                                                                                                                                          
ATACGGCTCCATGCCCTCTTCGTTGAAAGCATGAATGCGAATATTTTTCTGATATTTAAAAGTTTGCAATAACCAACGCAGACGATCCGGCACGGTCGGC                                                                                                        >  1:1891299/1‑100
  ACGTCTACATGCCCTCTTCGTTGAAAGCATTAACGCGAATATTTTTCTGATATTTAAAAGTTTGCAATAACCAACGCAGACGATCCGGCACGGTCGGC                                                                                                        <  2:86662/98‑1
   CGGCTCCATGCCCTCTTCGTTGAAAGCATGAATGCGAATATTTTTCTGATATTTAAAAGTTTGCAATAACCAACGCAGACGATCCGGCACGGTCGGCTGC                                                                                                     >  1:2665994/1‑100
   CGGCTCCATGCCCTCTTCGTTGAAAGCATGAATGCGAATATTTTTCTGATATTTAAAAGTTTGCAATAACCAACGCAGACGATCCGGCACGGTCGGCTGC                                                                                                     <  1:753972/100‑1
    GGCTCCATGCGCTCTTCGTTGAAAGCATGCATGCGAATATTGTTCTGATATTTAAAAGTTTGCAATAACCAACGCAGACGATCCGGCACGGTCGGCTGC                                                                                                     <  2:770937/99‑1
                                                                          AACGCAGACGATCCGGCACGGTCGGCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAG                              >  1:1457047/1‑100
                                                                           ACGCAGACGATCCGGCACGGTCGGCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAG                              >  1:2476057/1‑99
                                                                             GCAGACGATCCGGCACGGTCGGCTGCGTTGGTTATTGCAACCTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCG                           >  2:2440502/1‑100
                                                                             GCAGACGATCCGGCACGGTCGGCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCG                           >  1:2013157/1‑100
                                                                             GCAGACGATCCGGCACGGTCGGCTGCGTTGGTTATTGAAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAGGGAGCCG                           >  2:1474638/1‑100
                                                                               AGACGATCCGGCACGGTCGGCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGT                          >  1:705188/1‑99
                                                                                                    TGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTG    >  2:787388/1‑100
                                                                                                     GCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGG   <  1:33369/100‑1
                                                                                                     GCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGG   <  1:666903/100‑1
                                                                                                      CGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGA  >  1:2697069/1‑100
                                                                                                      CGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGA  <  1:427134/100‑1
                                                                                                                                                                                                          
ATACGGCTCCATGCCCTCTTCGTTGAAAGCATGAATGCGAATATTTTTCTGATATTTAAAAGTTTGCAATAACCAACGCAGACGATCCGGCACGGTCGGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4627764‑4627662
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgcGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGA  >  NC_000913/4627688‑4627786

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