New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	15 (0.510)	11 (0.430)	11/232	NT	44.7%	coding (290/630 nt)	ycfK	e14 prophage; protein StfP
	?	NC_000913	1209619 =	14 (0.550)					pseudogene (37/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment

TTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAA  >  NC_000913/1209489‑1209752                                                                                                                                   |                                                                                                                                      tctGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCCGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                                                                                           <  2:188726/141‑1 (MQ=255)    tcattccacccccacggaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                        <  1:485154‑M2/16‑1 (MQ=255)              cccacggaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                              <  2:80223‑M2/26‑1 (MQ=255)                      aatgattactgattgctcaatagcatcaatcataaccggtatagctatgaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                                                                                      <  2:436203‑M2/34‑1 (MQ=255)                       cagcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACGCCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGCTTATGAAAATGAAGAAATAACAAGCAAGTGGc                                                                                                     >  2:210427/1‑143 (MQ=255)                       cagcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                                                     >  1:271161/1‑143 (MQ=255)                             aaaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                                                               >  1:352045/1‑143 (MQ=255)                                attgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATACCAACCAAATCGCACCATTCCt                                                                                            >  2:135436‑M2/100‑143 (MQ=255)                                 aGCGACACCATTGCATTCCAGCCAGCCTGTTGGCGGAGTGCCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACGAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                           <  1:552165/143‑1 (MQ=255)                                        ataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                                                    <  2:113745‑M2/52‑1 (MQ=255)                                         cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc                                                                                   <  1:319729/143‑1 (MQ=255)                                           ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                 >  1:398606/1‑143 (MQ=255)                                               aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                             <  2:28195/143‑1 (MQ=255)                                                 ttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAACTGGCATCATTCGTGCTTTTACAAGGGCGAt                                                                           >  2:65481/1‑143 (MQ=255)                                                     taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                       <  1:117415‑M2/65‑1 (MQ=255)                                                         cggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgcgctttctcccaaaccaaCGTTTATGAATATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                                   <  2:287485‑M2/69‑1 (MQ=255)                                                                            gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg                                                >  2:556313‑M2/56‑143 (MQ=255)                                                                                  cACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                          <  1:140116/143‑1 (MQ=255)                                                                                  cACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                          <  1:515459/143‑1 (MQ=255)                                                                                      gAACCGGGAGACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                      <  1:224661/143‑1 (MQ=255)                                                                                              acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTACCTACGCCTAGCAGAAAAGCACCAGa                              >  2:61981/1‑143 (MQ=255)                                                                                                 ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGCTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                           <  1:153978/143‑1 (MQ=255)                                                                                                       ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCCGGGGGATTTAACATGCTTATTGGCTATGTACGCGGATCAACAAATGACCAGACCACAGATc                     >  1:426034‑M2/29‑143 (MQ=255)                                                                                                            ttttgctccttctcccacaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACaa                <  1:102342‑M2/120‑1 (MQ=255)                                                                                                             tttgctccttctcccaaaccaaCGGTTAGGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGACTTAACATGCTTACTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc               >  2:484038‑M2/23‑143 (MQ=255)                                                                                                                        tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTCCCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTg    >  1:162776‑M2/12‑143 (MQ=255)                                                                                                                        tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTg    <  2:271161‑M2/132‑1 (MQ=255)                                                                                                                          ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa  <  2:221007‑M2/134‑1 (MQ=255)                                                                                                                                   |                                                                                                                                      TTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAA  >  NC_000913/1209489‑1209752

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.