New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =3 (0.090)16 (0.540) 15/232 NT 80.7% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 5 (0.170)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAA  >  NC_000913/1209480‑1209752
                                                                                                                                           |                                                                                                                                     
gctgcttttgcggcattccaccccctctgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                                                                                    <  2:221475‑M2/4‑1 (MQ=255)
    aaGCCTGTGCCAGTGCCGGGGATGCTTCAGCAGAAAAAGCGTCACCAGTGCATTTCAGCCAGCCTGTTGGCCGAGCGGCTGAAGGCCACTGAACCGGAACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATg                                                                                                                                <  1:250991/143‑1 (MQ=255)
             cattccacccccactgaatgattactgattgctcaataccaccaagcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                                                                                       <  1:608170‑M2/17‑1 (MQ=255)
                tccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTACGAAAATGAAGaaa                                                                                                                    >  2:524598‑M2/124‑143 (MQ=255)
                 acacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgccccttctccccaacaaaCCATTATGAAAATGACGAAAt                                                                                                                   >  2:285528‑M2/123‑143 (MQ=255)
                    cccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                                                                <  2:388962‑M2/24‑1 (MQ=255)
                     ccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                                                                                               <  1:330001‑M2/25‑1 (MQ=255)
                       ccactgactgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                                                                                             <  2:29957‑M2/27‑1 (MQ=255)
                           tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagcagttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                                                                                         <  2:324302‑M2/31‑1 (MQ=255)
                                   actgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaaccaaCGTTTATGAAAATGAAGAAATAACCAGCAAAAGGcaaca                                                                                                 >  2:449828‑M2/105‑143 (MQ=255)
                                       aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACCGGTACTGCAGAGCCTTCTCCCAAACCAACGCTTATGAAAATGAAGACACCACAAGCAACTGGcatgagcca                                                                                             >  2:482337/1‑137 (MQ=255)
                                         ttgctcaataccaccaatcataaccggtatagctatgaacccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaatccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                           <  1:428654‑M2/45‑1 (MQ=255)
                                               aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                                                                                     <  1:324896‑M2/51‑1 (MQ=255)
                                                        aatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                                                                            >  1:488663‑M2/84‑143 (MQ=255)
                                                                    tatagccatgaacccggtgcgtcctaatgtagcccttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                                                                <  2:96339‑M2/72‑1 (MQ=255)
                                                                           atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                         <  1:380917‑M2/79‑1 (MQ=255)
                                                                                               aaCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                     <  1:520578/143‑1 (MQ=255)
                                                                                                       cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                             <  1:133324/143‑1 (MQ=255)
                                                                                                            aCAGGTAATGCAGGACCTTCTCCCAAAATAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAGCATGCTTATTGTCTATGTACGCGTATCAACAAATGACCAGAACACAg                        <  2:444862/143‑1 (MQ=255)
                                                                                                              attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                      <  1:371563‑M2/114‑1 (MQ=255)
                                                                                                                    ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAACAACAAATGACCAGAACACAGATCTACaa                >  2:506473‑M2/24‑143 (MQ=255)
                                                                                                                     tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAAGAAATGACCAGAACACGGATCTACAAc               >  2:123302‑M2/23‑143 (MQ=255)
                                                                                                                            cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATCCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAACTGACCAGAACACAGATCTACAacgtacacc        >  2:120343‑M2/16‑139 (MQ=255)
                                                                                                                                  ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa  >  2:295048‑M2/10‑143 (MQ=255)
                                                                                                                                  ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa  >  2:589209‑M2/10‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                     
GGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAA  >  NC_000913/1209480‑1209752

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.