New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5596092	6 (0.360)	7 (0.520)	5/230	1.0	61.3%	intergenic (‑66/+79)	PP_4921/thiC	NCS1 family transporter/phosphomethylpyrimidine synthase
	?	NC_002947	= 5596129	4 (0.300)					intergenic (‑103/+42)	PP_4921/thiC	NCS1 family transporter/phosphomethylpyrimidine synthase

GGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCGTAC  >  NC_002947/5595985‑5596276                                                                                                                                                 |                                                                                                                                                    ggggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCCCGAACACCGGCAAAGCCGGTGCCATTCACCCCAGTGCTTTCttcgc                                                                                                                                                <  1:120560‑M2/150‑6 (MQ=255)   ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgg                                                                                                                                              >  2:228120‑M2/1‑143 (MQ=255)         gATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGTCGGTGCCATGCACCGCAGTGCTTGCttc                                                                                                                                                  <  1:291190‑M2/140‑4 (MQ=255)         gATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGTCGGTGCCATGCACCGCAGTGCTTGCttc                                                                                                                                                  >  2:291190‑M2/1‑137 (MQ=255)                                                                                         ggAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactccagcggacatgaccatgacatcacccagccagttctctccc                                                        <  2:245455‑M2/150‑94 (MQ=255)                                                                                                     cccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaatc                                            <  2:260565‑M2/150‑106 (MQ=255)                                                                                                      ccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggaaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaatcc                                           <  2:110243‑M2/150‑107 (MQ=255)                                                                                                      ccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCtt                                           >  1:274275/1‑150 (MQ=255)                                                                                                       cgcgAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTg                                          <  2:264843/150‑1 (MQ=255)                                                                                                          gAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaatccccat                                       <  2:339026‑M2/150‑111 (MQ=255)                                                                                                          gAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaatccccat                                       <  2:339056‑M2/150‑111 (MQ=255)                                                                                                              acCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGAcc                                   >  2:182401/1‑150 (MQ=255)                                                                                                                                           tgcttgctTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGc      >  2:333612/1‑150 (MQ=255)                                                                                                                                               tgcttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaatccccatccagcagcgcattttcggtgctcgcgacctgttttca  >  1:339026‑M2/1‑3 (MQ=255)                                                                                                                                               tgcttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaatccccatccagcagcgcattttcggtgctcgcgacctgttttca  >  1:339056‑M2/1‑3 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCGTAC  >  NC_002947/5595985‑5596276

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