New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi 97574 =5 (0.330)3 (0.210) 3/276 1.8 38.1% coding (488/1317 nt) murD UDP‑N‑acetylmuramoyl‑L‑alanine:D‑glutamate ligase
?NC_000913_3_vch_pgi 97595 = 5 (0.350)coding (509/1317 nt) murD UDP‑N‑acetylmuramoyl‑L‑alanine:D‑glutamate ligase

CACGCTAGTGGGTGAAATGGCGAAAGCGGCGGGGGTTAACGTTGGTGTGGGTGGCAATATTGGCCTGCCTGCGTTGATGCTACTGGATGATGAGTGTGAACTGTACGTGCTGGAACTGTCGAGCTTCCAGCTGGAAACCACCTCCAGCTTACAGGCGGTAGCAGCGACCATTCTGAACGTGACTGAAGATCATATGGATCGCTATCCGTTTGGTTTACAACAGTATCGTGCAGCAAAACTGCGCATTTACGAAAACGCGAAAGTTTGCGTG  >  NC_000913_3_vch_pgi/97446‑97716
                                                                                                                                |                                                                                                                                              
attaaccacgcaaactttcgcgttttcgtaaatgcgcagttttgctgcacgatactgttgtaaaccaaacggatagcgatccatatgatcttcagtcacgttcagaatggtcgctgctaccgcctgtaAGCTGGAAACCACCTCCAGCt                                                                                                                            >  1:225178‑M1/129‑149 (MQ=255)
           gTGAAATGGCGAAAGCGGCGGGGGTTAACGTTGGTGTGGGTGGCAATATTGGCCTGCCTGCGTTGATGCTGCAGGATGATGAGTGAGAACTGTACGTGCTGGAACTGTCGAGCTACCAGCTGGAAACCACCTCCAGCTTACAGGCGGTa                                                                                                                 >  1:237103/1‑149 (MQ=255)
                                                                  aaacggatagcgatccatatgatcttcagtcacgttcagaatggtcgctgctaccgcctgtaAGCTGGAAACCACCTCCAGCTTACAGGCg                                                                                                                    <  1:68622‑M1/29‑1 (MQ=255)
                                                                  aaacggatagcgatccatatgatcttcagtcacgttcagaatggtcgctgctaccgcctgtaAGCTGGAAACCACCTCCAGCTTACAGGCg                                                                                                                    >  2:68622‑M1/63‑91 (MQ=255)
                                                                                                      ggaCGGGCTGGAACTGTCGAGCTTCCAGCTGGAAACCACCTCCAGCTTACAGGCGGTAGCAGCGACCATTCTGAACGTGACTGAAGATCATATGGATCGCTATCCGTTTGGTTTACAACAGTATCGTGCAGCAAAACTGCGCATTTACg                      <  1:25091/147‑1 (MQ=255)
                                                                                                        aCGTGCTGGAACGGTCGAGCTTCCAGCTGGAAACCACCTCCAGCTTACAGGCGGTAGCAGCGACCATTCTGAACGTGACTGAAGATCATATGGATCGCTATCCGTTTGGTTTACAACAGAATCGTTGAGCAAAACTGCGCATTTACGaa                    <  2:164193/149‑1 (MQ=255)
                                                                                                          gTGCTGGAACTGTCGAGCTTCCAGCTGGAAACCACCTCCAGCTTACAGGCGGTAGCAGCGACCATTCTGAACGTGACTGAAGATCATATGGATCGCTATCCGTTTGGTTTACAACAGTATCGTGCAGCAAAACTGCGCATTTACGaaaa                  >  1:143702/1‑149 (MQ=255)
                                                                                                                gaccTGTCGAGCTTCCAGCTGGAAACCACCTCCAGCTTACAGGCGGTAGCAGCGACCATTCTGAACGTGACTGAAGATCATATGGATCGCTATCCGTTTGGTTTACAACAGTATCGTGCAGCAAAACTGCGCATTTACGAAAACGCGaa            <  2:83839/146‑1 (MQ=255)
                                                                                                                          ggttctAGCTGGAAACCACCTCCAGCTTACAGGCGGTAGCAGCGGCCATTCTGAACGGGACTGAAGATCATATGGATCGCTATCCGTTTGGTTTACAACAGTATCGTGCAGCAAAACTGCGCATTTACGAAAACGCGAAAGTTTGCGTg  >  2:127208‑M1/7‑149 (MQ=255)
                                                                                                                                |                                                                                                                                              
CACGCTAGTGGGTGAAATGGCGAAAGCGGCGGGGGTTAACGTTGGTGTGGGTGGCAATATTGGCCTGCCTGCGTTGATGCTACTGGATGATGAGTGTGAACTGTACGTGCTGGAACTGTCGAGCTTCCAGCTGGAAACCACCTCCAGCTTACAGGCGGTAGCAGCGACCATTCTGAACGTGACTGAAGATCATATGGATCGCTATCCGTTTGGTTTACAACAGTATCGTGCAGCAAAACTGCGCATTTACGAAAACGCGAAAGTTTGCGTG  >  NC_000913_3_vch_pgi/97446‑97716

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