New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	1207790 =	19 (0.570)	11 (0.370)	11/248	1.8	46.0%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_hsa_tpiA	1209619 =	9 (0.300)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGG  >  NC_000913_3_hsa_tpiA/1207644‑1207937                                                                                                                                                   |                                                                                                                                                    catgctgcacagttctgcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatacccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTg                                                                                                                                                   >  2:438147‑M1/147‑149 (MQ=255)  atcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTgg                                                                                                                                                  >  2:184899‑M1/146‑149 (MQ=255)  atcctgcacaggtcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTgg                                                                                                                                                  >  2:60440‑M1/146‑149 (MQ=255)    cctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGtt                                                                                                                                                >  2:348154‑M1/144‑149 (MQ=255)      tacacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                                                                                                              >  2:390860‑M1/142‑149 (MQ=255)        gtggAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGG                                                                                                                                            <  2:384785‑M1/149‑1 (MQ=255)        aactgttcatcgtattacgttgtagatctgtgttctggttatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                                                                                                            >  1:497891‑M1/140‑149 (MQ=255)         tggAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGA                                                                                                                                           <  1:233691‑M1/149‑1 (MQ=255)         acagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtcaacgaagtgaggatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGga                                                                                                                                           >  1:315765‑M1/139‑149 (MQ=255)          ggAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTCTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACATTGGTTCGGGAG                                                                                                                                          >  2:221659‑M1/1‑149 (MQ=255)                aTCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAg                                                                                                                                    <  1:400395/149‑1 (MQ=255)                    cattacgttgtagatctgggttctggtcatttgttgatacgcgtacatagcaaataagcatgttaaatccccctggtaaaagtaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                                                                                                >  2:52351‑M1/128‑149 (MQ=255)                     gTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaaa                                                                                                                               <  2:557339/149‑1 (MQ=255)                        caaggctggagAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcatta                                                                                                                            <  2:368547‑M1/144‑11 (MQ=255)                          aGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGGACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacc                                                                                                                          <  2:241290‑M1/149‑13 (MQ=255)                                  tctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggtatgatgccatttgcttgttatttcttcattttcataaacgTGGGTTGGGGAGAAGGGGCAAGACTCAATGCAGCaa                                                                                                                  <  1:20191‑M1/36‑1 (MQ=255)                                             gtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                                                                       >  1:277716‑M1/103‑149 (MQ=255)                                             gtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGATTTGGGAAAAGGACCAAAACTCAATGCAGCAACGGCTAcataa                                                                                                       <  2:191562‑M1/47‑3 (MQ=255)                                               gACCAGTTACTTACCGCACTTAGCGCGCCGGAGTTAAGCCGGAAGAATCCGTCTGGCGATATCAAATCGGATGGCAGTGTGCAAACGGTTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                                                                                                     <  2:490959‑M1/149‑34 (MQ=255)                                                    gTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcctt                                                                                                >  1:63815‑M1/1‑111 (MQ=255)                                                     ttacttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAc                                                                                               >  1:161260/1‑149 (MQ=255)                                                     ttacctacCCCACTTCGCGCGCTGCTGTTAAGCCCCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGACCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAc                                                                                               <  2:275664/149‑1 (MQ=255)                                                             cGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCa                                                                                       >  1:44894/1‑149 (MQ=255)                                                              gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCAt                                                                                      <  2:427042/149‑1 (MQ=255)                                                               cACTTCGCGGGCTGCTGTTAAGCCGGAAGGATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                                                                     <  2:29937/149‑1 (MQ=255)                                                                       gcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAcc                                                                             >  1:357670/1‑149 (MQ=255)                                                                            gctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAACGGTTTAATAGCTATACCGGTTa                                                                        >  2:491948/1‑149 (MQ=255)                                                                                ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                                                                                                         <  1:463854‑M1/80‑1 (MQ=255)                                                                                ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                                                                                                         <  1:463954‑M1/80‑1 (MQ=255)                                                                                ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                                                                                                         >  2:463854‑M1/1‑80 (MQ=255)                                                                                ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                                                                                                         >  2:463954‑M1/1‑80 (MQ=255)                                                                                    aaatcccccgtgtaaaaggaggaatgatgccatttgcttggtatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                                                                >  2:20191‑M1/64‑149 (MQ=255)                                                                                     aatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                               >  2:426600‑M1/63‑149 (MQ=255)                                                                                         cccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAt                                                           >  1:385736‑M1/59‑149 (MQ=255)                                                                                           cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaaagTTGGTTTAGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGGGGTATTg                                                         <  1:198232‑M1/93‑1 (MQ=255)                                                                                            aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgc                                                        <  2:143952‑M1/149‑79 (MQ=255)                                                                                            aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgc                                                        <  2:143562‑M1/149‑79 (MQ=255)                                                                                                 aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                                   >  2:169783‑M1/51‑149 (MQ=255)                                                                                                         gtatCAAATCGGATGGCACTGTGCAAACGGCTCGCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATc                                           <  1:233373/148‑1 (MQ=255)                                                                                                          aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCa                                          <  1:447287‑M1/108‑1 (MQ=255)                                                                                                            tCAAATCGGATGGCAGTGGGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaa                                        <  2:524254‑M1/149‑95 (MQ=255)                                                                                                            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaa                                        >  1:454345‑M1/1‑55 (MQ=255)                                                                                                            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCCtt                                        >  2:234878/1‑149 (MQ=255)                                                                                                            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCAtt                                        >  1:50109/1‑149 (MQ=255)                                                                                                              aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaaga                                      >  2:179855‑M1/1‑53 (MQ=255)                                                                                                              aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                         >  2:60192/1‑98 (MQ=255)                                                                                                              aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                         <  1:60192/98‑1 (MQ=255)                                                                                                               aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgacatgcaatggtgccgctttttctgctgaagaa                                     >  2:214165‑M1/1‑52 (MQ=255)                                                                                                                gccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGt                                    >  1:460861‑M1/36‑149 (MQ=255)                                                                                                                 tCGGATGGCACTGTGCAAACGGCTCACGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTATGACGCACCGGTTTCATAGCTGTACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTg                                   <  1:527281/149‑1 (MQ=255)                                                                                                                        gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTgg                            <  2:44894/149‑1 (MQ=255)                                                                                                                                     ggCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                                                                     >  1:229683/1‑79 (MQ=255)                                                                                                                                     ggCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                                                                     <  2:229683/79‑1 (MQ=255)                                                                                                                                          tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccaccccgccaacaggctagctgaaatgccatggtcccgctttttctgctgaagaatacccggaactggcaaaggcttatccg          >  2:147836‑M1/1‑25 (MQ=255)                                                                                                                                            gAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgac        >  1:231017‑M1/1‑23 (MQ=255)                                                                                                                                            gAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgac        >  1:479058‑M1/1‑23 (MQ=255)                                                                                                                                                  cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg  >  1:142540/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGG  >  NC_000913_3_hsa_tpiA/1207644‑1207937

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