New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3946823	14 (0.210)	10 (0.150)	5/280	9.2	17.3%	intergenic (+5/‑48)	rrfC/aspT	5S ribosomal RNA of rrnC operon/tRNA‑Asp
	?	NC_000913	= 4040680	82 (1.230)					intergenic (+45/+225)	rrfA/mobB	5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAG  >  NC_000913/3946681‑3946965                                                                                                                                                |                                                                                                                                               gTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAAt                                                                                                                                                <  1:147442/143‑1 (MQ=32)  tCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGCAGTG‑TGGGGGCTCCCCATGCGAGAGGAGGGAGCTGCCAGGCACCAAAtt                                                                                                                                               >  2:773942/1‑143 (MQ=11)  tCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTAAGCAGAGACTCGCGGCAGAGATG‑GGATTCTGTTGGCGTCTCCCCATGATAGGGTAGGGAACGGCCAGGCATCAAAtt                                                                                                                                               >  2:247985/1‑143 (MQ=2)     gatgaACCAGAGTTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGATCCCATCCCGAACTCAGAAGTGAAACGCCGGAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAg                                                                                                                                            <  2:1017012/143‑1 (MQ=11)     gatgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAg                                                                                                                                            <  1:586107/143‑1 (MQ=32)                                                                                                                                aCTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCAcc                  <  1:826237/143‑1 (MQ=11)                                                                                                                                    ccAGGCATCAAATTAAGCAGTACGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTACAATACCTGCCTGTAACCACACAACTCGTGTGGCCGAGTCCCGTCCGTTCCGCCACCCTaa              >  2:242223/1‑143 (MQ=11)                                                                                                                                    ccAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCCGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTGCCGTCCGGTGCGCCCCCCtca              >  1:65196/1‑141 (MQ=11)                                                                                                                                       ggCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTa           <  1:77009/143‑1 (MQ=35)                                                                                                                                          aTCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAggg        <  2:598026/143‑1 (MQ=38)                                                                                                                                             aaaTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGCGGGGGGTCGGAGTCCCG‑CCGGTTCGCCAcctacttgggaggctc    >  2:991704/1‑129 (MQ=11)                                                                                                                                              gaTTACGCGGTAAGCCGGTCATACAACCGGTGGTTGTAAAAGAAGTCGGTGGAGCGGTAGTTCAGTCGGTTAGACTACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGCGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGt    <  1:676067/142‑1 (MQ=14)                                                                                                                                              aaTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCCCCCTAATCACGGGcct    >  1:328867/1‑141 (MQ=255)                                                                                                                                               aTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAACTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTa   <  1:916122/143‑1 (MQ=255)                                                                                                                                                ttAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCAAGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAg  >  1:205595/1‑143 (MQ=255)                                                                                                                                                |                                                                                                                                               GTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAG  >  NC_000913/3946681‑3946965

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