New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 12078057 (0.210)20 (0.660) 19/248 0.7 67.6% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA = 1209602 13 (0.430)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACA  >  NC_000913_3_hsa_tpiA/1207658‑1207953
                                                                                                                                                   |                                                                                                                                                    
cATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGa                                                                                                                                                     <  1:186159/149‑1 (MQ=255)
      attacggtgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaaa                                                                                                                                               >  1:393224‑M1/127‑149 (MQ=255)
          cAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                                                                                                           <  1:269530‑M1/149‑12 (MQ=255)
               tGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                                                                                                                      >  1:173358‑M1/1‑133 (MQ=255)
                 gagaAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttgg                                                                                                                                    <  2:592561‑M1/149‑19 (MQ=255)
                      ggggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcc                                                                                                                               >  1:552062‑M1/1‑126 (MQ=255)
                       gggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCACGCCTGCTTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                                                                                                                              <  1:441027‑M1/149‑25 (MQ=255)
                            tggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTTGGACAAGGAGCAACACTCAATGCAGCAACGGCTCCAt                                                                                                                         <  2:364518‑M1/45‑1 (MQ=255)
                               cGACCAGTTACTTACCGCACTTCGCGCGCAGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                                                                                                                      <  2:338786‑M1/149‑33 (MQ=255)
                                 atttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                                                                                                                    <  1:322464‑M1/50‑1 (MQ=255)
                                       tacttacCGCACGTCGCGTGTGGCTGGTAGGCCGTAAGAACCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAACCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttca                                                                                                              <  2:176544‑M1/149‑41 (MQ=255)
                                          ttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttgatgtcccggttccgtggccttcagcc                                                                                                           >  2:345556‑M1/1‑106 (MQ=255)
                                          ttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                                                                                                           >  2:484600/1‑149 (MQ=255)
                                           acgcgtatattggcaataatcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGtt                                                                                                          >  2:239506‑M1/90‑149 (MQ=255)
                                            acCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttccgcccc                                                                                                         >  1:218349‑M1/1‑104 (MQ=255)
                                                tacatagccaataagcatgttaaatccccctggtaaaagcaggcatgatgccatttgctagttatctcttcattttcataaaagTTGGTTTGGGAGAAGGAGCAAAACACAATGCAGCAGCGGCTAACTTAGGCCGCACCGGTTTCATa                                                                                                     <  2:402172‑M1/65‑1 (MQ=255)
                                                 cattcgGGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgc                                                                                                    <  2:173358‑M1/147‑51 (MQ=255)
                                                       cgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                                                                              <  2:579445/149‑1 (MQ=255)
                                                          ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg                                                                                           <  1:283271‑M1/75‑1 (MQ=255)
                                                              ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg                                                                                            >  2:212798/1‑144 (MQ=255)
                                                              ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg                                                                                            <  1:212798/144‑1 (MQ=255)
                                                              cggttAAGACGCAAGAATACGTTTGGCGATATCAAATCGTATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggct                                                                                       <  2:224157‑M1/147‑64 (MQ=255)
                                                                  tAAGCCGCAAGGATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaa                                                                                   <  1:508243‑M1/149‑68 (MQ=255)
                                                                  gttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                                                      >  1:125533‑M1/67‑98 (MQ=255)
                                                                  gttaaatccccctggtaaaaccaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                                                      <  2:125533‑M1/32‑1 (MQ=255)
                                                                     aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGCTATGCttggg                                                                                <  1:160148‑M1/86‑2 (MQ=255)
                                                                      ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcccctccgccaacaggctggctgaaatgca                                                                               >  2:210012‑M1/1‑78 (MQ=255)
                                                                          aaGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAt                                                                           <  1:484600/149‑1 (MQ=255)
                                                                           ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                                          <  1:229348‑M1/92‑1 (MQ=255)
                                                                                ccGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgc                                                                     <  1:402361‑M1/149‑82 (MQ=255)
                                                                                    ttGGCGATATCAAAGCGGAGGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATCCAGCAACGGCTACATTCGGACGCACCGGTTTCATAGCTATACCGGTTATGATTTGTGCTATTGAGCAATCa                                                                 <  2:144681/149‑1 (MQ=38)
                                                                                       caggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGAACCGGTTTCATAGCTCCACCGGTTATGATTGTAGGTATTGAGCAATGAGTa                                                              >  2:504906‑M1/46‑149 (MQ=255)
                                                                                         ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAAt                                                            >  1:205152‑M1/44‑149 (MQ=255)
                                                                                          gaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATc                                                           >  1:32718‑M1/43‑149 (MQ=255)
                                                                                               aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgcagaaga                                                      >  1:391824‑M1/1‑53 (MQ=255)
                                                                                                aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgatggtgtcccggttccgtggccttcagccactacgccaacaggctggctgaaatgcaatagtgccgccttttattctcaagaa                                                     >  2:137615‑M1/1‑52 (MQ=255)
                                                                                                   catttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgg                                                  <  1:53736‑M1/116‑1 (MQ=255)
                                                                                                     gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaaggcaatggtgccgctttttctgctgaagaataccc                                                <  1:591321‑M1/149‑103 (MQ=255)
                                                                                                     gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggcaggctgaaatgcaatggtgccgctttttctgctgaagaataccc                                                <  1:590905‑M1/149‑103 (MQ=255)
                                                                                                          ttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGGGCAATCAGTAATCATTCAGTGGGGGTGGa                                           <  1:306343‑M1/123‑1 (MQ=255)
                                                                                                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaact                                          >  1:519231‑M1/1‑41 (MQ=255)
                                                                                                             ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGCCGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTTAGCAATCAGTAATCATTCAGTGGGGGTAGAATg                                        <  1:192297‑M1/126‑1 (MQ=255)
                                                                                                               atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCCGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGcc                                      >  1:350952‑M1/22‑149 (MQ=255)
                                                                                                               atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGcc                                      >  1:376175‑M1/22‑149 (MQ=255)
                                                                                                                  aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggctccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaag                                   >  2:94153‑M1/1‑34 (MQ=255)
                                                                                                                     ttattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                                >  2:82990‑M1/16‑149 (MQ=255)
                                                                                                                           tgatgaacgTTGGTTTGGGAGAAGGAGTAAAACTCAATGGAGCAAGGGCTAAATTAGGACGCACCGGTGTCATAGCTTGACCGGTTATGATTGGTGGGATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCt                          >  2:471207‑M1/10‑149 (MQ=255)
                                                                                                                                    TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggg                 >  1:485534‑M1/1‑149 (MQ=255)
                                                                                                                                      ggTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataa                  >  2:373440‑M1/1‑14 (MQ=255)
                                                                                                                                      ggTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataa                  <  1:373440‑M1/146‑133 (MQ=255)
                                                                                                                                        tttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgcttgttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgc             <  1:398778‑M1/149‑138 (MQ=255)
                                                                                                                                              AGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGaa       <  2:184267‑M1/149‑1 (MQ=255)
                                                                                                                                               gaAGGctctgcattacctgttggtctcccggttccgtggccttcagccactccgccaacatgctggctgaaatgcaatggtgccgctttttctgctgcagaatacccggacctggcaacggcttatccgagaaataaattacctgattt      >  2:278021‑M1/1‑5 (MQ=255)
                                                                                                                                                   gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaagggcttatccgacaaataaattgactgatttacgt  >  1:150503‑M1/1‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACA  >  NC_000913_3_hsa_tpiA/1207658‑1207953

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.