New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4530003 =	42 (1.260)	5 (0.160)	5/250	3.3	11.7%	coding (241/1314 nt)	sgcC	putative PTS system EIIC permease component
	?	NC_000913	4530040 =	35 (1.090)					coding (204/1314 nt)	sgcC	putative PTS system EIIC permease component
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCGTTTTGGTCAGACGGGTGACCAGCATGGCGACATTGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGA  >  NC_000913/4529901‑4530188                                                                                                                                            |                                                                                                                                                     tCGTTTTGGTCAGACGGGTGACCAGCATGGCGACATTGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTc                                                                                                                                              >  6:112780/1‑148 (MQ=255)        ggTCAGACGGGTGACCAGCATGGCGACATTGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGa                                                                                                                                                     <  6:25576/134‑1 (MQ=255)        ggTCAGACGGGTGACCAGCATGGCGACATTGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGa                                                                                                                                                     >  5:25576/1‑134 (MQ=255)          tCAGACGGGTGACCAGCATGGCGACATTGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACt                                                                                                                                    >  1:120674/1‑149 (MQ=255)                      ccAGCATGGCGACATTGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGa                                                                                                                        >  1:163538/1‑149 (MQ=255)                            tGGCGACATTGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAgg                                                                                                                  <  2:120674/149‑1 (MQ=255)                                gACATTGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGt                                                                                                              >  2:8925/1‑149 (MQ=255)                                      gagaAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATc                                                                                                         <  5:27603/148‑1 (MQ=255)                                      gagaAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATc                                                                                                         <  1:8925/148‑1 (MQ=255)                                         aagaaAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAACCGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGcc                                                                                                      >  6:8345/1‑148 (MQ=255)                                         aagaaAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCaa                                                                                                                                            >  2:25845/1‑110 (MQ=255)                                         aagaaAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCa                                                                                                     >  1:88017/1‑149 (MQ=255)                                         aagaaAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCa                                                                                                     >  3:58122/1‑149 (MQ=255)                                           gaaAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATg                                                                                                   <  2:163538/149‑1 (MQ=255)                                                               aTGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATAc                                                                                <  3:5293/148‑1 (MQ=255)                                                               aTGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATAc                                                                                >  4:5293/1‑148 (MQ=255)                                                                tGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCAc                                                                                           >  5:130320/1‑136 (MQ=255)                                                                tGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCAc                                                                                           <  6:130320/136‑1 (MQ=255)                                                                                   aTAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCAc                                                                                           >  2:11095/1‑117 (MQ=255)                                                                                   aTAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCAc                                                                                           <  1:11095/117‑1 (MQ=255)                                                                                       ccAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGa                                                                                               >  6:95751/1‑109 (MQ=255)                                                                                       ccAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGa                                                                                               <  5:95751/109‑1 (MQ=255)                                                                                         aaCGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATc                                                                                             >  5:39860/1‑109 (MQ=255)                                                                                         aaCGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATc                                                                                             <  6:39860/109‑1 (MQ=255)                                                                                                   gCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAAt                                           <  5:8345/149‑1 (MQ=255)                                                                                                     aGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCg                                         >  1:43061/1‑149 (MQ=255)                                                                                                                cAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAg                              >  3:27233/1‑149 (MQ=255)                                                                                                                        aaGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCa                                                                                         <  3:86616/82‑1 (MQ=255)                                                                                                                        aaGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCa                                                                                         >  4:86616/1‑82 (MQ=255)                                                                                                                        aaGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCg                      <  2:88017/149‑1 (MQ=255)                                                                                                                             aTGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTc                                                                                                             <  4:55726/57‑1 (MQ=255)                                                                                                                             aTGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTc                                                                                                             >  3:55726/1‑57 (MQ=255)                                                                                                                               gCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAt                >  3:73395/1‑148 (MQ=255)                                                                                                                                     cTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGc           <  8:10453/147‑1 (MQ=255)                                                                                                                                     cTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCa         >  7:10453/1‑149 (MQ=255)                                                                                                                                          ccGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCAc    >  1:154250/1‑149 (MQ=255)                                                                                                                                           cGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCAc    >  3:35336/1‑148 (MQ=255)                                                                                                                                            gAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGa  <  4:27233/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                     TCGTTTTGGTCAGACGGGTGACCAGCATGGCGACATTGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGA  >  NC_000913/4529901‑4530188

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