New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4602596 =155 (1.040)8 (0.060) 4/168 NT 5.2% coding (263/771 nt) yjjP DUF1212 family inner membrane protein
?NC_000913 4602621 = 146 (1.040)coding (238/771 nt) yjjP DUF1212 family inner membrane protein
Rejected: Frequency below cutoff threshold.

GCCATAGTGCTGACAACTATTAAAGATGGGCAATGCCTGACATCGACACGTAAAAATCACGATCGCGGCATTAATATGCATGTGGTGACTGAAGTCCAGCACAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4602699‑4602596
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cacatGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTT  >  NC_000913/4602621‑4602719
                                                                                                                                                                                                           
GCCATAGTGCTGACAACTATTAAAGATGGGCAATGCCTGACATCGACACGTAAAAATCACGATCGCGGCATTAATATGCATGTGGTGACTGAAGTCCAGC                                                                                                         <  2:1166561/100‑1
 CCATAGTGCTGACAACTATTAAAGATGGGCAATGCCTGACATCGACACGTAAAAATCACGATCGCGGCATTAATATGCATGTGGTGACTGAAGTCCAGCA                                                                                                        >  1:1828108/1‑100
 CCATAGTGCTGACAACTATTAAAGATGGGCAATGCCTGACATCGACACGTAAAAATCACGATCGCGGCATTAATATGCATGTGGTGACTGAAGTCCAGCA                                                                                                        <  2:2494666/100‑1
                                                                   GCATTAATATGCATGTGGTGACTGAAGTCCAGCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAG                                      >  1:550282/1‑100
                                                                        AATATGCATGTGGTGACTGAAGTCCAGCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTC                                 >  1:499426/1‑100
                                                                          nATGCATGTGGTGACTGAAGTCCAGCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAG                               >  2:3268193/2‑100
                                                                            TGCATGTGGTGACTGAAGTCCAGCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCA                             >  1:1266607/1‑100
                                                                            TGCATGTGGTGACTGAAGTCCAGCACATGCATATTAATGCCGCGACCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCA                             >  2:1677240/1‑100
                                                                            TGCATGTGGTGACTGAAGTCCAGAACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCA                             >  2:4149418/1‑100
                                                                             GCATGTGGTGACTGAAGTCCAGCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCAC                            >  1:275612/1‑100
                                                                             GCATGTGGTGACTGAAGTCCAGCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCAC                            >  2:1408873/1‑100
                                                                                                   CACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAA      <  1:1654943/100‑1
                                                                                                     CATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACT    <  1:3063463/100‑1
                                                                                                       TGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTT  >  2:782411/1‑100
                                                                                                                                                                                                           
GCCATAGTGCTGACAACTATTAAAGATGGGCAATGCCTGACATCGACACGTAAAAATCACGATCGCGGCATTAATATGCATGTGGTGACTGAAGTCCAGCACAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4602699‑4602596
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cacatGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTT  >  NC_000913/4602621‑4602719

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