New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	3645087 =	17 (0.930)	4 (0.230)	4/274	1.9	19.9%	coding (96/2043 nt)	prlC	oligopeptidase A
	?	NC_000913_3_vch_pgi	3645120 =	16 (0.920)					coding (63/2043 nt)	prlC	oligopeptidase A
Rejected: Frequency below/above cutoff threshold.

GGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGC  >  NC_000913_3_vch_pgi/3644945‑3645215                                                                                                                                               |                                                                                                                                 ggCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCTGCTGGCAGAGCTTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTc                                                                                                                            >  1:288033/1‑149 (MQ=255) ggCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTc                                                                                                                            >  2:210843/1‑149 (MQ=255)   cTATTTTGCACCGATGTCCGGGGGCTGATCGGGGAGAAGATGCGACCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCaa                                                                                                                          <  2:103499/149‑1 (MQ=255)      tctttCACCGAGTTCAGGTGGCCGACCGGGGAGAAGATTCTTCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGc                                                                                                                       <  2:285532/147‑1 (MQ=255)            aCCGAGTTCAGGTGGCTGACCGGGGAGACGATACGCCCCAACACATCGTCAACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCACTACTACGCGCCCCACCTTTTCGCGGCAATCGTTCCATCCCCTatc                                                                                                                 >  2:104176/1‑147 (MQ=255)              cGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCa                                                                                                               <  1:210843/149‑1 (MQ=255)                aGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCCATGCCTTAGTCACg                                                                                                             >  1:179431/1‑149 (MQ=255)                  ttCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGc                                                                                                           >  1:214727/1‑149 (MQ=255)                                cGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATggt                                                                                             >  2:188822/1‑147 (MQ=255)                                       aagaTACACCCCAACACATCGTCCACTACCGCCAGCCGCCGGCACAGATAGTCCCAGGTGTACGGTGCCCCATACGCTACGACGCGCTCCACAATCTCCCTGCAGCCCGTCAACGCCTTAGCCCCGGCTGAACCGACAAGTCCCggaaa                                                                                      >  2:101805/1‑146 (MQ=255)                                                       caTCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGa                                                                      <  2:55415/149‑1 (MQ=255)                                                          ctctttaaccaggactgcgcgaatgacgaatccgttactgactccctttgaattgcctccgttttctaaaattctcccggaacatGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGga                                                                                       >  1:21726‑M1/86‑129 (MQ=255)                                                          ctctttaaccaggactgcgcgaatgacgaatccgttactgactccctttgaattgcctccgttttctaaaattctcccggaacatGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGga                                                                                       <  2:21726‑M1/44‑1 (MQ=255)                                                                                            ccAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTCTAAAAAACGGAGGCAAATCAAAGGGAGTCAGTAACGGCTTCGTCATTc                                 >  1:43444/1‑149 (MQ=255)                                                                                                       gTGCCCCTTGCGCTAGTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACTGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTg                      <  2:258154/149‑1 (MQ=255)                                                                                                                 cgcTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAgg            >  2:146155/1‑149 (MQ=255)                                                                                                                  gcTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGt           >  1:306730/1‑149 (MQ=255)                                                                                                                           cgcTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAAGTTAAGGAAgc  <  1:15261/149‑1 (MQ=255)                                                                                                                           cgcTCCACGTTTGCGCGGCAGTCGTTCAATCCCTTAGTCACGGCTGGAACGATATGTGCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACTGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAgc  <  1:14365/149‑1 (MQ=255)                                                                                                                                              tGTCGTTCAATGCCTTAGTCACGGCTGGAACGACa                                                                                                 >  1:48128‑M1/2‑35 (MQ=255)                                                                                                                                              tGTCGTTCAATGCCTTAGTCACGGCTGGAACGACa                                                                                                 <  2:48128‑M1/34‑1 (MQ=255)                                                                                                                                               |                                                                                                                                 GGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGC  >  NC_000913_3_vch_pgi/3644945‑3645215

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