New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =9 (0.170)31 (0.690) 27/166 0.9 80.0% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 8 (0.180)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

AGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATG  >  NC_000913/1209521‑1209718
                                                                                                  |                                                                                                   
tttctcaataccaccagtcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                                                     <  2:1023517‑M2/1‑1 (MQ=255)
    tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                                                                >  1:8305‑M2/95‑100 (MQ=255)
      aataccaccaatcatagccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATg                                                                                              <  2:1087763‑M2/8‑1 (MQ=255)
      aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATg                                                                                              <  2:1067295‑M2/8‑1 (MQ=255)
          ccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                                                          >  1:405997‑M2/89‑100 (MQ=255)
          ccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                                                          <  1:865841‑M2/12‑1 (MQ=255)
           caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAt                                                                                         <  2:490727‑M2/13‑1 (MQ=255)
           caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAt                                                                                         >  2:76609‑M2/88‑100 (MQ=255)
            accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                                                          >  2:352809‑M2/87‑98 (MQ=255)
                ttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                                                                    <  2:953202/100‑1 (MQ=255)
                  cataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAtgaagc                                                                                    <  2:1286427‑M2/18‑2 (MQ=255)
                  cataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                                                                  >  2:993906‑M2/81‑100 (MQ=255)
                               gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                     >  1:1239660‑M2/68‑100 (MQ=255)
                                 tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                   >  2:55579‑M2/66‑100 (MQ=255)
                                 tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                   >  2:221326‑M2/66‑100 (MQ=255)
                                  atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                  >  1:744764‑M2/65‑100 (MQ=255)
                                    gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                                                <  1:344190‑M2/38‑1 (MQ=255)
                                          gatgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                          <  2:791305‑M2/44‑1 (MQ=255)
                                           gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                         <  2:1195774/100‑1 (MQ=255)
                                               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                     >  2:1019831/1‑100 (MQ=255)
                                                 cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                                                   <  2:446960‑M2/51‑1 (MQ=255)
                                                   taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc                                                 <  2:1187800‑M2/53‑1 (MQ=255)
                                                   cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc                                                 <  2:1061890/100‑1 (MQ=255)
                                                    aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                >  2:702815‑M2/47‑100 (MQ=255)
                                                      tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATAAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACAAgg                                              >  2:1225674‑M2/45‑100 (MQ=255)
                                                       gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                             <  2:31544‑M2/57‑1 (MQ=255)
                                                        tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                            <  1:457934‑M2/58‑1 (MQ=255)
                                                        ccGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                            >  2:994887/1‑100 (MQ=255)
                                                               tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                     >  2:668154‑M2/36‑100 (MQ=255)
                                                                gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                    <  1:102443‑M2/66‑1 (MQ=255)
                                                                ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                    >  2:870963/1‑100 (MQ=255)
                                                                 ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                   <  2:744759‑M2/67‑1 (MQ=255)
                                                                   gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                 >  1:907695‑M2/32‑100 (MQ=255)
                                                                   gcattgagttttgctacttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                 >  2:997961‑M2/32‑100 (MQ=255)
                                                                    cattgagttttgctccttctcccaaaccaaCGTTGATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                                <  1:55580‑M2/70‑1 (MQ=255)
                                                                     attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                               >  1:116225‑M2/30‑100 (MQ=255)
                                                                       gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa                             >  2:15137/1‑100 (MQ=255)
                                                                             cAGAGCCTTCTACCAAACCAACGTTTATGAAACTGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGAGGATTTAAACTGATTATTGGc                        >  2:428813/1‑99 (MQ=255)
                                                                                 tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                   >  2:1028621‑M2/18‑100 (MQ=255)
                                                                                    ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcg                <  2:552346‑M2/86‑1 (MQ=255)
                                                                                    ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcg                <  1:951947‑M2/86‑1 (MQ=255)
                                                                                       tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt             >  1:80859‑M2/12‑100 (MQ=255)
                                                                                        cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa            >  2:153203‑M2/11‑100 (MQ=255)
                                                                                        cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa            >  1:991145‑M2/11‑100 (MQ=255)
                                                                                        cccAAATCTACGTCTACAATAATGAAGTAATAACAAGCAAATGGCATCATTCAACCTTTTACGAGGTGGATATAACATGCTTATTGGCTATGTACGTGTa            >  2:990257/1‑100 (MQ=255)
                                                                                               caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaaca      <  2:979854‑M2/96‑1 (MQ=255)
                                                                                                aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa    >  2:423853‑M2/3‑100 (MQ=255)
                                                                                                aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa    <  2:116223‑M2/98‑1 (MQ=255)
                                                                                                aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa    <  1:1170185‑M2/98‑1 (MQ=255)
                                                                                                 aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt   <  2:1239649‑M2/99‑1 (MQ=255)
                                                                                                  cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATg  <  2:390197/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
AGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATG  >  NC_000913/1209521‑1209718

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