New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3974772 =	10 (0.580)	3 (0.180)	3/268	1.9	23.4%	intergenic (‑256/‑85)	PP_3503/PP_3504	sigma‑54 dependent transcriptional regulator/hypothetical protein
	?	NC_002947	3974818 =	10 (0.600)					intergenic (‑302/‑39)	PP_3503/PP_3504	sigma‑54 dependent transcriptional regulator/hypothetical protein

GCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCGTGTCAACGAAGCACTGCTGATTGCCGACCACGCCTTCGAACCTTTCCAGTGCGTAGCCTGGGATGCACCCAACGGTACTGGTGAACTGAGCCT  >  NC_002947/3974681‑3974966                                                                                                                                          |                                                                                                                                                     aggcgtggtcggcaatcagcagtgcttcgttgacacggagcggggcattcatggttcggtctctctatggtcccaaacatcgaagtgaagtCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGa                                                                                                                                                      <  1:389125‑M1/47‑1 (MQ=255) aggcgtggtcggcaatcagcagtgcttcgttgacacggagcggggcattcatggttcggtctctctatggtcccaaacatcgaagtgaagtCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGa                                                                                                                                                      >  2:389125‑M1/92‑138 (MQ=255)                         aCCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGaa                                                                                                                  >  2:9913/1‑150 (MQ=255)                             tGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCAt                                                                                                              <  1:56870/150‑1 (MQ=255)                               aaaaaaaTTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGa                                                                                                            <  1:371777/150‑1 (MQ=255)                                aaaaaaTTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGaa                                                                                                           <  1:273352/150‑1 (MQ=255)                                             ccAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCgt                                                                                              <  1:9913/150‑1 (MQ=255)                                                               aCGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCGTGTCAACGAAGCa                                                                                  >  1:84683/1‑144 (MQ=255)                                                               aCGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCGTGTCAACGAAGCa                                                                                  <  2:84683/144‑1 (MQ=255)                                                               aCGAAAAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCGTGTCAACGAAGCActgctg                                                                            >  1:102664/1‑150 (MQ=255)                                                                        gAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAGTCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCGTGTCAACGAAGCACTGCTGATTGCCGAc                                                                   <  1:109840/150‑1 (MQ=255)                                                                                           tccgcgttgcttcaaagtcgcctcattgatcgcgaagtgcagcgcggACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCGTGTCAACGAAGCACTGCTGATTGCCGACCACGCCt                                                            <  1:389125‑M2/91‑1 (MQ=255)                                                                                           tccgcgttgcttcaaagtcgcctcattgatcgcgaagtgcagcgcggACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCGTGTCAACGAAGCACTGCTGATTGCCGACCACGCCt                                                            >  2:389125‑M2/48‑138 (MQ=255)                                                                                             cgcgttgcttcaaagtcgcctcattgatcgcgaagtgcagcgcggACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCGTGTCAACGAAGCACTGCTGATTGCCGACCACGCCTTCGAACCTTTCCAg                                              >  2:63406‑M2/46‑150 (MQ=255)                                                                                                                           tttGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCGTGTCAACGAAGCACTGCTGATTGCCGACCACGCCTTCGAACCTTTCCAGTGCGTAGCCTGCGCAGCCCGCAACGGAACt                >  2:116482/1‑150 (MQ=255)                                                                                                                                         gACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCGTGTCAACGAAGCACTGCTGATTGCCGACCAAACATTTCACGTTGGGCACTGCGTAGCCTGGGATGCACCCAACGGTACTGGTGAACTGAGCCt  <  2:102664‑M2/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                     GCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCATGAATGCCCCGCTCCGTGTCAACGAAGCACTGCTGATTGCCGACCACGCCTTCGAACCTTTCCAGTGCGTAGCCTGGGATGCACCCAACGGTACTGGTGAACTGAGCCT  >  NC_002947/3974681‑3974966

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