New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4638648 =116 (1.030)4 (0.040) 3/166 NT 3.6% coding (471/1353 nt) creD inner membrane protein
?NC_000913 4638675 = 106 (1.010)coding (498/1353 nt) creD inner membrane protein
Rejected: Frequency below cutoff threshold.

CGCCAGTCCCCTTCAGGTAAAGGGATATGCACGCCCTGCCCGCCTTGCTCTAACCCGGTGCCGGGTTCAATGGTCAGCGCCGTTCCGTTAACTTCAGGCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4638749‑4638648
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgctGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATA  >  NC_000913/4638675‑4638772
                                                                                                                                                                                                        
CGCCAGTCCCCTTCAGGTAAAGGGATATGCACGCCCTGCCCGCCTTGCTCTAACCCGGTGCCGGGTTCAATGGTCAGCGCCGTTCCGTTAACTTCAGGCG                                                                                                      <  1:2613297/100‑1
 GCCAGTCCCCTTCAGGTAAAGGGATATGCACGCCCTGCCCGCCTTGCTCTAACCCGGTGCCGGGTTCAATGGTCAGCGCCGTTCCGTTAACTTCAGGCGC                                                                                                     >  1:2834761/1‑100
 GCCAGTCCCCTTCAGGTAAAGGGATATGCACGCCCTGCCCGCCTTGCTCTAACCCGGTGCCGGGTTCAATGGTCAGCGCCGTTCCGTTAACTTCAGGCGC                                                                                                     >  2:798264/1‑100
                                                                    AATGGTCAGCGCCGTTCCGTTAACTTCAGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGG                                  >  1:2566737/1‑100
                                                                        GTCAGCGCCGTTCCGTTAACTTCAGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGAC                              >  1:12856/1‑100
                                                                        GTCAGCGCCGTTCCGTTAACTTCAGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGAC                              >  1:12876/1‑100
                                                                          CAGCGCCGTTCCGTTAACTTCAGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTG                            >  1:1384673/1‑100
                                                                                                GGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTG      >  1:7994/1‑100
                                                                                                GGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTG      >  2:20951/1‑100
                                                                                                 GCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGA     <  2:2302161/100‑1
                                                                                                 GCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGAATGGCGGAAGCAGAACCTGAAGCTGA     >  1:784331/1‑100
                                                                                                  CGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAA    <  1:2003730/100‑1
                                                                                                  CGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAA    <  2:1921068/100‑1
                                                                                                    CTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATA  <  2:2945812/100‑1
                                                                                                    CTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATA  <  2:3049493/100‑1
                                                                                                                                                                                                        
CGCCAGTCCCCTTCAGGTAAAGGGATATGCACGCCCTGCCCGCCTTGCTCTAACCCGGTGCCGGGTTCAATGGTCAGCGCCGTTCCGTTAACTTCAGGCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4638749‑4638648
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgctGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATA  >  NC_000913/4638675‑4638772

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