New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	1811415 =	18 (0.700)	3 (0.140)	3/272	NT	18.0%	intergenic (+22/+106)	ispD/PP_1615	2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase/Transcriptional regulator, LysR family
	?	AE015451	1811495 =	12 (0.540)					intergenic (+102/+26)	ispD/PP_1615	2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase/Transcriptional regulator, LysR family

TCAGGCGCCGCGGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTG  >  AE015451/1811294‑1811688                                                                                                                                                                                                          |                                                                                                                                                                                                  tCAGGCGCCGCGGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACGCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCCGCCGaa                                                                                                                                                                 >  1:236436/1‑236 (MQ=255)                              gTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                                                                                               >  1:458123/1‑177 (MQ=255)                                                  cGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACcagca                                                                                                                                                                                    >  2:483275/1‑168 (MQ=255)                                                                    tGGTTGCGCCAGCGCCGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCGTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTCGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGACACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATccac                                                                                                                    <  1:162627/214‑1 (MQ=255)                                                                                                                ggCTGGCCTTTTCGCGGGTAACCCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCaa                                                                                                                                         <  1:94640/149‑1 (MQ=255)                                                                                                                     gCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGtt                                             <  1:274548/236‑1 (MQ=255)                                                                                                                        ttttCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCtt                                                                                                                                >  2:316143/1‑150 (MQ=255)                                                                                                                        ttttCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAAACCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTAcc                                                                                                       >  1:91327/1‑175 (MQ=255)                                                                                                                             gcgGGTCAACCCGCTCCGACAATGGATTGCATCGGTGGTTGGACCCAGACCGTCGTAGGAGCGGGTTTACCCGCGGAGCGGCCAGCACCAGCAGCACAAccc                                                                                                                                                                           >  2:94636/1‑102 (MQ=25)                                                                                                                               gggtaaacccgctcctacgacggtctgtgtccatccaccgatgcaatccattgtaggagcgggtttacccgcgaaGAGGCCAGTa                                                                                                                                                                                          <  1:54844‑M2/10‑1 (MQ=255)                                                                                                                               GGGTAACCCCGCTCCTACAGTGGATTGCATCGGTGGTTGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAAcc                                                                                                                                       >  1:330124‑M1/1‑136 (MQ=255)                                                                                                                                              tacgacggtctgtgtccaaccaccgatgcaatccattgtaggagcgggtttacccgcgaaGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCcct                                                                                                                                    >  1:449614‑M2/61‑124 (MQ=255)                                                                                                                                              tacgacggtctgtgtccaaccaccgatgcaatccattgtaggagcgggtttacccgcgaaGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCcct                                                                                                                                    <  2:449610‑M2/64‑1 (MQ=255)                                                                                                                                                        tGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCAACCAGCCGCCGTCCCTTc                                                                                                    >  2:412638/1‑146 (MQ=255)                                                                                                                                                                tggtTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGgtg  <  1:418743/236‑1 (MQ=255)                                                                                                                                                                  gtTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTTCGCCTCCAGCAGCGACACCAGCGcacca             >  2:210724/1‑223 (MQ=255)                                                                                                                                                                         acaGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAg                                                                       <  1:403274/158‑1 (MQ=255)                                                                                                                                                                                  tgtaggagcgggtttacccgcgaaGAGGc                                                                                                                                                                                               >  2:84735‑M2/25‑29 (MQ=255)                                                                                                                                                                                   gtaggagcgggtttacccgcgaagag                                                                                                                                                                                                 >  2:52057‑M2/24‑26 (MQ=255)                                                                                                                                                                                               ttacccgcgaaGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCaa                                                                                                                                         <  2:119575‑M2/59‑1 (MQ=255)                                                                                                                                                                                               ttacccgcgaaGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCaa                                                                                                                                         >  1:119579‑M2/12‑70 (MQ=255)                                                                                                                                                                                                   ccgcgaaGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCg                                                                                                                                             <  2:49567‑M2/55‑1 (MQ=255)                                                                                                                                                                                                          |                                                                                                                                                                                                  TCAGGCGCCGCGGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTG  >  AE015451/1811294‑1811688

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