New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4040789	31 (0.690)	6 (0.140)	6/266	4.7	14.8%	intergenic (+153/+117)	rrfA/mobB	5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
	?	NC_000913	= 4171888	40 (0.950)					intergenic (+132/‑169)	rrfB/murB	5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTC  >  NC_000913/4171741‑4172033                                                                                                                                                     |                                                                                                                                                  gggAACTGCCAGGCATCAAATAAGACGAAAGGCTCAGTCGAAAGCCTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCGCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGt                                                                                                                                                  <  2:633016/149‑1 (MQ=12)  ggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTg                                                                                                                                                 >  1:582067/1‑149 (MQ=32)  ggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTg                                                                                                                                                 <  2:641801/149‑1 (MQ=32)    aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGc                                                                                                                                               >  1:203386/1‑149 (MQ=32)          cAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCggg                                                                                                                                            >  2:542460/1‑146 (MQ=32)          cAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCggg                                                                                                                                            <  1:542460/146‑1 (MQ=32)          cAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAg                                                                                                                                         <  1:763317/149‑1 (MQ=32)          cAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAg                                                                                                                                         <  1:259503/149‑1 (MQ=32)               aTCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CACCGCCCCGCCGAGTGGCGGGCCGGAccc                                                                                                                                    >  1:348063/1‑147 (MQ=11)                tCAAATACAACGAAAGACTCAGTCGAAAGACTGGGCCTTTCGTATTATCTCTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAGAGGCCCGGGGGGTGGCGGGCAGGAcgcc                                                                                                                                   >  2:85059/1‑149 (MQ=11)                 cAAATAAAACGAAAGGCTTAGTCGAAAGACTGGGCCGTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACACT‑CGCAGTGATCGGCCCGGAGGGTGGCGGGCAGGAcgccc                                                                                                                                  <  2:752399/149‑1 (MQ=11)                  aaaTAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAACTCCGCCGGGAGCGGATTTGACCGTTGCGAAG‑CAACGGCCCGGAGGGGGGCGGGCAGGCcgcccg                                                                                                                                 >  2:625795/1‑149 (MQ=255)                           gAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGCACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACt                                                                                                                        >  2:540010/1‑149 (MQ=17)                               ggCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGGTGATTTTAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCa                                                                                                                    <  2:324233/149‑1 (MQ=12)                                gCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAg                                                                                                                   <  2:70575/149‑1 (MQ=32)                                  tCAGTCGAAAGACTGGGCCTTACGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGc                                                                                                                 >  2:268796/1‑149 (MQ=17)                                      tCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGc                                                                                                                                               <  1:234099/115‑1 (MQ=32)                                      tCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGc                                                                                                                                               >  2:234099/1‑115 (MQ=32)                                        gAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaac                                                                                                           >  2:676843/1‑148 (MQ=14)                                         aaaGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAt                                                                                                          >  1:135447/1‑149 (MQ=11)                                                                               gAACGCGCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGt                                                                    <  2:607566/149‑1 (MQ=14)                                                                                 aCGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGttt                                                                  <  2:135447/149‑1 (MQ=32)                                                                                    ctctCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTa                                                               <  2:254694/149‑1 (MQ=35)                                                                                       tcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACaa                                                            >  2:92957/1‑149 (MQ=38)                                                                                        ccTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACaaa                                                           <  2:636829/149‑1 (MQ=38)                                                                                        ccTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACaaa                                                           <  1:531267/149‑1 (MQ=38)                                                                                        ccTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACaaa                                                           <  2:582067/149‑1 (MQ=38)                                                                                              tAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CACCGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTc                                                     >  2:220053/1‑149 (MQ=255)                                                                                                ggACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCt                                                   <  1:268796/149‑1 (MQ=255)                                                                                                  aCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgt                                                 >  2:348179/1‑149 (MQ=255)                                                                                                  aCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgt                                                 >  1:388534/1‑149 (MQ=255)                                                                                                         cgccgGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCata                                          >  2:276523/1‑149 (MQ=255)                                                                                                             gggAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTATCGTCATATCTa                                      >  1:68333/1‑149 (MQ=255)                                                                                                                ggCGGATGTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACCCCGGCCAGTAACTGCCAGGCATCAAATCAAGCAGAAGGCCATCCTGACGGATGGCCTCTTTGCGTTTCTACAACCTCTTCCTGTCGTCATATCTTCaa                                   <  2:718842/148‑1 (MQ=17)                                                                                                                 gCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAg                                  >  1:320564/1‑149 (MQ=255)                                                                                                                       ttGAACGGTGCTAAG‑TAACGGGCGGGGGGGTGGCGGGCAGGACTCTCGCCATAAACTGGCAGGCATCAAATTAAGCAGAAGGGCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATcc                            <  2:82960/149‑1 (MQ=18)                                                                                                                        tGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCGGTCTCCCTCTTCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTGGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccc                           <  1:540010/149‑1 (MQ=18)                                                                                                                         gAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATcccc                          <  1:6930/149‑1 (MQ=255)                                                                                                                          aaCGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccccc                         <  1:276523/149‑1 (MQ=255)                                                                                                                           aCGGTGCGAGG‑CACGGGCCCGGCGGGCGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGAGGGTCTTTTTGCTTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATcccccc                        <  2:714457/149‑1 (MQ=18)                                                                                                                             gTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcac                      >  1:684306/1‑149 (MQ=255)                                                                                                                                   aaG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATAc                >  1:374262/1‑149 (MQ=255)                                                                                                                                       cAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGt             <  1:182903/149‑1 (MQ=255)                                                                                                                                       cAACGGCCCGGAGGGTGCCAGGCAGCAACCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATAGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGt             <  2:480164/149‑1 (MQ=25)                                                                                                                                             cccGGAGGGGGGCGGGCAGGACGCACGCCAACACCTGCGCGGACTCAATTTAAGCACAAGACGATCTTCAGATATGCTCTGTTGGCGGTGCGACATACTCTGCCTGGCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTa       <  2:382328/149‑1 (MQ=255)                                                                                                                                              ccGGAGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCAGCAAATTAAGCAGAAGGCCATCCTGACGGATGGTCTTGTTGCGTTTCTACCAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAg      <  2:3541/149‑1 (MQ=25)                                                                                                                                               cGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGc     <  1:451018/149‑1 (MQ=255)                                                                                                                                                  aGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTc  <  2:60557/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTC  >  NC_000913/4171741‑4172033

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.