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.

TATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCGGGGGGAATCACCAAACCTGTACTGTTACTCGAAGGCTTTTTTGAT  >  NC_000913/4265767‑4266030
                                                                                                                                                |                                                                                                                       
tATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGagac                                                                                                                     <  4:75356‑M2/149‑5 (MQ=255)
    aaGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGa                                                                                                                        >  6:37766‑M2/1‑141 (MQ=255)
    aaGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGa                                                                                                                        <  5:37766‑M2/142‑2 (MQ=255)
                               ttgtGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCtt                                                                                       <  5:91869/148‑1 (MQ=255)
                               ttgtGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCtt                                                                                       >  6:91869/1‑148 (MQ=255)
                                     ttAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTc                                                                                                        <  8:26616/125‑1 (MQ=255)
                                     ttAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTc                                                                                                        >  7:26616/1‑125 (MQ=255)
                                      tAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGc                                                                                   >  8:110634/1‑145 (MQ=255)
                                      tAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGc                                                                                   <  7:110634/145‑1 (MQ=255)
                                                                 gCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACtgc                                                    >  2:134697/1‑149 (MQ=255)
                                                                  cAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAAc                                                                                                            >  1:159513/1‑92 (MQ=255)
                                                                  cAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAAc                                                                                                            <  2:159513/92‑1 (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)
                                                                                           ccAGTAAAATGGTTGCGGTGGTGAAAGCGATCGCTTATGGTCACGGTCTGCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGGCTTTGGCGTAGCCCGGCTCGAAGAAGCTCTGCGTCTGCGTGCGGGGGGAATCACCAAACCTGTa                          <  5:83042/149‑1 (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)
                                                                                                                                                gagaCCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCGGGGGGAATCACCAAAc                               >  6:71815/1‑91 (MQ=255)
                                                                                                                                                gagaCCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCGGGGGGAATCACCAAAc                               <  5:71815/91‑1 (MQ=255)
                                                                                                                                                |                                                                                                                       
TATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAGACCGCGCGAACGCTCCCCGATGCTGACGCCTTTGGCGTAGCCCGTCTCGAAGAAGCTCTGCGACTGCGTGCGGGGGGAATCACCAAACCTGTACTGTTACTCGAAGGCTTTTTTGAT  >  NC_000913/4265767‑4266030

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.