New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA 1207790 =8 (0.130)45 (0.840) 38/250 0.5 88.2% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA 1209619 = 5 (0.090)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_hsa_tpiA/1209475‑1209767
                                                                                                                                                |                                                                                                                                                    
aggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                                                                                                                  >  1:961186‑M2/145‑149 (MQ=255)
  gcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaacggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCATTTAt                                                                                                                                                >  1:47677‑M2/143‑149 (MQ=255)
     gatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccacCGTTTACGaa                                                                                                                                             >  2:161909‑M2/140‑149 (MQ=255)
       tgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                                                                                                           <  2:510573‑M2/12‑1 (MQ=255)
                       gaagcccactgaaagattactgattgcgcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgcgtcatggagttttgctccttctccaacaccaaCGTTCATTAAAATGAAGAAATAACCAGc                                                                                                                           >  2:967316‑M2/122‑149 (MQ=255)
                       cacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                                  <  1:834123‑M2/21‑1 (MQ=255)
                       cacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                                  >  2:834123‑M2/122‑142 (MQ=255)
                           cccactgaatgattactgattgctcaataccaccactcataacagaaatagcaatgaaaccggtgcgtcctaatgtaaccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAACAAATAACAAGCAAAt                                                                                                                       >  2:440924‑M2/118‑149 (MQ=255)
                               ctgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                                                   <  2:842875‑M2/36‑1 (MQ=255)
                                tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                                  <  1:691442‑M2/37‑1 (MQ=255)
                                 gaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                                                                                                 <  1:576605‑M2/38‑1 (MQ=255)
                                    cagcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCATACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                                                                              >  1:1047637/1‑149 (MQ=255)
                                     agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                             <  2:1047637/149‑1 (MQ=255)
                                            gattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGACATCATTCCTGCttt                                                                                                      >  2:720565‑M2/101‑149 (MQ=255)
                                             attgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                                                                                     >  2:347558‑M2/100‑149 (MQ=255)
                                                gctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc                                                                                                  >  1:552018‑M2/97‑149 (MQ=255)
                                                  tcaataccaccaatcataaccggtatagctctgacaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTtacccg                                                                                                >  1:115025‑M2/95‑147 (MQ=255)
                                                      taccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagcctgtgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                                            >  1:731231‑M2/91‑149 (MQ=255)
                                                                      cggtatagctatgaagccgttgcgtcctaatgtagccgttgctgcagtgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                                                                            <  2:308990‑M2/75‑1 (MQ=255)
                                                                          atagcttggaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                                        <  1:99491‑M2/79‑1 (MQ=255)
                                                                            agatatgaaaccggttcgtcctaaggtagccgttgctgcatcgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                                                      <  2:383084‑M2/81‑1 (MQ=255)
                                                                             gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGCTTTAACATGCTTCTTGGCTATg                                                                     >  2:162506‑M2/68‑149 (MQ=255)
                                                                             gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATg                                                                     >  2:118461‑M2/68‑149 (MQ=255)
                                                                              ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                                  <  1:546104‑M2/21‑1 (MQ=255)
                                                                              ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                                  >  2:546104‑M2/67‑87 (MQ=255)
                                                                               tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGCCACCATTCCTGCTTGTACCAGGGGGCTGTAAAAGACTTATTGGCGCTGTa                                                                   >  2:407344‑M2/66‑149 (MQ=255)
                                                                                atgataccggtgggtcctaatgtagccgttgctgcattgagttttgctccgtctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTGAGTAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                                  <  2:227343‑M2/85‑1 (MQ=255)
                                                                                atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                                  >  1:170465‑M2/65‑149 (MQ=255)
                                                                                         gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGATGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                         <  1:967085‑M2/94‑1 (MQ=255)
                                                                                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                                     <  2:835996/149‑1 (MQ=255)
                                                                                                ctaatgtagccgttgctgcattgagtttggccccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGCCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa                                                  <  1:486987‑M2/101‑1 (MQ=255)
                                                                                                  aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAcc                                                <  1:169833‑M2/103‑1 (MQ=255)
                                                                                                   atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGGTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAATCGTATCAACAAATAACCa                                               >  2:730780‑M2/46‑149 (MQ=255)
                                                                                                     gtagccgttgctgcattgagttttggtccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                        >  1:950483‑M2/44‑138 (MQ=255)
                                                                                                     gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                        <  2:950483‑M2/95‑1 (MQ=255)
                                                                                                        gggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                                          >  1:126813/1‑149 (MQ=255)
                                                                                                         ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                    <  2:713146‑M2/67‑1 (MQ=255)
                                                                                                         ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                    >  1:713146‑M2/40‑106 (MQ=255)
                                                                                                             tggtgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                                     <  1:658431‑M2/114‑1 (MQ=255)
                                                                                                             tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                                     <  1:382410‑M2/114‑1 (MQ=255)
                                                                                                             tgctgcatggagttttgctccttctcccaaaccaaCGTTTATCAAAATGAAGAAATAACAATCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                                     <  1:720565‑M2/114‑1 (MQ=255)
                                                                                                                tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTa                                  >  2:437309‑M2/33‑149 (MQ=255)
                                                                                                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa                          >  1:137867/1‑149 (MQ=255)
                                                                                                                         ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt                         >  2:587595‑M2/24‑149 (MQ=255)
                                                                                                                         ctttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGAGGATTTAACATGCTCATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt                         >  2:796514‑M2/24‑149 (MQ=255)
                                                                                                                          tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATg                        <  2:552018‑M2/127‑1 (MQ=255)
                                                                                                                              ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTGCAACGTAAt                         <  2:10304‑M2/126‑1 (MQ=255)
                                                                                                                              ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt                         >  1:10304‑M2/19‑144 (MQ=255)
                                                                                                                              ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt                    >  1:1014680‑M2/19‑149 (MQ=255)
                                                                                                                               tccttctcccaaaccaaCGTTTATGAAAATGAAGGAATAACAAGCAAATGGCATCATTCCTGCCTTTCCCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTg                   >  1:296092‑M2/18‑149 (MQ=255)
                                                                                                                                ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGa                  >  1:143214‑M2/17‑149 (MQ=255)
                                                                                                                                 cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa                 >  1:32646‑M2/16‑149 (MQ=255)
                                                                                                                                  ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAc                >  1:133448‑M2/15‑149 (MQ=255)
                                                                                                                                      cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCCCTGAACtgtg            >  2:579374‑M2/11‑149 (MQ=255)
                                                                                                                                       ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATAAACAAATGACCATAACACAGATCTACAACGTAATGCGCTAACCTGTGc           >  2:187368‑M2/10‑149 (MQ=255)
                                                                                                                                        caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAACAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGATGAACTGTGCa          >  1:369881‑M2/9‑149 (MQ=255)
                                                                                                                                              aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGc    >  1:905098‑M2/3‑149 (MQ=255)
                                                                                                                                                cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGa  <  1:118461/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_hsa_tpiA/1209475‑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.