New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 176399	NA (NA)	15 (0.400) +CGAAA	9/256	3.1	25.3%	intergenic (+43/‑420)	PP_5SA/PP_16SB	5S ribosomal RNA/16S ribosomal RNA
	?	NC_002947	= 182013	46 (1.170)					intergenic (+2/+63)	PP_5SB/PP_0160	5S ribosomal RNA/ferrioxamine receptor
Rejected: Coverage evenness skew score above cutoff.

CACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCACCTCCCGCGACGAGAGATCGAAGCGAGTTAA  >  NC_002947/176264‑176529                                                                                                                                        |                                                                                                                                   caccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAg                                                                                                                      <  1:175710/150‑1 (MQ=35)          ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCtt                                                                                                             >  1:43178/1‑150 (MQ=255)                     tCAGTAGTGAAACGATGCATCGCAGACCGGAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTaa                                                                                                  <  1:495222/150‑1 (MQ=255)                               aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGt                                                                                        <  2:220391/150‑1 (MQ=255)                                       aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGt                                                                                >  2:907701/1‑150 (MQ=255)                                          gCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                             >  2:561346/1‑102 (MQ=31)                                          gCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                             <  1:561346/102‑1 (MQ=31)                                           ccGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCa                                                                            >  1:825320/1‑150 (MQ=255)                                             gATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGt                                                                          >  1:66380/1‑150 (MQ=255)                                                            tttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGAc                                                           >  1:55456/1‑150 (MQ=255)                                                                cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAgtg                                                       >  2:379990/1‑150 (MQ=255)                                                                      tgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTg                                                                                                      <  1:559549/97‑1 (MQ=255)                                                                      tgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTg                                                                                                      >  2:559549/1‑97 (MQ=255)                                                                                tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGtat                                       <  2:66380/150‑1 (MQ=255)                                                                                   tCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGtattat                                    <  1:379990/150‑1 (MQ=255)                                                                                     gTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTc                                  <  1:907701/150‑1 (MQ=255)                                                                                      tCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCa                                 <  2:55456/150‑1 (MQ=255)                                                                                         aGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACa                                                          >  2:178413/1‑122 (MQ=255)                                                                                         aGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACa                                                          <  1:178413/122‑1 (MQ=255)                                                                                           aTTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCACCTcc                            >  1:891777/1‑150 (MQ=255)                                                                                                 ttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAggtgg                                                                >  2:26938/1‑108 (MQ=255)                                                                                                 ttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAggtgg                                                                <  1:26938/108‑1 (MQ=255)                                                                                                            ccTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCACCTCCCGCGACGAGAGATCGaa           <  2:43178/150‑1 (MQ=255)                                                                                                                     gcATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCACCTCCCGCGACGAGAGATCGAAGCGAGTTaa  >  2:869435/1‑150 (MQ=255)                                                                                                                             gTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTg                                          <  2:235694/102‑1 (MQ=255)                                                                                                                             gTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTg                                          >  1:235694/1‑102 (MQ=255)                                                                                                                                        |                                                                                                                                   CACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCACCTCCCGCGACGAGAGATCGAAGCGAGTTAA  >  NC_002947/176264‑176529

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