New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078054 (0.100)18 (0.620) 12/98 0.7 87.8% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 2 (0.070)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA  >  NC_000913/1209535‑1209669
                                                                   |                                                                   
aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                     >  2:867495/1‑68 (MQ=255)
aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                     >  2:15363/1‑68 (MQ=255)
             tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                        <  1:814010‑M2/68‑14 (MQ=255)
                tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaa                                                     <  1:242886‑M2/68‑17 (MQ=255)
                      aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtttt                                               >  2:143002‑M2/1‑46 (MQ=255)
                               aGGCCACGGAACCGGGACACCAACAGTTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                      >  1:305789‑M2/1‑37 (MQ=255)
                               aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                      >  2:439736‑M2/1‑37 (MQ=255)
                               aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                      >  1:464227‑M2/1‑37 (MQ=255)
                               aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                      >  2:11972‑M2/1‑37 (MQ=255)
                                       gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcaca                              <  1:664588‑M2/68‑40 (MQ=255)
                                         aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                            >  1:1162324‑M2/1‑27 (MQ=255)
                                              gACACCAACAGGTAATGCAGAGCCATCTCACAAACAAACGTTTATGAAAATGAAGAAATAACAAGCaa                       >  2:1021981/1‑68 (MQ=255)
                                                  ccAACTGGTAATGCAGTGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccga                   >  2:1134241‑M2/1‑18 (MQ=255)
                                                    aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                 <  1:1288779‑M2/68‑53 (MQ=255)
                                                         gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata            <  1:441976‑M2/68‑58 (MQ=255)
                                                         gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata            <  2:477570‑M2/68‑58 (MQ=255)
                                                         gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata            <  1:978142‑M2/68‑58 (MQ=255)
                                                             tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc        >  1:443825‑M2/1‑7 (MQ=255)
                                                              gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcc       <  1:711392‑M2/68‑63 (MQ=255)
                                                              gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc       >  2:736987/1‑68 (MQ=255)
                                                               cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcca      <  1:147322‑M2/68‑64 (MQ=255)
                                                                agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa     <  1:662122‑M2/68‑65 (MQ=255)
                                                                agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa     <  1:1408297‑M2/68‑65 (MQ=255)
                                                                   gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacg  <  1:1269041‑M2/68‑68 (MQ=255)
                                                                   |                                                                   
ATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA  >  NC_000913/1209535‑1209669

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