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

GCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913/1209486‑1209760
                                                                                                                                     |                                                                                                                                             
tttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaa                                                                                                                                      >  1:74629‑M2/134‑143 (MQ=255)
   gcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAt                                                                                                                                   <  1:102481‑M2/13‑1 (MQ=255)
     tGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                                                                                                                 <  1:307373/143‑1 (MQ=255)
           ccaccgccactgcatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctcatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                           <  2:158469‑M2/21‑1 (MQ=255)
               gggTATCCGTCAGCAGAAAAAGCGGCACCATTGCAATTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                                                                                                                       <  1:48726/143‑1 (MQ=255)
                cccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                      <  1:272061‑M2/26‑1 (MQ=255)
                      cctcACCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCATACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                                                                                <  1:297996/141‑1 (MQ=255)
                         tgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGCTTCTGAAAATGAAGACAGAACacacaaacgcc                                                                                                             >  2:245922‑M2/109‑139 (MQ=255)
                          gattactgattgctcaataccaccaatcataaccggtatagctatgagaccggtgcgtcctaatgtagccgttgctgcattccgttttgctccttctcccaaccccaCGTTTATGAAAATAAAGAAATAACAAACAAATGCca                                                                                                            >  1:537676‑M2/108‑143 (MQ=255)
                             tactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                                                                         >  1:182277‑M2/105‑143 (MQ=255)
                                tgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCCTTc                                                                                                      >  1:576354‑M2/102‑143 (MQ=255)
                                tgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                      >  1:566371‑M2/102‑143 (MQ=255)
                                 aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                                                                     >  1:703961/1‑143 (MQ=255)
                                       tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccatctcccaaaccaaCGTTTAAGAAAAGTAGGTAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                               >  1:362767‑M2/95‑143 (MQ=255)
                                        caataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                                                                              <  2:254305‑M2/50‑1 (MQ=255)
                                             ccaccaaccataaccggtatagctatgacaccggtgcgtcctaacgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                         <  1:23576‑M2/55‑1 (MQ=255)
                                               accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                       <  2:333716‑M2/57‑1 (MQ=255)
                                                ccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaccccaCCTTTATGAAAATTAAGAAATAACAATCAAATTGCCTCATTACTCCTTCtaccaaaga                                                                                      >  1:309462‑M2/86‑139 (MQ=255)
                                                       taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgaggtttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                               <  1:20775‑M2/65‑1 (MQ=255)
                                                          ccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                                                                            <  2:157498‑M2/68‑1 (MQ=255)
                                                               atagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGCTTTAACATGCTTa                                                                       >  1:137789‑M2/71‑143 (MQ=255)
                                                                         aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                             <  2:62468‑M2/83‑1 (MQ=255)
                                                                             ggtgcgtcctgatgttgccggtgttgcattgagttttgctccttctacgaagctaaCGTTTATGAAAATGAAGAAATTACAAGTAAATGGTATCATTTCTGTTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                                                         <  2:230379‑M2/87‑1 (MQ=255)
                                                                                 cgtccgagtgcagccgttgcggcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAACTAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                     <  2:407166‑M2/91‑1 (MQ=255)
                                                                                     ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGCAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGCGATTTACCATGCTTATTGGCTCTGCCCGCGTctcagc                                                 >  1:186622‑M2/49‑141 (MQ=255)
                                                                                       agtgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACCAGCAAATGGCATCATTCCTGCTTTTACCAGGGGTATTTAACATGCTTATTGACTATGTACGAGTATcaacaa                                               >  2:322498‑M2/47‑143 (MQ=255)
                                                                                                 ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                     <  2:323606‑M2/107‑1 (MQ=255)
                                                                                                 tggctggagtgagttttgctccttcccccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                     <  1:110871‑M2/107‑1 (MQ=255)
                                                                                                     tgcattgaattgtgctccttctcccaaaccacCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAAcaca                                 <  2:233134‑M2/111‑1 (MQ=255)
                                                                                                        attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                              <  2:566371‑M2/114‑1 (MQ=255)
                                                                                                         ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATc                             <  2:300384‑M2/115‑1 (MQ=255)
                                                                                                         ttgaggtgtgctcctgcgccaaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATc                             <  1:191636‑M2/115‑1 (MQ=255)
                                                                                                              ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACAAAAATATACaa                        >  2:400696‑M2/24‑143 (MQ=255)
                                                                                                                  gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTCCAACGTa                    >  1:22036‑M2/20‑143 (MQ=255)
                                                                                                                    tccttctcccaaaccaaCGTTTCTGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGCCCAGCACGCAGAGCTGCAACGTAAt                  >  1:279089‑M2/18‑143 (MQ=255)
                                                                                                                    tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGTGGATTTAACATGCTTATTGGCTTTGTACGCGTATCAACAAATTACCCGAGCACGGATCTACAACgaaaa                  >  1:28116‑M2/18‑142 (MQ=255)
                                                                                                                    tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGTTTTTTTAAGCTTTCAACCGAAGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt                  >  1:711526‑M2/18‑143 (MQ=255)
                                                                                                                      cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACTAGAACACAGATCTACAACGTAATgc                >  2:328780‑M2/16‑143 (MQ=255)
                                                                                                                         ctcccaacccaaCTGTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt             <  1:245922‑M2/131‑1 (MQ=255)
                                                                                                                            ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa          >  2:160656‑M2/10‑143 (MQ=255)
                                                                                                                                    aCGTTTATGAAAATGAAGAAATAACAAGCAAATGTCATCATTCCTGCTTTTACCAGGGGGATTTAACATTCTTATGAGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAg  <  2:244808‑M2/142‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                             
GCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913/1209486‑1209760

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