New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3946824	37 (0.280)	22 (0.170)	12/256	10.1	23.9%	intergenic (+5/‑48)	rrfC/aspT	5S ribosomal RNA of rrnC operon/tRNA‑Asp
	?	NC_000913	= 4040681	104 (0.800)					intergenic (+45/+225)	rrfA/mobB	5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
Rejected: Coverage evenness skew score above cutoff.

AAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAA  >  NC_000913/3946677‑3946971                                                                                                                                                    |                                                                                                                                                    aaGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAAt                                                                                                                                                     >  1:494772/1‑148 (MQ=32)  aGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAAtt                                                                                                                                                    >  8:456892/1‑148 (MQ=17)    cGGTCTGATGAAACAGAATTTGCCTGGCTGCCGTAGCGCGGGGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAg                                                                                                                                                 <  3:262801/149‑1 (MQ=17)    cGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGGGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAg                                                                                                                                                 <  5:423291/149‑1 (MQ=17)      gTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCa                                                                                                                                               <  5:204931/149‑1 (MQ=32)      gTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCa                                                                                                                                               <  6:394284/149‑1 (MQ=32)       tCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAg                                                                                                                                              >  6:11028/1‑149 (MQ=32)           atgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTaa                                                                                                                                           <  7:310955/148‑1 (MQ=32)            tgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCa                                                                                                                                               >  2:163839/1‑143 (MQ=31)            tgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCa                                                                                                                                               <  1:163839/143‑1 (MQ=31)                                                                                                                                gggAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCAccc                      >  8:407226/1‑148 (MQ=31)                                                                                                                                gggAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCAccc                      >  8:171653/1‑148 (MQ=31)                                                                                                                                 ggAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGATCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCt                     >  6:163544/1‑148 (MQ=17)                                                                                                                                 ggAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCt                     >  8:459950/1‑148 (MQ=32)                                                                                                                                  gAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTaa                   >  4:124239/1‑149 (MQ=33)                                                                                                                                  gAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTaa                   >  1:478650/1‑149 (MQ=33)                                                                                                                                        ccAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAggg             <  7:129734/149‑1 (MQ=38)                                                                                                                                        ccAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAggg             <  7:621905/149‑1 (MQ=38)                                                                                                                                        ccAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAggg             >  6:181860/1‑149 (MQ=38)                                                                                                                                        ccAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAggg             >  5:517947/1‑149 (MQ=38)                                                                                                                                        ccAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAggg             <  4:145563/149‑1 (MQ=38)                                                                                                                                        ccAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAgg              <  8:229612/148‑1 (MQ=38)                                                                                                                                         cAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAgggg            >  4:341804/1‑149 (MQ=39)                                                                                                                                         cAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAgggg            <  8:256145/149‑1 (MQ=39)                                                                                                                                             cATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTa        <  2:601626/149‑1 (MQ=255)                                                                                                                                             cATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTa        <  7:199555/149‑1 (MQ=255)                                                                                                                                              aTCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAg       >  7:518020/1‑149 (MQ=255)                                                                                                                                              aTCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAg       >  4:138217/1‑149 (MQ=255)                                                                                                                                               tCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGt      >  8:339864/1‑149 (MQ=255)                                                                                                                                               tCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCCTCCGTTCCGCCACCCTAATTAGGTGCGTAGt      >  7:497128/1‑149 (MQ=255)                                                                                                                                                 aaaTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGtt     <  2:306391/148‑1 (MQ=255)                                                                                                                                                 aaaTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGtt     >  4:251654/1‑148 (MQ=255)                                                                                                                                                 aaaTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGtt     <  3:251654/148‑1 (MQ=255)                                                                                                                                                 aaaTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCAGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGtt     >  8:93432/1‑148 (MQ=255)                                                                                                                                                  aaTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTc    >  8:259489/1‑148 (MQ=255)                                                                                                                                                   aTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCaa  <  4:327737/149‑1 (MQ=255)                                                                                                                                                   aTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCaa  >  6:542915/1‑149 (MQ=255)                                                                                                                                                   aTTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCaa  >  6:313762/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAA  >  NC_000913/3946677‑3946971

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