New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 1207805	10 (0.210)	36 (0.840)	32/250	0.4	82.9%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_hsa_tpiA	= 1209602	6 (0.140)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGT  >  NC_000913_3_hsa_tpiA/1209454‑1209741                                                                                                                                                     |                                                                                                                                            cacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGATATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCCGAGGGGCTTAAGGCCACGGCACCGAGACACCAACAGGTAATGCagag                                                                                                                                             >  2:40761/1‑149 (MQ=255)  aCGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                                                                                            <  1:706068‑M2/149‑2 (MQ=255)   cGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGGAGAGAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGGGTGGCTGAAAGCCACGGAACCGGGACACCAACAGGTATTGCCGAGcc                                                                                                                                           <  2:396891‑M2/149‑3 (MQ=255)   cGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcc                                                                                                                                           <  2:734162‑M2/149‑3 (MQ=255)     tAAATCAGGCAATTTATGTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcctt                                                                                                                                         >  2:584152‑M2/1‑145 (MQ=255)          cAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctccc                                                                                                                                    >  1:734162‑M2/1‑140 (MQ=255)            ggCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaa                                                                                                                                  <  1:106765‑M2/149‑12 (MQ=255)             ggtatttaattgtaGGATAGGCGTTTGCAAGTGCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATCTCATCCAGCCTGTTGGCGGAGTGCCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                                                                                                 >  2:376314‑M2/4‑137 (MQ=255)             gCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                                                                                                 >  2:106765‑M2/1‑137 (MQ=255)                 tttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaa                                                                                                                             >  2:447682‑M2/1‑133 (MQ=255)                  ttatttGTCGGATAAGCATTTGCCAGTTCCGTGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCCGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                                                                                            >  2:14995‑M2/1‑132 (MQ=255)                        gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                                                                                      >  1:809674‑M2/1‑126 (MQ=255)                         tCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcg                                                                                                                     <  1:336673‑M2/149‑25 (MQ=255)                            gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                  >  1:364331‑M2/1‑122 (MQ=255)                            gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                  >  1:353713‑M2/1‑122 (MQ=255)                             aTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCGCAATTGCTTTTCACGAAGCCTGTTGGAAGTGTGCCCGCCTCCTTCAGACCCGCCACATCAAGAGGCAAAAGAGAACCTTCTCCCAAACCAAGGttttcgagag                                                                                                                 >  2:201689/1‑142 (MQ=255)                                 gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                                                                             <  1:363039‑M2/149‑33 (MQ=255)                                   cTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                                                                                           >  2:446744‑M2/1‑115 (MQ=255)                                   cTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                                                                                           >  1:579645‑M2/1‑115 (MQ=255)                                      tGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt                                                                                                        <  1:50050/149‑1 (MQ=255)                                       gCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                                                                                           >  2:323259/1‑129 (MQ=255)                                       gCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                                                                                           <  1:323259/129‑1 (MQ=255)                                               cGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccat                                                                                               <  2:57412‑M2/149‑47 (MQ=255)                                               cGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccat                                                                                               >  2:58476‑M2/1‑103 (MQ=255)                                                gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACAGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatc                                                                                              <  2:666765‑M2/149‑48 (MQ=255)                                                  gTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGCGTGAACGCCACGGAACCGGGACACCAACAGCTAATGCATAGccttctcccaaacccacactttcgagagccgcttgcccagtgccatccg                                                                                            <  2:412590‑M2/149‑50 (MQ=255)                                                        ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgaacagtgccatccgatttga                                                                                      <  2:330723‑M2/149‑56 (MQ=255)                                                               gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcc                                                                               >  1:780583‑M2/1‑87 (MQ=255)                                                                            cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                  >  1:227765/1‑149 (MQ=255)                                                                              ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgctatccccaaacgggttcttgcg                                                                <  2:63206‑M2/149‑78 (MQ=255)                                                                                 cttttcaGCCAGCTTGTTGGGGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggct                                                             >  1:324324‑M2/3‑69 (MQ=255)                                                                                       agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACAACAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgctctcgaatcgcgggtgcttgcggcttagcag                                                       <  2:389961‑M2/149‑87 (MQ=255)                                                                                        gccagccTGGTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGccttctcccaaaccaaggttgtcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagc                                                      >  1:63206‑M2/1‑62 (MQ=255)                                                                                         ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg                                                                                 <  2:277891‑M2/121‑61 (MQ=255)                                                                                         ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg                                                                                 >  1:277891‑M2/1‑61 (MQ=255)                                                                                         ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                                                                             <  1:126784‑M2/93‑33 (MQ=255)                                                                                         ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                                                                             >  2:126784‑M2/1‑61 (MQ=255)                                                                                          cagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcag                                                    >  1:348663‑M2/1‑60 (MQ=255)                                                                                              cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcg                                                <  2:577432‑M2/149‑94 (MQ=255)                                                                                                    gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                          <  2:87207/149‑1 (MQ=255)                                                                                                      ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgc                                        >  2:706068‑M2/1‑48 (MQ=255)                                                                                                       gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcg                                       >  1:40761‑M2/1‑47 (MQ=255)                                                                                                              tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagta                                <  2:542701‑M2/149‑110 (MQ=255)                                                                                                                     cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtc                         >  1:49824‑M2/1‑33 (MQ=255)                                                                                                                     cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggctgaacagcagcgcgcgaagtgcggtaagtaactggtc                         <  1:63452‑M2/149‑117 (MQ=255)                                                                                                                       cGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgt                       >  2:363039‑M2/1‑31 (MQ=255)                                                                                                                       cGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgt                       >  1:57412‑M2/1‑31 (MQ=255)                                                                                                                        ggAACCGGGACACCAACATGTAATGCAGAGcctgcgcccaaaccaatgttttcgagagccgtttgcacagtgccatccgatttgatatctccaaacggattcttgcggcttaacagcagctcgcgaagtgcggtaagtaactggtcgtg                      >  1:389961‑M2/1‑30 (MQ=255)                                                                                                                              gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc                <  2:77700‑M2/149‑126 (MQ=255)                                                                                                                               ggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttgacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccc               <  1:641740‑M2/149‑127 (MQ=255)                                                                                                                                            tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctgg  >  1:77700‑M2/1‑10 (MQ=255)                                                                                                                                            tAATGCAGAGccttctcccaaaccaaggttttcgagagcagtttgcacagtgccatacgatttgatatctccaaacggattcttgcggattaacagcagagcgcgaagtgaggtaagtaactggtcgtgccgacccttctcaaggatgg  <  2:372419‑M2/149‑140 (MQ=255)                                                                                                                                                     |                                                                                                                                            CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGT  >  NC_000913_3_hsa_tpiA/1209454‑1209741

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.