New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 2958779	15 (0.470)	5 (0.160)	5/278	3.6	24.8%	intergenic (+45/+97)	PP_2589/PP_2590	aldehyde dehydrogenase family protein/ferric siderophore receptor
	?	NC_002947	2958841 =	16 (0.520)					intergenic (+107/+35)	PP_2589/PP_2590	aldehyde dehydrogenase family protein/ferric siderophore receptor
Rejected: Coverage evenness skew score above cutoff.

AAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTGTCCCCCAGGTTGGCGATGTTCGC  >  NC_002947/2958705‑2958986                                                                                                                                         |                                                                                                                                                  aaGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGCGGCCGTa                                                                                                                                        >  5:125487/1‑148 (MQ=255)                cAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCCGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggcg                                                                                                                         >  1:62032/1‑148 (MQ=255)                           ttgttGTTTTTATCATTGACCGGTTCTTCTGCTTCGCGGGCTTGCCCGCTCCCACTGTTTCATCACATGTCTGTGGTTTGTGTTGGACCTTTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTa                                                                                                              <  6:61289/144‑1 (MQ=25)                                gTATTCAGCATTGCCCGCGTCGGCTGCTGCGCGGGCATGCCCGCGCCCACAGGTTCAGCACAAGTCGGAGGGTGGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGt                                                                                                         <  3:75585/148‑1 (MQ=38)                                               cgcgTCGTCTGCTTCGCGGGCATGCCCGCTCCCACCGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGaaa                                                                                          <  3:45849/148‑1 (MQ=255)                                                        tGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAgg                                                                                 >  5:155873/1‑148 (MQ=255)                                                                gggCATGCCCGCGCCCCCCGGGCCCGCACTAGTAGTAGGGTTGCGCTGGCCCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGTTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCgcggc                                                                         <  3:138244/148‑1 (MQ=37)                                                                            tCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGCATCTGGCGCCGCGGCGGGCCCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAAaccc                                                             >  3:32789/1‑148 (MQ=255)                                                                               cacaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACGTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg                                                          <  4:88235/148‑1 (MQ=255)                                                                                    gTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTg                                                     >  1:133927/1‑148 (MQ=255)                                                                                         gCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCa                                                >  4:138871/1‑148 (MQ=255)                                                                                          cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAg                                               >  7:89509/1‑148 (MQ=255)                                                                                                    aGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACAcc                                      >  3:78551/1‑147 (MQ=255)                                                                                                         tgtgCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCagta                                <  8:89509/148‑1 (MQ=255)                                                                                                                           gCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTGTCCCCCAGGtt              >  5:54531/1‑148 (MQ=255)                                                                                                                                cgcgcCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTGTCCCCCAGGTTGGCGa         <  2:62032/148‑1 (MQ=255)                                                                                                                                       gcgAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTGTCCCCCAGGTTGGCGATGTtcgc  >  1:35974/1‑148 (MQ=255)                                                                                                                                         |                                                                                                                                                  AAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTGTCCCCCAGGTTGGCGATGTTCGC  >  NC_002947/2958705‑2958986

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