New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780530 (1.060)9 (0.420) 6/100 1.1 34.9% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 11 (0.510)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCT  >  NC_000913/1209538‑1209665
                                                                |                                                               
tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcct                                                              <  1:854668‑M2/68‑4 (MQ=255)
        ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaa                                                      <  1:988695‑M2/68‑12 (MQ=255)
               gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                               <  1:664880/68‑1 (MQ=255)
                  gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                            <  2:188799/68‑1 (MQ=255)
                     tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcg                                         <  1:97562‑M2/68‑25 (MQ=255)
                            agtccACGGAACCGTGACACCATCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                  <  2:176879/65‑1 (MQ=255)
                             ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                 >  1:432005/1‑68 (MQ=255)
                              gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                >  2:208614/1‑68 (MQ=255)
                                      aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                        >  1:245801‑M2/1‑27 (MQ=255)
                                      aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                        >  2:243145‑M2/1‑27 (MQ=255)
                                       ccGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtg                       >  2:757243‑M2/1‑26 (MQ=255)
                                           gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                        >  2:241931‑M2/1‑22 (MQ=255)
                                           gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                        <  1:241931‑M2/63‑42 (MQ=255)
                                           gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                   >  1:488196/1‑68 (MQ=255)
                                                 aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt             <  2:176605‑M2/68‑53 (MQ=255)
                                                    aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca          >  1:997163/1‑68 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata        <  2:234877‑M2/68‑58 (MQ=255)
                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata        <  1:855102‑M2/68‑58 (MQ=255)
                                                      gTAATGCAGAGCCTTCTCCCAATCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt        <  1:837368/68‑1 (MQ=255)
                                                          tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg    >  1:734483/1‑68 (MQ=255)
                                                           gCAGAGCCGTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc   <  2:591702/68‑1 (MQ=255)
                                                            cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt  >  2:929712/1‑68 (MQ=255)
                                                                |                                                               
TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCT  >  NC_000913/1209538‑1209665

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.