New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA 1207790 =4 (0.140)23 (0.920) 19/248 0.4 92.8% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA 1209619 = 0 (0.000)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCC  >  NC_000913_3_hsa_tpiA/1207644‑1207932
                                                                                                                                                  |                                                                                                                                              
cAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTG                                                                                                                                              <  1:40898‑M1/149‑1 (MQ=255)
     tgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                                                                                                         >  1:24005‑M1/142‑149 (MQ=255)
       caaagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatcccccttgtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                                                                                                       >  1:295911‑M1/140‑149 (MQ=255)
                 ccGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                                                                                             <  2:87458‑M1/149‑5 (MQ=255)
                        cgttgtagatctgtgttctggtcattggttgatacgcgtacatagcaaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                                                                                                      >  1:317913‑M1/123‑149 (MQ=255)
                        cAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                                                                                      >  1:456392‑M1/1‑138 (MQ=255)
                                    gtgttctggacgtttgttgtgacgcgtacatggccaagaagaatgttaaattccccgggtataagcaggaatgatgccatttgcgtgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                          >  2:488055‑M1/111‑149 (MQ=255)
                                     gggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                                                                                                         <  2:284775‑M1/149‑25 (MQ=255)
                                       gcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccgc                                                                                                       >  1:470118‑M1/1‑123 (MQ=255)
                                            gtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCCATGCAGCAACGGCTACATTa                                                                                                  <  1:280407‑M1/47‑1 (MQ=255)
                                             cggccagtTACTGACCGCACTTCGCGCGCGGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGATAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                                                                                                 <  2:470118‑M1/146‑33 (MQ=255)
                                                ccAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                                                                                              <  1:177572‑M1/149‑36 (MQ=255)
                                                       cgtaccgCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagc                                                                                       <  1:78571‑M1/147‑43 (MQ=255)
                                                          acCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGATTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGTTTTCGGAGAAGGctcggcattacctgttggtgtcccggttccgtggccttcagccac                                                                                    >  1:195142‑M1/1‑104 (MQ=255)
                                                              cACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccg                                                                                >  1:73071‑M1/1‑100 (MQ=255)
                                                                      caataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAcc                                                                        <  1:8319‑M1/73‑1 (MQ=255)
                                                                         ctagcatgttaaatccccctggtacacgcaggaattatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt                                                                     >  2:230831‑M1/74‑149 (MQ=255)
                                                                            gcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCCAGATTAGGCCACACGGGCTTACTagcgatacaggctag                                                                  <  1:230831‑M1/79‑8 (MQ=255)
                                                                              atgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGa                                                                <  2:465807‑M1/81‑1 (MQ=255)
                                                                               ttAAGCCGCAACAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAGGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaa                                                               >  1:307572‑M1/1‑83 (MQ=255)
                                                                                gttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcgtcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt                                                              >  2:312810‑M1/67‑149 (MQ=255)
                                                                                 ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAGTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATCGCTATACCGGTTATGTTtg                                                             <  1:191502‑M1/84‑1 (MQ=255)
                                                                                 ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtg                                                             <  1:361143‑M1/84‑1 (MQ=255)
                                                                                  taaatccacctgttaaaagcaggaatgatgctatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                                                            >  2:320866‑M1/65‑149 (MQ=255)
                                                                                    ccGAAAGAATGCGGTGGGCGATATCAAATCGCATGGCACTGGGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                          <  2:308872/149‑1 (MQ=38)
                                                                                    aacccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                          >  1:296553‑M1/63‑149 (MQ=255)
                                                                                         ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                     >  2:361143‑M1/58‑149 (MQ=255)
                                                                                                     gCGATATCAAATCGGATGGCACTGGGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcctgtttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttct                                         <  2:181675‑M1/149‑89 (MQ=255)
                                                                                                       gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgc                                       >  2:370686‑M1/1‑59 (MQ=255)
                                                                                                        gaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATc                                      >  1:465807‑M1/43‑149 (MQ=255)
                                                                                                         tatCAAATTGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCa                                     <  2:287197/149‑1 (MQ=255)
                                                                                                          atCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCAt                                    <  2:366073/149‑1 (MQ=255)
                                                                                                             aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaaga                                 >  1:344349‑M1/1‑53 (MQ=255)
                                                                                                                   gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                           >  2:19927/1‑149 (MQ=255)
                                                                                                                      gcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATGAGTAGTCATTCTGTGGGGGTg                        <  1:180449‑M1/121‑1 (MQ=255)
                                                                                                                      gcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTAATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTg                        <  1:180625‑M1/121‑1 (MQ=255)
                                                                                                                         tgttacttcttcatttacataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCTGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGaa                     >  1:44154‑M1/26‑149 (MQ=255)
                                                                                                                          cTGTGCAAACGGCTCTCGAAAACCTCGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggccgaaatgcaatggtgccgctttttctgctgaagaatacccggaactg                    <  1:451337‑M1/149‑110 (MQ=255)
                                                                                                                                tcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCa              <  1:59231‑M1/131‑1 (MQ=255)
                                                                                                                                aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcactggtgccgctttttctgctgaagaatacccggaactggcaaag              >  1:150252‑M1/1‑34 (MQ=255)
                                                                                                                                aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagtatacccggaactgtcaaag              >  1:150250‑M1/1‑34 (MQ=255)
                                                                                                                                aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaag              >  1:420031‑M1/1‑34 (MQ=255)
                                                                                                                                 cttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCAt                                                                                 <  2:206332‑M1/64‑1 (MQ=255)
                                                                                                                                 cttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCAt                                                                                 >  1:206332‑M1/18‑81 (MQ=255)
                                                                                                                                 cttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTCTTGAGCCATCAGTAATCATTCAGTGGGGGTGGAATGCCGcat             <  2:438225‑M1/132‑2 (MQ=255)
                                                                                                                                      ctctCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttat        <  1:370686‑M1/149‑122 (MQ=255)
                                                                                                                                       tttcctaaacgGTCGCTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCAt       >  1:278827‑M1/12‑149 (MQ=255)
                                                                                                                                       tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgcaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatc       >  2:150768‑M1/1‑27 (MQ=255)
                                                                                                                                       tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatc       >  1:59045‑M1/1‑27 (MQ=255)
                                                                                                                                       tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaaga                                 >  1:178240‑M1/1‑27 (MQ=255)
                                                                                                                                       tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaaga                                 <  2:178240‑M1/123‑97 (MQ=255)
                                                                                                                                         tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccg     >  1:204466‑M1/1‑25 (MQ=255)
                                                                                                                                          cGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccga    <  1:485379‑M1/149‑126 (MQ=255)
                                                                                                                                           atgaaggttGGTTTGGGGGAAGGTTGAAGACTTGATGGAGGAACGGGTGTTTTAGGATGGATCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGc   <  2:178968/142‑1 (MQ=255)
                                                                                                                                            taaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGcc  <  1:258006‑M1/143‑1 (MQ=255)
                                                                                                                                            aaaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgaca  >  2:417397‑M1/1‑22 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCC  >  NC_000913_3_hsa_tpiA/1207644‑1207932

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.