New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 1207805	19 (0.570)	13 (0.440)	11/248	1.8	50.2%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_hsa_tpiA	= 1209602	9 (0.300)					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                                                                                                                                                     |                                                                                                                                            cacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                                                                                             >  2:63815/1‑149 (MQ=255)  aCGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                                                                                            <  1:365567‑M2/149‑2 (MQ=255)  aCGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCAAGGAAACGGGACACACACAGGGAATGCAGAGc                                                                                                                                            <  2:574191‑M2/149‑2 (MQ=255)   cGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcc                                                                                                                                           >  1:61968‑M2/1‑147 (MQ=255)      aaaTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttc                                                                                                                                        >  2:49220‑M2/1‑144 (MQ=255)         tccgccaatctATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCGTCAGCAGAAAAAGCGGCCCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCACAGccttctcc                                                                                                                                     >  1:344284‑M2/6‑141 (MQ=255)            ggCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGCAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGGAccttctcccta                                                                                                                                  <  1:89475‑M2/149‑12 (MQ=255)                    atttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGt                                                                                                                          >  1:199986/1‑149 (MQ=255)                      ttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                                                                                        <  1:163412/149‑1 (MQ=255)                        gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                                                                                      >  1:231017‑M2/1‑126 (MQ=255)                        gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                                                                                      >  1:479058‑M2/1‑126 (MQ=255)                          cGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGGACCATGGCATTTCAGCTAGCCTGTTGGCGGGGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                                                                    >  2:147836‑M2/1‑124 (MQ=255)                                         cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGACAGACGCAAACTGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAc                                                                                                     <  1:573100/149‑1 (MQ=255)                                                     ttcttcAGCAGAAAAAGCGGCACCATTGCATGTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                                                                                         >  2:214165‑M2/1‑97 (MQ=255)                                                      tcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                                                                                        >  2:179855‑M2/1‑96 (MQ=255)                                                        ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcccactgccatccgatttga                                                                                      <  2:524254‑M2/149‑56 (MQ=255)                                                        ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga                                                                                      >  1:454345‑M2/1‑94 (MQ=255)                                                        ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCACCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAACAACACGCAAATgcactca                                                                                      >  2:135497/1‑143 (MQ=255)                                                             cagAAAAAGCGGCTCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                 <  1:568626/149‑1 (MQ=255)                                                                       ggCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGACATCATTCCTGCTTTTAcct                                                                       >  2:43569/1‑148 (MQ=255)                                                                        gCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggatt                                                                      <  2:143952‑M2/149‑72 (MQ=255)                                                                        gCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggatt                                                                      <  2:143562‑M2/149‑72 (MQ=255)                                                                                         ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTATCATGCtt                                                     <  2:365567/149‑1 (MQ=255)                                                                                                              tGAAGGCCACGGAACCGGGACACCCACAGGTAATGCAGAAGATTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATca                                >  1:84686/1‑149 (MQ=255)                                                                                                                aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaac                              >  1:63815‑M2/1‑38 (MQ=255)                                                                                                                     cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagaaccgtttgcacactgccatccgatttgatatcgccagacggattcttccggcttaactccggcgcgctaagtgcggtaagtaactggtc                         <  2:490959‑M2/149‑117 (MQ=255)                                                                                                                      aCGAAAACGGGACCCCACACGGTAATGCATAAGCTTCTCACAAACCACCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATTAATCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa                        <  2:190787/149‑1 (MQ=255)                                                                                                                                          ggTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtccggtaagtaactggtcgtgccgccccttctccaggct    <  2:241290‑M2/149‑138 (MQ=255)                                                                                                                                            tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccagccttg  <  2:368547‑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.