New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 250236	36 (0.800)	5 (0.120)	4/228	NT	13.3%	coding (967/1347 nt)	nfeD	membrane bound peptidase NefD
	?	NC_002947	= 250260	32 (0.770)					coding (943/1347 nt)	nfeD	membrane bound peptidase NefD

TTGAGAGCCATGCCCAGCACGCCACCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/250088‑250236 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgaagccTTCTTGCCCAGTTTCGGCGTGGTCGGCTTCGGCGGCATCGTTGCCTTTGTGGTCGGCGCATTGATCCTGATCGACACCGACGCTCCCGGTTTCGGCATCCCCCTGGCCTTGATCGGCACCTTGGCGCTGCTCTCGGCACTCCTGATCGG  <  NC_002947/250260‑250112                                                                                                                                                                                                                                                                                                            TTGAGAGCCATGCCCAGCACGCCACCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCT                                                                                                                                                      >  2:676336/1‑150              CCAGCACGCCACCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCTTCTTGCCCAGTTT                                                                                                                                         >  1:646404/1‑150                               AGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCTTCTTGCCCAGTTTCGGCGTGGTCGGCTTCG                                                                                                                        <  1:676336/150‑1                                                 GCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCTTCTTGCCCAGTTTCGGCGTGGTCGGCTTCGGCGGCATCGTTGCCTTTG                                                                                                      >  2:146530/1‑150                                                   CAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCTTCTTGCCCAGTTTCGGCGTGGTCGGCTTCGGCGGCATCGTTGCCTTTGTG                                                                                                    <  1:146530/150‑1                                                   CAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCTTCTTGCCCAGTTTCGGCGTGGTCGGCTTCGGCGGCATCGTTGCCTTTGTG                                                                                                    <  2:646404/150‑1                                                                                                                                              CCGAATCCTTCTTACCCAGTTTCGGCGTGGTCGGCTTCGGCGGCATCGTTGCCTTTGTGGTCGGCGCATTGATCCTGATCGACACCGACGCTCCCGGTTTCGGCATCCCCCTGGCCTTGATCGGCACCTTGGCGCTGCTCTCGGCACTCC         <  1:849751/150‑1                                                                                                                                              CCGAAGCCTTCTTGCCCAGTTTCGGCGTGGTCGGCTTCGGCGGCATCGTTGCCTTTGTGGTCGGCGCATTGATCCTGATCGACACCGACGCTCCCGGTTTCGGCATCCCCCTGGCCTTGATCGGCACCTTGGCGCTGCTCTCGGCACTCC         >  2:455144/1‑150                                                                                                                                                 AAGCCTTCTTGCCCAGTTTCGGCGTGGTCGGCTTCGGCGCCATCGTTGCCCTTGTGGTTGGAGCATTGATGCTGATCGACACCGCCGCTCCCGGGTTCGGCATCCCCCTGGCCTTGATCGGCACCTTGGCGCTGCTCTCGGCACTCCTGA      >  1:929425/1‑150                                                                                                                                                     CTTCTTGCCCAGTTTCGGCGTGGTCGGCTTCGGCGGCATCGTTGCCTTTGTGGTCGGCGCATTGATCCTGATCGACACCGACGCTCCCGGTTTCGGCATCCCCCTGGCCTTGATCGGCACCTTGGCGCTGCTCTCGGCACTCCTGATCGG  >  1:24549/1‑150                                                                                                                                                                                                                                                                                                            TTGAGAGCCATGCCCAGCACGCCACCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/250088‑250236 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgaagccTTCTTGCCCAGTTTCGGCGTGGTCGGCTTCGGCGGCATCGTTGCCTTTGTGGTCGGCGCATTGATCCTGATCGACACCGACGCTCCCGGTTTCGGCATCCCCCTGGCCTTGATCGGCACCTTGGCGCTGCTCTCGGCACTCCTGATCGG  <  NC_002947/250260‑250112

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