New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi 1207790 =4 (0.100)27 (0.770) 23/250 0.6 87.7% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_pae_pgi 1209619 = 4 (0.110)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

AGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_pae_pgi/1209485‑1209764
                                                                                                                                      |                                                                                                                                                 
ttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaacccaaCGTTTATGAAAATga                                                                                                                                     <  1:494701‑M2/15‑1 (MQ=255)
ttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                                                                                                     <  2:302797‑M2/15‑1 (MQ=255)
 tttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcggggccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                                    <  1:650048‑M2/16‑1 (MQ=255)
    gcggcattccacccctactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaagaa                                                                                                                                 <  2:371849‑M2/19‑1 (MQ=255)
      ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                                                                                                              <  1:212967‑M2/6‑1 (MQ=255)
      ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                                                                                                              >  2:212967‑M2/129‑134 (MQ=255)
      ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                               >  2:512521‑M2/129‑149 (MQ=255)
                     ctgcatgattactgattgctcaataccaccaatcacaaccggtatagctatgaaaccggcgcgtcctaatctagccgttgctgcattgagttttgctccttctccccaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                                                <  2:74740‑M2/36‑1 (MQ=255)
                               gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                      >  1:251851/1‑149 (MQ=255)
                                        tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgatgcactgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                             >  1:479209‑M2/95‑149 (MQ=255)
                                          aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctactgtagccgttgctgcattgagttttgctccttcgcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                           <  1:58730‑M2/57‑1 (MQ=255)
                                          aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                           <  2:439359‑M2/57‑1 (MQ=255)
                                           tatctaccaatcataaccggtatagctatgaaaccgctgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAACGATGAAATAACACGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                          <  1:548567‑M2/58‑1 (MQ=255)
                                           ataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                          >  1:236601‑M2/92‑149 (MQ=255)
                                            cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                                         >  1:368832/1‑149 (MQ=255)
                                             accaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                                                                                        >  2:444600‑M2/90‑149 (MQ=255)
                                              ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                                                                                       >  1:539965/1‑149 (MQ=255)
                                                 ccaatcatatccggtataacgatgaaaccggtgcgtcctaatgtagccgttgctgaattgagttttgctccttctcccaaaccgaCGTGTATGAAAATGAAGAACTGACAAGCTAATGGCATCATTCCTGCTTTTAGGAGggggagatc                                                                                    >  2:548567‑M2/86‑145 (MQ=255)
                                                    atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctacttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCATGCTTTTACCATGGTGGTTTAACa                                                                                 >  1:126489‑M2/83‑149 (MQ=255)
                                                     tcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                                                                                >  2:517439‑M2/82‑149 (MQ=255)
                                                      cataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                                               <  2:300458‑M2/69‑1 (MQ=255)
                                                              ggatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATCAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                       >  2:640742‑M2/73‑149 (MQ=255)
                                                                  tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                                                   <  2:539965/149‑1 (MQ=255)
                                                                    ctatgaaaccggtgcgtcctaatgtagccgctgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                                 <  1:531450‑M2/83‑1 (MQ=255)
                                                                           accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa                                                          >  1:436125‑M2/60‑149 (MQ=255)
                                                                               gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                      >  1:358473‑M2/56‑149 (MQ=255)
                                                                                      ctaatgtagccgttgctgcattgagttttgctccttctcccaaacaaaCGTTTATGAAGATGAAGAAATAACAAGCAAATGGCATCATTTCTGCTTTTACGAGGGGGATTTAACATGCTTGTTGGCTTTGTTCGCGTTTCATCACATGa                                               >  2:14470‑M2/49‑149 (MQ=255)
                                                                                        aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCATGCTTTTACCAGGGGGATGTAACTTGATTATTGGCTATGTACGCGTATGAACAAATGAcc                                             >  2:77790‑M2/47‑149 (MQ=255)
                                                                                                        cattgagttttgctccgtctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACa                             >  1:604071‑M2/31‑149 (MQ=255)
                                                                                                         attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACaa                            >  1:168528‑M2/30‑149 (MQ=255)
                                                                                                                  tgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACACGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcg                   <  1:659279‑M2/129‑1 (MQ=255)
                                                                                                                         tctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACt            <  2:446941‑M2/136‑1 (MQ=255)
                                                                                                                           tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt          >  1:325760‑M2/12‑149 (MQ=255)
                                                                                                                            cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgtg         <  2:251851‑M2/139‑1 (MQ=255)
                                                                                                                                   caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATg  <  1:447787‑M2/146‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                 
AGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_pae_pgi/1209485‑1209764

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