New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4171896	43 (0.960)	6 (0.140)	6/266	4.7	12.6%	intergenic (+140/‑161)	rrfB/murB	5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding
	?	NC_000913	= 4171901	43 (1.020)					intergenic (+145/‑156)	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.

CAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGC  >  NC_000913/4171750‑4172041                                                                                                                                                    |                                                                                                                                                  cAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAg                                                                                                                                                 <  1:259503/149‑1 (MQ=32) cAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAg                                                                                                                                                 <  1:763317/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)                               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                                                                   <  2:582067/149‑1 (MQ=38)                                                                               ccTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACaaa                                                                   <  1:531267/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)                                                                                                                                              ggCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCCACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGtttt     <  1:92957/149‑1 (MQ=255)                                                                                                                                               gCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCCCAGATACGGTAAACTAGCCTCGttttt    >  1:698338/1‑149 (MQ=255)                                                                                                                                                 gggCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGc  <  1:651722/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG‑CAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCGTTTTTGC  >  NC_000913/4171750‑4172041

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.