New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 1207805	13 (0.380)	17 (0.570)	16/248	1.0	64.8%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_hsa_tpiA	= 1209602	7 (0.230)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTA  >  NC_000913_3_hsa_tpiA/1209454‑1209735                                                                                                                                                     |                                                                                                                                      cacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                                                                                       <  2:177203/149‑1 (MQ=255) cacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                                                                                       <  1:515163/149‑1 (MQ=255)       aaTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttct                                                                                                                                 >  1:157115‑M2/1‑143 (MQ=255)        aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                                                                                                <  1:326484‑M2/149‑8 (MQ=255)               aaTTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacc                                                                                                                         <  2:360269‑M2/149‑15 (MQ=255)               aaTTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacc                                                                                                                         <  2:252892‑M2/149‑15 (MQ=255)                           ggATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGAAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgag                                                                                                             >  1:252892‑M2/1‑123 (MQ=255)                                aGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                                                                                        <  2:272421/149‑1 (MQ=255)                                         cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAc                                                                                               >  2:480909/1‑149 (MQ=255)                                            ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttcgcccaaaccaaggttttcgagagccgtttgcacagtgc                                                                                            <  1:15312‑M2/149‑44 (MQ=255)                                                     ttcttcAGCAGAATAAGCGNCNCNANTGCAGTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGTGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                                                                                   >  1:335640‑M2/1‑97 (MQ=255)                                                                          aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattct                                                              <  2:181592‑M2/149‑74 (MQ=255)                                                                                  aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggctt                                                      <  1:530393‑M2/149‑82 (MQ=255)                                                                                    ttCATCGATCCTGTTGGGGTAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaatcttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgccgcttaa                                                    >  2:576526‑M2/1‑66 (MQ=255)                                                                                               tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgc                                         <  1:127115‑M2/149‑95 (MQ=255)                                                                                                    gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat                                                                             >  2:108114‑M2/1‑50 (MQ=255)                                                                                                    gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat                                                                             <  1:108114‑M2/108‑59 (MQ=255)                                                                                                           ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCATAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttggtatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaa                             >  2:251257‑M2/1‑43 (MQ=255)                                                                                                           ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa                             <  1:480909/149‑1 (MQ=255)                                                                                                            gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaag                            <  1:274607‑M2/149‑108 (MQ=255)                                                                                                            gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACGAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                           <  1:394278/134‑1 (MQ=255)                                                                                                            gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACGAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                           >  2:394278/1‑134 (MQ=255)                                                                                                               gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaa                         >  2:291623‑M2/1‑39 (MQ=255)                                                                                                                      aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa                  >  1:197188/1‑149 (MQ=255)                                                                                                                      aCGGAACCGGGACAACAACAGGTAATGCCGCGCCTACTCCCAAGCCCAAGTTAACGCAAGTGAATAAATAACAATCAAATGGGATCATTCCTTCGGTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa                  <  2:502337/149‑1 (MQ=255)                                                                                                                                     cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc          >  1:570204‑M2/1‑17 (MQ=255)                                                                                                                                     cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc          <  2:570204‑M2/142‑126 (MQ=255)                                                                                                                                     cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca                                               <  1:375130‑M2/105‑89 (MQ=255)                                                                                                                                     cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca                                               >  2:375130‑M2/1‑17 (MQ=255)                                                                                                                                      aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcca  >  1:181592‑M2/1‑16 (MQ=255)                                                                                                                                                     |                                                                                                                                      CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTA  >  NC_000913_3_hsa_tpiA/1209454‑1209735

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