New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 857864	35 (1.410)	6 (0.250)	6/260	2.2	16.6%	coding (963/1443 nt)	phr	deoxyribodipyrimidine photolyase
	?	NC_002947	= 857920	26 (1.060)					coding (907/1443 nt)	phr	deoxyribodipyrimidine photolyase
Rejected: Frequency below/above cutoff threshold.

GCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCACTGTCGAACTCGCCCCGGTTGCCGGC  >  NC_002947/857783‑858052                                                                                                                                          |                                                                                                                                     gcgcAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGcgcg                                                                                                                          <  1:581423/150‑1 (MQ=255)   gcAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGt                                                                                                                        <  2:525725/150‑1 (MQ=255)   gcAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGt                                                                                                                        >  1:525725/1‑150 (MQ=255)                 aTCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCtgt                                                                                                          <  2:162849/150‑1 (MQ=255)                 aTCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCtgt                                                                                                          >  1:565255/1‑150 (MQ=255)                       ccGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGgcgc                                                                                                                   >  1:148487/1‑135 (MQ=255)                       ccGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGgcgc                                                                                                                   <  2:148487/135‑1 (MQ=255)                        cGGTGTGCAGCAATTGGCGCATGGCGGCTTCAATGCTCGGCATGCCGGTGCTGCCCTGCTCCCAGGCCGCCAGATCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAAc                                                                                                   <  1:141123/150‑1 (MQ=255)                              gcagcaATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCa                                                                                             <  1:381924/150‑1 (MQ=255)                               cagcaATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCga                                                                                                                 <  2:586343/129‑1 (MQ=255)                               cagcaATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCga                                                                                                                 >  1:586343/1‑129 (MQ=255)                                              ggcggcATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGaa                                                                             <  1:400467/150‑1 (MQ=255)                                                           aTCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGt                                                                >  2:286316/1‑150 (MQ=255)                                                               ggAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCgg                                                                                                                         >  1:43114/1‑89 (MQ=255)                                                               ggAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCgg                                                                                                                         <  2:43114/89‑1 (MQ=255)                                                                  aGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCAt                                                         >  1:39787/1‑150 (MQ=255)                                                                       gTGCGGCCCTGCTCCCAAGCCGCCAGGTCTGCCGAGGCGTCACGCCAGTGCAGGGCTTCGGTCTGGGCGCATCAGGCGCTGTGGCGCAAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATc                                                    <  2:256785/150‑1 (MQ=255)                                                                          cGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAg                                                 <  2:595268/150‑1 (MQ=255)                                                                            gCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGAgcagggccgcaccggcttcccgatcattgatgccgccatgcgccaattgctg                                                                                  >  2:177439‑M2/1‑63 (MQ=255)                                                                            gCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGAgcagggccgcaccggcttcccgatcattgatgccgccatgcgccaattgctg                                                                                  <  1:177439‑M2/115‑53 (MQ=255)                                                                              ccTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCgcagggccgcaccggcttcccgatcattgatgccgccatgcgccaattgctgcacaccgggtggatgaacaaccgcttgccgatggtcg                                             <  2:207223‑M2/150‑90 (MQ=255)                                                                               cTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCgcagggccgcaccggcttcccgatcattgatgccgccatgcgccaattgctgcacaccgggtggatgcacaaccgcttgcg                                                     >  1:156995‑M2/1‑60 (MQ=255)                                                                               cTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCgcagggccgcaccggcttcccgatcattgatgccgccatgcgccaattgctgcacaccgggtggatgcacaaccgcttgcg                                                     <  2:156995‑M2/141‑82 (MQ=255)                                                                               cTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCtgtg                                                                                                         <  2:39658/89‑1 (MQ=255)                                                                               cTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCtgtg                                                                                                         >  1:39658/1‑89 (MQ=255)                                                                                     ccAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCgcagggccgcaccggcttcccgatcattgatgccgccatgcgccaattgctgcacaccgggtggatgcacaaccgcttgcgcatggtcgtggccat                                      >  1:207223‑M2/1‑54 (MQ=255)                                                                                      cAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCt                                     <  2:39787/150‑1 (MQ=255)                                                                                          ccgccAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGtgc                                 >  2:601504/1‑150 (MQ=255)                                                                                                 gTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCtgt                                                                                                          >  1:622800/1‑70 (MQ=255)                                                                                                 gTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCtgt                                                                                                          <  2:622800/70‑1 (MQ=255)                                                                                                                      gggcagggcTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTACACGCCTGTGCTGCCACTGTCGAACTCGCCCCGGTTGcc     >  2:305414/1‑150 (MQ=255)                                                                                                                         cagggcTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCACTGTCGAACTCGCCCCGGTTGCCGGc  >  1:467655/1‑150 (MQ=255)                                                                                                                                          |                                                                                                                                     GCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCACTGTCGAACTCGCCCCGGTTGCCGGC  >  NC_002947/857783‑858052

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