New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1604785 =34 (1.370)4 (0.170) 4/256 2.6 10.6% coding (739/993 nt) lsrD autoinducer 2 import system permease protein
?NC_000913 1604815 = 35 (1.490)coding (769/993 nt) lsrD autoinducer 2 import system permease protein
Rejected: Frequency below/above cutoff threshold.

ACCTGATTTGGCACTCCTGCCATTTGCAAACCTTGTTGCAAATATCCCACTAATAAAACCGCAATGGCGGTGCCGATAATGGAACCGGAACCACCATAAATATTGGCCCCGCCAAGCACCACGGCGGTGATGGCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1604921‑1604785
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcgggGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCGGTGCGCTACTTATCG  >  NC_000913/1604815‑1604953
                                                                                                                                                                                                                                                                                     
ACCTGATTTGGCACTCCTGCCATTTGCAAACCTTGTTGCAAATATCCCACTAATAAAACCGCAATGGCGGTGCCGATAATGGAACCGGAACCACCATAAATATTGGCCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTAT                                                                                                                                 <  1:107479/149‑1
                                                                                                           CCCGCCAAGCACCACGGCGGTGATGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCG                                                                                             >  1:3604/1‑78
                                                                                                           CCCGCCAAGCACCACGGCGGTGATGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCG                                                                                             <  2:3604/78‑1
                                                                                                                  AGCACCACGGCGGTGATGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCGGT               <  2:107479/149‑1
                                                                                                                                   TGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCGGTGCGCTAgatcggg  >  3:43093/1‑139
                                                                                                                                   TGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCGGTGCGC           <  4:43093/137‑1
                                                                                                                                                                                                                                                                                     
ACCTGATTTGGCACTCCTGCCATTTGCAAACCTTGTTGCAAATATCCCACTAATAAAACCGCAATGGCGGTGCCGATAATGGAACCGGAACCACCATAAATATTGGCCCCGCCAAGCACCACGGCGGTGATGGCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1604921‑1604785
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcgggGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCGGTGCGCTACTTATCG  >  NC_000913/1604815‑1604953

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