New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 1357753	NA (NA)	4 (0.300)	4/242	1.3	NA	intergenic (‑85/+49)	ymjA/puuP	DUF2543 family protein/putrescine importer
	?	NC_000913_3_bme_pgi	= 1357792	NA (NA)					noncoding (23/30 nt)	REP111 (repetitive extragenic palindromic) element; contains 1 REP sequences	REP111 (repetitive extragenic palindromic) element; contains 1 REP sequences

TCATTATTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCATCCGACTATTACGTTTCACTCACCGGCGTTCTGTCACCATCATACAACGGCACTTTGCGATAGCGGCGGATCAGATACCATAAATATGCGCCGCCTAACGAAGCCCATACCAGACCGAGTGTCAGTGACGTTGACTCGAGGTTAACC  >  NC_000913_3_bme_pgi/1357665‑1357940                                                                                                                                |                                                                                                                                                     tCATTATTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCcattaactataaaaagatgtat                                                                                                                                <  4:85155‑M2/150‑23 (MQ=255)          gCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCcattaactataaaaagatgtatgacccggat                                                                                                                       >  3:85155‑M2/1‑119 (MQ=255)                 gggTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCcattaactataaaaagatgtatgacccggatcacgagt                                                                                                                >  1:158974‑M2/1‑112 (MQ=255)                     tAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCATCCGACTATTACGTTTCACTCACCGGCGTTCTGTCACcatc                                                                                                            <  3:59361/150‑1 (MQ=255)                                                                 gTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCcattaactataaaaagatgtatgacccggatcacgagtgctatcatccgcccttatgcttagcgcattaaccccaggggcaaataa                                                                >  1:96030‑M2/1‑64 (MQ=255)                                                                         tCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCATCCGACTATTACGTTTCACTCACCGGCGTTCTGTCACCATCATACAACGGCACTTTGTGATAGCGGCGGATCAGATACCATAAat                                                                <  1:36901/142‑1 (MQ=255)                                                                         tCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCATCCGACTATTACGTTTCACTCACCGGCGTTCTGTCACCATCATACAACGGCACTTTGTGATAGCGGCGGATCAGATACCATAAat                                                                >  2:36901/1‑142 (MQ=255)                                                                                        tGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCATCCGACTATTACGTTTCACTCACCGGCGTTCTGTCACCATCATACAACGGCACTTTGCGATAGCGGCGGATCAGATACCATAAATATGCGCCGCCTAACGAAGCCCAt                                         <  2:70416/150‑1 (MQ=255)                                                                                                             ggcaggATAAGGCGCCGCGCATCCGACTATTACGTTTCACTCATCGGCGTTCTGTCACCAGCATACAACGGCACTTTGCGATAGCGGCGATTCAGATACCATAAACATGCGCCGCCGCACGAAGCCTATACCAGCCCGACTGTCCGTGAc                    >  2:132085/1‑150 (MQ=255)                                                                                                              gcaggATAATGCGCCGCGCATCCGACTATTACGTTTCACTCACCGGCGTTCTGTCACCATCATACAACGGCACTTTGCGATAGCGGCGGATCAGATACCATAAATATGCGCCGCCTAACGAAGCCCATACCAGACCGAGTGTCAGTGACg                   >  2:52699/1‑150 (MQ=255)                                                                                                                               gcATCCGACTATTACGTTTCACTCACCGGCGTTCTGTCACCATCATACAACGGCACTTTGCGATAGCGGCGGATCAGATACCATAAATATGCGCCGCCTAACGAAGCCCATACCAGACCGAGTGTCAGTGACGTTGACTCGAGGTTAAcc  <  1:52699/150‑1 (MQ=255)                                                                                                                                |                                                                                                                                                     TCATTATTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCATCCGACTATTACGTTTCACTCACCGGCGTTCTGTCACCATCATACAACGGCACTTTGCGATAGCGGCGGATCAGATACCATAAATATGCGCCGCCTAACGAAGCCCATACCAGACCGAGTGTCAGTGACGTTGACTCGAGGTTAACC  >  NC_000913_3_bme_pgi/1357665‑1357940

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