New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA 1207790 =3 (0.080)26 (0.760) 26/256 0.4 81.7% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA 1209619 = 9 (0.260)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913_3_hsa_tpiA/1209473‑1209759
                                                                                                                                                  |                                                                                                                                            
atttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGAGACCAACAGGTAATGCAGAGCCTTCTCCCCAACCAACGt                                                                                                                                            >  1:354957/1‑149 (MQ=255)
         tgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTTTGaaaa                                                                                                                                   >  2:73185‑M2/138‑149 (MQ=255)
              ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                                                                                              >  2:434984/1‑149 (MQ=255)
                  ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGGTTATTAAAATGAGGAAAt                                                                                                                          >  1:287398‑M2/129‑149 (MQ=255)
                   gcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                                                                                                         <  2:440864‑M2/22‑1 (MQ=255)
                         cacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                                                                                   >  1:451939‑M2/122‑149 (MQ=255)
                          acccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAAAAGCa                                                                                                                  >  1:312773‑M2/121‑149 (MQ=255)
                              gTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                              <  1:533348/149‑1 (MQ=255)
                                  tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                          >  1:282309‑M2/113‑149 (MQ=255)
                                  tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                          <  1:486092‑M2/37‑1 (MQ=255)
                                          agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                                  <  1:157347/149‑1 (MQ=255)
                                           gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                 >  2:465521/1‑149 (MQ=255)
                                            ctgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCCTTCCTGCt                                                                                                >  1:470923‑M2/103‑149 (MQ=255)
                                             tgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                               >  1:71066‑M2/102‑149 (MQ=255)
                                              gattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                              <  2:585544‑M2/49‑1 (MQ=255)
                                                    tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAACTAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                        >  1:354049‑M2/95‑149 (MQ=255)
                                                          gcaccaaacaaaaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                                                                                  <  2:287398‑M2/61‑1 (MQ=255)
                                                                 tcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttctgctccttctcccaaaccaaCTTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                                                                           <  2:37320‑M2/68‑1 (MQ=255)
                                                                        cggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAACCAAATGGCATCATTCCTGCTTTTACAAGGGGGATTTAACATGCTTAtt                                                                    >  2:220590‑M2/75‑149 (MQ=255)
                                                                         ggtatagctatgaaaccggggcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                                   <  2:1500‑M2/76‑1 (MQ=255)
                                                                                ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                            >  2:621325‑M2/67‑149 (MQ=255)
                                                                                 tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgacccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGTTTATTGGCTATGTa                                                           >  1:76568‑M2/66‑149 (MQ=255)
                                                                                      aaccggtgggtactaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGACTATGTACGCGt                                                      >  2:68096‑M2/61‑149 (MQ=255)
                                                                                          tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCCACGTTTATGAAAATGAAGAAGTAACAAGCAAATAGCATCATTCCTGCTTTTACCAAGGGGATTTAACATACTTATTGGCTATCTCGGCGTATca                                                  >  2:331350/1‑149 (MQ=255)
                                                                                                 cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaagaa                                                                                                                            >  1:537632‑M2/50‑68 (MQ=255)
                                                                                                 cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaagaa                                                                                                                            <  2:537632‑M2/19‑1 (MQ=255)
                                                                                                  ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGTGATTTAACATGCTTATTGTCTATGTACGCGTATCAACAAATGa                                          >  2:46693‑M2/49‑149 (MQ=255)
                                                                                                     atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGAGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATCaccc                                       >  2:58078‑M2/46‑148 (MQ=255)
                                                                                                         cGGGACACCAACAGGTAATGCAGAGCCTTCTCCCACACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAac                                   <  1:315258/149‑1 (MQ=255)
                                                                                                         agccgttgctgcattgagttttgctccttcccccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAac                                   <  1:43693‑M2/108‑1 (MQ=255)
                                                                                                                ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATc                            >  1:622085/1‑149 (MQ=255)
                                                                                                                      ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                      >  1:63777‑M2/29‑149 (MQ=255)
                                                                                                                       gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACTAGCACCCAGATCTACAACg                     <  2:565931/149‑1 (MQ=255)
                                                                                                                          gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGACCACAGATCTACAACGTaa                  <  2:63777‑M2/125‑1 (MQ=255)
                                                                                                                                ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt            >  1:401823‑M2/19‑149 (MQ=255)
                                                                                                                                         ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAGGACCAGAACACAGATCTACAACGTAATGCGCTGGACTGTGc   >  1:28558‑M2/10‑149 (MQ=255)
                                                                                                                                          caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa  <  2:354957‑M2/141‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
ATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913_3_hsa_tpiA/1209473‑1209759

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.