New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078057 (0.250)21 (0.870) 20/192 0.2 80.8% coding (305/630 nt) stfP e14 prophage, uncharacterized protein
?NC_000913 = 1209602 4 (0.170)pseudogene (18/501 nt) stfE pseudogene, e14 prophage, side tail fiber protein fragment family,Phage or Prophage Related

TGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATC  >  NC_000913/1209476‑1209733
                                                                                                                              |                                                                                                                                   
tGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                                                                                         <  1:599018‑M2/139‑13 (MQ=255)
               tGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                                                                          >  2:222080/1‑139 (MQ=255)
                     ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                                                                                    <  1:458618‑M2/139‑34 (MQ=255)
                        cGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                                                                                                 >  2:458618‑M2/1‑103 (MQ=255)
                         gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgca                                                                                                >  2:430155‑M2/1‑102 (MQ=255)
                               tcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCTTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                                                                                          <  1:41228‑M2/139‑44 (MQ=255)
                                   cagcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatc                                                                                      >  2:599018‑M2/1‑92 (MQ=255)
                                           aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                                              >  2:366792/1‑139 (MQ=255)
                                                        tGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgg                                                                 >  2:49671‑M2/1‑71 (MQ=255)
                                                               cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaacccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttg                                                          >  2:45170‑M2/1‑64 (MQ=255)
                                                               cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttg                                                          >  2:247801‑M2/1‑64 (MQ=255)
                                                                      ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattct                                                            <  2:271176‑M2/130‑74 (MQ=255)
                                                                      ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattct                                                            >  1:271176‑M2/1‑57 (MQ=255)
                                                                                 aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTTCTTTTACCAGGGGGATTTAACATGCTTATTgg                                        <  1:302578/139‑1 (MQ=255)
                                                                                    ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcgg                                                       >  2:420530‑M2/1‑43 (MQ=255)
                                                                                    ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcgg                                                       <  1:420530‑M2/121‑79 (MQ=255)
                                                                                           ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcgg                              >  2:598773‑M2/1‑36 (MQ=255)
                                                                                           ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcgg                              <  1:598773‑M2/139‑104 (MQ=255)
                                                                                                cGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagt                         <  2:57741‑M2/139‑109 (MQ=255)
                                                                                                             ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCATGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAg            <  2:395695/139‑1 (MQ=255)
                                                                                                                  aGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccc       >  2:41228‑M2/1‑13 (MQ=255)
                                                                                                                      aaTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgc                                                         <  1:502866‑M2/85‑77 (MQ=255)
                                                                                                                      aaTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgc                                                         >  2:502866‑M2/1‑9 (MQ=255)
                                                                                                                       aTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctc  <  2:323219‑M2/139‑132 (MQ=255)
                                                                                                                           agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata                                                                            >  2:505538‑M2/1‑4 (MQ=255)
                                                                                                                           agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata                                                                            <  1:505538‑M2/61‑58 (MQ=255)
                                                                                                                              |                                                                                                                                   
TGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATC  >  NC_000913/1209476‑1209733

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