New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA 1207790 =7 (0.150)42 (1.030) 37/248 0.2 85.5% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA 1209619 = 8 (0.200)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_hsa_tpiA/1209481‑1209767
                                                                                                                                          |                                                                                                                                                    
gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                                                                                                            >  1:581843/1‑149 (MQ=255)
          ggcattccaccccgactcaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                                  <  1:791619‑M2/21‑1 (MQ=255)
          ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                                                                                                                   >  1:594118‑M2/129‑148 (MQ=255)
          ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                                                                                                                   <  2:594118‑M2/20‑1 (MQ=255)
             atcccaaccccactgaaggattactgattgctcaagtccacccatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                                                                               <  1:731065‑M2/24‑1 (MQ=255)
                     cccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAGGAAGAAATAACAAGCAAAt                                                                                                                       >  1:43311‑M2/118‑149 (MQ=255)
                      ccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                                      >  1:109289‑M2/117‑149 (MQ=255)
                      ccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                                      >  2:107604‑M2/117‑149 (MQ=255)
                             atgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                                                                               >  2:639786‑M2/110‑133 (MQ=255)
                             atgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                                                                               <  1:639786‑M2/24‑1 (MQ=255)
                               gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                                                                               >  1:406126‑M2/108‑131 (MQ=255)
                               gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                                                                               <  2:406126‑M2/24‑1 (MQ=255)
                                attactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                                                                            >  1:237484‑M2/107‑149 (MQ=255)
                                 ttactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGTAGAAATAACAAGCAAATGGCATCATTCCt                                                                                                           >  1:675506‑M2/106‑149 (MQ=255)
                                    ctgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                                                                                                        <  1:66531‑M2/47‑1 (MQ=255)
                                       attgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                                                                                     >  1:389444‑M2/100‑149 (MQ=255)
                                           ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                                                                                                      >  1:228742‑M2/96‑112 (MQ=255)
                                           ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                                                                                                      <  2:228742‑M2/17‑1 (MQ=255)
                                           ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                                 >  1:645075‑M2/96‑149 (MQ=255)
                                             caataccagcaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                               <  2:150609‑M2/56‑1 (MQ=255)
                                             cACCATTGCAGTTCAGCCAGCCTGTTGGCGTAGTGGCTGAAGGCCACCGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                               <  2:417236/149‑1 (MQ=255)
                                               ccAGTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCATCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                             <  2:476514/149‑1 (MQ=255)
                                                  ccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATACCAAGCAAATGGCATCATTCCTTCTTTTACCAGGGGGAt                                                                                          >  2:703651‑M2/89‑149 (MQ=255)
                                                   caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt                                                                                         >  1:222732‑M2/88‑149 (MQ=255)
                                                      caatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                      >  1:509303‑M2/85‑149 (MQ=255)
                                                       tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACAGTCGCTCCAAGAGGTACTACAGCGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTGTACCAGGGGGATTTAAc                                                                                     >  1:603880/1‑149 (MQ=255)
                                                        atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgcgctttctcccaaaccaaCGTTTATGAAAATTAAGGAATAACAAGCAAATGGCATCATTCCTGCTTTTGCGAGGGGGTTTTAACa                                                                                    >  1:136482‑M2/83‑149 (MQ=255)
                                                           ataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                                                                                 <  2:765539‑M2/70‑1 (MQ=255)
                                                            taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                                                                                                     >  1:704157‑M2/79‑112 (MQ=255)
                                                            taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                                                                                                     <  2:704157‑M2/34‑1 (MQ=255)
                                                               ccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                                                                             >  1:401965‑M2/76‑149 (MQ=255)
                                                                cggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGTTTTAACATGcttctt                                                                            >  1:717002‑M2/75‑149 (MQ=255)
                                                                                ccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                            <  2:401965‑M2/91‑1 (MQ=255)
                                                                                   gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                         <  1:107604‑M2/94‑1 (MQ=255)
                                                                                    tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                        <  2:222732‑M2/95‑1 (MQ=255)
                                                                                        tcctaatgtcgccgttgctgcattgagttttgctccttctcccaaaccacCGTTTATGAAAATGAAGAACTAACAAGCAACTGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAAGCTCAATGGCTTCGTGCGCGTATCAAccaac                                                    >  2:435411‑M2/51‑148 (MQ=255)
                                                                                        tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGTAAATGGCATCATTCCTGCTTTTAACAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                    <  1:68752‑M2/99‑1 (MQ=255)
                                                                                        tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                    <  2:228013‑M2/99‑1 (MQ=255)
                                                                                          ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGGAGAAATAACAAGCAAATGGCATCATTCCTGCGTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACCAGCATCAACAaatta                                                  >  1:702349‑M2/49‑147 (MQ=255)
                                                                                          aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCACCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa                                                  <  2:581843/149‑1 (MQ=255)
                                                                                                      ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                            <  2:255918‑M2/11‑1 (MQ=255)
                                                                                                      ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                            >  1:255918‑M2/37‑47 (MQ=255)
                                                                                                       acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                                     >  1:650451/1‑149 (MQ=255)
                                                                                                        gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATc                                    <  2:675506‑M2/115‑1 (MQ=255)
                                                                                                         ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCCTCATTCCTGCTTTTACCAGGGAGATTTAACATGCTTGTTGTCTATGTACGCGTATCAACAAATGCCAAGAACACAGctct                                   >  2:616746‑M2/34‑149 (MQ=255)
                                                                                                           aCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTAc                                 >  2:483612/1‑149 (MQ=255)
                                                                                                             attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACaa                               <  2:434456‑M2/120‑1 (MQ=255)
                                                                                                              ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAACCGTATAAGAATCTAACTATAATACAGATCTACAAc                              >  2:206234‑M2/29‑149 (MQ=255)
                                                                                                              ggTAATGCAGAGCCTTTTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                              <  1:682051/149‑1 (MQ=255)
                                                                                                               tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCATGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACg                             >  1:303830‑M2/28‑149 (MQ=255)
                                                                                                                          ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACACGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGa                  >  2:425168‑M2/17‑149 (MQ=255)
                                                                                                                                 ccaaaccaaCGCTTATGAAAATGCAGAAATCACAAGCAATCGGCATCTTTCCTCCTTTTACCAGGGGGATTTAACATGCTTATTCGCTTGGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc           <  1:481917‑M2/140‑1 (MQ=255)
                                                                                                                                      ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGAt      >  2:530240‑M2/5‑149 (MQ=255)
                                                                                                                                       caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTATACTGTGCAGGATg     <  2:509303‑M2/146‑1 (MQ=255)
                                                                                                                                          cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGa  <  1:565754/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_hsa_tpiA/1209481‑1209767

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.