New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3955391 =97 (1.070)4 (0.050) 4/232 NT 4.2% coding (61/1545 nt) ilvA l‑threonine dehydratase, biosynthetic; also known as threonine deaminase
?NC_000913 3955421 = 92 (1.070)coding (91/1545 nt) ilvA l‑threonine dehydratase, biosynthetic; also known as threonine deaminase
Rejected: Frequency below/above cutoff threshold.

CATGGCGTATGCGCCGCGCAGCTTAAAGCTGTGCACTGGCTGGCGATCTTCGCGCTTCACCAGAATGACGTTATCAAGACGCGACGACAGTTTTTCCATTTTTTGTAGCGGCGTAACCTGCGCCGCCTCGTAAACCGGCGCGCGCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3955537‑3955391
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgcagGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCG  >  NC_000913/3955421‑3955567
                                                                                                                                                                                                                                                                                                      
CATGGCGTATGCGCCGCGCAGCTTAAAGCTGTGCACTGGCTGGCGATCTTCGCGCTTCACCAGAATGACGTTATCAAGACGCGACGACAGTTTTTCCATTTTTTGTAGCGGCGTAACCTGCGCCGCCTCGTAAACCGGCGCGCGCAG                                                                                                                                                     >  2:299481/1‑147
       TATGCGCCGCGCAGCTTAAAGCTGTGCACTGGCTGGCGATCTTCGCGCTTCACCAGAATGACGTTATCAAGACGCGACGACAGTTTTTCCATTTTTTGTAGCGGCGTAACCTGCGCCGCCTCGTAAACCGGCGCGCGCA                                                                                                                                                      <  7:282353/139‑1
       TATGCGCCGCGCAGCTTAAAGCTGTGCACTGGCTGGCGATCTTCGCGCTTCACCAGAATGACGTTATCAAGACGCGACGACAGTTTTTCCATTTTTTGTAGCGGCGTAACCTGCGCCGCCTCGTAAACCGGCGCGCGCA                                                                                                                                                      >  8:282353/1‑139
                                                                                                               CGTAACCTGCGCCGCCTCGTAAACCGGCGCGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACG                                     >  6:31137/1‑148
                                                                                                                    CCTGCGCCGCCTCGTAAACCGGCGCGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTA                                                       <  3:146114/125‑1
                                                                                                                    CCTGCGCCGCCTCGTAAACCGGCGCGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTA                                                       >  4:146114/1‑125
                                                                                                                            GCCTCGTAAACCGGCGCGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGG                        <  5:31137/148‑1
                                                                                                                                             GCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAG        >  6:9897/1‑147
                                                                                                                                               GCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATTGCGGGCCTGACGGAAGAACAGAAA     <  3:286805/148‑1
                                                                                                                                                CAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAG    <  6:122774/148‑1
                                                                                                                                                  GGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCG  <  5:317288/148‑1
                                                                                                                                                                                                                                                                                                      
CATGGCGTATGCGCCGCGCAGCTTAAAGCTGTGCACTGGCTGGCGATCTTCGCGCTTCACCAGAATGACGTTATCAAGACGCGACGACAGTTTTTCCATTTTTTGTAGCGGCGTAACCTGCGCCGCCTCGTAAACCGGCGCGCGCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3955537‑3955391
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgcagGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCG  >  NC_000913/3955421‑3955567

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.