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.

CGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCGGTGCGCTACTTATCGTC  >  NC_000913/1604671‑1604955
                                                                                                                                                 |                                                                                                                                            
cGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCaa                                                                                                                                          >  3:153145/1‑149 (MQ=255)
        acctgatttggcactcctgccatttgcaaaccttgttgcaaatatcccactaataaaaccgcaatggcggtgccgataatggaaccggaaccaccataaatattggccccgccatcaccgccgtggtgcttggcgggGCCAATATTTAt                                                                                                                                   <  1:107479‑M2/12‑1 (MQ=255)
        aTGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGTGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTAt                                                                                                                                   <  5:7161/148‑1 (MQ=255)
                      ccATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTggttccggttc                                                                                                                     <  4:149519/148‑1 (MQ=255)
                            cGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATc                                                                                                              <  1:4893/149‑1 (MQ=255)
                                  tGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCAcc                                                                                                        >  5:134280/1‑149 (MQ=255)
                                       tCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATAGTGGTTCCGGTTCCATTATCGGCACCGCCAt                                                                                                   >  5:110753/1‑149 (MQ=255)
                                         tGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCAcc                                                                                                        <  4:158942/142‑1 (MQ=255)
                                         tGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCAcc                                                                                                        >  3:158942/1‑142 (MQ=255)
                                              tCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATAGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGtt                                                                                            <  6:110753/149‑1 (MQ=255)
                                                          tgctgGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTggttc                                                                                                                           <  4:79047/106‑1 (MQ=255)
                                                          tgctgGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTggttc                                                                                                                           >  3:79047/1‑106 (MQ=255)
                                                            ctgGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTttatta                                                                                      >  2:61781/1‑141 (MQ=255)
                                                            ctgGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTttatta                                                                                      <  1:61781/141‑1 (MQ=255)
                                                                    gtgttttGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCAGTTTTATTAGTGGGATATTTGcaa                                                                       >  1:5684/4‑148 (MQ=255)
                                                                        tttGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAgg                                                                  >  4:60453/1‑149 (MQ=255)
                                                                           ggTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGata                                                                              <  1:10381/134‑1 (MQ=255)
                                                                           ggTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGata                                                                              >  2:10381/1‑134 (MQ=255)
                                                                            gTTCAGCACGTTCCGATCTCGGTGCGTCGATTCTGATGACCGCAATCACAGCCGTGGTGCTTGGCGGGGCCAATATTAATGGAGGTTCCGCTTACATTATCGGCACCGCCATTGCGGTTGTATAAGTCGGATATTTGCAACAAGGttag                                                              >  2:123389/1‑147 (MQ=255)
                                                                                 gCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCa                                                                                                                   <  2:133909/91‑1 (MQ=255)
                                                                                 gCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCAGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAAt                                                         <  2:5684/149‑1 (MQ=255)
                                                                                  cACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTg                                                                                   >  1:26150/1‑122 (MQ=255)
                                                                                  cACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTg                                                                                   <  2:26150/122‑1 (MQ=255)
                                                                                   aCGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATgg                                                       <  3:60453/149‑1 (MQ=255)
                                                                                   aCGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATgg                                                       >  3:69348/1‑149 (MQ=255)
                                                                                     gTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCa                                                     >  6:59202/1‑149 (MQ=255)
                                                                                          gATCTCGGTGCGTCGTTTCTGATTCCCGACAAAACCGCCGTGGTGCTTGGCGGGGCAAATATTTATGGTGGTTCCGGTGCCATTATCGGCACCGCCATAGCGGTTTTATAAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAgg                                                >  6:106213/1‑148 (MQ=255)
                                                                                            tctcGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGc                                                                                                      >  3:100330/1‑93 (MQ=255)
                                                                                            tctcGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGc                                                                                                      <  4:100330/93‑1 (MQ=255)
                                                                                                ggTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGc                                              <  4:153145/145‑1 (MQ=255)
                                                                                                      tcgtTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAg                                      <  1:5968/147‑1 (MQ=255)
                                                                                                      tcgtTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAg                                      >  2:5968/1‑147 (MQ=255)
                                                                                                      tcgtTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGgt                                    <  6:134280/149‑1 (MQ=255)
                                                                                                      ataaatattggcCCCGCCAAGCACCACGGCGGTGATGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGgt                                    <  1:107479‑M1/137‑1 (MQ=255)
                                                                                                                   cccgccaagcaccacggcggtgatggcgggGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCg                                                                                               <  2:3604‑M2/48‑1 (MQ=255)
                                                                                                                   cccgccaagcaccacggcggtgatggcgggGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCg                                                                                               >  1:3604‑M2/31‑78 (MQ=255)
                                                                                                                          agcaccacggcggtgatggcgggGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCGGt                 <  2:107479‑M2/126‑1 (MQ=255)
                                                                                                                             ccgccgTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGtt                                                                                            >  5:60266/1‑71 (MQ=255)
                                                                                                                             ccgccgTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGtt                                                                                            <  6:60266/71‑1 (MQ=255)
                                                                                                                                         ttGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCGGTGCGCTACTTAtcgtc  >  4:144890/1‑149 (MQ=255)
                                                                                                                                          tggcgggGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCGgtgcgctagatcggg    >  3:43093‑M2/8‑139 (MQ=255)
                                                                                                                                          tggcgggGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCGGTgcgc             <  4:43093‑M2/130‑1 (MQ=255)
                                                                                                                                                 gCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCgg                  >  5:96753/1‑125 (MQ=255)
                                                                                                                                                 gCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCgg                  <  6:96753/125‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
CGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCATTATCGGCACCGCCATTGCGGTTTTATTAGTGGGATATTTGCAACAAGGTTTGCAAATGGCAGGAGTGCCAAATCAGGTGTCCAGCGCCCTTTCCGGTGCGCTACTTATCGTC  >  NC_000913/1604671‑1604955

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