New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 12078055 (0.210)22 (1.040) 20/254 0.2 83.9% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA = 1209602 4 (0.190)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913_3_hsa_tpiA/1209454‑1209748
                                                                                                                                                    |                                                                                                                                                  
cacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGGTCCGGGTATTCTTCAGCAGAAACAGCGGAACCATTGTCTTTAAGTTGGCCTTTTGGCGGAGTGGCTGAAGTCCACGGCACCGGCACACGGACCGGTAATGCagag                                                                                                                                                    >  1:98099/1‑149 (MQ=255)
 aCGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                                                                                                   >  1:216366‑M2/1‑148 (MQ=255)
             cAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                                                                                                       >  1:181256‑M2/1‑136 (MQ=255)
                tttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaa                                                                                                                                    >  2:389283‑M2/1‑133 (MQ=255)
                         cccacaagccttAGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                                                                           >  2:400676‑M2/6‑124 (MQ=255)
                         cGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                                                                           >  1:166750‑M2/1‑124 (MQ=255)
                           gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                         <  2:74829‑M2/149‑28 (MQ=255)
                                   cttgCCAGGTCTGGGTAGTCTTCAGCAGAAAAAGCGGCACCATTTCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttg                                                                                                                 >  2:168101‑M2/2‑114 (MQ=255)
                                                        tcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgat                                                                                            <  1:386099‑M2/149‑57 (MQ=255)
                                                                            aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgc                                                                        >  2:386099‑M2/1‑73 (MQ=255)
                                                                             ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcccagtgccatccgatttgatatcgccaaacggattcttgcg                                                                       <  1:359592‑M2/149‑78 (MQ=255)
                                                                                        ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGACTGCTCCCGAAACAACGATTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAACAGGGGGATTAAACATGcgt                                                            >  2:335984/1‑147 (MQ=255)
                                                                                                      gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcg                                              >  2:343494‑M2/1‑47 (MQ=255)
                                                                                                       aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcgg                                             >  1:74829‑M2/1‑46 (MQ=255)
                                                                                                             tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaa                                                                 <  1:218358‑M2/123‑84 (MQ=255)
                                                                                                             tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaa                                                                 >  2:218358‑M2/1‑40 (MQ=255)
                                                                                                                 ggCCACGGAACCGGGACACCAACCGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactg                                   >  1:286049‑M2/1‑36 (MQ=255)
                                                                                                                    cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgggttcttgcggcgtaacagcagcccgccaagtccggtaagtaagtggtc                                <  1:150086‑M2/149‑117 (MQ=255)
                                                                                                                    cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtc                                >  2:162616‑M2/1‑33 (MQ=255)
                                                                                                                         aaCCGGGACACCAACAGGTAATGCAGAGccttctcccaaacaaaggttgtcgagagccgtttgcacagtgccatccgatttgatatcgcccaacattgacttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgcc                           <  1:412240‑M2/149‑122 (MQ=255)
                                                                                                                             gggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAAAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                       >  1:365950/1‑149 (MQ=255)
                                                                                                                                   ccAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctc                 <  1:367512‑M2/149‑132 (MQ=255)
                                                                                                                                    cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaactacctggtcgtgccgccccttctcc                <  1:140723‑M2/149‑133 (MQ=255)
                                                                                                                                          gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcgccttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctg          <  2:254739‑M2/149‑139 (MQ=255)
                                                                                                                                           tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctgg         >  1:417184‑M2/1‑10 (MQ=255)
                                                                                                                                           tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctgg         <  1:395372‑M2/149‑140 (MQ=255)
                                                                                                                                           tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctgg         <  2:416880‑M2/149‑140 (MQ=255)
                                                                                                                                                 agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                    <  2:315213/132‑1 (MQ=255)
                                                                                                                                                 agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                    >  1:315213/1‑132 (MQ=255)
                                                                                                                                                  gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccgga  <  1:213843‑M2/149‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913_3_hsa_tpiA/1209454‑1209748

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.