New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4627665128 (0.920)6 (0.050) 6/170 NT 4.6% coding (351/1233 nt) nadR trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
?NC_000913 = 4627683 127 (0.950)coding (369/1233 nt) nadR trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Rejected: Frequency below cutoff threshold.

TATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4627563‑4627665
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagACGATCCGGCACGGTCGGCTGCTGCGACATGGCACTGTCTTCGAACAACGCGCGGTCACGGGTATCGTCAAAACCCATAATGATATGCAGCTCGTCAAC  <  NC_000913/4627683‑4627585
                                                                                                                                                                                                          
TATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAG                                                                                                        <  2:656039/100‑1
TATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAG                                                                                                        <  2:265787/100‑1
 ATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGC                                                                                                       <  1:1175724/100‑1
 ATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGC                                                                                                       <  1:1354004/100‑1
 ATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGC                                                                                                       <  2:2920276/100‑1
 ATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGC                                                                                                       >  2:1085553/1‑100
 ATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGC                                                                                                       <  2:2053152/100‑1
  TCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCA                                                                                                      >  2:600065/1‑100
  TCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCA                                                                                                      >  1:3783568/1‑100
  TCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCA                                                                                                      >  2:2104695/1‑100
                           CGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGACGATCCGGCACGGTCGGCTGCTG                                                                             <  1:1176971/100‑1
                             AGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGACGATCCGGCACGGTCGGCTGCTGCG                                                                           <  1:637085/100‑1
                                 GCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGACGATCCGGCACGGTCGGCTGCTGC                                                                            <  1:2990566/95‑1
                                   ATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGACGATCCGGCACGGTCGGCTGCTGCGACATGG                                                                     <  2:121603/100‑1
                                    TATCATTATGTGTTTTGACATTACCCGTGACCGCGCATTGTGCGAAGACAGTGCCATGTCGCAGCAGACGATACGGCACGGTCGGCTGCTGCGACATG                                                                      <  1:2200733/98‑1
                                      TCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGACGATCCGGCACGGTCGGCTGCTGCGACATGGCAC                                                                  <  2:1190053/100‑1
                                                                                                   GCAGACGATCCGGCACGGTCGGCTGCTGCGACATGGCACTGTCTTCGAACAACGCGCGGTCACGGGTATCGTCAAAACCCATAATGATATGCAGCTCGTC     <  2:1321900/100‑1
                                                                                                    CAGACGATCCGGGACGGTCGGCTGCTGCGACATGGCACTGTCTTCGAACAACGCCCGGTCACGGGTAGCGTGAACACCCATAATGATATGCAGGTCGTC     >  2:3782510/1‑99
                                                                                                    CAGACGATCCGGCACGGTCGGCTGCTGCGACATGGCACTGTCTTCGAACAACGCGCGGTCACGGGTATCGTCAAAACCCATAATGATATGCAGCTCGTCA    >  1:74718/1‑100
                                                                                                    CAGACGATCCGGCACGGTCGGCTGCTGCGACATGGCACTGTCTTCGAACAACGCGCGGTCACGGGTATCGTCAAAACCCATAATGATATGCAGCTCGTCA    >  2:350595/1‑100
                                                                                                     AGACGATCCGGCACGGTCGGCTGCTGCGACATGGCACTGTCTTCGAACAACGCGCGGTCACGGGTATCGTCAAAACCCATAATGATATGCAGCTCGTCAA   >  1:3713426/1‑100
                                                                                                      GACGATCCGGCACGGTCGGCTGCTGCGACATGGCACTGTCTTCGAACAACGCGCGGTCACGGGTATCGTCAAAACCCATAATGATATGCAGCTCGTCAAC  <  2:2633764/100‑1
                                                                                                                                                                                                          
TATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4627563‑4627665
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagACGATCCGGCACGGTCGGCTGCTGCGACATGGCACTGTCTTCGAACAACGCGCGGTCACGGGTATCGTCAAAACCCATAATGATATGCAGCTCGTCAAC  <  NC_000913/4627683‑4627585

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.