New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =8 (0.130)44 (0.780) 35/252 0.7 87.1% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 6 (0.110)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACT  >  NC_000913/1209491‑1209754
                                                                                                                                |                                                                                                                                       
gccattccccctccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctccccatccatCGTTTATGAAAATga                                                                                                                           <  2:85685‑M2/15‑1 (MQ=255)
 gcattccacccccactgaatgagtactgattgctcaataccaccaatcctaaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccacCGTGTACGAAAATgaa                                                                                                                          >  2:197959‑M2/128‑143 (MQ=255)
    ttccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaagaa                                                                                                                       <  2:509134‑M2/19‑1 (MQ=255)
    ttccacccccactgaatgattactgattgctcaataccaccaatcagaaccgttttagctatgcaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaagaa                                                                                                                       <  2:510262‑M2/19‑1 (MQ=255)
           cccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                <  1:231835‑M2/26‑1 (MQ=255)
            ccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgagcttatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                                                                                               >  1:432741‑M2/117‑143 (MQ=255)
                  aatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAACATGAAGAAATAACAAGCAAATg                                                                                                         >  1:1108369‑M2/111‑143 (MQ=255)
                     gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgggttcagcccccgctccccaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                                      >  1:839449‑M2/108‑143 (MQ=255)
                       ttactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGCcatc                                                                                                    >  2:842037‑M2/106‑143 (MQ=255)
                          ctgattgctcactaccaccattcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctccccaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                                                                 <  2:54753‑M2/41‑1 (MQ=255)
                                gcgcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                                                                                           <  1:496819‑M2/47‑1 (MQ=255)
                                         cactcatcataagtggtatagctatgaaacccgtgcgtcctaatgtagccgtggctgcattgagttttgctccttctcccaaaccaaTGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAccacg                                                                                  <  2:207968‑M2/56‑3 (MQ=255)
                                         caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTGTGGAAATGAGGAAATAACAAGCAAATGGGATCATTCCTGGTTTTACAAgg                                                                                  >  1:960985‑M2/88‑143 (MQ=255)
                                          accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                 >  1:979291‑M2/87‑143 (MQ=255)
                                            caatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                               >  1:316006‑M2/85‑143 (MQ=255)
                                              ttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCCACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                                                                             >  1:962954/1‑143 (MQ=255)
                                                 ataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa                                                                          >  1:656584‑M2/80‑143 (MQ=255)
                                                  taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGGAGAAATAACAAGCAACTGGCATCATTCATGCATTGACACAGGGGAGTTaa                                                                         >  2:527606‑M2/79‑143 (MQ=255)
                                                     ccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                                                                      >  2:718060‑M2/76‑143 (MQ=255)
                                                      gccTGTTGGCGGAGTGGCTGAAGGCGACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAGAAGGGGGCTTTAACagg                                                                     >  2:548874/1‑141 (MQ=255)
                                                      cggtatagccatgaaaccggtgcgtcctaatgttgccccggcttcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                                     <  2:505921‑M2/69‑1 (MQ=255)
                                                               tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                            >  2:765483‑M2/66‑143 (MQ=255)
                                                               tatgaaaccggtgcgtcctaatgtagcagttgctgcattgagttttgcgccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                            <  1:591511‑M2/78‑1 (MQ=255)
                                                                  gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCAACATTCCTACCTTTCCCAGGGGGATTTAACATGCTTATTGGCTAt                                                         >  2:282870‑M2/63‑143 (MQ=255)
                                                                   aaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATg                                                        >  1:462503‑M2/62‑143 (MQ=255)
                                                                     accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                                                      >  2:471235‑M2/60‑143 (MQ=255)
                                                                         gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTGACCAGGGGGATTTAAAATGCTTATTAGCCGTGTAcgcg                                                  >  2:717625‑M2/56‑143 (MQ=255)
                                                                            ggCCATGGAACCGGGACACCAAAAGGAAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                               <  2:115576/143‑1 (MQ=255)
                                                                            ctccctaacgtaacggttgctgcattgagttttgctccttatcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                               <  2:734138‑M2/91‑1 (MQ=255)
                                                                            cgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                               >  1:85036‑M2/53‑143 (MQ=255)
                                                                              tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGATATGTACGCGTCTca                                             >  1:470113‑M2/51‑143 (MQ=255)
                                                                               cctaatgtagccgtagctgcatttagcattgccacttcgcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATCCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCCATCTCCGCGTCTcaa                                            >  2:518246‑M2/50‑143 (MQ=255)
                                                                                 taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAGAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGACTATGTACGAGTATcaaca                                          >  2:1043689‑M2/48‑143 (MQ=255)
                                                                                     aCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTTGCTATGTACGCGTATCAACAAATg                                      >  2:323758/1‑143 (MQ=255)
                                                                                            ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                               <  2:299775‑M2/107‑1 (MQ=255)
                                                                                            ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTGATTGGCTATGTACGCGTATCAACAAATGACCAGaa                               <  2:168585‑M2/107‑1 (MQ=255)
                                                                                            cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                               <  1:637265/143‑1 (MQ=255)
                                                                                                tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAAcaca                           <  2:114932‑M2/111‑1 (MQ=255)
                                                                                                 gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                          >  2:975192‑M2/32‑143 (MQ=255)
                                                                                                 aCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAacatgg                          >  1:72211/1‑140 (MQ=255)
                                                                                                   attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                        <  1:334660‑M2/114‑1 (MQ=255)
                                                                                                   attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                        <  2:878033‑M2/114‑1 (MQ=255)
                                                                                                   attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                        <  1:471235‑M2/114‑1 (MQ=255)
                                                                                                       agttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCGGGGGGCTCTCACATGCTTATTGGCTATGTACGCGTATCAACAAATGCCCAGAACACAGATCTAc                    >  1:839514‑M2/26‑143 (MQ=255)
                                                                                                       agttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTAc                    >  1:823378‑M2/26‑143 (MQ=255)
                                                                                                       agttttgctccccctcccaaaccacCTTTGATGAAAATGAAGAAATAACACGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTAc                    <  1:625384‑M2/118‑1 (MQ=255)
                                                                                                          tttgctccttctcccaaaccaaCGTTTATGAAAATATAGAAATAACAAGAAAATGGCATCATTCCTGCTTTTACCAGGTGGATATAACATGCGTATTGGCTATGTACGCGTATCAACAAATGACCAGAACCCAGATCTACAAc                 >  2:1068770‑M2/23‑143 (MQ=255)
                                                                                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTa              >  2:161891‑M2/20‑143 (MQ=255)
                                                                                                              ctccttctcccaaaccacCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa             <  1:891694‑M2/125‑1 (MQ=255)
                                                                                                               tccttctcccaatccaaCTTTTATTAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt            <  2:644963‑M2/126‑1 (MQ=255)
                                                                                                                   tctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcgc        <  1:752733‑M2/130‑1 (MQ=255)
                                                                                                                     tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGCTTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTg      >  1:508354‑M2/12‑143 (MQ=255)
                                                                                                                      cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGa     >  1:555795‑M2/11‑143 (MQ=255)
                                                                                                                      cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGa     >  1:405557‑M2/11‑143 (MQ=255)
                                                                                                                       ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa    >  2:1037195‑M2/10‑143 (MQ=255)
                                                                                                                       ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCAGaa    >  2:1043588‑M2/10‑143 (MQ=255)
                                                                                                                         aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACt  <  1:513053‑M2/136‑1 (MQ=255)
                                                                                                                                |                                                                                                                                       
TGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACT  >  NC_000913/1209491‑1209754

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