New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 6171392 =75 (1.150)3 (0.050) 3/254 NT 3.5% coding (910/1836 nt) glmS L‑glutamine/D‑fructose‑6‑phosphate aminotransferase
?NC_002947 6171430 = 94 (1.500)coding (872/1836 nt) glmS L‑glutamine/D‑fructose‑6‑phosphate aminotransferase
Rejected: Frequency below/above cutoff threshold.

ACCTGTTTGCCAAGGTGCGCAACGTGCAGATCGTTGCCTGTGGCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/6171436‑6171392
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcaAACAGGTCGGCAGCTTGCGGGCCAAAGGCCTGGACCAGCACGTTGTCCTTGCCCAGGCGGCCTTCAAGCGTGCGCTGAACAACGCTTGGCTGCTCGTGAATTTCCTTGAGCATGAAGTGACGGTAAGCGCCCTTGTCGGCAGCCTCG  >  NC_002947/6171430‑6171576
                                                                                                                                                                                                
ACCTGTTTGCCAAGGTGCGCAACGTGCAGATCGTTGCCTGTGGCAAACAGGTCGGCAGCTTGCGGGCCAAAGGCCTGGACCAGCACGTTGTCCTTGCCCAGGCGGCCTTCAAGCGTGCGCTGAACAACGCTTGGCTGCTCGTGAATTTC                                             >  5:176550/1‑149
           AAGGTGCGCAACGTGCAGATCGTTGCCTGTGGCAAACAGGTCGGCAGCTTGCGGGCCAAAGGCCTGGACCAGCACGTTGTCCTTGCCCAGGCGGCCTTCAAGCGTGCGCTGAACAACGCTTGGCTGCTCGTGAATTTCC                                            <  7:229658/139‑1
           AAGGTGCGCAACGTGCAGATCGTTGCCTGTGGCAAACAGGTCGGCAGCTTGCGGGCCAAAGGCCTGGACCAGCACGTTGTCCTTGCCCAGGCGGCCTTCAAGCGTGCGCTGAACAACGCTTGGCTGCTCGTGAATTTCC                                            >  8:229658/1‑139
                                           CAAACAGGTCGGCAGCTTGCGGGCCAAAGGCCTGGACCAGCACGTTGTCCTTGCCCAGGCGGCCTTCAAGCGTGCGCTGAACAACGCTTGGCTGCTCGTGAATTTCCTTGAGCATGAAGTGACGGTAAGCGCCCTTGTCGGCAGCCTCG  <  1:182025/149‑1
                                           CAAACAGGTCGGCAGCTTGCGGGCCAAAGGCCTGGACCAGCACGTTGTCCTTGCCCAGGCGGCCTTCAAGCGTGCGCTGAACAACGCTTGGCTGCTCGTGAATTTCCTTGAGCATGAAGTGACGGTAAGCGCCCTTGTCGGCAGCCTCG  <  2:268533/149‑1
                                           CAAACAGGTCGGCAGCTTGCGGGCCAAAGGCCTGGACCAGCACGTTGTCCTTGCCCAGGCGGCCTTCAAGCGTGCGCTGAACAACGCTTGGCTGCTCGTGAATTTCCTTGAGCATGAAGTGACGGTAAGCGCCCTTGTCGGCAGCCTCG  >  6:149101/1‑149
                                                                                                                                                                                                
ACCTGTTTGCCAAGGTGCGCAACGTGCAGATCGTTGCCTGTGGCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/6171436‑6171392
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcaAACAGGTCGGCAGCTTGCGGGCCAAAGGCCTGGACCAGCACGTTGTCCTTGCCCAGGCGGCCTTCAAGCGTGCGCTGAACAACGCTTGGCTGCTCGTGAATTTCCTTGAGCATGAAGTGACGGTAAGCGCCCTTGTCGGCAGCCTCG  >  NC_002947/6171430‑6171576

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.