New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5552330 =42 (0.760)3 (0.060) 3/242 NT 7.5% coding (1559/1566 nt) PP_4882 iron ABC transporter permease
?NC_002947 5552377 = 34 (0.650)intergenic (+40/+621) PP_4882/PP_5717 iron ABC transporter permease/hypothetical protein

ACTCAGGCTTGTGCCCAGGCCACATGGCGGCCTGGGCAAGCTTTCAAGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5552380‑5552330
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caagcctGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAGACCTTGGGAAACCCTCGAGATCTGAAA  >  NC_002947/5552377‑5552522
                                                                                                                                                                                                     
ACTCAGGCTTGTGCCCAGGCCACATGGCGGCCTGGGCAAGCTTTCAAGCCTGAGTGCCTAGACTTCATCACAAGTTCAGTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGG                                                  >  3:264591/1‑149
  TCAGGCTTGTGCCCAGGCCACATGGCGGCCTGGGCAAGCTTTCAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGT                                                <  3:213406/149‑1
  TCAGGCTTGTGCCCAGGCCACATGGCGGCCTGGGCAAGCTTTCAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGG                                                 >  4:213406/1‑148
                                            CAAGCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAGACCTTGGGAAACCCTCGAGATCT      <  5:95609/149‑1
                                               GCCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCA                                >  2:301536/1‑120
                                                CCTGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAACACTCTTTTCTTAGACCTTGGGACACCCTCGAGATCTGAAA  >  4:36834/1‑149
                                                                                                                                                                                                     
ACTCAGGCTTGTGCCCAGGCCACATGGCGGCCTGGGCAAGCTTTCAAGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5552380‑5552330
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caagcctGAGTGCCTAGACTTCATCACAAGTTCACTTTAAAACTCTCACTAGCCAGTTGAAAGCAAAGAACAAGGTCGTAGATACTACAAGCCCCACCCAAGAGGGTAACAAAACTCTTTTCATAGACCTTGGGAAACCCTCGAGATCTGAAA  >  NC_002947/5552377‑5552522

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