New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP000730 = 21219785 (0.060)6 (0.080) 6/262 12.0 55.8% coding (132/177 nt) USA300HOU_1999 hypothetical bacteriophage protein
?CP000730 = 2122008 NA (NA)coding (102/177 nt) USA300HOU_1999 hypothetical bacteriophage protein
Rejected: Coverage evenness skew score above cutoff.

TTATCCAAAAGATAGAGTTGTCAGTTTAGATCATGTTTTGTATTTCGACATTAACGTTAAAGGTGAAAAGAATAACGATAGTCCTTATCCATCTTAAAATCCATAGTGCTTATAGTTATTTGCTATTTGTTGACAGTTATCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTATTC  >  CP000730/2121840‑2122118
                                                                                                                                          |                                                                                                                                            
ttATCCAAAAGATAGAGTTGTCAGTTTAGATCATGTTTTGTATTTCGACATTAACGTTAAAGGTGAAAAGAATAACGATAGTCCTTATCCATCTTAAAATCCATAGTGCTTATAGTTATTTGCTATTTGTTGACAGTTATc                                                                                                                                            >  1:867937/1‑141 (MQ=31)
                                                                                                                                       gTTATCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTa     <  1:411581/141‑1 (MQ=11)
                                                                                                                                          aTCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTATtc  >  2:209691/1‑141 (MQ=32)
                                                                                                                                          aTCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTATtc  >  2:666518/1‑141 (MQ=32)
                                                                                                                                          aTCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTATtc  <  2:674997/141‑1 (MQ=32)
                                                                                                                                          |                                                                                                                                            
TTATCCAAAAGATAGAGTTGTCAGTTTAGATCATGTTTTGTATTTCGACATTAACGTTAAAGGTGAAAAGAATAACGATAGTCCTTATCCATCTTAAAATCCATAGTGCTTATAGTTATTTGCTATTTGTTGACAGTTATCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTATTC  >  CP000730/2121840‑2122118

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