New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 1758576 =8 (0.490)3 (0.200) 3/238 NT 24.6% coding (1590/1686 nt) PP_1563 terminase large subunit
?NC_002947 1758613 = 11 (0.730)coding (1627/1686 nt) PP_1563 terminase large subunit

CCGGTACTGACCTGGATGGCCTCGAATGTCGTGGCGAAGCTGGACGTGAAGGACAACATCTACCCCAACAAGGAGCGCCCGGAGAACAAGATCGACGGCATGGTGAGCTTGATCACCGGCTGCGCCGTGGCCATCAAGCTCGGTATCGACGACTCCGGCCACTTCGATGACTTTCTTGCCAGCCCGATCGTGGTTGGCTAACGGGAAAACCTATGAAAACTGGCCTGATCATCTTTCT  >  NC_002947/1758472‑1758709
                                                                                                        |                                                                                                                                      
ccGGTACTGAGCTGAATGGCCTCGAATGTCGTGGCGAAGCTGGACGTGAAGGACAACATCTACCCCAACAAGGAGCGCCGGGAGGACAAGATCGACGGCATGGTGAGCTTGATCACCGGCTGCGCCATGCCCATCACGCTCGCTATcgaa                                                                                          >  2:172021/1‑149 (MQ=255)
ccGGTACTGACCTGGATGGCCTCGAATGTCGTGGCGAAGCTGGACGTGAAGGACAACATCTACCCCAACAAGGAGCGCCCGGAGAACAAGATCGACGGCATGGTGAGCTTGATCACCGGCTGCGCCGTGGCCATCAAGCTCGGTATcgac                                                                                          <  1:198208/150‑1 (MQ=255)
aagcaccagaaagatgatcaggccagttttcataggttttcccgttagccaaccacgatcgggctggcaagaaagtcatcgaagtggccggagtcgtcgataccGAGCTTGATCACCGGCTGTGCCGTGGCCATCAAGCTCg                                                                                                  <  1:113025‑M1/38‑1 (MQ=255)
aagcaccagaaagatgatcaggccagttttcataggttttcccgttagccaaccacgatcgggctggcaagaaagtcatcgaagtggccggagtcgtcgataccGAGCTTGATCACCGGCTGCGCCGTGGCCATCAAGCTCg                                                                                                  >  2:113025‑M1/105‑142 (MQ=255)
                                             gTGAAGGACAACATCTACCCCAACAAGGAGCGCCCGGAGAACAAGATCGACGGCATGGTGAGCTTGATCACCGGCTGCGCCGTGGCCATCAAGCTCGGTATCGACGACTCCGGCCACTTCGATGACTTTCTTGCCAGCCCGATCGTGGtt                                             >  1:157945/1‑150 (MQ=255)
                                                    acaacaTCTACCCCAACAAGGAGCGCCCGGAGAACAAGATCGACGGCATGGTGAGCTTGATCACCGGCTGCGCCGTGGCCATCAAGCTCGGTATCGACGACTCCGGCCACTTCGATGACTTTCTTGCCAGCCCGATCGTGGTTGGCTAAc                                      >  2:60505/1‑150 (MQ=255)
                                                     caacaTCTACCCCAACAAGGAGCGCCCGGAGAACAAGATCGACGGCATGGTGAGCTTGATCACCGGCTGCGCCGTGGCCATCAAGCTCGGTATCGACGACTCCGGCCACTTCGATGACTTTCTTGCCAGCCCGATCGTGGTTGGCTAACg                                     <  1:60505/150‑1 (MQ=255)
                                                            tACCCCAACAAGGAGCGCCCGGAGAACAAGATCGACGGCATGGTGAGCTTGATCACCGGCTGCGCCGTGGCCATCAAGCTCGGTATCGACGACTCCGGCCACTTCGATGACTTTCTTGCCAGCCCGATCGTg                                                >  1:437402/1‑132 (MQ=255)
                                                            tACCCCAACAAGGAGCGCCCGGAGAACAAGATCGACGGCATGGTGAGCTTGATCACCGGCTGCGCCGTGGCCATCAAGCTCGGTATCGACGACTCCGGCCACTTCGATGACTTTCTTGCCAGCCCGATCGTg                                                <  2:437402/132‑1 (MQ=255)
                                                                                        cggagtcgtcgataccGAGCTTGATGGCCACGGCGCAGCCGGTGATCAAGCTCGGTATCGACGACTCCGGCCACTTCGATGACTTTCTTGCCAGCCCGATCGTGGTTGGCTAACGGGAAAACCTATGAAAACTGGCCTGATCATCTTTCt  <  1:236404‑M1/134‑1 (MQ=255)
                                                                                        aGATCGACGGCATGGTGAGCTTGATCACCGGCTGCGCCGTGGCCATCAAGCTCGGTATCGACGACTCCGGCCACTTCGATGACTTTCTTGCCAGCCCGATCGTGGTTGGCTAACGGGAAAACCTATGAAAACTGGCCTGATCATCTTTCt  >  2:389736/1‑150 (MQ=255)
                                                                                                        |                                                                                                                                      
CCGGTACTGACCTGGATGGCCTCGAATGTCGTGGCGAAGCTGGACGTGAAGGACAACATCTACCCCAACAAGGAGCGCCCGGAGAACAAGATCGACGGCATGGTGAGCTTGATCACCGGCTGCGCCGTGGCCATCAAGCTCGGTATCGACGACTCCGGCCACTTCGATGACTTTCTTGCCAGCCCGATCGTGGTTGGCTAACGGGAAAACCTATGAAAACTGGCCTGATCATCTTTCT  >  NC_002947/1758472‑1758709

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.