New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 4040689	32 (1.020)	6 (0.190)	6/284	3.2	18.9%	intergenic (+53/+217)	rrfA/mobB	5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
	?	NC_000913_3_bme_pgi	= 4213156	20 (0.650)					noncoding (117/120 nt)	rrfE	5S ribosomal RNA of rrnE operon
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCA  >  NC_000913_3_bme_pgi/4040542‑4040799                                                                                                                                                    |                                                                                                               cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTGTTATCtgttgt                                                                                                               >  1:517325/1‑149 (MQ=21)         aCCCCATGCCGGACTGAGAATTTAAACGCCATAGCGGAGATTGTAGTATTGTGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                                                       <  2:57504/149‑1 (MQ=11)           cccATGCCGAACTCAGAAGTGAAACGCCGTCGAGCCGGGAATAGGTTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGCCTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGa                                                                                                     >  1:404305/1‑149 (MQ=11)                   gAACTCAGAAGGGAAACGCCGTAGCGCAGATGGTAGTGTGGGGTCTCCGCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCGCAGTCGGAAGACTGGGCCTTTCGTTCTATCTGTTGTTTGTCGGTGAACGCtctc                                                                                             <  2:150196/149‑1 (MQ=11)                    aaCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGACCGCTCTcc                                                                                            >  1:27513/1‑149 (MQ=21)                     aCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCt                                                                                           <  1:115355/149‑1 (MQ=35)                       tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGa                                                                                         <  1:389973/149‑1 (MQ=35)                       caaGAAGTGAAACGCCGTAGCGCCGATGGGAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGa                                                                                         <  2:383553/147‑1 (MQ=21)                        cAGAAGTCGAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCACGCGAGAGTAGGGAACTGCCGGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                                                        <  2:10889/149‑1 (MQ=255)                         aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGCACGCTCTCCTGAGt                                                                                       >  2:285794/1‑149 (MQ=21)                           aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGCACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATATGTTGTTTGTCGGTGCACGCTCTCATGAGtgg                                                                                     >  1:522451/1‑147 (MQ=255)                                         aGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCCTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCgg                                                                       <  1:63871/149‑1 (MQ=21)                                            gcCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAg                                                                    >  2:511880/1‑149 (MQ=35)                                                aTGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGa                                                                <  1:327582/149‑1 (MQ=35)                                                      gtgtGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAACACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGaa                                                          >  1:416184/1‑149 (MQ=25)                                                       tgtgGGGTCTCCTCATGCGAGGGTAGGGAACAGGAAGGAAACTAATAAAACGAAAAGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTAGTCGGTGAACGCTCTCCGGAGTCAGACAAATTCGCCGGGGGCGGATTTGAAc                                                         >  2:366147/1‑149 (MQ=11)                                                                  ctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                         >  2:111886/1‑138 (MQ=34)                                                                  cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                         <  1:111886/136‑1 (MQ=14)                                                                   ccgtGCGAGGTTAGGGAACTGCCAGGAATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCaa                                             <  2:260201/146‑1 (MQ=255)                                                                     aTGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACg                                           >  1:76312/1‑149 (MQ=32)                                                                     aTGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACg                                           <  2:4808/149‑1 (MQ=32)                                                                         gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                       >  1:109966/1‑149 (MQ=32)                                                                             gTAGGGAACTGCCAGGCATCAAAGAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAACTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAg                                   <  1:297774/149‑1 (MQ=17)                                                                                      tGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACg                                           >  1:53152/1‑132 (MQ=33)                                                                                      tGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACg                                           <  2:53152/132‑1 (MQ=33)                                                                                         cAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAg                       <  2:109966/149‑1 (MQ=32)                                                                                             cATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcg                   >  1:373230/1‑149 (MQ=32)                                                                                               tCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGGGCGGATTTGAACGTTGCGAAGCAACGGCCCGGCGGGTGGCGGGCAGGAcgcc                 >  1:154140/1‑149 (MQ=17)                                                                                                 aaaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccg               >  2:360455/1‑149 (MQ=32)                                                                                                       aaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaa         >  1:306466/1‑149 (MQ=32)                                                                                                         cGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAAc       >  2:314565/1‑149 (MQ=32)                                                                                                              ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCCCCCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGCCACCCCGCCATACACCGCCa  >  2:420428/1‑149 (MQ=11)                                                                                                                                                    |                                                                                                               CCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCA  >  NC_000913_3_bme_pgi/4040542‑4040799

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