New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1849289 =	21 (0.620)	5 (0.160)	5/228	3.4	21.9%	coding (453/1857 nt)	sppA	protease IV (signal peptide peptidase)
	?	NC_000913	1849327 =	16 (0.500)					coding (491/1857 nt)	sppA	protease IV (signal peptide peptidase)
Rejected: Coverage evenness skew score above cutoff.

GACGACCGCAATATCACCGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAA  >  NC_000913/1849161‑1849427                                                                                                                                 |                                                                                                                                           gacgacCGCAATATCACCGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAggg                                                                                                                        >  2:12047/1‑149 (MQ=255)             atCACCGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGGTTATGCCGTTGGCGAGAACTACAGCCACGGGCAATATTAtct                                                                                                            >  8:141999/1‑149 (MQ=255)              tCACCGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGtt                                                                                                                                            >  2:166957/1‑116 (MQ=255)                 ccGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCgaga                                                                                                                                     <  1:155541/120‑1 (MQ=255)                 ccGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCgaga                                                                                                                                     >  2:155541/1‑120 (MQ=255)                    gTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCa                                                                                                       >  6:95301/1‑147 (MQ=255)                    gTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCa                                                                                                       <  5:95301/147‑1 (MQ=255)                          tgATGGAGCTGAAAAACTTCGCAGGCGGGGACCAACCGTCTATGCAGTACATCGGCAATGCTCTGAAAGAGTTTCGTGTCAGCGGGAAACCGGTTTATGCCGttt                                                                                                                                           <  4:152440/105‑2 (MQ=255)                          tgATCGATCTGAAAAACTTCCCAGCCCGCGACCAACCGTCTAAGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGCATTATGCCGTTg                                                                                                                                           >  3:152440/1‑105 (MQ=255)                                     aaaaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCaa                                                                                             <  1:109579/140‑1 (MQ=255)                                     aaaaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACAGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCaa                                                                                             >  2:109579/1‑140 (MQ=255)                                                                     gCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTg                                                                                    >  5:26480/1‑117 (MQ=255)                                                                     gCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTg                                                                                    <  6:26480/117‑1 (MQ=255)                                                                                  aaaGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATaaaa                                                                                        >  2:53839/1‑100 (MQ=255)                                                                                  aaaGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATaaaa                                                                                        <  1:53839/100‑1 (MQ=255)                                                                                               agTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGtactac                          <  1:12047/149‑1 (MQ=255)                                                                                                               gAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGc                                                                                                                      <  2:103695/41‑1 (MQ=255)                                                                                                               gAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGc                                                                                                                      >  1:103695/1‑41 (MQ=255)                                                                                                                 aaCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCttt                                               >  3:63823/1‑110 (MQ=255)                                                                                                                 aaCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCttt                                               <  4:63823/110‑1 (MQ=255)                                                                                                                       tttATGCAGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCATGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGaaa  <  7:141999/149‑1 (MQ=255)                                                                                                                                 |                                                                                                                                           GACGACCGCAATATCACCGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAA  >  NC_000913/1849161‑1849427

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