New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	171284 =	NA (NA)	3 (0.060)	3/238	NT	NA	intergenic (+465/‑105)	PP_0159/PP_16SA	family III CoA‑transferase/16S ribosomal RNA
	?	NC_002947	171337 =	NA (NA)					intergenic (+518/‑52)	PP_0159/PP_16SA	family III CoA‑transferase/16S ribosomal RNA

GACGAAAAGCTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAGCAGTATTGAACTGAAGAGTTTGATCATGGCTCAGATT  >  NC_002947/171191‑171403                                                                                                                                                   |                                                                   gaCGAAAAGCTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCa                                                                 >  2:227007/1‑150 (MQ=32)    gAAAAGCTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGt                                                              >  2:1116586/1‑150 (MQ=32)      aaaGCTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGttt                                                            >  2:286118/1‑150 (MQ=32)             ttAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGttt                                                     >  1:342071/1‑150 (MQ=32)                           tttAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAg                                       >  1:322199/1‑150 (MQ=32)                           tttAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAg                                       <  2:840521/150‑1 (MQ=32)                           tttAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAg                                       <  2:322199/150‑1 (MQ=32)                            ttAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAGc                                      >  1:886543/1‑150 (MQ=32)                              aaCAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTAAGCTACAAAAAGATTATCAGCATCACAAGTGACCGTGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAGCAg                                    >  1:32860/1‑150 (MQ=11)                               aCAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAGCAGt                                   <  2:1249851/150‑1 (MQ=32)                                 aaaTCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAGCAGTAt                                 >  1:466717/1‑150 (MQ=32)                                     cGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAGCAGTATTGaa                             <  2:466717/150‑1 (MQ=32)                                         tCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAGCAGTATTGAACTGa                         >  1:1177172/1‑150 (MQ=32)                                              cAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAGCAGTATTGAACTGAAGAGt                    <  1:286118/150‑1 (MQ=32)                                                aTTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAg                                                          >  2:843661/1‑110 (MQ=33)                                                aTTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAg                                                          <  1:843661/110‑1 (MQ=33)                                                               tgtgAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAGCAGTATTGAACTGAAGAGTTTGATCATGGCTCAGAt   >  1:1270357/1‑150 (MQ=32)                                                                gtgAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTAATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAGCAGTATTGAACTGAAGAGTTTGATCATGGCTCAGAtt  <  2:512662/150‑1 (MQ=17)                                                                                                                                                   |                                                                   GACGAAAAGCTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTATGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGACCATGCGAGAAATCACATAGTCATTTGAGATTGCTGAGCCAAGTTTAGGGTTTCTTAAAAACCCAAGCAGTATTGAACTGAAGAGTTTGATCATGGCTCAGATT  >  NC_002947/171191‑171403

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