New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780543 (1.210)12 (0.450) 9/98 0.9 33.1% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 16 (0.600)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA  >  NC_000913/1209535‑1209669
                                                                   |                                                                   
aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACTCCAACAGGTAATGCagag                                                                     >  1:44644/1‑68 (MQ=255)
aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                     <  2:238890/68‑1 (MQ=255)
 tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                    >  1:171310‑M2/1‑67 (MQ=255)
   tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGCGcct                                                                  <  2:238663‑M2/68‑4 (MQ=255)
   tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcct                                                                  <  2:480730‑M2/68‑4 (MQ=255)
       ccagccTGTTGTCGGAGTTGCTGAATGCAACGGAACTGTGACATCAACAGGTATTGCATATCCTtctc                                                              <  2:1038042/68‑1 (MQ=255)
       ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                              >  2:522621‑M2/1‑61 (MQ=255)
           ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaa                                                          <  1:1265571‑M2/68‑12 (MQ=255)
                     gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                >  2:1008275/1‑68 (MQ=255)
                                 gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                         <  1:330107‑M2/47‑13 (MQ=255)
                                 gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                         >  2:330105‑M2/1‑35 (MQ=255)
                                      ggAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATa                               <  2:1231442/68‑1 (MQ=255)
                                       gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                      <  1:658872‑M2/60‑32 (MQ=255)
                                       gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                      >  2:658870‑M2/1‑29 (MQ=255)
                                         aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgattgcacagt                            <  1:1101007‑M2/68‑42 (MQ=255)
                                           cGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                          >  2:1044185‑M2/1‑25 (MQ=255)
                                            gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcc                         <  1:920947‑M2/68‑45 (MQ=255)
                                               acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                      >  1:388773/1‑68 (MQ=255)
                                                       aGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga              <  1:320945‑M2/68‑56 (MQ=255)
                                                         gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata            <  1:1096757‑M2/68‑58 (MQ=255)
                                                         gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccga                   >  2:93718‑M2/1‑11 (MQ=255)
                                                         gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccga                   <  1:93718‑M2/61‑51 (MQ=255)
                                                         gTAATGCAGAGCCTTCTCCCAAACCAACTTTTATGAAAATGAAGAAATAACAATCAAATGGCATCttt            <  1:1299167/68‑1 (MQ=255)
                                                         gTAATGCAGAGCCTTCTCCCAAACCAACGTTAATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt            >  2:1283414/1‑68 (MQ=255)
                                                             tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg        >  2:133665/1‑68 (MQ=255)
                                                             tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg        >  2:331083/1‑68 (MQ=255)
                                                              gCAGAGCCTTCTCCCAAACCAACGTTTATGAATATGATGAATTAACTAGCAAATGGCATCATTCCTGc       <  2:771552/68‑1 (MQ=255)
                                                              gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc       <  2:853796/68‑1 (MQ=255)
                                                                agcgCCTTCTCCCAAACCACCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCAGCATTCCTGCtt     >  2:77433/4‑68 (MQ=255)
                                                                agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa     <  1:823559‑M2/68‑65 (MQ=255)
                                                                agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa     <  2:63690‑M2/68‑65 (MQ=255)
                                                                agagccttctcccaaaccaaggatttcgagagccgtttgcacagtgccaaccgatatgatatcgccaa     <  1:827496‑M2/68‑65 (MQ=255)
                                                                  agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAAAACAAGCAAATGGCATCATTCCTGCtttt   <  2:177906/68‑1 (MQ=255)
                                                                   gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                          <  1:1163347/44‑1 (MQ=255)
                                                                   gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                          >  2:1163346/1‑44 (MQ=255)
                                                                   gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa  <  2:200712/68‑1 (MQ=255)
                                                                   gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa  >  1:1000909/1‑68 (MQ=255)
                                                                   gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCATGCTTTTa  >  2:525750/1‑68 (MQ=255)
                                                                   |                                                                   
ATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA  >  NC_000913/1209535‑1209669

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