New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 12078059 (0.260)22 (0.720) 20/244 NT 71.0% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA = 1209602 10 (0.330)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGG  >  NC_000913_3_hsa_tpiA/1207658‑1207946
                                                                                                                                                   |                                                                                                                                             
ccgcgcatgatggtgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                                                                                                              >  1:118371‑M1/133‑149 (MQ=255)
cagcgagttacgttgtagatctgggttctggtcatttgttcatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                                                                                                              >  2:145281‑M1/133‑149 (MQ=255)
 agcgcattacgttgtagatctgtgttctggtcctttgttgatacgcgtacatagccaataagcatgttaaatcccccttgtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAg                                                                                                                                             >  1:75611‑M1/132‑149 (MQ=255)
       ttacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAc                                                                                                                                       <  1:199307‑M1/24‑1 (MQ=255)
       tGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTGGGGAGAAGGctctgcat                                                                                                                                       >  1:439591‑M1/1‑141 (MQ=255)
               tatatctgtgttctggtcatttattgatacgcgtatatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                                               >  1:577208‑M1/118‑149 (MQ=255)
                 gagaAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCCGTTAAGCCGCAAGAATCCGACTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCGCCCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctttttg                                                                                                                             >  2:465272‑M1/1‑131 (MQ=255)
                    aaggtcaggcacgaccaGGTACTTACCGCACTTGTCACGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCTAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgt                                                                                                                          <  2:398828‑M1/142‑22 (MQ=255)
                    aaGGGGCGGCACGACCAGTTACTTCCCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgt                                                                                                                          <  2:40726‑M1/149‑22 (MQ=255)
                    aaGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgt                                                                                                                          >  1:620176‑M1/1‑128 (MQ=255)
                     tgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAGACTCAATGCAGAAACg                                                                                                                         <  1:334698‑M1/38‑1 (MQ=255)
                         gcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTa                                                                                                                     <  1:83496/149‑1 (MQ=255)
                               tcatttgttgatacgcgtacatagccaataagaatgctaacaacccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGGGAAGCAGCAAAACTCAATGCAGCAACGGCTACATTAg                                                                                                               <  1:156762‑M1/48‑1 (MQ=255)
                                       tacttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttca                                                                                                       <  1:163936‑M1/149‑41 (MQ=255)
                                         atacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCgg                                                                                                     <  1:110370‑M1/58‑1 (MQ=255)
                                              cGCACTTCGCGCGCTGCTGTTAAGCCGGAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactc                                                                                                <  2:196896‑M1/149‑48 (MQ=255)
                                              cGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                                                                                                                  <  1:329105/131‑1 (MQ=255)
                                              cGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                                                                                                                  >  2:329105/1‑131 (MQ=255)
                                                tacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaaagTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                                                                              <  2:594762‑M1/65‑1 (MQ=255)
                                                cACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccg                                                                                              >  2:448936‑M1/1‑100 (MQ=255)
                                                   ttCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcccctccgcca                                                                                           >  1:45536‑M1/1‑97 (MQ=255)
                                                   ttCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcca                                                                                           >  1:231705‑M1/1‑97 (MQ=255)
                                                       cgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacag                                                                                       >  2:455772‑M1/1‑93 (MQ=255)
                                                          gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggct                                                                                    >  1:369051‑M1/1‑90 (MQ=255)
                                                            tgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCAGTTTCATAGCTATACCGGtt                                                                                  >  2:52909/1‑149 (MQ=255)
                                                                  tAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt                                                                            <  2:250393/149‑1 (MQ=255)
                                                                         ccccctggtaaaagcaggaatgatgccagttgcttgttatttcttcagtttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAGTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTa                                                                     >  2:29973‑M1/60‑149 (MQ=255)
                                                                         aaaaaatccgTGGGGCGATATCAAATCGGATGGCACTGGGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatg                                                                     <  1:613836‑M1/145‑75 (MQ=255)
                                                                          cccctggtaaaagcaggaatgatgccctgtgcttgttatttcttcagtttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAAGGCTACATTAGGACGCACCGGGTTCATAGCTGTACCGGTTATGATTGGTGGTAt                                                                    <  2:509668‑M1/91‑1 (MQ=255)
                                                                             ctggtaaaagcaggaatgatgccctttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                                                                 <  1:557086‑M1/94‑1 (MQ=255)
                                                                             aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAACCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgc                                                                 <  1:182513‑M1/149‑79 (MQ=255)
                                                                                    aagcaggaatgctgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCa                                                          >  2:26597‑M1/49‑149 (MQ=255)
                                                                                      ggCGATATCAAATCTGATGGCACTGTGCAAACGGCTTTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttc                                                        <  2:159535‑M1/149‑88 (MQ=255)
                                                                                          atatCAAATTGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgct                                                    <  2:589235‑M1/149‑92 (MQ=255)
                                                                                               atgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACAGGTTTCATAGCTATAACGGGTATGATTGGTGGTATTGATCAATCAGTAATcatttg                                               <  2:299170‑M1/112‑3 (MQ=255)
                                                                                                 aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaat                                             >  1:263800‑M1/1‑51 (MQ=255)
                                                                                                     tttgcttgttatttcttcattttcataaacgTTTGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                         >  2:339517‑M1/32‑149 (MQ=255)
                                                                                                     gATGGCAGTGTGCATACGGCTCTCGAAAGCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatgttgccgctttttctgctgaagaataccc                                         <  2:650392‑M1/149‑103 (MQ=255)
                                                                                                     gATGGCACTGTACAAACGGCTCTCGAAAACCTTGGTTTGGGAGACGGctctgaattaactgttggtgtcccggttccgtggccttcagccactccgccaacaggctgggtgaaatgcattggtgccgctttttttgatgaaaaataccc                                         >  2:251322‑M1/1‑47 (MQ=255)
                                                                                                       tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccgg                                       >  1:364619‑M1/1‑45 (MQ=255)
                                                                                                          ttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAAATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGATATTGAGCAATCAGTAATCATTCAGGTGCGGTGGa                                    <  2:289927‑M1/123‑1 (MQ=255)
                                                                                                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAGACTCAATCCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATCCCGGTTATTAATGATGGTATTGAGCAATCAGTAATCATTCAGGAGGGGTGGaa                                   >  2:206855/1‑149 (MQ=255)
                                                                                                                 cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGATATGATTGGTGGTATGGAGCAATCAGTAATCATTGAGTGACGGTGGAATGCCGc                             >  2:374176/1‑149 (MQ=255)
                                                                                                                     tcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                         <  2:147565‑M1/134‑1 (MQ=255)
                                                                                                                     tcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATAAGAAATCAGTCAGTGGGGGTGGAAGGCCGCaaaa                         <  1:172668‑M1/134‑1 (MQ=255)
                                                                                                                      ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttaagccactcgcccaacagcctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactgggaacggctt                        >  2:237057‑M1/1‑30 (MQ=255)
                                                                                                                        ttttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATACTGCAGAAGCa                      <  2:477713‑M1/137‑1 (MQ=255)
                                                                                                                         tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatc                     <  1:448936‑M1/149‑123 (MQ=255)
                                                                                                                                   cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTCTTGAGCAATCAGTAATCATTCAGTGGGGGTGGCATGCCGCAAAAGCATCTGCCTCTgg           >  2:625865/1‑149 (MQ=255)
                                                                                                                                      tgtTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaatt        <  1:48542‑M1/148‑136 (MQ=255)
                                                                                                                                      tgcttGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTCGAGTAGGACGCGCCGGTGTCCAAGAAATACCCGTCACGAGTGTGGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggg        <  2:40/146‑1 (MQ=255)
                                                                                                                                      GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggg        >  2:226381‑M1/1‑149 (MQ=255)
                                                                                                                                       gTTTGGGAGAAGGctctgcatcacctgggggtgtcccgattccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattg       >  2:334228‑M1/1‑13 (MQ=255)
                                                                                                                                       GTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgggggg       >  1:562010‑M1/1‑149 (MQ=255)
                                                                                                                                          TGGGGGAAGGACCAAAACCCAATGCAGCAACGGCTACTTTAGGACGCACTGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATTATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGAt    >  2:26338‑M1/1‑149 (MQ=255)
                                                                                                                                           GGGGGAAGGAGCAAAACTAAATCCAGCAACGGCTAAATTAGGACGCACCGGTTTCATAGCTATACCGGTGATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATg   >  2:373700‑M1/1‑149 (MQ=255)
                                                                                                                                            GGAGAAGGAGCAAAACTCAGTTCAGCAACGGCTAGAATAGGAGTCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCCATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGTAGCTGCCTCTGGGGGGGATgg  >  2:81037‑M1/1‑149 (MQ=255)
                                                                                                                                            GGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCATCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCACAAGCATCTGCCTCTGGGGGGGATgg  <  1:596359‑M1/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGG  >  NC_000913_3_hsa_tpiA/1207658‑1207946

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