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.

TGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCTGTAAAAACTCTGAATCTCACTACGCCCGGTTTTTAAACTCTGTAGCCTCACCCCCGCCAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCGGCACCATTCCACAGCCGTTCCCCGGCCTTGCTGTGGCTGGGCGTACGCGGCACCACCTCCATTTCTCGGCGCTTGGTAAACCAGTTGAGCGGCGAGCGCGGTGCATACAACCAGCGGTTGAGCCGGTCGAGCACATCCAGCAAGCGCCGGGGGAAC  >  AE015451/2472106‑2472503                                                                                                                                                                                                                |                                                                                                                                                                                               tgtgCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCTGTAAAAACTCTGAATCTCACTACGCCCGGTTTTTAAACTCTGTAGCCTCACCCCCGCCAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCACAAAACACCTGTCCCa                                                                                                                                                                   >  1:132572/1‑237 (MQ=255)                                             tgcctgaccgacttgccgtcgcaggaccagaatgaagatgtggtgcgctccaaggtgctggaggcgcacaagatcctcatgggcctgaatgactccaaccgcgaagccttccatgatctggtcaatgcgctggagcggagctgatctatagcctgtgccggccGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCGGCACCATTCCACAGCCGTTCCCCGGCCTTGCTGTGGCTg                                                                                                                       <  1:426880‑M2/74‑1 (MQ=255)                                                           tgtgtgTCCACCCCCAATCAGGTGTAGGTGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCTGTAAAAACTCTGAATCTCACTACGCCCGGTTTTTAAACTCTGTAGCCTCACCCCCGCCAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCGGCACCATTCCACAGCCGtt                                                                                                                                           <  1:162667/203‑1 (MQ=255)                                                                  ccACCCCCAATCAGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCTGTAAAAACTCTGAATCTCACTACGCCCGGTTTTTAAACTCTGTAGCCTCACCCCCGCCAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGcc                                                                                                                                                                <  2:258661/175‑1 (MQ=255)                                                                   taCCCCCAATCAGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCTGTAAAAACTCTGAATCTCACTACGCCCGGTTTTTAAACTCTGTAGCCTCACCCCCGCCAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCGGCACCATTCCACAGCCGTTCCCCGGCCTTGCTGTGGCTGGGCGTACGCGGCACCACCTCCa                                                                                                 >  1:415818/2‑237 (MQ=255)                                                                                        tgcgctccaaggtgctggaggcgcacaagatcctcatgggcctgaatgactccaaccgcgaagccttccatgatctggtcaatgcgctggagcggagctgatctatagcctgtgccggccGCCCCTGTTCATCAAAAAACACCTGTCCCAGcc                                                                                                                                                                >  1:242578‑M2/121‑153 (MQ=255)                                                                                        tgcgctccaaggtgctggaggcgcacaagatcctcatgggcctgaatgactccaaccgcgaagccttccatgatctggtcaatgcgctggagcggagctgatctatagcctgtgccggccGCCCCTGTTCATCAAAAAACACCTGTCCCAGcc                                                                                                                                                                <  2:242576‑M2/33‑1 (MQ=255)                                                                                                                                   tGAATCTCACTACGCCCGGTTTTTAAACTCTGTAGCCTCACCCCCGCCAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCGGCACCATTCCACAGCCGTTCCCCGGCCTTGCTGTGGCTGGGCGTACGCGGCACCACCTCCATTTCTCGGCGCTTGGTAAACCAGTTGAGCGGCGAGCGCGGTGc                                                      >  2:192247/1‑216 (MQ=255)                                                                                                                                                                   tAGCCTCACCCCCGCCAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCGGCACCATTCCACAGCCGTTCCCCGGCCTTGCTGTGGCTGGGCGTACGCGGCACCACCTCCATTTCTCGGCGCTTGGTAAACCAGTTGAGCGGCGAGCGCGGTGCATACAACCAGCGGTTGAGCCGGTCGAGCACATCCAGCAAGCGCCGGGGCAAc  >  1:229484/1‑236 (MQ=255)                                                                                                                                                                                  cAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTg                                                                                                                                                                        <  1:215695/55‑1 (MQ=255)                                                                                                                                                                                  cAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTg                                                                                                                                                                        >  2:215694/1‑55 (MQ=255)                                                                                                                                                                                   aGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCGGCACCATTCCACAGCCGTTCCCCGGCCTTGCTGTGGCTGGGCGTACGCGGCACCACCTCCATTTCTCg                                                                                          <  1:401421/132‑1 (MQ=255)                                                                                                                                                                                   aGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCGGCACCATTCCACAGCCGTTCCCCGGCCTTGCTGTGGCTGGGCGTACGCGGCACCACCTCCATTTCTCg                                                                                          >  2:401419/1‑132 (MQ=255)                                                                                                                                                                                               gTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGt                                                                                                                                                                       <  2:195750/43‑1 (MQ=255)                                                                                                                                                                                                  cacCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCGGCACCATTCCACAGCCGTTCCCCGGCCTTGCTGTGGCCGGGCGTACGCCGCACCACCTCCAt                                                                                                >  2:17828/1‑111 (MQ=255)                                                                                                                                                                                                        cgcCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCgg                                                                                                                                                             >  1:209402/1‑44 (MQ=255)                                                                                                                                                                                                        cgcCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCgg                                                                                                                                                             <  2:209401/44‑1 (MQ=255)                                                                                                                                                                                                             gccGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCGGCAcc                                                                                                                                                         >  1:195751‑M2/4‑43 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                               TGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGCCTGCGCCCAAACACCAATGTGTGTCCACCCCCAATCAGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCTGTAAAAACTCTGAATCTCACTACGCCCGGTTTTTAAACTCTGTAGCCTCACCCCCGCCAGCATCCAGCTGGTACACCCGCGCCGGCCGCCCCTGTTCATCAAAAAACACCTGTCCCAGCCCGGCACCATTCCACAGCCGTTCCCCGGCCTTGCTGTGGCTGGGCGTACGCGGCACCACCTCCATTTCTCGGCGCTTGGTAAACCAGTTGAGCGGCGAGCGCGGTGCATACAACCAGCGGTTGAGCCGGTCGAGCACATCCAGCAAGCGCCGGGGGAAC  >  AE015451/2472106‑2472503

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