New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 250473	60 (1.190)	4 (0.080)	4/246	NT	7.1%	coding (730/1347 nt)	nfeD	membrane bound peptidase NefD
	?	NC_002947	= 250522	48 (1.000)					coding (681/1347 nt)	nfeD	membrane bound peptidase NefD

CGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAATTCGAACAACAGGCCGTAAACGCCGATCATGATCAGGATCAGCGCTACGCTGGGGTTGGTGATGACCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/250321‑250473 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaccgcAGGCGCCAGCGTGGTCGAGCATCTGCCGGACTGGCGCACACGGTTTCTGGCGGTCATCA  <  NC_002947/250522‑250464                                                                                                                                                                                                                      CGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAATTCGAACAACAGGCCGTAAACGCCGATCATGATCAGGATCAGCGCTACGCTGGGGTTGGTGATGA                                                                 >  2:137481/1‑149     AACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAATTCGAACAACAGGCCGTAAACGCCGATCATGATCAGGATCAGCGCTACGCTGGGGTTGGTGATGACCGC                                                             <  6:251580/149‑1                                     GCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAATTCGAACAACAGGCCGTAAACGCCGATCATGATCAGGATCAGCGCTACGCTGGGGTTGGTGATGACCGCAGGCGCCAGCGTGGTCGAGCATCTGCCGGACT                             >  8:45439/1‑149                                                  ACCACACCGCCCACCGCAGAACCTGGGTTCATGAATTCGAACAACAGGCCGTAAACGCCGATCATGATCAGGATCAGCGCTACGCTGGGGTTGGTGATGACCGCAGGCGCCAGCGTGGTCGAGCATCTGCCGGACTGGCGCACACGGTT                >  2:22748/1‑149                                                               CCGCAGAACCTGGGTTCATGAATTCGAACAACAGGCCGTAAACGCCGATCATGATCAGGATCAGCGCTACGCTGGGGTTGGTGATGACCGCAGGCGCCAGCGTGGTCGAGCATCTGCCGGACTGGCGCACACGGTTTCTGGCGGTCATC   <  7:45439/149‑1                                                                CGCAGAACCTGGGTTCATGAATTCGAACAACAGGCCGTAAACGCCGATCATGATCAGGATCAGCGCTACGCTGGGGTTGGTGATGACCGCAGGCGCCAGCGTGGTCGAGCATCTGCCGGACTGGCGCACACGGTTTCTGGCGGTCATCA  <  1:22748/149‑1                                                                                                                                                                                                                      CGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAATTCGAACAACAGGCCGTAAACGCCGATCATGATCAGGATCAGCGCTACGCTGGGGTTGGTGATGACCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/250321‑250473 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaccgcAGGCGCCAGCGTGGTCGAGCATCTGCCGGACTGGCGCACACGGTTTCTGGCGGTCATCA  <  NC_002947/250522‑250464

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