New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	1224228 =	16 (0.680)	4 (0.180)	4/264	2.6	23.0%	intergenic (+321/+51)	ycgI/minE	pseudogene/cell division topological specificity factor
	?	NC_000913_3_pae_pgi	1224254 =	12 (0.550)					intergenic (+347/+25)	ycgI/minE	pseudogene/cell division topological specificity factor

AAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATC  >  NC_000913_3_pae_pgi/1224080‑1224341                                                                                                                                                     |                                                                                                                  aacgTTCTTTACCCAACCAGGCCTAAACATGTTGTTGTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACt                                                                                                                   <  2:120938/146‑1 (MQ=255)   ggTTATTTCCCCAAACAGGCGTAACCATGTTGTTTTGTGCGTGGGAACAGACTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACACGGTAGATAAAAGAGATAAAAAACTCa                                                                                                                 >  2:256846/1‑149 (MQ=255)     ttATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGAGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGAAGAGACAGATCAAAAAATCAAg                                                                                                               >  2:177693/1‑149 (MQ=255)                 ccAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCAACACGTGGATGGGAGGGAGGAAAAACTCAAGGCAGAGATaaat                                                                                                   >  2:125150/1‑147 (MQ=255)                   aGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActct                                                                                                 >  2:360235/1‑149 (MQ=255)                              ttttgttttgTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGGGGGATGAAAACCTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                                                      >  1:166303/1‑149 (MQ=255)                                           gTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATACGCTTTTa                                                                         >  1:179141/1‑149 (MQ=255)                                                       ttAAGGTGTAAAGGGGGAGGGGGAAATAGCAATGAGGAGTATCAGCAAGAAGAATCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAt                                                             <  2:23718/149‑1 (MQ=255)                                                          aGGTGTAAAGGGGGAGGTGGAAATAGCAATGCGGAGTAGCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGACAGAGATAACTCTGCCTTGAAGATAAATGCGCTGTTACAGCGGGCTTATTTc                                                          >  2:180192/1‑149 (MQ=255)                                                            gagtAAAGGGGGAGGTGGAAAAAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGTCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                        <  2:410100/147‑1 (MQ=255)                                                             tgtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                       >  2:300788/1‑149 (MQ=255)                                                                               aaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg                                     >  1:162809/1‑149 (MQ=255)                                                                                      tgaacgtgaccttaccggaagcagaagagctgaaataagcccgctgtaaaagcgcatttatctTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                                  <  1:34109‑M1/50‑1 (MQ=255)                                                                                      tgaacgtgaccttaccggaagcagaagagctgaaataagcccgctgtaaaagcgcatttatctTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                                  >  2:34109‑M1/64‑113 (MQ=255)                                                                                             gTATCATCAAGAATACTCGCCGCTTTACCAAAACCTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTc                       <  1:320205/149‑1 (MQ=255)                                                                                              tATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCa                      >  1:271628/1‑149 (MQ=255)                                                                                                     cggaagcagaagagctgaaataagcccgctgtaaaagcgcatttatctTCAAGGCAGAGATaa                                                                                                     >  2:170413‑M1/49‑63 (MQ=255)                                                                                                     cggaagcagaagagctgaaataagcccgctgtaaaagcgcatttatctTCAAGGCAGAGATaa                                                                                                     <  1:170413‑M1/15‑1 (MQ=255)                                                                                                               gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGcc     >  2:257734/1‑149 (MQ=255)                                                                                                                  gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATc  >  1:220932/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                  AAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATC  >  NC_000913_3_pae_pgi/1224080‑1224341

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