New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 421449263 (1.740)5 (0.150) 5/224 2.6 7.4% coding (213/930 nt) metA homoserine O‑transsuccinylase
?NC_000913 = 4214509 67 (1.980)coding (230/930 nt) metA homoserine O‑transsuccinylase
Rejected: Frequency below/above cutoff threshold.

TTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCG‑‑‑C‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4214340‑4214492
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgttgcg‑‑‑cGATTCACGGGAATCGATGCGCAACAGCTGAATATCGACCTGCAAAGGTGAGTTTGAAAGCAGGCGCAGAAACTGA  <  NC_000913/4214509‑4214435
                                                                                                                                                                                                                                       
TTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTT                                                                                    >  1:66214/1‑149
    TGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCG                                                                                 >  5:129758/1‑148
     GATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCG‑‑‑C                                                                             >  6:85531/1‑148
       TGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCG‑‑‑CGAT                                                                          >  6:123892/1‑149
                             TCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCG‑‑‑CGATTCACGGGAATCGATGCGCAACA                                                    <  5:123892/149‑1
                                  AAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCG‑‑‑CGATTCACGGGAATCGATGCGCAACAGCTGA                                               <  5:85531/149‑1
                                                                                                                          ACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATC‑‑‑GCGCAACAGCTGAATATCGACCTGCAAAGGTGAGTTTGAAAGCAGGCGCAGAAACTGA  <  5:169354/106‑1
                                                                                                                          ACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATC‑‑‑GCGCAACAGCTGAATATCGACCTGCAAAGGTGAGTTTGAAAGCAGGCGCAGAAACTGA  >  6:169354/1‑106
                                                                                                                                                                                                                                       
TTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCG‑‑‑C‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4214340‑4214492
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgttgcg‑‑‑cGATTCACGGGAATCGATGCGCAACAGCTGAATATCGACCTGCAAAGGTGAGTTTGAAAGCAGGCGCAGAAACTGA  <  NC_000913/4214509‑4214435

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