New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 43055372 (0.040)6 (0.130)
+CCGCTTCA		 6/472 2.4 86.1% coding (1061/1986 nt) yjcS linear primary‑alkylsulfatase
?NC_000913 4324315 = 0 (0.000)coding (38/1017 nt) phnD phosphonate/phosphate ABC transporter periplasmic binding protein

CGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGT  >  NC_000913/4305295‑4305750
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                     
cgtAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGctgaagc                                                                                                                                                                                                               <  1:247941/251‑7 (MQ=255)
  tAATAGCTAAGATAGAAGTTATACACCGCGGGGGCGTTGTGGCTGACAGAACCGTAATATCCGCGGCTTGCCCAGTTATTGGCAATTTCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTATCCCTTGTTCGCCATGTGCAGGGTCTGGTCGTGAAGGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTAGTTCCCCATACCGGCCAGTTGTGCggtccaagcg                                                                                                                                                                                                              <  2:303075/250‑9 (MQ=255)
                                                                              tATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGccgcttcactggtgaaggccagcgaggcaattatcttagcgttcattcgttttcctccatgatggggatgtcaggcagctttcgcg                                                                                                                                 >  1:139791‑M1/1‑165 (MQ=255)
                                                                                                                                                                              tGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGccgcttcactggtgaaggccagcgaggcaattatcttagcgttcattcgttttcctcaatgatggggatgtcaggcagctttcgcgttctcttcgacgcggttaatgctgcggtagagatggtcaaaacgttcgttatcaaactgttggctgctgccgtcgtagaagacgtgcccctggc                                   <  1:222906‑M1/249‑181 (MQ=255)
                                                                                                                                                                                          gggatttACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGccgcttcactggtgaaggccagcgaggcaattatcttagcgttcattcgttgtcctcaatgatggggatgtcaggcagctttcgcgttctcttcgacgcggttaatgctgcggtagagatgggcaaaacggtcgttatcaaactgtgggcggctgccggcgtagaagactggcccctggcgcagggc                            <  2:232092‑M1/241‑188 (MQ=255)
                                                                                                                                                                                          ggTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACAGGCCAGGTGTGCGccgcttcactggtgaaggccagcgcggcaattatcttagcgttcattcgttttcctcaatgatggggatgtcaggcagctttcgcgttctcttcgacgcggttaatgctgcggtagagatggtcaaaacgttcgttatcaaactgttggctgctgccgtcgtagaagacgtgccactggcgcagggc                            >  1:232092‑M1/1‑57 (MQ=255)
                                                                                                                                                                                                        taATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGccgcttcactggtgaaggccagcgaggcaattatcttagcgttcattcgttttcctcaatgatggggatgtcaggcagctttcgcgttctcttcgacgcggttaatgctgcggtagagatggtcaaaacgttcgttatcaaactgttggctgctgccgtcgtagaagacgtgcccctggcgcagggcgacgatgcgttcgcagtagcg       >  1:299506‑M1/1‑43 (MQ=255)
                                                                                                                                                                                                             aTTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGccgcttcactggtgaaggccagcgaggcaattatcttagcgttcattcgttttcctcaatgatggggatgtcaggcagctttcgcgttctcttcgacgcggttaatgctgcggtagagatggtcaaaacgttcgttatcaaactgttggctgctgccgtcgtagaagacgtgcccctggcgcagggcgacgatgcgttcgcagtagcgcaggg  >  1:412404‑M1/1‑38 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                     
CGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAGTTGTGCAGGGTATGCGTGGCGTTCTCGGCGGTACACAGGGCTTTCAGGGCCGGAATATAGAAGTGCATTTCGGCTGGCGCTTCGCTACCTGGGGTCATCAGGAAGT  >  NC_000913/4305295‑4305750

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