New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =10 (0.290)15 (0.540) 14/128 1.1 60.0% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 12 (0.430)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGC  >  NC_000913/1209541‑1209696
                                                                              |                                                                             
gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGt                                                                             >  2:266210/1‑81 (MQ=255)
   agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTa                                                                          <  1:266169/81‑1 (MQ=255)
     ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATg                                                                        <  1:266170/81‑1 (MQ=255)
           gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAt                                                                   >  2:1014509‑M2/68‑80 (MQ=255)
           gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                                  <  2:265866‑M2/14‑1 (MQ=255)
            gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                                  <  1:266172/80‑1 (MQ=255)
               gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                               >  2:266215/1‑80 (MQ=255)
                   accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                           >  2:266218‑M2/60‑80 (MQ=255)
                   accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                           >  1:266205‑M2/60‑80 (MQ=255)
                          cgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                    <  1:265865‑M2/28‑1 (MQ=255)
                                  tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                           <  2:265858‑M2/37‑1 (MQ=255)
                                  aaCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCCAATGGcat                                           <  1:266175/81‑1 (MQ=255)
                                  aaCCGGGACACCAACAGGTAATGCAGAGCCTACTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                           <  2:266174/81‑1 (MQ=255)
                                    tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                          <  2:265864‑M2/38‑1 (MQ=255)
                                     cGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                        <  1:266176/81‑1 (MQ=255)
                                         gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                    <  1:265873‑M2/44‑1 (MQ=255)
                                         acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                    >  1:266207/1‑81 (MQ=255)
                                          cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                   <  1:266178/81‑1 (MQ=255)
                                                 attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                            >  2:266220‑M2/30‑81 (MQ=255)
                                                 attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                            <  1:265855‑M2/52‑1 (MQ=255)
                                                    gagttttgctccttctccaaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                          <  1:265857‑M2/54‑1 (MQ=255)
                                                     agttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                        <  1:265860‑M2/56‑1 (MQ=255)
                                                     aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                        >  1:266221/1‑81 (MQ=255)
                                                         ttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                    >  2:266213‑M2/22‑81 (MQ=255)
                                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                   <  2:266181‑M2/61‑1 (MQ=255)
                                                            ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt                  <  2:266182‑M2/62‑1 (MQ=255)
                                                            agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt                  >  1:266209/1‑80 (MQ=255)
                                                              ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa               >  1:266202‑M2/17‑81 (MQ=255)
                                                                     ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt        <  1:266183‑M2/72‑1 (MQ=255)
                                                                      caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa       >  1:266214‑M2/9‑81 (MQ=255)
                                                                           caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc  <  1:266184‑M2/78‑1 (MQ=255)
                                                                            aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc  >  2:266226‑M2/3‑80 (MQ=255)
                                                                              |                                                                             
GCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGC  >  NC_000913/1209541‑1209696

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