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

CCGGGCCAATGTATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGG  >  NC_002947/4253455‑4253746                                                                                                                                                  |                                                                                                                                                   ccGGGCCAATGTATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCtgattga                                                                                                                                                <  2:170699/150‑1 (MQ=255)            tATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGccc                                                                                                                                     <  2:99258/150‑1 (MQ=255)                aatccagctcacccgcgccctgaaaagccggcacatcttcatgctgtcgctgggtggggttattggaaccggcctgttcatgggctcgggcgtaaccatcaatCAGGATCGCGCCCACCGGCCCGCCCtgattgat                                                                                                                                               <  2:593275‑M1/33‑1 (MQ=255)                aatccagctcacccgcgccctgaaaagccggcacatcttcatgctgtcgctgggtggggttattggaaccggcctgttcatgggctcgggcgtaaccatcaatCAGGATCGCGCCCACCGGCCCGCCCtgattgat                                                                                                                                               >  1:593275‑M1/104‑136 (MQ=255)                       gTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGaa                                                                                                                          <  2:88525/150‑1 (MQ=255)                                  cctgaaaagccggcacatcttcatgctgtcgctgggtggggttattggaaccggcctgttcatgggctcgggcgtaaccatcaatCAGGATCGCGCCCACCGGCCCGCCCTGATTGATgg                                                                                                                                             >  2:29674‑M1/86‑120 (MQ=255)                                  cctgaaaagccggcacatcttcatgctgtcgctgggtggggttattggaaccggcctgttcatgggctcgggcgtaaccatcaatCAGGATCGCGCCCACCGGCCCGCCCTGATTGATgg                                                                                                                                             <  1:29674‑M1/35‑1 (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)                                                                                                               ccatcaatcagggcgggccggtgggcgcgatcctgATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAgg                                                                >  2:29674‑M2/36‑120 (MQ=255)                                                                                                               ccatcaatcagggcgggccggtgggcgcgatcctgATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAgg                                                                <  1:29674‑M2/85‑1 (MQ=255)                                                                                                                 gTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCAGGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTgtc                                >  2:396611/1‑150 (MQ=255)                                                                                                                 atcaatcagggcgggccggtgggcgcgatcctgATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGAtt                                              >  1:593275‑M2/34‑136 (MQ=255)                                                                                                                 atcaatcagggcgggccggtgggcgcgatcctgATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGAtt                                              <  2:593275‑M2/103‑1 (MQ=255)                                                                                                                              gcgcCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGACCTGGATTTTTTCAACGCTGTCAGTCATGGCacac                   <  1:340451/150‑1 (MQ=255)                                                                                                                                               ctgATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCgg  >  1:400077‑M2/4‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CCGGGCCAATGTATTTGGTGGCGTGGGCCTGAAACGAACCGGACACCGGCATCTGCACCGACAGCTCGCCCAGGCAGACCATCACCAGGTACATCAGCAAGCCTGCGACCAGGTAAGCCAGGATCGCGCCCACCGGCCCGCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGG  >  NC_002947/4253455‑4253746

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.