New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 185887328 (0.880)7 (0.230) 6/242 NT 20.5% coding (2010/2259 nt) PP_1661 dehydrogenase
?NC_002947 = 1858902 27 (0.870)coding (2039/2259 nt) PP_1661 dehydrogenase

TGCTCAAGCGCACCGTGCAGGTGCAGGCCAACCAGGCATGGCAGGGCACGGGCATGCAGGTCAATGGGCAGAGCTGGCAGTTGGCGGTATGCACGGGCGGGCAGTGGACGGCGAACCCGGCCACCGGCATGGTTGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/1858737‑1858873
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcTTGGCGATCAACCTTGGGTTGCCCGCCGCGCCAACCATGCCGGTG  <  NC_002947/1858902‑1858858
                                                                                                                                                                                       
TGCTCAAGCGCACCGTGCAGGTGCAGGCCAACCAGGCATGGCAGGGCACGGGCATGCAGGTCAATGGGCAGAGCTGGCAGTTGGCGGTATGCACGGGCGGGCAGTGGACGGCGAACCCGGCCACCGGCATAGCCGGGCTTGGCGATCAAC                                   >  1:716284/1‑150
                                CAGGCATGGCAGGGCACGGGCATGCAGGTCAATGGGCAGAGCTGGCAGTTGGCGGTATGCACGGGCGGGCAGTGGACGGCGAACCCGGCCACCGGCATAGCCGGGCTTGGCGATCAACCTTGGGTTGCCCGCCGCGCCAACCATGCCGG  <  2:716284/150‑2
                                 AGGCATGGCAGGGCACGGGCATGCAGGTCAATGGGCAGAGCTGGCAGTTGGCGGTATGCACGGGCGGGCAGTGGACGGCGAACCCGGCCACCGGCATAGCCGGGCTTGGCGATCAACCTTGGGTTGCCCGCCGCGCCAACCATGCCGGTG  <  1:68276/150‑1
                                          AGGGCACGGGCATGCAGGTCAATGGGCAGAGCTGGCAGTTGGCGGTATGCACGGGCGGGCAGTGGACGGCGAACCCGGCCACCGGCATGGTTGGCTTGGCGATCAACCTTGGGTTGCCCGCCGC                  >  1:548620/1‑124
                                          AGGGCACGGGCATGCAGGTCAATGGGCAGAGCTGGCAGTTGGCGGTATGCACGGGCGGGCAGTGGACGGCGAACCCGGCCACCGGCATGGTTGGCTTGGCGATCAACCTTGGGTTGCCCGCCGC                  <  2:548620/124‑1
                                                         AGGTCAATGGGCAGAGCTGGCAGTTGGCGGTATGCACGGGCGGGCAGTGGACGGCGAACCCGGCCACCGGCATGGTTGGCTTGGCGATCAACCTTGGGTTGCCCGCCGCGCCAACCATG        >  1:292551/1‑119
                                                         AGGTCAATGGGCAGAGCTGGCAGTTGGCGGTATGCACGGGCGGGCAGTGGACGGCGAACCCGGCCACCGGCATGGTTGGCTTGGCGATCAACCTTGGGTTGCCCGCCGCGCCAACCATG        <  2:292551/119‑1
                                                                                                                                                                                       
TGCTCAAGCGCACCGTGCAGGTGCAGGCCAACCAGGCATGGCAGGGCACGGGCATGCAGGTCAATGGGCAGAGCTGGCAGTTGGCGGTATGCACGGGCGGGCAGTGGACGGCGAACCCGGCCACCGGCATGGTTGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/1858737‑1858873
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcTTGGCGATCAACCTTGGGTTGCCCGCCGCGCCAACCATGCCGGTG  <  NC_002947/1858902‑1858858

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