New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_pgi	= 1207805	4 (0.110)	23 (0.740)	18/242	NT	80.0%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_hsa_pgi	= 1209602	8 (0.260)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGC  >  NC_000913_3_hsa_pgi/1209454‑1209746                                                                                                                                                     |                                                                                                                                                 cacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                                                                                                  >  1:219845/1‑149 (MQ=255)         tCAGGCAATCTATTTTGCGGATAAGCCTTTGCCAGTTCCGGCTATTCTTGAGCAGAAAAAGCGGCACCATTGAATTTCACCCAGCCACTTGGCTGAGCGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                                                                                                          >  2:204420‑M2/1‑141 (MQ=255)            ggCAATTTATTTGTCGGATAAGGCTTTGCCAGTTCCGGGTGTTCTTCAGCTGAAAAAGCGGCACCATTGGATTTCAGCCAGCCTGTTGGCGGAGTGGCGGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGgcttctcccag                                                                                                                                       <  2:625149‑M2/149‑12 (MQ=255)              cAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGGATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGATCGAACACAGTTAAGGCAGACGCTTCTccacaac                                                                                                                                     >  2:546332/1‑144 (MQ=255)                 tctacttgttgGATAAGCCGTCGCCAGTGCCGGGCATTCTTCAGCGGAAAAAACGGCACCATTGCAGTTCAGACAGCCTGTTGGCGGGGTGGCTTGAGGCCACGGCAACGGGAAACCAACCGGTAATGCAGAGCCTTCTCCCAAACCaa                                                                                                                                  <  2:444601/144‑1 (MQ=255)                  ttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                                                                                                 >  1:33776‑M2/1‑132 (MQ=255)                  ttaattGTCGGAAAAGCCTTTTCCATTGCCGGGTATTCTTCTGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGTAGTGGCTGAAGGCCACGGAACCGTGACACCAACAGGTAATTCAGAGCCTTCTCCCAAACCAAc                                                                                                                                 <  2:331816/149‑1 (MQ=255)                           ggATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGCCACCGCCAGGTAATGCAGAGcctgctcccaaaccaaggtcttcaag                                                                                                                        <  1:124319‑M2/149‑27 (MQ=255)                            gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                       <  1:127483‑M2/149‑28 (MQ=255)                            gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                                                                                       >  1:540865/1‑149 (MQ=255)                             aTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTATTTGCCGAGTCGCTAAAGGCCACGGAACCGGGACACCAACAAGTAATGCAGGGCCTTCTCCCCAACCCACATTTACGaaaa                                                                                                                      >  2:609702/1‑149 (MQ=255)                                 gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctccccaaccaaggttttcgagagtcgt                                                                                                                  <  1:296994‑M2/149‑33 (MQ=255)                                 gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctaccaaaccaagagtttcgagagccgt                                                                                                                  <  2:483022‑M2/149‑33 (MQ=255)                                      tGCCAGTGCCGGGGAGTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgca                                                                                                             >  2:649877‑M2/1‑112 (MQ=255)                                        ccAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATaa                                                                                                           >  1:631891/1‑149 (MQ=255)                                        ccAGGTGCGGGTATTGTGCATCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGATGGCGTAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGACTTCTCCCAATCCAACGTTTATGAAAATGAAGAAATaa                                                                                                           <  2:308854/149‑1 (MQ=255)                                         cAGTTCGGGGTATTCTTCATCAGAAAAAGCGGCACCATTGCATTTCAGACAGCCTGTTGGCGGAGGTGCTTAAGGCCACGGAACCGGTACATCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacag                                                                                                          >  2:178195‑M2/1‑109 (MQ=255)                                                     ttcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                                                                                              >  1:544465‑M2/1‑97 (MQ=255)                                                     ttcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                                                                                              >  1:287568‑M2/1‑97 (MQ=255)                                                        ttcAGCAGAAAAAGCGGCACCATTGCATTTCCGCCCGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga                                                                                           <  1:508392‑M2/149‑56 (MQ=255)                                                        ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTACTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcctca                                                                                           >  1:245959/1‑149 (MQ=255)                                                             cagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg                                                                                      <  1:283755‑M2/149‑61 (MQ=255)                                                              agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc                                                                                     >  2:283755‑M2/1‑88 (MQ=255)                                                               gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                    <  1:28120/149‑1 (MQ=255)                                                                aaaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGCAGGGGCTGAGGACCACGAAACCGCGACACCAACCGGTAATGCCGAAccttctaccaaaccaaggttttcgagagccgtttgcacagttgcatccgatttgatatcccca                                                                                   <  2:327043‑M2/149‑64 (MQ=255)                                                                          aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgggttct                                                                         <  1:305743‑M2/149‑74 (MQ=255)                                                                          aCCATTGCATTTCAGCCAGCCTCTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatctcgccaaacgggttct                                                                         <  2:493430‑M2/149‑74 (MQ=255)                                                                            cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGGGccttctccccaaccaaggttttcgagacccgtttgcacgctgccatccgatttgatatcgcaaaacggtatcttc                                                                       <  2:412992‑M2/149‑76 (MQ=255)                                                                            cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttg                                                                       <  1:414054‑M2/149‑76 (MQ=255)                                                                            cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttg                                                                       <  1:414759‑M2/149‑76 (MQ=255)                                                                                              cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcg                                                     <  1:523456‑M2/149‑94 (MQ=255)                                                                                              cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcg                                                     >  1:203704‑M2/1‑56 (MQ=255)                                                                                                  tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaa                                                 >  1:575631‑M2/1‑52 (MQ=255)                                                                                                            gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgctttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaag                                       >  1:167161‑M2/1‑42 (MQ=255)                                                                                                                                   atcAACAGGTAATGCAGAGCCTTCTTCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                <  1:106640/147‑1 (MQ=255)                                                                                                                                   acCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttct                <  1:116364‑M2/149‑131 (MQ=255)                                                                                                                                           gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtggtgccgccccttttccaggctg        <  1:565797‑M2/149‑139 (MQ=255)                                                                                                                                                gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAATACAGATCTACAACGGAATg   >  1:345507/1‑149 (MQ=255)                                                                                                                                                 cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcagcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccg  >  2:546261‑M2/1‑5 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGC  >  NC_000913_3_hsa_pgi/1209454‑1209746

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.