New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4378781	30 (0.900)	6 (0.190)	5/278	3.1	17.4%	intergenic (‑69/‑20)	PP_3854/PP_5639	lysozyme/hypothetical protein
	?	NC_002947	= 4378793	28 (0.870)					intergenic (‑81/‑8)	PP_3854/PP_5639	lysozyme/hypothetical protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTCGGACCAGGCCAGCATGTCTAGGAAGGCAATCACGTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCGTGTCGAGCGAATCAAAGGCCGTGGGCGACCTGGCGGAATGGAATTCCTAC  >  NC_002947/4378638‑4378917                                                                                                                                                |                                                                                                                                         ttCGGACCAGGCCAGCATGTCTAGGAAGGCAATCACGTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATc                                                                                                                                        >  7:8125/1‑146 (MQ=255)     gACCAGGCCAGCATGTCTAGGAAGGCAATCACGTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACagag                                                                                                                                  <  5:14869/148‑1 (MQ=255)             cAGCATGTCTAGGAAGGCAATCACGTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGTCCATCACAGAGGGAGTGGc                                                                                                                          >  7:167606/1‑148 (MQ=255)              aGCATGTCTAGGAAGGCAATCACGTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAACAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTCTGATCTGGCCATCACAGAGGGAGTGGcc                                                                                                                         <  7:57973/148‑1 (MQ=255)                cATGTCTAGGAAGGCAATCACGTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCtg                                                                                                                       <  4:162358/148‑1 (MQ=255)                  tGTCTAGGAAGGCAATCACGTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCtgat                                                                                                                     <  8:2898/148‑1 (MQ=255)                         gAAGGCAATCACGTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAag                                                                                                                 <  2:174485/145‑1 (MQ=255)                             gCAATCACGTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTat                                                                                                          <  3:5109/148‑1 (MQ=255)                                    cGTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTatggatgg                                                                                                    <  5:58804/147‑1 (MQ=255)                                     gTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGa                                                                                                  <  1:141785/148‑1 (MQ=255)                                     gTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGa                                                                                                  <  4:143001/148‑1 (MQ=255)                                       ttACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGt                                                                                                <  1:18536/148‑1 (MQ=255)                                           gccgccGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGAc                                                                                             <  6:169468/147‑1 (MQ=255)                                             cgccgGCATCTGCGGCAGAAATACGGGCCATGTTTTTCTCCAGGCAAAAAGAAGACAGATCTAGGAGGGCCTGATTATTTAATACGCTATGATCTGTCCATCACAGAGGGATTGGCATGATTAAGAAGCTATGGATGGGAGTGACGtt                                                                                          >  8:141917/1‑148 (MQ=255)                                                cgGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTgg                                                                                        >  7:71608/1‑147 (MQ=255)                                                    aTCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATg                                                                                   <  2:6530/148‑1 (MQ=255)                                                                  tACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTg                                                                     <  1:37988/148‑1 (MQ=255)                                                                    cGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCg                                                                   >  4:98021/1‑148 (MQ=255)                                                                         cATGGTTTTCTCCAGGCAAAAAGAAGCCGGATCTAGGCGGGCCTGTTTATTTAATACGCTATGATGTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTAAGGAGGGGAGTGACGTGGGTCATGGCGCTAGGCGGTTGCGGGCCg                                                              >  2:150155/1‑148 (MQ=255)                                                                             gTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTg                                                           >  4:56567/1‑147 (MQ=255)                                                                                        gCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCgtgt                                                 >  3:152164/1‑146 (MQ=255)                                                                                          aaaaaGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCGTGTCGAg                                             >  3:82912/1‑148 (MQ=255)                                                                                                                ggCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCGTGTCGAGCGAATCAAAGGCCGTGGGCGAc                       <  5:154397/148‑1 (MQ=255)                                                                                                                   cTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCGTGTCGAGCGAATCAAAGGCCGTGGGCGACCTg                    >  6:94276/1‑148 (MQ=255)                                                                                                                    tGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCGTGTCGAGCGAATCAAAGGCCGTGGGCGACCTgg                   <  3:98021/148‑1 (MQ=255)                                                                                                                        tattTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCGTGTCGAGCGAATCAAAGGCCGTGGGCGACCTGGCgga               <  4:172734/148‑1 (MQ=255)                                                                                                                               tACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCGTGTCGAGCGAATCAAAGGCCGTGGGCGACCTGGCggaatggaa         <  3:56567/147‑1 (MQ=255)                                                                                                                                   cTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCGTGTCGAGCGAATCAAAGGCCGTGGGCGACCTGGCGGAATGGAATTCCt    >  4:118390/1‑148 (MQ=255)                                                                                                                                    tATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCGTGTCGAGCGAATCAAAGGCCGTGGGCGACCTGGCGGAATGGAATTCCTa   <  5:25592/148‑1 (MQ=255)                                                                                                                                     aTGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCGTGTCGAGCGAATCAAAGGCCGTGGGCGACCTGGCGGAATGGAATTCCTAc  >  7:139469/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                         TTCGGACCAGGCCAGCATGTCTAGGAAGGCAATCACGTTTACGCCGCCGGCATCTGCGGCAGAAATACGGGCCATGGTTTTCTCCAGGCAAAAAGAAGCCCGATCTAGGCGGGCCTGATTATTTAATACGCTATGATCTGGCCATCACAGAGGGAGTGGCCTGATGAAGCAGCTATGGATGGGAGTGACGTTGGTCATGGCGCTAGGCGGTTGCGGGACGGTGCGAACCGTGTCGAGCGAATCAAAGGCCGTGGGCGACCTGGCGGAATGGAATTCCTAC  >  NC_002947/4378638‑4378917

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.