New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	1154932 =	20 (0.740)	6 (0.240)	5/270	2.4	29.7%	intergenic (‑224/+33)	PP_1008/gapA	RNA polymerase subunit sigma‑24/glyceraldehyde‑3‑phosphate dehydrogenase
	?	NC_002947	= 2076208	10 (0.400)					intergenic (+168/+59)	PP_1846/PP_1847	group II intron‑encoding maturase/TonB‑dependent ferric siderophore receptor

CGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCCACATCGCGCCA  >  NC_002947/2076066‑2076356                                                                                                                                               |                                                                                                                                                     cGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgg                                                                                                                                               >  1:497854‑M2/1‑143 (MQ=255) cGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgg                                                                                                                                               >  2:539582‑M2/1‑143 (MQ=255)                                          tCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCgg                                                                                                      >  2:229535/1‑150 (MQ=255)                                            gaaattaaGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTAGACCGAGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGACCGTGTTAGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGcc                                                                                                    >  2:465648/8‑150 (MQ=25)                                                            ggggCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGc                                                                                    <  2:178648/150‑1 (MQ=255)                                                                    tgtgCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaa                                                                            >  2:225641/1‑150 (MQ=255)                                                                                   gACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCg                                                             >  2:400695/1‑150 (MQ=255)                                                                                               cTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGC‑‑‑‑ggcaaggccgctcccacacgggtcgcgtcagccagcccggtgcagcgccaggcagttgtccagcatccggttgctgaacccccactcgttgtcataccaggccag                                              <  2:92212‑M2/150‑106 (MQ=255)                                                                                                ttGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgggcaaggccgctcccacacgggtcgcgtcagccagcccggtgcagcgccaggcagttgtccagcatccggttgctgaacccccactcgttgtcata                                                <  2:597373‑M2/150‑103 (MQ=255)                                                                                                ttGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgggcaaggccgctcccacacgggtcgcgtcagccagcccggtgcagcgccaggcagttgtccagcatccggttgctgaacccccactcgttgtcata                                                <  2:597369‑M2/150‑103 (MQ=255)                                                                                                    ccGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTg                                            >  2:336748/1‑150 (MQ=255)                                                                                                     cgcgAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgggcaaggccgctcccacacgggtcgcgtcagccagcccggtgcagcgccaggcagttgtccagcatccggttgctgaacccccactcgttgtcataccagg                                           <  2:157279‑M2/150‑108 (MQ=255)                                                                                                        gAAGCCAACACCGTGGTGGATGGCATGGGCTTCGCCCGTGttcgcgggcaagcccgctcccccagggatcgcgccagccagcccggtgcagcgccaggcagttgtccagcatccggttgctgaacccccactcgttgtcataccaggcca                                        <  2:74892‑M2/150‑111 (MQ=255)                                                                                                                 acCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgggcaaggccgctcccacacgggtcgcgtcagccagcccggtgcagcgccaggcagttgtccagcatccggttgctgaacccccactcgttgtcataccaggccagcaccttga                               >  1:157279‑M2/1‑31 (MQ=255)                                                                                                                   cGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATc                             <  2:323755/150‑1 (MQ=255)                                                                                                                                          cccGTGTTGGGGGGCAAGCCCGCTCCCACAGGTCGCTAAAGGGCCTCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCCACATcgc      <  1:465648/150‑1 (MQ=255)                                                                                                                                           ccGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCCACATcgcg     >  2:194868/1‑150 (MQ=255)                                                                                                                                              tgtTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCCACATCGCGCCa  <  2:32195/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     CGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCCACATCGCGCCA  >  NC_002947/2076066‑2076356

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.