New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2720039	36 (1.060)	4 (0.130)	4/252	3.7	11.2%	coding (87/2661 nt)	pka	protein lysine acetyltransferase
	?	NC_000913	= 2720049	30 (0.950)					coding (97/2661 nt)	pka	protein lysine acetyltransferase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAA  >  NC_000913/2719901‑2720185                                                                                                                                                     |                                                                                                                                         gAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGc                                                                                                                                          <  2:134344/149‑1 (MQ=255)      gTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAAc                                                                                                                                     <  4:7227/149‑1 (MQ=255)          gAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGTAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCTTCGATGAAACCCAATCTCGCAGGTTACCTGATGATGCGTAACatgc                                                                                                                                 >  6:17157/1‑149 (MQ=255)               cGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCg                                                                                                                            >  8:27451/1‑149 (MQ=255)                gTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCgg                                                                                                                           >  5:38535/1‑149 (MQ=255)                        tCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCAAAAACGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGTTTACCAGATGATGCGTAACCTCCTGGCGGGAGGCtt                                                                                                                    >  6:29905/1‑148 (MQ=255)                                                   gTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGCCGc                                                                                          >  7:93230/1‑147 (MQ=255)                                                        aGTCAGCGAGGACTGGAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTg                                                                                   >  5:122933/1‑149 (MQ=255)                                                               gAGGACTGGAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCg                                                                            >  2:109400/1‑149 (MQ=255)                                                                aGGACTGGAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCgg                                                                           <  6:38535/149‑1 (MQ=255)                                                                      ggAAGCACTACTGCGACCAAAATCGATAGAGGTAATTGGCGCGTCGATGAAACCCAATCGCTCAGGTTACCTGATGATTCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCAGGTGACGCCAGCCTGGAAAGCGGTGTTgg                                                                     <  4:18876/149‑1 (MQ=255)                                                                       gAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTggg                                                                    >  2:158453/1‑149 (MQ=255)                                                                       gAAGCACTACTGAGACCAAAATCGATAGCGTTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGTCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTgg                                                                     <  5:17157/148‑1 (MQ=255)                                                                          gCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGgt                                                                   >  3:15462/1‑147 (MQ=255)                                                                          gCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGgt                                                                   <  4:15462/147‑1 (MQ=255)                                                                            actactGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGgtgtgt                                                               >  1:140618/1‑149 (MQ=255)                                                                                ctGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTtggc                                                           <  7:27451/149‑1 (MQ=255)                                                                                       cAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCg                                                    >  5:68513/1‑149 (MQ=255)                                                                                         aaaTCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGTCTTTAACGGACCGGTACTCCCGGTGACGc                                                                                          <  2:114586/109‑1 (MQ=255)                                                                                         aaaTCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGc                                                                                          >  1:114586/1‑109 (MQ=255)                                                                                           aTCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGata                                                >  1:115448/1‑149 (MQ=255)                                                                                               aTAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGc                                            >  5:65334/1‑149 (MQ=255)                                                                                                    ggTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGgt                                                                          <  6:101983/114‑1 (MQ=255)                                                                                                    ggTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGgt                                                                          >  5:101983/1‑114 (MQ=255)                                                                                                    ggTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCt                                       <  6:68513/149‑1 (MQ=255)                                                                                                        aTTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGcc                                   <  8:93230/149‑1 (MQ=255)                                                                                                             cgcgTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTa                              <  2:115448/149‑1 (MQ=255)                                                                                                                    aTGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTa                              <  5:75882/142‑1 (MQ=255)                                                                                                                    aTGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTa                              >  6:75882/1‑142 (MQ=255)                                                                                                                               tCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGttt            <  6:122933/149‑1 (MQ=255)                                                                                                                                 gcgcAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCt                                       <  8:15167/120‑1 (MQ=255)                                                                                                                                    cAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGt       <  1:109400/149‑1 (MQ=255)                                                                                                                                          aCCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGGGTGTTAGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCaa  <  1:158453/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                         GAACAGTTAGAAAGCGTTTAAAATCATTCGGTCACTTCTGCGGGAGACCGGTATGAGTCAGCGAGGACTGGAAGCACTACTGCGACCAAAATCGATAGCGGTAATTGGCGCGTCGATGAAACCCAATCGCGCAGGTTACCTGATGATGCGTAACCTGCTGGCGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAA  >  NC_000913/2719901‑2720185

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