New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 = 291210665 (1.300)3 (0.060) 3/96 NT 4.1% coding (202/2235 nt) relA (p)ppGpp synthetase I/GTP pyrophosphokinase
?AP009048 = 2912128 78 (1.660)coding (180/2235 nt) relA (p)ppGpp synthetase I/GTP pyrophosphokinase
Rejected: Frequency below/above cutoff threshold.

GCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/2912035‑2912106
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgtTGAGATGGTGGAGATCCTCTCGACATTAAGTATGGACATTGACACGCTGCGGGCGGCGCTGCTTTTCCC  <  AP009048/2912128‑2912060
                                                                                                                                             
GCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGTGTGTCAATGTCCATACTTAATG                                                                        >  1:2328948/1‑71
 CTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGT                                                                       <  1:1741808/71‑1
                  CGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAAT                                                                         >  1:2914116/1‑52
                          GGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTTGAGATGGTGGAGATCCTCTCGACA                                              >  1:3985996/1‑71
                          GGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTTGAGATGGTGGAGATCCTCTCGACA                                              <  1:885763/71‑1
                                GCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTTGAGATGGTGGAGATCCTCTCGACATTAAGT                                        <  1:2962009/71‑1
                                                                     TGTTGAGATGGTGGAGATCCTCTCGACATTAAGTATGGACATTGACACGCTGCGGGCGGCGCTGCTTTT     >  1:365522/1‑69
                                                                     TGTTGAGATGGTGGAGATCCTCTCGACATTAAGTATG                                     <  1:3201696/37‑1
                                                                      GTTGAGATGGTGGAGATCCTCTCGACATTAAGTATGGACATTGACACGCTGCGGGCG                >  1:3831124/1‑57
                                                                       TTGAGATGGTGGAGATCCTCTCGACATTAAGTATGGACATTGACACGCTGCGGGCGGCGCTGCTTTTCCC  <  1:1810574/70‑1
                                                                                                                                             
GCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/2912035‑2912106
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgtTGAGATGGTGGAGATCCTCTCGACATTAAGTATGGACATTGACACGCTGCGGGCGGCGCTGCTTTTCCC  <  AP009048/2912128‑2912060

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