New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 2472115	18 (0.800)	3 (0.140)	3/302	2.6	16.1%	intergenic (+19/+130)	PP_2166/PP_2167	Anti‑anti‑sigma factor/putative PhoD‑like phosphatase
	?	AE015451	2472313 =	14 (0.650)					coding (1783/1851 nt)	PP_2167	putative PhoD‑like phosphatase
Rejected: Frequency below/above cutoff threshold.

CCCTGCACCGATTGCCTGACCGACTTGCCGTCGCAGGACCAGAATGAAGATGTGGTGCGCTCCAAGGTGCTGGAGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCTGTAAAAACTCTGAATCTCACTACGCCCGGTTTTTAAACTCTGTAGCCTCACCCCCGCCAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCA  >  AE015451/2471941‑2472342                                                                                                                                                                               |                                                                                                                                                                                                                                    cccTGCACCGATTGCCTGCCCGACTTGCCGTCGCAGGACCAGAATGAAGATGTGGTGCGCCCCAAGGTGCTGGAGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGc                                                                                                                                                                                                      <  1:83627/206‑1 (MQ=255) cccTGCACCGATTGCCTGACCGACTTGCCGTCGCAGGACCAGAATGAAGATGTGGTGCGCTCCAAGGTGCTGGAGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGACCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGt                                                                                                                                                                                                                          >  2:228943/1‑186 (MQ=255)             tGCCTGACCGACTTGCCGTCGCAGGACCAGAATGAAGATGTGGTGCGCTCCAAGGTGCTGGAGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCgcccctgttcatcaaaaaacacctgtcccagcccggcaccattccacagccgttccccggccttgctgtggctg                                                                                                                                                           <  1:426880‑M1/237‑75 (MQ=255)              gCCTGACCGACTTGCCGTCGCAGGACCAGAATGAAGATGTGGTGCGCTCCAAGGTGCTGGAGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTAg                                                                                                                                                          >  1:395409/1‑237 (MQ=255)                                cGCAGGACCAGAATGAAGATGTGGTGCGCTCCAAGGTGCTGGAGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTAGGAGCGGGTTTACCCGCGa                                                                                                                                        <  1:396353/237‑1 (MQ=255)                                             tGAAGATGTGGTGCGCTCCAAGGTGCTGGAGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCa                                                                                                                                                                     <  1:29815/195‑1 (MQ=255)                                                        tgCGCTCCAAGGTGCTGGAGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCgcccctgttcatcaaaaaacacctgtcccagcc                                                                                                                                                                                                    <  2:242576‑M1/153‑34 (MQ=255)                                                        tgCGCTCCAAGGTGCTGGAGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCgcccctgttcatcaaaaaacacctgtcccagcc                                                                                                                                                                                                    >  1:242578‑M1/1‑120 (MQ=255)                                                                         gAGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCTGTAAAAACTCTGAATCTCACTACg                                                                                               <  1:190375/237‑1 (MQ=255)                                                                          aGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTcgc                                                                                                                                                                                                                              >  2:169761/1‑109 (MQ=255)                                                                          aGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACACGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTAGGAGCggg                                                                                                                                                   >  1:258663/1‑184 (MQ=255)                                                                                                       gACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAAtgtgt                                                                                                                                                                                >  2:263073/1‑126 (MQ=255)                                                                                                       gACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAAtgtgt                                                                                                                                                                                <  1:263075/126‑1 (MQ=255)                                                                                                            cAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACGACGGCCTGCGCCCAAACACGaa                                                                                                                                                                                     >  2:162666/1‑116 (MQ=255)                                                                                                                                    ggTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTa                                                                                                                                                           >  2:16434/1‑118 (MQ=255)                                                                                                                                    ggTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTa                                                                                                                                                           <  1:16434/118‑1 (MQ=255)                                                                                                                                        aaTGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTc                                                                                                                                                                             <  1:141581/96‑1 (MQ=255)                                                                                                                                        aaTGCGCTGGAGCGGAGCTGATCTATAGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTc                                                                                                                                                                             >  2:141580/1‑96 (MQ=255)                                                                                                                                                               tataGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCcaca                                                                                                                                                                     <  2:139827/81‑1 (MQ=255)                                                                                                                                                               tataGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCa                                                                                                                                                                     >  1:139828/1‑81 (MQ=255)                                                                                                                                                                      tgtgCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCTGTAAAAACTCTGAATCTCACTACGCCCGGTTTTTAAACTCTGTAGCCTCACCCCCGCCAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCACAAAACACCTGTCCCa  >  1:132572/1‑237 (MQ=255)                                                                                                                                                                             gccgcccctgttcatcaaaaaacacctgtcccagcccggcacc                                                                                                                                                                                             >  1:195751‑M1/1‑3 (MQ=255)                                                                                                                                                                               |                                                                                                                                                                                                                                    CCCTGCACCGATTGCCTGACCGACTTGCCGTCGCAGGACCAGAATGAAGATGTGGTGCGCTCCAAGGTGCTGGAGGCGCACAAGATCCTCATGGGCCTGAATGACTCCAACCGCGAAGCCTTCCATGATCTGGTCAATGCGCTGGAGCGGAGCTGATCTATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCTGTAAAAACTCTGAATCTCACTACGCCCGGTTTTTAAACTCTGTAGCCTCACCCCCGCCAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCA  >  AE015451/2471941‑2472342

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