New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1647054 =22 (0.840)3 (0.120) 3/276 2.8 14.2% intergenic (+317/+120) flxA/intK Qin prophage; uncharacterized protein/pseudogene, integrase fragment, Qin prophage;Phage or Prophage Related
?NC_000913 1647120 = 16 (0.660)intergenic (+383/+54) flxA/intK Qin prophage; uncharacterized protein/pseudogene, integrase fragment, Qin prophage;Phage or Prophage Related
Rejected: Frequency below cutoff threshold.

TTACCTAAAAAGAAGCCTGATCCCGATAAAAAAGATTGAAACCAAGATCGATTCGGTTGAGTGCATATCCATTCATAGGGTAGATTCTTAAGTCGCGTTTCTGGTGTTCATTTTCGGGTGGTTTGTTACTTGTTTTACCGGGGATATGCCAGAAACGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1647212‑1647054
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccagaaacgcgACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGA  >  NC_000913/1647120‑1647264
                                                                                                                                                                                                                                                                                                                
TTACCTAAAAAGAAGCCTGATCCCGATAAAAAAGATTGAAACCAAGATCGATTCGGTTGAGTGCATATCCATTCATAGGGTAGATTCTTAAGTCGCGTTTCTGGTGTTCATTTTCGGGTGGTTTGTTACTTGTTTTACCGGGGATATGCC                                                                                                                                                            <  2:122391/150‑1
                                   TTGAAACCAAGATCGATTCGGTTGAGTGCATATCCATTCATAGGGTAGATTCTTAAGTCGCGTTTCTGGTGTTCATTTTCGGGTGGTTTGTTACTTGTTTTACCGGGGATATGCCAGAAACGCGACTTAAGAATCTACCCTATGAATGGAT                                                                                                                        <  1:410852/151‑1
                                         CCAAGATCGATTCGGTTGAGTGCATATCCATTCATAGGGTAGATTCTTAAGTCGCGTTTCTGGTGTTAATTTTCGGGTGGTTTGTTACTTGTTTTACCGGGGATATGCCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCAC                                                                                                                  <  1:157974/151‑1
                                                   TTCGGTTGAGTGCATATCCATTCATAGGGTAGATTCTTAAGTCGCGTTTCTGGTGTTCATTTTCGGGTGGTTTGTTACTTGTTTTACCGGGGATATGCCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAAT                                                                                                        <  1:317220/151‑1
                                                                                                                                                       GAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGG    >  1:139370/1‑151
                                                                                                                                                         AACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGA  <  1:142869/151‑1
                                                                                                                                                         AACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGA  <  2:139403/151‑1
                                                                                                                                                                                                                                                                                                                
TTACCTAAAAAGAAGCCTGATCCCGATAAAAAAGATTGAAACCAAGATCGATTCGGTTGAGTGCATATCCATTCATAGGGTAGATTCTTAAGTCGCGTTTCTGGTGTTCATTTTCGGGTGGTTTGTTACTTGTTTTACCGGGGATATGCCAGAAACGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1647212‑1647054
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccagaaacgcgACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGA  >  NC_000913/1647120‑1647264

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