New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4129541 =78 (0.640)3 (0.030) 3/240 NT 3.9% coding (870/1161 nt) metB cystathionine gamma‑synthase, PLP‑dependent
?NC_000913 4129565 = 72 (0.620)coding (894/1161 nt) metB cystathionine gamma‑synthase, PLP‑dependent
Rejected: Frequency below/above cutoff threshold.

GTCATGGTTGCGGCGTGAGAGATTAAACTTTCCACTCCCCCTAATGATTCCGCCAGCGTAAACAACGACAGCCCGCCCAGGAAACGACGCAGCGTCTGCTCATCGCCATCCAGTTCAAAACTCAACATTGCGCCAAAGCCTTTTTGCTGGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4129693‑4129541
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcgcAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCA  >  NC_000913/4129565‑4129709
                                                                                                                                                                                                                                                                                                          
GTCATGGTTGCGGCGTGAGAGATTAAACTTTCCACTCCCCCTAATGATTCCGCCAGCGTAAACAACGACAGCCCGCCCAGGAAACGACGCAGCGTCTGCTCATCGCCATCCAGTTCAAAACTCAACATTGCGCCAAAGCCTTTTTGCTG                                                                                                                                                       >  4:293994/1‑149
 TCATGGTTGCGGCGTGAGAGATTAAACTTTCCACTCCCCCTAATGATTCCGCCAGCGTAAACAACGACAGCCCGCCCAGGAAACGACGCAGCGTCTGCTCATCGCCATCCAGTTCAAAACTCAACATTGCGCCAAAGCCTTTTTGCTGG                                                                                                                                                      >  4:987177/1‑149
   ATGGTTGCGGCGTGAGAGATTAAACTTTCCACTCCCCCTAATGATTCCGCCAGCGTAAACAACGACAGCCCGCCCAGGAAACGACGCAGCGTCTGCTCATCGCCATCCAGTTCAAAACTCAACATTGCGCCAAAGCCTTTTTGCTGGCG                                                                                                                                                    >  2:167489/1‑149
   ATGGTTGCGGCGTGAGAGATTAAACTTTCCACTCCCCCTAATGATTCCGCCAGCGTAAACAACGACAGCCCGCCCAGGAAACGACGCAGCGTCTGCTCATCGCCATCCAGTTCAAAACTCAACATTGCGCCAAAGCCTTTTTGCTGGCG                                                                                                                                                    >  4:690337/1‑149
                         AACTTTCCACTCCCCCTAATGATTCCGCCAGCGTAAACAACGACAGCCCGCCCAGGAAACGACGCAGCGTCTGCTCATCGCCATCCAGTTCAAAACTCAACATTGCGCCAAAGCCTTTTTGCTG                                                                                                                                                       <  3:816402/124‑1
                         AACTTTCCACTCCCCCTAATGATTCCGCCAGCGTAAACAACGACAGCCCGCCCAGGAAACGACGCAGCGTCTGCTCATCGCCATCCAGTTCAAAACTCAACATTGCGCCAAAGCCTTTTTGCTG                                                                                                                                                       >  4:816402/1‑124
                                                                                                                     AAACTCAACATTGCGCCAAAGCCTTTTTGCTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTC                                  >  1:1077725/1‑149
                                                                                                                        CTCAACATTGCGCCAAAGCCTTTTTGCTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGT                                                   <  3:500448/129‑1
                                                                                                                        CTCAACATTGCGCCAAAGCCTTTTTGCTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGT                                                   >  4:500448/1‑129
                                                                                                                                                    GGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGC   <  4:165438/149‑1
                                                                                                                                                     GCGCAATGTTGAGTTTTGAACTGGATGACGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCA  <  2:515105/149‑1
                                                                                                                                                                                                                                                                                                          
GTCATGGTTGCGGCGTGAGAGATTAAACTTTCCACTCCCCCTAATGATTCCGCCAGCGTAAACAACGACAGCCCGCCCAGGAAACGACGCAGCGTCTGCTCATCGCCATCCAGTTCAAAACTCAACATTGCGCCAAAGCCTTTTTGCTGGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4129693‑4129541
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcgcAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCA  >  NC_000913/4129565‑4129709

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