New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 396770354 (1.330)6 (0.160) 8/236 2.7 10.3% intergenic (+27/‑213) rho/wecA transcription termination factor/UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
?NC_000913 = 3967707 55 (1.490)intergenic (+31/‑209) rho/wecA transcription termination factor/UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
Rejected: Frequency below/above cutoff threshold.

CAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3967548‑3967703
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaacgccacgtAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCA  <  NC_000913/3967707‑3967564
                                                                                                                                                                                                                                                                                                            
CAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTACTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAAC                                                                                                                                                         >  4:126902/1‑149
             TCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACAC                                                                                                                                            >  1:179431/1‑149
                  CGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGC                                                                                                                                       <  2:179431/149‑1
                                                                                AATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATC                                                                                                   >  7:156655/1‑123
                                                                                AATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATC                                                                                                   <  8:156655/123‑1
                                                                                                             AAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGA                                                  >  3:156292/1‑143
                                                                                                             AAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGA                                                  <  4:156292/143‑1
                                                                                                                                                AAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATT         <  2:126175/149‑1
                                                                                                                                                  AACGCCACGTGTTTACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATC                     >  3:92354/1‑135
                                                                                                                                                  AACGCCACGTGTTTACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATC                     <  4:92354/135‑1
                                                                                                                                                  AACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCA                       >  1:206782/1‑133
                                                                                                                                                  AACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCA                       <  2:206782/133‑1
                                                                                                                                                       CACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCA  <  3:42316/149‑1
                                                                                                                                                                                                                                                                                                            
CAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3967548‑3967703
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaacgccacgtAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCA  <  NC_000913/3967707‑3967564

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