New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 296246916 (0.820)4 (0.240) 3/114 1.5 23.7% coding (296/324 nt) ppdC putative prepilin peptidase‑dependent protein
?NC_000913 = 2962525 12 (0.720)coding (240/324 nt) ppdC putative prepilin peptidase‑dependent protein

CGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTG  >  NC_000913/2962403‑2962533
                                                                  |                                                                
cGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGt                                                                 >  2:625147/1‑68 (MQ=255)
                cccTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTg                                                 >  1:618360/1‑68 (MQ=255)
                      aCATAGCGGCTCCTGACAACAGACGAATCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCg                                           >  1:71494/1‑68 (MQ=255)
                       cATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                          >  1:221035/1‑68 (MQ=255)
                       cATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                          >  1:390207/1‑68 (MQ=255)
                       cATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                          >  1:322580/1‑68 (MQ=255)
                       cATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                          >  2:64852/1‑68 (MQ=255)
                       cATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                          >  1:209369/1‑68 (MQ=255)
                        aTAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCggg                                         <  2:311101/68‑1 (MQ=255)
                              gCTCCTGACTACTGACGATTAGGGCAATGCAGGAGGGTCATCTCGCCTTCTCTGCACCCGGGTGaaaa                                   >  2:613158/1‑67 (MQ=255)
                              gCTCCTGACTACTGACGAATCGGGCAATGCAGGCAGGTCATATCGCCTTCTCTGCCCACGGGTGatca                                   >  2:613801/1‑65 (MQ=255)
                                        actGACGATTCGGGCAATGCAGGCGGGatgtgtcagcatcagcgttacgctagtttcacccgggggca                         >  2:530930‑M1/1‑27 (MQ=255)
                                          tGACGATTCGGGCAATGCAGGCGGGatgtgtcagcatcagcgttacgctagtttcacccgggggcaga                       <  1:530930‑M1/68‑44 (MQ=255)
                                            aCGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGa                     <  2:109979/68‑1 (MQ=255)
                                             cGATTCGGGCAATGCAGGCGGGatgtgtcagcatcagcgttacgctagtttcacccgggggcagagaa                    >  2:390207‑M1/1‑22 (MQ=255)
                                             cGATTCGGGCAATGCAGGCGGGatgtgtcagcatcagcgttacgctagtttcacacgggggcagagag                    >  2:351584‑M1/1‑22 (MQ=255)
                                              gATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATg                   <  1:323748/68‑1 (MQ=255)
                                                    ggCAATGCTGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGaaa                                    <  2:623090/45‑1 (MQ=255)
                                                    ggCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGaaa                                    >  1:623089/1‑45 (MQ=255)
                                                      cAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGAcaca           <  2:694296/68‑1 (MQ=255)
                                                           cAGGCGGGatgtgtcagcatcagcgttacgctagtttcacccggtggcagagaaggcgagatgacccg      <  2:317101‑M1/68‑61 (MQ=255)
                                                           cAGGCGGGatgtgtcagcatcagcgttacgctagattcacccgggggcagagaaggcgagatgacccg      >  2:221035‑M1/1‑8 (MQ=255)
                                                             ggCGGGatgtgtcagcatcagcgttacgctagtttcacccgggggcagagaaggcgagatgacccgcc    >  1:254026‑M1/1‑6 (MQ=255)
                                                             ggCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATcccgcc    >  1:254026‑M2/1‑62 (MQ=255)
                                                               cGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAATCTAGCGTAACGCTGATGCTGACACATcccgcctg  >  2:221035‑M2/1‑60 (MQ=255)
                                                               cGGGTCATCTCGCCTTCTCTGCCACCGGGTGAAACTAGCGTAACGCTGATGCTGACACATcccgcctg  <  2:317101‑M2/68‑9 (MQ=255)
                                                                  |                                                                
CGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTG  >  NC_000913/2962403‑2962533

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