New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	1207790 =	6 (0.170)	17 (0.530)	16/244	NT	71.9%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_hsa_tpiA	1209619 =	8 (0.250)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_hsa_tpiA/1209479‑1209767                                                                                                                                             |                                                                                                                                                     agatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagtttttctccttctcccaaaccaaCGTTTATGa                                                                                                                                              >  1:28892‑M2/141‑149 (MQ=255)     gcttttgcggcattccacccccactgaatgattactgattgctcaataccaccgatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccacaccaaCGTTTATGAAAAt                                                                                                                                          <  2:667164‑M2/13‑1 (MQ=255)        gCCTTGGCAAGTTACGGGTATTCTGCAGAAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                                                                                                                       <  2:253174/149‑1 (MQ=255)            ttGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaagaac                                                                                                                                   >  2:474116/1‑148 (MQ=255)                             tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                                  <  1:392533‑M2/37‑1 (MQ=255)                                            gCGGCACCATTGCATTTCAGGCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCCGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTGATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                                                                   <  1:504250/149‑1 (MQ=255)                                                   cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                                            <  1:248382/149‑1 (MQ=255)                                                       accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                        >  2:535573‑M2/87‑149 (MQ=255)                                                       accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgccccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                        >  1:367687‑M2/87‑149 (MQ=255)                                                         caatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                      <  1:688886‑M2/65‑1 (MQ=255)                                                                aaccggtgtagctatgaaaccggtgggtcctaatgtagccgttgctgcattgagttttgctccttatcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                                                                               <  2:537289‑M2/72‑1 (MQ=255)                                                                    cgtGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAGGTTTATGAAAATGAAGAAATAACAAGCAAATGGTATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                                           <  2:16137/147‑1 (MQ=255)                                                                       atagctatgaaaccggtgcgtcctaatgtagccgttgctgcatagagttttactgcttcgcccaaaccaaCGTTTATGAAAATTAATAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                                        <  2:487398‑M2/79‑1 (MQ=255)                                                                                 aactggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt                                                              <  2:610089‑M2/89‑1 (MQ=255)                                                                                       tgcgtcctattgtagtctttgctgtattgagttttgctccttctccaaaactaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATGCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                        <  1:70234‑M2/95‑1 (MQ=255)                                                                                          gtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                                     <  1:66079‑M2/98‑1 (MQ=255)                                                                                            cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAGTGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATg                                                   >  1:504010‑M2/50‑149 (MQ=255)                                                                                             ctaatgtagccgtagctgcattgagtcttgctccttctcccaaaccaaCGATTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa                                                  >  1:563987‑M2/49‑149 (MQ=255)                                                                                                   gagccgttggcgcattgacttttgctccttctcccaacccaaAGTTTATGAAAATGAATAAATAACAAGCAAAGGGCATCATTCCTGCTCTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                            <  2:80490‑M2/107‑1 (MQ=255)                                                                                                   ccGGGACACCAACAGGTAATGCAGGGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                            <  1:193951/149‑1 (MQ=255)                                                                                                      ggACACCAACAGGTAATGCAGAGCCGTCTCCGAAACCAACGTCTATGAAAATGACGAAATAACAACAAAATGGCATCATTCCTGCTTTTACCAGGGGTATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAacac                                         <  2:170886/149‑1 (MQ=255)                                                                                                         cacCAACTGGGAATTCAGGTGCGTCTCCGGAACCCTCCTTGAGGAAAATGAAGGAATAAGCAGCAAATGGCATCCTTCCTGCTTTTACCAGGGGGATTTAACCTGCTTTTTGGCTATGTACGCGTATCTACAAATGACCATAACACAGa                                      <  2:247758/149‑1 (MQ=255)                                                                                                               cattgagttttgctccttctcccacaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCGTTTACCAGGGGGATTTAACATGCTTATTGGCTCTGTACGCGTAGCAACCAATGCCCCGAACACAAATCTACa                                >  2:171277‑M2/31‑149 (MQ=255)                                                                                                                     gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa                          <  1:354980‑M2/125‑1 (MQ=255)                                                                                                                      ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACCAATCTACACCGTAAt                         >  2:76507‑M2/24‑149 (MQ=255)                                                                                                                           ctccttctcccaaaccaaCGTTTATGAAACTGAAGAAATAACAAGCAAATGGCATCATTCCTGCTCGTACCAGGGGAATTTAACATGCCTATTGGCTATGTACGCGCATCAACAAATGACCAGAACACAGCTCTACAACGTaacgcgat                    >  2:461384‑M2/19‑147 (MQ=255)                                                                                                                           ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCCGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACACATCTCCAACCTAATGCCCt                    >  1:167707‑M2/19‑149 (MQ=255)                                                                                                                                  gcgtaaaaCAACGTTTATGAAAATGAAGATATAACAAGCAAATGGCGTCATTCCTGCTTTTACCAGGGGGATTTAACATACTTACTGGCTATGTACGCGTATCAACAAAAGGCCAGAACACAGATCTACAGCGTAATGCTCTGATCtgt             <  2:122540/145‑1 (MQ=255)                                                                                                                                             cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGa  <  1:475962/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                     CGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_hsa_tpiA/1209479‑1209767

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