New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 172788563 (1.200)5 (0.100) 4/236 NT 7.4% coding (321/2487 nt) plsB glycerol‑3‑phosphate acyltransferase
?NC_002947 = 1727900 65 (1.290)coding (306/2487 nt) plsB glycerol‑3‑phosphate acyltransferase

CAGCCGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/1727732‑1727885
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcCCCGCCCACGCTCGAGCGCCTGGTG  <  NC_002947/1727900‑1727876
                                                                                                                                                                                   
CAGCCGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAG                               <  1:1254945/150‑1
  GCCGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGC                             <  2:395801/150‑1
                             AGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCCCCGCCCACGCTCGAGCGCCTGGTG  <  2:1397622/150‑1
                                                                                    CCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCCCCGCCCACGCTCGAGCGCCTGG    >  1:962961/1‑93
                                                                                    CCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCCCCGCCCACGCTCGAGCGCCTGG    <  2:962961/93‑1
                                                                                                      AATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCCCCGCCCACGCTCGAGCGCCTGG    >  1:849246/1‑75
                                                                                                      AATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCCCCGCCCACGCTCGAGCGCCTGG    <  2:849246/75‑1
                                                                                                                                                                                   
CAGCCGGCGCAGGCGCCCGGTTACCGCCCAGCTGTCGGCGAACAGCAGCTTCCAGGGGCTGGACTCGCTGGCCGGGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/1727732‑1727885
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcCCCGCCCACGCTCGAGCGCCTGGTG  <  NC_002947/1727900‑1727876

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