New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA 1207790 =7 (0.200)16 (0.520) 15/244 NT 66.4% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA 1209619 = 10 (0.330)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGC  >  NC_000913_3_hsa_tpiA/1209482‑1209765
                                                                                                                                         |                                                                                                                                                  
tgcttctgcagtattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                                                                                                         <  2:477713‑M2/12‑1 (MQ=255)
   ttttgcggccttccacccccactgactgatttcttattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                                                                                                      <  1:172668‑M2/15‑1 (MQ=255)
   ttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                                                                                                      <  2:147565‑M2/15‑1 (MQ=255)
         tGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt                                                                                                                                <  1:553791/149‑1 (MQ=255)
         tGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt                                                                                                                                <  1:552760/149‑1 (MQ=255)
              tccaccgcacctgaatgattactgattgctcaatatcaccaatcataaccggtatagctatgaaaccggtgcgtcctaatttagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                           <  2:289927‑M2/26‑1 (MQ=255)
                   ccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaacaaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                                                                                                      >  2:339517‑M2/119‑149 (MQ=255)
                         caaatgattactgattgatcaataccaccaatcatacccgttatagctatgaaacctgtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                                <  2:299170‑M2/37‑1 (MQ=255)
                           ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                                                                              >  1:602492/1‑149 (MQ=255)
                              agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATGATTc                                                                                                           >  1:10150/1‑149 (MQ=255)
                                    tgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCAGCATTCCTGCtt                                                                                                     >  2:26597‑M2/102‑149 (MQ=255)
                                           tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAGGGCATCATTCCTGCTTTTACCAg                                                                                              <  1:557086‑M2/55‑1 (MQ=255)
                                              ataccaccaatcataaccggtacagctatgaacccggtgcgtcctaatgtagcctttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAACTGAAGAAATAACAAGCACAGGGCATCATTCCTGCTTTTACCAgggg                                                                                           <  2:509668‑M2/58‑1 (MQ=255)
                                               taccaccaatcataaccggtatagctatgaaaccggtgcgtcctactgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAACTGAAGAAATAACAAGCAACTGGCATCATTCCTGCTTTTACCAggggg                                                                                          >  2:29973‑M2/91‑149 (MQ=255)
                                                        gcagccagcCTGTTGGCGGAGTGACTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                                                                                 <  2:135872/148‑1 (MQ=255)
                                                           gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCCGGGGGATTTAACATGCt                                                                              >  1:81185/1‑149 (MQ=255)
                                                           gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                                                                              >  1:81047/1‑149 (MQ=255)
                                                                        tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCTTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                                                                 <  2:594762‑M2/84‑1 (MQ=255)
                                                                               ccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                          <  1:110370‑M2/91‑1 (MQ=255)
                                                                                      gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                                   <  2:629278/149‑1 (MQ=255)
                                                                                        cACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATg                                                 >  1:505905/1‑149 (MQ=255)
                                                                                         ctaatgtagccgttgctgcattgagttttgctgcttcccccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGTTGTTAGCATTCTTATTGGCTATGTACGCGTATCAACAAATGa                                                <  1:156762‑M2/101‑1 (MQ=255)
                                                                                                  ggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTGGACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCACCAAATGACCAGAacac                                       >  1:620086/1‑149 (MQ=255)
                                                                                                   cgtttctgcattgagtcttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAAcaca                                      <  1:334698‑M2/111‑1 (MQ=255)
                                                                                                         tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATATACGCGTATCAATAAATGACCAGAACACAGAtata                                >  1:577208‑M2/33‑149 (MQ=255)
                                                                                                                 gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa                        <  1:199307‑M2/125‑1 (MQ=255)
                                                                                                                       ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACAAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAGGACCAGAACACAGATCTACAACGTAATGCGCt                  >  1:75611‑M2/19‑149 (MQ=255)
                                                                                                                        tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATGAACAAATGACCAGAACCCAGATCTACAACGTAACTCGCTg                 >  2:145281‑M2/18‑149 (MQ=255)
                                                                                                                        tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACACCATCatgcgcgg                 >  1:118371‑M2/18‑147 (MQ=255)
                                                                                                                          cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGACATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACCACGTAATGCGCTGaa               >  1:188236‑M2/16‑149 (MQ=255)
                                                                                                                           ttctcccaaaccaaCGTCTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAc              <  1:466829‑M2/135‑1 (MQ=255)
                                                                                                                                ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGCTCTACAACGTAATGGGCTCAACTGTGc         >  2:376706‑M2/10‑149 (MQ=255)
                                                                                                                                    accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGAATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGa     <  2:362743‑M2/144‑1 (MQ=255)
                                                                                                                                       aaCGTTTATGAAAATGAAGAAATAACAAGCAAATTGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGc  <  2:483102‑M2/147‑1 (MQ=255)
                                                                                                                                       aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATTGACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGACCTgtgcagagtgc  >  2:405427‑M2/3‑144 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
ATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGC  >  NC_000913_3_hsa_tpiA/1209482‑1209765

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