New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =8 (0.190)30 (0.810) 27/252 0.3 86.8% coding (290/630 nt) stfP e14 prophage, uncharacterized protein
?NC_000913 1209619 = 2 (0.050)pseudogene (1/501 nt) stfE pseudogene, e14 prophage, side tail fiber protein fragment family,Phage or Prophage Related

GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913/1209477‑1209760
                                                                                                                                              |                                                                                                                                              
gcagatgcttttgcggcattccacccccactgaatgattactgattgctaaatacatctaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                                                                                                >  1:613230‑M2/143‑143 (MQ=255)
      gcttttgcggcattccactcccattgaatgattactgattgctcaatatcaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                                                                                          <  2:38994‑M2/7‑1 (MQ=255)
        ttttgcggcattccacccccactgaatgattactgattgctcaataccaccaagcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctcaatgcaggtatgctccttctccccgaccgacgttgctca                                                                                                                                        >  1:195136‑M2/135‑138 (MQ=255)
            gcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcatttagttttgctccttctcccaaaccaaCGTTTATGAAAAt                                                                                                                                    <  1:534489‑M2/13‑1 (MQ=255)
                             ctgaatgattactgattgctcaataccaccaatcataaccggtatagccctacggtcgaatccgcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCaa                                                                                                                  <  1:530264‑M2/30‑1 (MQ=255)
                              tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgccttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCaaa                                                                                                                 >  2:492731‑M2/113‑143 (MQ=255)
                                aatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATg                                                                                                               >  2:461303‑M2/111‑143 (MQ=255)
                                 ttgattactgattgcgcaataccacgaatcataaccgatatagctatgaaacctttgcttcctaatgtatccgttgattcattaaattttactccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATgg                                                                                                              <  2:643775‑M2/34‑1 (MQ=255)
                                               ctcaataccaccaatcataaccggtatagctatgaaaccagtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaacccgCGGATATGAAAATGAAGAAATC‑agaaccaaatggcaccgatcatgctg                                                                                                >  2:211271‑M2/96‑131 (MQ=255)
                                                  aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTa                                                                                             >  1:708879‑M2/93‑143 (MQ=255)
                                                          caatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaacgaaCGTTTATGAAAATGAAGAAATC‑ACAAGCAACTGGCATCATTCCTGCTTATACCCggggg                                                                                     >  2:102806‑M2/85‑143 (MQ=255)
                                                                taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                               >  2:675549‑M2/79‑143 (MQ=255)
                                                                taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATAATTCCTGCTGTTACAAGGgggatctca                                                                               >  1:76334‑M2/79‑139 (MQ=255)
                                                                    cggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAGTA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                                           <  1:412167‑M2/69‑1 (MQ=255)
                                                                      gtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATG‑GGAAATATACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                                                                         >  1:529955‑M2/73‑143 (MQ=255)
                                                                       tatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                                                                        >  1:174649‑M2/72‑143 (MQ=255)
                                                                          tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                                                                     <  1:299948/143‑1 (MQ=255)
                                                                            ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACCAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                   >  1:632767‑M2/67‑143 (MQ=255)
                                                                              atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                                 <  1:542567‑M2/79‑1 (MQ=255)
                                                                                gaaaccggtgcttcctaatgtcgccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTCGCTAt                                                               <  2:604316‑M2/81‑1 (MQ=255)
                                                                                  aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                             >  1:77364‑M2/61‑143 (MQ=255)
                                                                                             cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTCTGTACGCGTATcaa                                                  >  1:275428‑M2/50‑143 (MQ=255)
                                                                                              ctaacggagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                 <  1:635092‑M2/95‑1 (MQ=255)
                                                                                                 atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                              <  1:688389‑M2/98‑1 (MQ=255)
                                                                                                    cagcagttgctgcattgagttttgttgcttcgcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAAGGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa                                           <  2:449222‑M2/101‑1 (MQ=255)
                                                                                                     agccgttgcggcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGCAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                          <  1:609834‑M2/102‑1 (MQ=255)
                                                                                                        cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGGATCAACAAATGACCAg                                       >  1:665281‑M2/39‑143 (MQ=255)
                                                                                                          ttgctgcattcagttttgctccttctcccaaaccaaCGTTCATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCGTTCCTGCTTTTACCAGAGGGATTTAACATGCTGATTCCCTATGTACGATAATCACcaaatgaagtgac                                     >  2:207104‑M2/37‑137 (MQ=255)
                                                                                                            gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                                   >  1:175901‑M2/35‑143 (MQ=255)
                                                                                                               aCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                                >  1:600122/1‑143 (MQ=255)
                                                                                                                  ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATc                             >  1:679595‑M2/29‑143 (MQ=255)
                                                                                                                  ttgagttttgatccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCCTCATTCCTGCATTTAACAGAGGGATTTAACATGCTTATTGGCTATGTCCGCGTATCAACAAACGACCACAACACAGTTc                             >  2:100842‑M2/29‑143 (MQ=255)
                                                                                                                              ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTCTTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACACCGTAata                 >  2:180404‑M2/17‑142 (MQ=255)
                                                                                                                                    cccaaaccaaCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACACCGTCATGCGCTGa           >  1:196967‑M2/11‑143 (MQ=255)
                                                                                                                                            tgcGTCTATTAAAATGAAGAAACA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAGTTAACATGCGTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAGCTACAACGTAATGCGCTGAACTGTGCa   <  2:611415/141‑1 (MQ=255)
                                                                                                                                             aCGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAg  <  2:174649‑M2/142‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATA‑ACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913/1209477‑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.