New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4860181	24 (0.850)	3 (0.120)	3/224	NT	11.1%	coding (914/2331 nt)	ligA	DNA ligase
	?	NC_002947	= 4860205	26 (1.010)					coding (890/2331 nt)	ligA	DNA ligase

ACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGCCCGCCGCATGTTTGAGCTCGCGGCTGATCGGC  >  NC_002947/4860058‑4860348                                                                                                                                                    |                                                                                                                                                aCCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTcg                                                                                                                                               <  2:147461‑M2/150‑3 (MQ=255)   cTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTAAGGTCGTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTcgcg                                                                                                                                             <  1:46443‑M2/150‑5 (MQ=255)    ttGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTcgcgc                                                                                                                                            >  2:382725‑M2/1‑145 (MQ=255)    ttGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTcgcgc                                                                                                                                            <  1:546135‑M2/150‑6 (MQ=255)         cGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTcgcgctgggc                                                                                                                                       >  2:622189‑M2/1‑140 (MQ=255)                     cgcCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGCGGCTCTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGtt                                                                                                                           >  1:13403/1‑150 (MQ=255)                         acaGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGAcc                                                                                                                       <  2:439216/150‑1 (MQ=255)                            ggCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTg                                                                                                                    <  2:397043/150‑1 (MQ=255)                                                      cTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACgg                                                                                          <  1:147461/150‑1 (MQ=255)                                                         gAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAg                                                                                       >  1:193735/1‑150 (MQ=255)                                                                    ccAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTAcgcc                                                                            >  1:305042/1‑150 (MQ=255)                                                                          ccTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGa                                                                                                                         >  2:3591/1‑99 (MQ=255)                                                                          ccTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGa                                                                                                                         <  1:3591/99‑1 (MQ=255)                                                                               gTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATg                                                                 >  2:394174/1‑150 (MQ=255)                                                                                            tGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAgg                                                    >  1:194442/1‑150 (MQ=255)                                                                                               ccGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCa                                                 >  2:444880/1‑150 (MQ=255)                                                                                                cGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAg                                                <  2:320076/150‑1 (MQ=255)                                                                                                       ttGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGa                                                                                                                   >  1:450611/1‑76 (MQ=255)                                                                                                       ttGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGa                                                                                                                   <  2:450611/76‑1 (MQ=255)                                                                                                       ttGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCa                                                                                        <  1:520363/103‑1 (MQ=255)                                                                                                       ttGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCa                                                                                        >  2:520363/1‑103 (MQ=255)                                                                                                             gCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACACCGTAAATCTCATACAGAAGGCTGGTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGagt                                   >  1:189605/1‑148 (MQ=255)                                                                                                              cAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGc                                  <  2:46443/150‑1 (MQ=255)                                                                                                               aaTCGCCCAGCGCGGCTGACGGGCGCGGAAGCCCAGCTcgcgctgggcgattgcccacaagttcccggccatggaagagctgaccgaggtgctggatgtcgagttccaggttggccgtacgggtgcagaaacgcctgtggcgcggctgaa                                 <  2:386077‑M2/150‑113 (MQ=255)                                                                                                                     ccAGCGCGGCTCACGGGCGCGGAAGCCCAGCTcgcgctgggcgattgcccacaagttcccggccatggaagagctgaccgaggtgctggatgtcgagttccag                                                                          <  1:94496‑M2/103‑72 (MQ=255)                                                                                                                     ccAGCGCGGCTCACGGGCGCGGAAGCCCAGCTcgcgctgggcgattgcccacaagttcccggccatggaagagctgaccgaggtgctggatgtcgagttccag                                                                          >  2:94496‑M2/1‑32 (MQ=255)                                                                                                                       aGCGCGACTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGCCCGCCGCAt                         >  2:165952/1‑150 (MQ=255)                                                                                                                           cgGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGCCCGCCGCATGttt                     <  2:305042/150‑1 (MQ=255)                                                                                                                              cTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGCCCGCCGCATGTTTGAg                  >  2:718978/1‑150 (MQ=255)                                                                                                                                 aCGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATa                                                                                             >  2:137861/1‑72 (MQ=255)                                                                                                                                 aCGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATa                                                                                             <  1:137861/72‑1 (MQ=255)                                                                                                                                     gcgcGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGCCCGCCGCATGTTTGAGCTCGCgg           <  2:492666/150‑1 (MQ=255)                                                                                                                                           ccgCCCAGCGCGAGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGc                                  >  2:432377/3‑121 (MQ=255)                                                                                                                                           ccgCCCAGCGCGAGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGc                                  <  1:432377/119‑1 (MQ=255)                                                                                                                                              cctaGCTCGCGCTGGGCGGCCAGGCTGTTCACCTTGACTACTACCCCGTCGATCTCATACGGCAGTCTGTTACGCCGCGCGCCAATGTCTCGATAGTAGGCCGGGCCCTCTTCGATGCCCGCCGCATGTTTTCGCTCGCGGCTGATCTGc  <  2:100016/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                ACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGCCCGCCGCATGTTTGAGCTCGCGGCTGATCGGC  >  NC_002947/4860058‑4860348

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