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.

TGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCT  >  NC_002947/857721‑858009                                                                                                                                                |                                                                                                                                                  tGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAg                                                                                                                                             <  2:616568/150‑1 (MQ=255)             ttGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCt                                                                                                                                 >  1:256785/1‑150 (MQ=255)                  ccAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCgg                                                                                                                            >  1:348999/1‑150 (MQ=255)                   cAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCggg                                                                                                                           >  2:381924/1‑150 (MQ=255)                         aTCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCa                                                                                                                     >  1:162849/1‑150 (MQ=255)                          tCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCAc                                                                                                                    <  2:348999/150‑1 (MQ=255)                          tCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCAc                                                                                                                    <  1:190665/150‑1 (MQ=255)                                           aGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCa                                                                                                                                      <  2:311382/115‑1 (MQ=255)                                           aGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCa                                                                                                                                      >  1:311382/1‑115 (MQ=255)                                           aGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCg                                                                                                   >  1:595268/1‑150 (MQ=255)                                             gAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGt                                                                                                 >  2:400467/1‑150 (MQ=255)                                                 aTGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCgcccagaccgaagccctgccctggcgtgacgccccggcagacctggcggcctgg                                                                                             >  1:207223‑M1/1‑96 (MQ=255)                                                        cGACCATCGGCAAGCGGTTGTTCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCgcccagaccgaagccctgccctggcgtgacgccccggcagacctggcggcctgggagcagg                                                                                      <  2:207223‑M1/150‑62 (MQ=255)                                                               gcgcAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGcgcg                                                                               <  1:581423/150‑1 (MQ=255)                                                                cgcAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCgcccagaccgaagccctgccctggcgtgacgccccggcagacctggcggcctgggagcag                                                                                       >  1:156995‑M1/1‑81 (MQ=255)                                                                cgcAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCgcccagaccgaagccctgccctggcgtgacgccccggcagacctggcggcctgggagcag                                                                                       <  2:156995‑M1/141‑61 (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)                                                                                             cagcaATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCtcccagaccgaagccctgccctggcgtgacgccccggcagacctggcggcctgggagcagggc                                                                                    <  1:177439‑M1/115‑64 (MQ=255)                                                                                             cagcaATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCtcccagaccgaagccctgccctggcgtgacgccccggcagacctggcggcctgggagcagggc                                                                                    >  2:177439‑M1/1‑52 (MQ=255)                                                                                             cagcaATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCga                                                                      >  1:586343/1‑129 (MQ=255)                                                                                             cagcaATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCga                                                                      <  2:586343/129‑1 (MQ=255)                                                                                                caATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTgg                                                                                             <  1:212376/103‑1 (MQ=255)                                                                                                caATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTgg                                                                                             >  2:212376/1‑103 (MQ=255)                                                                                                            ggcggcATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGaa                                  <  1:400467/150‑1 (MQ=255)                                                                                                                         aTCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCgg                                                                                                                            <  2:531340/47‑1 (MQ=255)                                                                                                                         aTCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCgg                                                                                                                            >  1:531340/1‑47 (MQ=255)                                                                                                                         aTCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGt                     >  2:286316/1‑150 (MQ=255)                                                                                                                             ggAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCgg                                                                              <  2:43114/89‑1 (MQ=255)                                                                                                                             ggAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCgg                                                                              >  1:43114/1‑89 (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                                                              >  1:39658/1‑89 (MQ=255)                                                                                                                                             cTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCtgtg                                                              <  2:39658/89‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  TGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCT  >  NC_002947/857721‑858009

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