New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078054 (0.110)25 (0.780) 22/232 0.4 90.1% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 2 (0.060)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCT  >  NC_000913/1209472‑1209734
                                                                                                                                  |                                                                                                                                    
tATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcc                                                                                                                                    >  2:158545‑M2/1‑131 (MQ=255)
        ggATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctccca                                                                                                                            <  2:601461‑M2/133‑11 (MQ=255)
                       aGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                                                             >  1:5010‑M2/1‑108 (MQ=255)
                         ttCCGGGTATCCTTCAGCAGAAAACGCGCCACTATTGCATTGCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                           >  2:208608‑M2/1‑106 (MQ=255)
                           ccGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                                                                         <  1:48643/133‑1 (MQ=255)
                               gTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                                                                                     >  1:572743‑M2/1‑100 (MQ=255)
                                     ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttgtctagagccatttgcaca                                                                                               <  1:296256‑M2/133‑40 (MQ=255)
                                       cagcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccataccaaggtttttgatagccgtgtgcacagt                                                                                             <  1:183373‑M2/133‑42 (MQ=255)
                                                           ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc                                                                         >  1:345540‑M2/1‑72 (MQ=255)
                                                                tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaac                                                                    <  2:464760‑M2/133‑67 (MQ=255)
                                                                   cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccttccgatttgatatagccaaacgga                                                                 <  2:262898‑M2/133‑70 (MQ=255)
                                                                   cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgga                                                                 <  1:577584‑M2/133‑70 (MQ=255)
                                                                         gccTGTTGGCTGAGTCGCTGAAGGCCACTGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttg                                                           >  2:183373‑M2/1‑58 (MQ=255)
                                                                         gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttg                                                           <  1:466598‑M2/133‑76 (MQ=255)
                                                                            tGTTGGCGGAGTGGCTGAAGGCCACGGACCCGCGACACCAACAGGTAATTCACAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcgg                                                        >  2:236670‑M2/1‑55 (MQ=255)
                                                                            tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcgg                                                        >  2:425030‑M2/1‑55 (MQ=255)
                                                                              ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttg                                                           >  1:494758‑M2/1‑53 (MQ=255)
                                                                              ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttg                                                           <  2:494758‑M2/128‑76 (MQ=255)
                                                                                   ggAGTGGCTGAAGGCCACTGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaaca                                                 >  1:202940‑M2/1‑48 (MQ=255)
                                                                                        ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatagcgcaaaagggattcttccggcttaacagcagc                                            <  2:202940‑M2/133‑91 (MQ=255)
                                                                                          cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgc                                          >  2:466598‑M2/1‑41 (MQ=255)
                                                                                              aGGCCACGGACCCGGAACACCAACCGGTCATGAAGAGccttatcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcga                                      >  2:172176‑M2/1‑37 (MQ=255)
                                                                                               ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggcttcttgcggcttaacagcagcgcgcgaa                                     <  2:280629‑M2/133‑98 (MQ=255)
                                                                                                  cACGGAACCGGGACACCAACAGGTGATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                                  <  2:474079/133‑1 (MQ=255)
                                                                                                   aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgc                                 <  1:401797‑M2/133‑102 (MQ=255)
                                                                                                    cGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcg                                >  2:527106‑M2/1‑31 (MQ=255)
                                                                                                      gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggt                              >  1:149592‑M2/1‑29 (MQ=255)
                                                                                                                           aTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc         >  1:601461‑M2/1‑8 (MQ=255)
                                                                                                                                  gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcc  <  1:130146‑M2/133‑133 (MQ=255)
                                                                                                                                  |                                                                                                                                    
TATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCT  >  NC_000913/1209472‑1209734

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