New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5441919 =	22 (0.810)	4 (0.160)	4/280	2.9	17.4%	coding (637/837 nt)	PP_4775	hypothetical protein
	?	NC_002947	5441973 =	17 (0.660)					coding (583/837 nt)	PP_4775	hypothetical protein
Rejected: Frequency below/above cutoff threshold.

CCGACTGGCTCACCGCCCAGCGATGGCACTGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAA  >  NC_002947/5441826‑5442121                                                                                                                                                    |                                                                                                                                                     cgcgggaagtgtggttgggtggtgcgctgttctttctggtggccgatacctattaccagtacctggcagacagccaggaatatgtggtggtcaACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGa                                                                                                                                                    >  2:581681‑M1/94‑150 (MQ=255) ccGACTGGCTCACCGCCCAGCGATGGCACTGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACCTCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGa                                                                                                                                                    <  2:412469/150‑1 (MQ=255)                cccAGCGATGGCACTGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCt                                                                                                                                     <  1:312786/150‑1 (MQ=255)                         cgctgttctttctggtggccgatacctattaccagtacctggcagacagccaggaatatgtggtggtcaACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTg                                                                                                                                                     >  2:474993‑M1/70‑125 (MQ=255)                         cgctgttctttctggtggccgatacctattaccagtacctggcagacagccaggaatatgtggtggtcaACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTg                                                                                                                                                     <  1:474993‑M1/56‑1 (MQ=255)                             cTGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCAGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAg                                                                                                                        <  1:109786/150‑1 (MQ=255)                             cTGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCAGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAg                                                                                                                        <  1:63761/150‑1 (MQ=255)                              tGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAgg                                                                                                                       <  2:378086/150‑1 (MQ=255)                                   gCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTg                                                                                                                  <  2:407219/150‑1 (MQ=255)                                       tcctgctcctgcTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTaa                                                                                                              >  1:515190/1‑150 (MQ=255)                                       tcctgctcctgcTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTaa                                                                                                              >  1:396698/1‑150 (MQ=255)                                           gctcctgctGTTCCTGACCCCGGCCATACATCGGATAAGCCACCACGTCATCCCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCCCATATTCCTGGCTGTCTGCCGGGTACTGGTAATAgg                                                                                                          >  2:45166/1‑150 (MQ=255)                                                tgctgTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCg                                                                                                     >  2:560134/1‑150 (MQ=255)                                                tgctgTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCg                                                                                                     >  2:252714/1‑150 (MQ=255)                                                       ccTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGccacca                                                                                              >  1:68414/1‑150 (MQ=255)                                                         tGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCaga                                                                                            >  2:360011/1‑150 (MQ=255)                                                         tGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCaga                                                                                            >  1:120119/1‑150 (MQ=255)                                                                        aTCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGcacca                                                                             >  1:237397/1‑150 (MQ=255)                                                                                           tcaacccgccggtcggcgtgccgtcgcccgtgccggtggccccggcaaatggcgggtTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCcgcg                                                          >  2:581681‑M2/58‑150 (MQ=255)                                                                                            caacccgccggtcggcgtgccgtcgcccgtgccggtggccccggcaaatggcgggtTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAgcg                                                                                  <  1:474993‑M2/69‑1 (MQ=255)                                                                                            caacccgccggtcggcgtgccgtcgcccgtgccggtggccccggcaaatggcgggtTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAgcg                                                                                  >  2:474993‑M2/57‑125 (MQ=255)                                                                                              acccgccggtcggcgtgccgtcgcccgtgccggtggccccggcaaatggcgggtTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGt                                                       >  1:118316‑M2/55‑150 (MQ=255)                                                                                                   CCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGc                                                  >  2:341457‑M1/1‑150 (MQ=255)                                                                                                   CCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGc                                                  <  2:298266‑M1/150‑1 (MQ=255)                                                                                                                          gggCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGcc                           >  1:367301/1‑150 (MQ=255)                                                                                                                          gggCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACACGCACGcc                           >  1:571413/1‑150 (MQ=255)                                                                                                                                               cgggtTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGAt      <  2:192827‑M2/145‑1 (MQ=255)                                                                                                                                               cgggtTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGAt      <  2:192820‑M2/145‑1 (MQ=255)                                                                                                                                                   tTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTaa  >  1:535982‑M2/2‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CCGACTGGCTCACCGCCCAGCGATGGCACTGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAA  >  NC_002947/5441826‑5442121

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