New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 426590621 (0.710)5 (0.180) 5/248 2.6 19.6% coding (125/1080 nt) alr alanine racemase, biosynthetic, PLP‑binding
?NC_000913 = 4265911 21 (0.750)coding (130/1080 nt) alr alanine racemase, biosynthetic, PLP‑binding
Rejected: Frequency below/above cutoff threshold.

TATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCGGGGGGAATCACCAAACCTGTACTGTTACTCGAAGGCTTTTTTGATGCCAGAGATCTGCCGACGATT  >  NC_000913/4265767‑4266051
                                                                                                                                           |                                                                                                                                                 
tATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGagac                                                                                                                                          <  4:75356‑M2/149‑5 (MQ=255)
    aaGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGa                                                                                                                                             >  6:37766‑M2/1‑141 (MQ=255)
    aaGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGa                                                                                                                                             <  5:37766‑M2/142‑2 (MQ=255)
         acacAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTcttctt                                                                                                                                               >  6:129548/1‑135 (MQ=255)
         acacAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTcttctt                                                                                                                                               <  5:129548/135‑1 (MQ=255)
                               ttgtGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCtt                                                                                                            >  6:91869/1‑148 (MQ=255)
                               ttgtGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCtt                                                                                                            <  5:91869/148‑1 (MQ=255)
                                    aTTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTcaagaagaccgtga                                                                                                                                     <  8:99717‑M1/118‑15 (MQ=255)
                                    aTTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTcaagaagaccgtga                                                                                                                                     >  7:99717‑M1/1‑104 (MQ=255)
                                     ttAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTc                                                                                                                             >  7:26616/1‑125 (MQ=255)
                                     ttAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTc                                                                                                                             <  8:26616/125‑1 (MQ=255)
                                      tAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGc                                                                                                        <  7:110634/145‑1 (MQ=255)
                                      tAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGc                                                                                                        >  8:110634/1‑145 (MQ=255)
                                                                 gCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACtgc                                                                         >  2:134697/1‑149 (MQ=255)
                                                                  cAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAAc                                                                                                                                 <  2:159513/92‑1 (MQ=255)
                                                                  cAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAAc                                                                                                                                 >  1:159513/1‑92 (MQ=255)
                                                                       tCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCt                                                                                                                              <  8:45471/90‑1 (MQ=255)
                                                                       tCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCt                                                                                                                              >  7:45471/1‑90 (MQ=255)
                                                                       tCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCgg                                                                   <  3:75356/149‑1 (MQ=255)
                                                                         ccattttaCTGGCAGCGGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCACTGCGACTGCGAGCgggg                                                                 >  6:83042/8‑149 (MQ=255)
                                                                           cGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCgggg                                                                 <  1:134697/147‑1 (MQ=255)
                                                                                         tGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTcaagaagaccgtga                                                                                                                                     >  2:46867‑M1/1‑51 (MQ=255)
                                                                                           ccAGTAAAATGGTTGCGGTGGTGAAAGCGATCGCTTATGGTCACGGTCTGCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGGCTTTGGCGTAGCCCGGCTCGAAGAAGCTCTGCGTCTGCGTGCGGGGGGAATCACCAAACCTGTa                                               <  5:83042/149‑1 (MQ=255)
                                                                                           ccAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTcaagaagaccgtgaccataagcgttcgctttca                                                                                                                  >  2:122471‑M1/1‑49 (MQ=255)
                                                                                           ccAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTcaagaagaccgtgaccataagcgttcgctttca                                                                                                                  <  1:122471‑M1/82‑34 (MQ=255)
                                                                                           ccAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCGGGGGGAATCACCAAACCTGTa                                               <  5:80992/149‑1 (MQ=255)
                                                                                             aGTAAAATGGTTGCGGTGGTGAAAGCGAACGCCTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCGGGGGGAATCACCAAACCTGTACt                                             >  3:106653/1‑149 (MQ=255)
                                                                                                                tgAAAGCGAACGCTTATGGTCACGGTCTTCTTGagaccgtgaccataagcgttcgctttcaccaccgcaaccattttactgg                                                                                             <  1:122471‑M2/82‑50 (MQ=255)
                                                                                                                tgAAAGCGAACGCTTATGGTCACGGTCTTCTTGagaccgtgaccataagcgttcgctttcaccaccgcaaccattttactgg                                                                                             >  2:122471‑M2/1‑33 (MQ=255)
                                                                                                                tgAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCGGGGGGAATCACCAAACCTGTACTGTTACTCGAAGGCtttttt                          >  6:149326/1‑149 (MQ=255)
                                                                                                                   aaGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGTGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCGGGGGGAATCACCAAACCTGTACTGTTACTCGAAGGCTTTTTTGAt                       <  1:135905/149‑1 (MQ=255)
                                                                                                                                   tCACGGTCTTCTTGagaccgtgaccataagcgttcgctttcaccaccgcaaccattttactggcaggggccagttcacgaagacgttgcaggttgtgtcgcagagcgcggcggttaat                                      >  7:99717‑M2/1‑14 (MQ=255)
                                                                                                                                   tCACGGTCTTCTTGagaccgtgaccataagcgttcgctttcaccaccgcaaccattttactggcaggggccagttcacgaagacgttgcaggttgtgtcgcagagcgcggcggttaat                                      <  8:99717‑M2/118‑105 (MQ=255)
                                                                                                                                   tCACGGTCTTCTTGagaccgtgaccataagcgttcgctttcaccaccgcaaccattttactggca                                                                                           >  2:46867‑M2/1‑14 (MQ=255)
                                                                                                                                        gTCTcaagaagaccgtgaccataagcgttcgctttcaccaccgcaaccattttactggcaggggccagttcacgaagacgttgcaggttgtgtcgcagagcgcggcggttaatcacaacagttgccgcttgcatttgtgttccttgata  <  4:75356‑M1/149‑146 (MQ=255)
                                                                                                                                           tcaagaagaccgtgaccataagcgttcgctttcaccaccgcaaccattttactggcaggggccagttcacgaagacgttgcaggttgtgtcgcagagcgcggcggttaatcacaacagttgccgcttgcatttgtgttcctt      >  6:37766‑M1/1‑1 (MQ=255)
                                                                                                                                           tcaagaagaccgtgaccataagcgttcgctttcaccaccgcaaccattttactggcaggggccagttcacgaagacgttgcaggttgtgtcgcagagcgcggcggttaatcacaacagttgccgcttgcatttgtgttcctt      <  5:37766‑M1/142‑142 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
TATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCGGGGGGAATCACCAAACCTGTACTGTTACTCGAAGGCTTTTTTGATGCCAGAGATCTGCCGACGATT  >  NC_000913/4265767‑4266051

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