New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4253573 =	14 (0.550)	4 (0.160)	4/264	2.7	24.6%	coding (150/1407 nt)	rocE	amino acid permease RocE
	?	NC_002947	4253600 =	11 (0.450)					coding (123/1407 nt)	rocE	amino acid permease RocE

GTACACCCAGCCGATCATGAACCCGGTGGCCGGGCCAATGTATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTC  >  NC_002947/4253426‑4253716                                                                                                                                                    |                                                                                                                                                gTACACCCAGCCGATCATGAACCCGGTGGCCGGGCCAATGTATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAg                                                                                                                                               <  1:310442/150‑1 (MQ=255) ccggccgcaatggcgaggtgtgccatgactgacagcgttgaaaaaatccagctcacccgcgccctgaaaagccggcacatcttcatgctgtcgctgggtggggttattggaaccggcctgttcatgggctcgggcgtaaccatcaatCAg                                                                                                                                               >  1:400077‑M1/148‑150 (MQ=255)            cGATCATGAACCCGGTGGCCGGGCCAATGTATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCa                                                                                                                                    >  1:88525/1‑150 (MQ=255)                     aCCCGGTGGCCGGGCCAATGTATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGgcccgc                                                                                                                           >  1:99258/1‑150 (MQ=255)                         ggTGGCCGGGCCAATGTATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCtg                                                                                                                       >  1:66246/1‑150 (MQ=255)                              ccGGGCCAATGTATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCtgattga                                                                                                                  <  2:170699/150‑1 (MQ=255)                                         tATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGccc                                                                                                       <  2:99258/150‑1 (MQ=255)                                             aatccagctcacccgcgccctgaaaagccggcacatcttcatgctgtcgctgggtggggttattggaaccggcctgttcatgggctcgggcgtaaccatcaatCAGGATCGCGCCCACCGGCCCGCCCtgattgat                                                                                                                 >  1:593275‑M1/104‑136 (MQ=255)                                             aatccagctcacccgcgccctgaaaagccggcacatcttcatgctgtcgctgggtggggttattggaaccggcctgttcatgggctcgggcgtaaccatcaatCAGGATCGCGCCCACCGGCCCGCCCtgattgat                                                                                                                 <  2:593275‑M1/33‑1 (MQ=255)                                                    gTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGaa                                                                                            <  2:88525/150‑1 (MQ=255)                                                               cctgaaaagccggcacatcttcatgctgtcgctgggtggggttattggaaccggcctgttcatgggctcgggcgtaaccatcaatCAGGATCGCGCCCACCGGCCCGCCCTGATTGATgg                                                                                                               <  1:29674‑M1/35‑1 (MQ=255)                                                               cctgaaaagccggcacatcttcatgctgtcgctgggtggggttattggaaccggcctgttcatgggctcgggcgtaaccatcaatCAGGATCGCGCCCACCGGCCCGCCCTGATTGATgg                                                                                                               >  2:29674‑M1/86‑120 (MQ=255)                                                                  aaCCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCaa                                                                              <  1:117640/150‑1 (MQ=255)                                                                                   tGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGAc                                                             >  1:78857/1‑150 (MQ=255)                                                                                        cGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCAt                                                        >  2:340451/1‑150 (MQ=255)                                                                                            aGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACACGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAg                                                    >  1:34506/1‑150 (MQ=255)                                                                                                           aCCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGc                                              >  2:313335/1‑141 (MQ=255)                                                                                                           aCCATCACCAGGTACATCACCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGc                                              <  1:313335/141‑1 (MQ=255)                                                                                                                                              gTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCAGGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTgtc  >  2:396611/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                GTACACCCAGCCGATCATGAACCCGGTGGCCGGGCCAATGTATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTC  >  NC_002947/4253426‑4253716

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.