New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4305533	2 (0.040)	19 (0.420) +CGCCGCTTCA	15/468	1.2	95.2%	coding (1065/1986 nt)	yjcS	linear primary‑alkylsulfatase
	?	NC_000913	4324315 =	0 (0.000)					coding (38/1017 nt)	phnD	phosphonate/phosphate ABC transporter periplasmic binding protein

CATCGCCGACAGATTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATCCCCTGGATAATGCCCGCGTAGTCCGGGGCAAAGAAGGCGTTGACCTTCACGCCCAGCTTCTTCTCCATATCCTGCAAGAACGGCGTCCATTGCGGTTTCAGGTTTTGCTGTGATTCCGTTGAAATAATGCCGAAATTCAACGCCTTTTCCTGCTCTTCGGCGTGCGCCGGGCTTAACAGGGTGCTGAGGCTGAACATGCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGCCGTC  >  NC_000913/4324065‑4324562                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        tcgtaatagccaagatagaagttatacaccgcgcgggcgttgtggctgacagaaccgtaatagccgcggctggcccagttattggcaagtgcaggcggcagcttaatcatgtcgccgatttcattcatggtgtagccctggttcgccaggtgcagggtctggtcgtgaatgtacttgatgggatcgcggtatttaccaatataatcattgatatgcttattgccccagaccggccaggtgcgccgcttcac                                                                                                                                                                                                                                                         >  1:353708‑M2/251‑251 (MQ=255)               atagacgttatacaccgcgcgggcgttgtggctgacagaaccgtaatagcagaggctggcccagttattggcaagtgcaggcggcagcttaatcatgtcgccggtttcattcatggtgtagccctggttcgcccggtgcaaggtctggtcgtgaatgtacttgattgtatctcggtatttaccaatataatcattgatatgtttattaccccagaccgaccaggagcgacgctacaATGGTGACGGCCAGc                                                                                                                                                                                                                                           >  2:225408‑M2/237‑251 (MQ=255)                                      cgttgtggctgacagaaccgtaatagccgcggctggcccagttattggcaagtgcaggcggcagcttaatcatgtcgccgatttcattcatggtgtagccctggttcgccaggtgcagggtctggtcgtgaatgtacttgatggtatcgcggtatttaccaatataatcattgatatgcttattgccccagaccggccaggtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTCGCGTTCAt                                                                                                                                                                                                                    >  2:7395‑M2/214‑251 (MQ=255)                                                         gtaatagccgcggctggcccagttattggcaagtgcaggcggcagcttaatcatgtcgccgatttcattcatggtgtagccctggttcgccaggtgcagggtctggtcgtgaatgtacttgatggtatcgcggtatttaccaatataatcattgatatgcttattgccccagaccggccaggtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATgg                                                                                                                                                                                                 >  1:136619‑M2/195‑251 (MQ=255)                                                                    ggctggcccagttattggcaagtgcaggcggcagctttatcatgtcgccgatttcattcatggtgtagccctggttcgccaggtgcagggtctggtcgtgaatgtacttgatggtatcgcggtatttaccaatataatcattgatatgcttattgccccagaccggccaggtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGc                                                                                                                                                                                      >  1:182782‑M2/184‑251 (MQ=255)                                                                           ccagttattggcaagtgcaggcggcagcttaatcatgtcgccgatttcattcatggtgtagccctggttcgccaggtgcagggtctggtcgtgaatgtacttgatggtatcgcggtatttaccaatataatcattgatatgcttattgccccagaccggccaggtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTc                                                                                                                                                                               >  1:93899‑M2/177‑251 (MQ=255)                                                                              tattggcaagtgcaggcggcagcttaatcatgtcgccgatttcattcatggtgtagccctggttcgccaggtgcagggtctggtcgtgaatgtacttgatggtatcgcggtatttaccaatataatcattgatatgcttattgccccagaccggccaggtgtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCCATGATGGGGATGTCAGGCAGCTTTcgcg                                                                                                                                                                            >  1:139791‑M2/174‑251 (MQ=255)                                                                                                        taatcatgtcgccgatttcattcatggtgtagccctggttcgccaggtgcagggtctggtcgtgaatgtacttgatggtatcgcggtatttaccaatataatcattgatatgcttattgccccagaccggccaggtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAAtgctg                                                                                                                                                    <  1:7395‑M2/102‑1 (MQ=255)                                                                                                                                 ggtgtagccctggttcgccaggtgcagggtctggtcgtgaatgtacttgatggtagcgcggtatttaccaatataatcattgatatgcttattgccccaggccggccaggtgcgccgcttcaCTGGTGAAGGGCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGGTCGTTa                                                                                                                         <  2:353708‑M2/129‑1 (MQ=255)                                                                                                                                  gtgtagccctggttcgccaggtgcagggtctggtcgtgaatgtacttgatggtatcgcggtatttaccaatataatcattgatatgcttattgccccagaccggccaggtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTAt                                                                                                                        >  1:69321‑M2/122‑251 (MQ=255)                                                                                                                                                          agggtctggtcgtgaatgtacttgatggtatcgcggtatttaccaatataatcattgatatgcttattgccccagaccggccaggtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCtt                                                                                                                                                                     >  1:220094‑M2/98‑182 (MQ=255)                                                                                                                                                          agggtctggtcgtgaatgtacttgatggtatcgcggtatttaccaatataatcattgatatgcttattgccccagaccggccaggtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCtt                                                                                                                                                                     <  2:220094‑M2/85‑1 (MQ=255)                                                                                                                                                                              tgatggtatcgcggtatttaccaatataatcattgatatgcttattgccccagaccggccaggtgtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGgc                                                                              <  1:222906‑M2/172‑1 (MQ=255)                                                                                                                                                                                          ggtatttaccaatataatcattgatatgcttattgccccagacaggccaggtgtgcgccgcttcaCTGGTGAAGGCCAGCGCGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCACTGGCGCAGGGc                                                                       >  1:232092‑M2/66‑244 (MQ=255)                                                                                                                                                                                          gggatttaccaatataatcattgatatgcttattgccccagaccggccaggtgtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTGTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGGCAAAACGGTCGTTATCAAACTGTGGGCGGCTGCCGGCGTAGAAGACTGGCCCCTGGCGCAGGGc                                                                       <  2:232092‑M2/179‑1 (MQ=255)                                                                                                                                                                                                     caatataatcattgatatgcttattgccccagaccggccaggtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGAGGTCAGGCAGCGTTCGCGTTCTCTTCTACGCGGTGAATGCTGCGGTAGAGATGGGCAAAACGTTCGTTATCAAACTTTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTa                                                     <  1:225408‑M2/197‑1 (MQ=255)                                                                                                                                                                                                        taatcattgatatgcttattgccccagaccggccaggtgtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAgcg                                                  >  1:299506‑M2/52‑251 (MQ=255)                                                                                                                                                                                                             attgatatgcttattgccccagaccggccaggtgtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAggg                                             >  1:412404‑M2/47‑251 (MQ=255)                                                                                                                                                                                                                                      cgtttttcgtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTgac                    >  1:27997‑M2/22‑251 (MQ=255)                                                                                                                                                                                                                                      cgtttttcgtgcgccgcttcaCTGGGGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTgac                    >  1:344935‑M2/22‑251 (MQ=255)                                                                                                                                                                                                                                       gtttttcgtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTTCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTgacg                   >  2:339760‑M2/21‑251 (MQ=255)                                                                                                                                                                                                                                       gtttttcgtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGAgacg                   >  2:255391‑M2/21‑251 (MQ=255)                                                                                                                                                                                                                                         gccaggtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTAGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTgacgac                 >  1:383317‑M2/19‑251 (MQ=255)                                                                                                                                                                                                                                              gtgcgccgcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACgg            >  1:210892‑M2/14‑251 (MQ=255)                                                                                                                                                                                                                                                gcgccgtttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTg          >  1:135831‑M2/12‑251 (MQ=255)                                                                                                                                                                                                                                                     gcttcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGcc     >  1:333678‑M2/7‑251 (MQ=255)                                                                                                                                                                                                                                                        tcaCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCACTGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGCcgtc  >  1:208976‑M2/4‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        CATCGCCGACAGATTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATCCCCTGGATAATGCCCGCGTAGTCCGGGGCAAAGAAGGCGTTGACCTTCACGCCCAGCTTCTTCTCCATATCCTGCAAGAACGGCGTCCATTGCGGTTTCAGGTTTTGCTGTGATTCCGTTGAAATAATGCCGAAATTCAACGCCTTTTCCTGCTCTTCGGCGTGCGCCGGGCTTAACAGGGTGCTGAGGCTGAACATGCTGGTGAAGGCCAGCGAGGCAATTATCTTAGCGTTCATTCGTTTTCCTCAATGATGGGGATGTCAGGCAGCTTTCGCGTTCTCTTCGACGCGGTTAATGCTGCGGTAGAGATGGTCAAAACGTTCGTTATCAAACTGTTGGCTGCTGCCGTCGTAGAAGACGTGCCCCTGGCGCAGGGCGACGATGCGTTCGCAGTAGCGCAGGGCGTAATCCACCTGATGCAGCGTGACGACCACGGTGATGCCGTC  >  NC_000913/4324065‑4324562

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