New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 12078057 (0.220)12 (0.420) 11/244 NT 58.3% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA = 1209602 11 (0.390)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATAC  >  NC_000913_3_hsa_tpiA/1207657‑1207952
                                                                                                                                                    |                                                                                                                                                   
gCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                                                                                                     >  1:52672‑M1/1‑149 (MQ=255)
   tCCGGTGCCAGGCTGGAGGAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCGGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctc                                                                                                                                                  <  2:387987‑M1/149‑4 (MQ=255)
       attacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCcaaa                                                                                                                                              <  1:46684‑M1/23‑1 (MQ=255)
                gggggagggggcggcaggGCCATTTAGTTACGGCAGTTCGCGCGCTTCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattagctgtt                                                                                                                                     <  2:38056‑M1/142‑17 (MQ=255)
                          ttctggtcatttgttgatacgcggacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTa                                                                                                                           >  1:355347‑M1/108‑149 (MQ=255)
                                tcatttgttgatacgcgtacatagcgaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAg                                                                                                                     >  1:353606‑M1/102‑149 (MQ=255)
                                   tttggtgatacgcctacatagccaacaagaatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaccgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAc                                                                                                                  >  1:607198‑M1/99‑149 (MQ=255)
                                      gttcatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                                                                                                               >  2:504407‑M1/96‑149 (MQ=255)
                                               cGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCACATCGGATGGACCTGTGCAAACGGCTCTCGAAACCCTTGGTTTGGGAGAAGActctgcattacctgttggtgtcccgattccgtgcccttcagccactc                                                                                                      >  2:373572‑M1/1‑102 (MQ=255)
                                                gCACTTCACGCGCTGCTGATAAGGCGCAAGAATCCGTTTGGCGATATCAAATCGGCTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGATAAGGctctgcattagctgttggtgtccaggttccgtggccttcagccactcc                                                                                                     >  1:569855‑M1/1‑101 (MQ=255)
                                                 tacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                                                                                    >  1:536993‑M1/85‑149 (MQ=255)
                                                    ttCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAACCCTTGGTTTGGGAGAAGCAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt                                                                                                 >  1:570341/1‑149 (MQ=255)
                                                                 atgttaaatccccatggtaaaagcaggaatgatgccatttgattgttatgtattcattttcacaaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGa                                                                                    >  1:579284‑M1/69‑149 (MQ=255)
                                                                   gaagCCGCAAGAATCCGTTTGGGGATATCAAATCGGATGGCAATGGGCAAACGGCTCTCGAAAACGTTGGTTTGGGCGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaa                                                                                  <  2:145373‑M1/148‑68 (MQ=255)
                                                                    ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAGAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTGTAACGGTTATGATtg                                                                                 <  2:268119‑M1/84‑1 (MQ=255)
                                                                       ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggggtcccggttccgtggccttcagccactccgccaacaggctgcctgaaatgca                                                                              >  2:473918‑M1/1‑78 (MQ=255)
                                                                            aGAAGCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                                         <  2:192956/149‑1 (MQ=255)
                                                                                     aagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCa                                                                >  1:443799‑M1/49‑149 (MQ=255)
                                                                                       cacggaatggagccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGAACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGt                                                              >  1:53876‑M1/47‑149 (MQ=255)
                                                                                          gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgc                                                           >  1:200130‑M1/1‑59 (MQ=255)
                                                                                            tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAAACTTGGTTTGAGAGAAGGAGCAAAACACAATGCAGCAACGGCTAAACTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCCATCAGtcatca                                                         >  2:328829/1‑149 (MQ=255)
                                                                                              tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaa                                                       <  2:508479‑M1/149‑95 (MQ=255)
                                                                                                         ggCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatgcccgga                                            >  1:559545‑M1/1‑44 (MQ=255)
                                                                                                           cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctcaagaatacccggaac                                          >  1:203319‑M1/1‑42 (MQ=255)
                                                                                                           cACTGTGCAAACGGCTCTCGAAAACCTTGGATTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaccaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaac                                          >  2:604380‑M1/1‑42 (MQ=255)
                                                                                                            aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTTGaa                                         <  1:458912/149‑1 (MQ=255)
                                                                                                               gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcatgacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggc                                      <  1:255733‑M1/149‑112 (MQ=255)
                                                                                                                    cttcattttcataaacgTTGGTATGGCAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACAGGTTATGGTTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGAGGTGGAGTGCCGcac                                 <  2:38967‑M1/132‑2 (MQ=255)
                                                                                                                        gCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCGGTAATCATTCAGTGGGGGGGGAATGCCGCAAAAga                             >  1:293938/1‑148 (MQ=255)
                                                                                                                          tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatc                           >  1:28136‑M1/1‑27 (MQ=255)
                                                                                                                           ttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATc                          <  1:208785‑M1/139‑1 (MQ=255)
                                                                                                                                    cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                 >  1:565272/1‑149 (MQ=255)
                                                                                                                                    cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGATATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGACGCAAAAGCATCTGCCTCTgg                 >  2:480193/1‑149 (MQ=255)
                                                                                                                                         tttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgc            >  1:595617‑M1/1‑12 (MQ=255)
                                                                                                                                           tGGGAGTAGGctcttccttaccggttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaagacccggaactggcaaaggcttatccgacaaataaattgcct          <  2:259756‑M1/149‑140 (MQ=255)
                                                                                                                                            gggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggccaaggcttatccgacaaataaattgcctg         >  1:148662‑M1/1‑9 (MQ=255)
                                                                                                                                                   ggctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacacataaattgcctgatttacg  >  1:231573‑M1/1‑2 (MQ=255)
                                                                                                                                                   ggctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacg  >  1:465908‑M1/1‑2 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATAC  >  NC_000913_3_hsa_tpiA/1207657‑1207952

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