New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4159083 =	20 (0.370)	7 (0.150)	7/232	NT	35.4%	intergenic (+31/‑178)	PP_3661/PP_3662	membrane protein/hypothetical protein
	?	NC_002947	4159143 =	8 (0.170)					intergenic (+91/‑118)	PP_3661/PP_3662	membrane protein/hypothetical protein

CCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAA  >  NC_002947/4158941‑4159226                                                                                                                                               |                                                                                                                                                ccAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGc                                                                                                                                           >  1:166327/1‑149 (MQ=255)   aGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTc                                                                                                                                         >  3:3925/1‑149 (MQ=255)      ttCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                      >  1:22320/1‑149 (MQ=255)       gacatggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctgcttcgcgggtgAACCCGCTCCTa                                                                                                                                      <  2:195328‑M1/12‑1 (MQ=255)          atggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctgcttcgcgggtgAACCCGCTCCTACAgg                                                                                                                                  >  1:195328‑M1/134‑149 (MQ=255)          atggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctgcttcgcgggtaAACCCGCTCCTACAgg                                                                                                                                  <  6:304222‑M1/16‑1 (MQ=255)                    ccACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGtgctgc                                                                                                                        >  5:194505/1‑149 (MQ=255)                     gcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctgcttcgcgggtgAACCCGCTCCTAc                                                                                                                                     <  1:65892‑M1/13‑1 (MQ=255)                     gcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctgcttcgcgggtgAACCCGCTCCTAc                                                                                                                                     >  2:65892‑M1/123‑135 (MQ=255)                     gcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctgcttcgcgggtaAACCCGCTCCTACAGGTGGGTGCTGcc                                                                                                                       <  3:21869‑M1/27‑1 (MQ=255)                       cGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCt                                                                                                                     <  2:103420/149‑1 (MQ=255)                       cGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCt                                                                                                                     <  6:254933/149‑1 (MQ=255)                        gggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCt                                                                                                                     >  8:5429/1‑148 (MQ=255)                        gggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTg                                                                                                                    <  2:22320/149‑1 (MQ=255)                          gcgcTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGaa                                                                                                                  <  3:129214/149‑1 (MQ=255)                          gcgcTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGaa                                                                                                                  <  2:166327/149‑1 (MQ=255)                          cccgtGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGaa                                                                                                                  <  2:99265/145‑1 (MQ=255)                                  cACCGGTGCGATCCTCACCCTGGGGCTGGGTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCaa                                                                                                          <  6:194505/149‑1 (MQ=255)                                  cACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCaa                                                                                                          <  5:270693/149‑1 (MQ=255)                                  cACCGGTGCGATCCTCAACCTGTGGCTGGGTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCaa                                                                                                          <  7:5429/149‑1 (MQ=255)                                     cGGTGCGATCCTCAACCTGTGGCTGGTTGTCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGc                                                                                                       <  1:237936/149‑1 (MQ=255)                                             taaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctgcttcgcgggtaAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAg                                                                                               >  4:21869‑M1/99‑149 (MQ=255)                                              ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGTGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGTGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAg                                                                                               <  7:303409/148‑1 (MQ=255)                                               cTCCACCTGGGGCTGGTCGGCCGGGGGCGGGCGATTCCCCTGGGGGCGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAgg                                                                                              <  3:215285/148‑1 (MQ=38)                                                cgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgactggagaccgagtcggctgcttcgcgggtaAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                                             >  7:15341‑M1/96‑148 (MQ=255)                                                 acatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctgcttcgcgggtgAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGc                                                                                           >  3:119810‑M1/95‑149 (MQ=255)                                                   atccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctgcttcgcgggtaAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCgg                                                                                         >  5:304222‑M1/93‑149 (MQ=255)                                                    ccTGGGGCTGGGTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCGGTTCGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                                        <  1:150424/149‑1 (MQ=255)                                                                          cgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctgcttcgcgggtgAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                                        >  5:11332‑M1/70‑127 (MQ=255)                                                                          cgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctgcttcgcgggtgAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                                        <  6:11332‑M1/58‑1 (MQ=255)                                                                                                                                  tgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtggt          >  6:106902/1‑149 (MQ=255)                                                                                                                                     ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATggtggt       >  6:87728/1‑149 (MQ=255)                                                                                                                                          gggTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCaa  <  5:106902/149‑1 (MQ=255)                                                                                                                                           ggTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCa   >  1:316636/1‑147 (MQ=255)                                                                                                                                            gTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACAATTCTCGGTATGGTGGTGTCa   >  5:188577/1‑146 (MQ=255)                                                                                                                                               |                                                                                                                                                CCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAA  >  NC_002947/4158941‑4159226

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