New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA 1207790 =12 (0.320)16 (0.480) 16/248 1.2 68.6% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA 1209619 = 4 (0.120)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGC  >  NC_000913_3_hsa_tpiA/1209470‑1209765
                                                                                                                                                     |                                                                                                                                                  
tttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCACACCAAc                                                                                                                                                    >  2:123572/1‑150 (MQ=255)
               aGCCTTTGCCAGTTCCGTGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                                                                                                                     <  2:185075/150‑1 (MQ=255)
                     ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                                                                                                               >  1:207663‑M2/129‑150 (MQ=255)
                      gcattccacccccactgaatgattactgattgctcaataccaccaataataaccggtatagctatgaaaccggtgcgtcctaatgtagcctttgctgcattacgtctttctccttctcccaaaacaaCCTTTATAAAAACGAACAACTaa                                                                                                                              >  2:116084‑M2/128‑150 (MQ=255)
                                         tgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                           >  2:1106‑M2/109‑150 (MQ=255)
                                                  attgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                                                                                                  <  1:10562‑M2/51‑1 (MQ=255)
                                                         aacaccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                           <  2:207663‑M2/58‑1 (MQ=255)
                                                           gaccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttcccccaaaccaaCGTTTATGAAAATTAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                                                                                         <  2:227259‑M2/60‑1 (MQ=255)
                                                              caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                      >  2:478385‑M2/88‑150 (MQ=255)
                                                                      ataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                                                                              <  2:185718‑M2/71‑1 (MQ=255)
                                                                      agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                                                                              >  1:339753/1‑150 (MQ=255)
                                                                          ccggtatagctatggaagcggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                                                                          <  2:91719‑M2/75‑1 (MQ=255)
                                                                                tagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttactccttctcccaaacaaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGCCATCATTCCTGCTTTTACCAGGGGGATTTAACATCCTTATTGGCTAt                                                                    >  1:9790‑M2/70‑150 (MQ=255)
                                                                                  gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAACTGGCATCATTCCTGCTTTTACCAGGGGGAGTTAACATGCTTATTGCCCATGt                                                                  >  2:573846‑M2/68‑150 (MQ=255)
                                                                                         aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGATTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTGATTGGCTATGTACGCGTa                                                           >  1:403286‑M2/61‑150 (MQ=255)
                                                                                            cggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATca                                                        >  1:454420‑M2/58‑150 (MQ=255)
                                                                                                  gtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                  <  2:88561‑M2/99‑1 (MQ=255)
                                                                                                     ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAGCATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                               >  1:351323‑M2/49‑150 (MQ=255)
                                                                                                         aaCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGATACGTACGCTTATCAACAAATGACCAGa                                           >  2:275055/1‑150 (MQ=255)
                                                                                                          gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                          <  2:351323‑M2/107‑1 (MQ=255)
                                                                                                                                      cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAACAGGATTCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAAAACAGATCTCCAAGGTAATGCGCTGAAc              >  1:625298‑M2/16‑150 (MQ=255)
                                                                                                                                      cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCCGGGGGGTTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACAAAGATCTACAACGTAAGCCGCTgaaa              >  2:188526‑M2/16‑149 (MQ=255)
                                                                                                                                           cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                    >  2:125428‑M2/11‑75 (MQ=255)
                                                                                                                                           cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                    <  1:125428‑M2/65‑1 (MQ=255)
                                                                                                                                            ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTCCAACGTAATGCGCTGAACTGTgcg        >  1:221826‑M2/10‑149 (MQ=255)
                                                                                                                                            ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa        <  1:275055‑M2/141‑1 (MQ=255)
                                                                                                                                                  taaCTTTTATTAAAATGAAGAAATAACAAGCACATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGc  <  1:75250‑M2/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGC  >  NC_000913_3_hsa_tpiA/1209470‑1209765

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.