New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780539 (0.720)5 (0.110) 5/178 2.6 14.5% coding (305/630 nt) stfP e14 prophage, uncharacterized protein
?NC_000913 = 1209602 26 (0.570)pseudogene (18/501 nt) stfE pseudogene, e14 prophage, side tail fiber protein fragment family,Phage or Prophage Related
Rejected: Frequency below/above cutoff threshold.

CACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1207706‑1207805
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttggtttgggagaaggCTCTGCATTACCTGTTGGTGTCCCGGTTCCGTGGCCTTCAGCCACTCCGCCAACAGGCTGGCTGAAATGCAATGGTGCCGCTTTTTCTGCTGAAGAATACCCGGAACTGGCAAAGGCTTATCC  <  NC_000913/1209602‑1209480
                                                                                                                                                                                                                               
CACTTCGCGCGCTGCTGTTAAGCCGCCAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAagag                                                                                                                           <  2:833710/102‑5
      GCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGCTCTGCATTACCTGTTGGTGTCCCGGTTCCGTGGCCTTCAGCCAC                                                                                >  2:263150/1‑139
                  TAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGCTCTGCATTACCTGTTGGTGTCCCGGTTCCGTGGCCTTCAGCCACTCCGCCAACAGG                                                                    <  2:118831/139‑1
                         CAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGCTCTGCATTACCTGTTGGTGTCCCGGTTCCGTGGCCTTCAGCCACTCCGCCAACAGGCTGGCTG                                                             <  1:263150/139‑1
                                                      ATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTGTGGGAGAAGGCTCTGCATTACCTGTTGGTGTCCCGGTTCCGTGGCCTTCAGCCACTCCGCCAACAGGCTGGCTGAAATGCAATGGTGCCGCTTTTTCTGCTGA                                >  2:1194095/1‑139
                                                            CTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                                                                             >  1:974501/1‑40
                                                            CTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                                                                             <  2:974501/40‑1
                                                                         TCTCGAAAACCTTGGTTTGGGAGAAGGCTCTGCATTACCTGTTGGTGTCCCGGTTCCGTGGCCTTCAGCCACTCCGCCAACAGGCTGGCTGAAATGCAATGGTGCCGCTTTTTCTGCTGAAGAATACCCGGAACTGGCA             >  2:1188203/1‑139
                                                                                    TTGGTTTGGGAGAAGGCTCTGCATTACCTGTTGGTGTCCCGGTTCCGTGGCCTTCAGCCACTCCGCCAACAGGCTGGCTGAAATGCAATGGTGCCGCTTTTTCTGCTGAAGAATACCCGGAACTGGCAAAGGCTTATCC  <  1:922481/139‑1
                                                                                                  GGCTCTGCATTACCTGTTGGTGTCCCGGTTCCGTGGCCTTCAGCCACTCCGCCAACAGGCTGGCTGAAATGCAATGGTGCCGCTTTTTCTGCTGAAGAATACCCGGA                    <  1:1164830/107‑1
                                                                                                  GGCTCTGCATTACCTGTTGGTGTCCCGGTTCCGTGGCCTTCAGCCACTCCGCCAACAGGCTGGCTGAAATGCAATGGTGCCGCTTTTTCTGCTGAAGAATACCCGGA                    >  2:1164830/1‑107
                                                                                                                                                                                                                               
CACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1207706‑1207805
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttggtttgggagaaggCTCTGCATTACCTGTTGGTGTCCCGGTTCCGTGGCCTTCAGCCACTCCGCCAACAGGCTGGCTGAAATGCAATGGTGCCGCTTTTTCTGCTGAAGAATACCCGGAACTGGCAAAGGCTTATCC  <  NC_000913/1209602‑1209480

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