New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA 1207790 =7 (0.210)20 (0.660) 19/248 0.7 67.6% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA 1209619 = 13 (0.430)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913_3_hsa_tpiA/1209475‑1209760
                                                                                                                                                |                                                                                                                                             
ttGTCGGATAAGCCTTTGCCAGTGCCGGGTATACTTCAGCAGAAAAAGCGGCACCATTGTATTTCAGCCAGACTGTTGGCGGAGTGGCTAAAGGCCACTGAACCGGGACATCAAGAGGTAATGCAGAGCCTTCTCCCAAACCAACCttt                                                                                                                                           <  2:174963/149‑1 (MQ=255)
    agatgcttttgcggcattccacccccactgaatgattactgattgctcaatcccaccaatcataaccggtcaagctatgacaccggtgcgtcctaatttagcccttgctccattgagttttactccttctcccaaaccaaCGTTcatca                                                                                                                                       >  2:471207‑M2/141‑147 (MQ=255)
      gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGGTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                                                                                                     >  1:443292/1‑149 (MQ=255)
          ttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAaataa                                                                                                                                 >  2:82990‑M2/135‑147 (MQ=255)
                ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                           >  1:376175‑M2/129‑149 (MQ=255)
                ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgcggcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                           >  1:350952‑M2/129‑149 (MQ=255)
                  cattctacccccactgaatgattactgattgctaaataccaccaatcataaccggtatagctatgaaaccggtgcggcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaa                                                                                                                         <  1:192297‑M2/23‑1 (MQ=255)
                     tccacccccactgaatgattactgattgcccaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                      <  1:306343‑M2/26‑1 (MQ=255)
                            ccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                               <  1:53736‑M2/33‑1 (MQ=255)
                                  ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                                                                         >  1:122678/1‑149 (MQ=255)
                                  ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAATGACTGAAGGACACGGAACCGGGACACCAACCGGTAATTCATAGCCTTCTCCAAAACCCACGTTACTGAAAATTAAGAAATATATAGCAAAAGGcatcc                                                                                                         >  2:62264/1‑148 (MQ=255)
                                     gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                      >  1:32718‑M2/108‑149 (MQ=255)
                                      attactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                                                                     >  1:205152‑M2/107‑149 (MQ=255)
                                        tactcattgctcaatacctacaatcataaccggtggagctatgaaaccggttcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                                   >  2:504906‑M2/105‑149 (MQ=255)
                                                    aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                       <  1:229348‑M2/57‑1 (MQ=255)
                                                          cccaagcatagccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                 <  1:160148‑M2/63‑1 (MQ=255)
                                                              ttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                             <  2:122678/149‑1 (MQ=255)
                                                                     ccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                                                                      <  1:283271‑M2/74‑1 (MQ=255)
                                                                         tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                  >  1:305204/1‑149 (MQ=255)
                                                                               tatgaaaccggtgcggcctaagttagccgctgctgcattgtgttttgctccttctcccaaaccaaCTTTTATGAAAATGAAGAGATAACTAGCAAATGGCATCATGCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                                                            <  2:402172‑M2/84‑1 (MQ=255)
                                                                                    aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGATTATTGCCAATATACGCGt                                                       >  2:239506‑M2/61‑149 (MQ=255)
                                                                                     ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCCTATTCCATGGGGATTTAACATGCTTATTGGCTATGTACGCGTa                                                      <  1:296828/149‑1 (MQ=255)
                                                                                              tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                             <  1:322464‑M2/99‑1 (MQ=255)
                                                                                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                             >  1:315547/1‑149 (MQ=255)
                                                                                                aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                                                                                                                    >  1:265706/1‑76 (MQ=255)
                                                                                                aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                                                                                                                    <  2:265706/76‑1 (MQ=255)
                                                                                                aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa                                           >  1:448194/1‑149 (MQ=255)
                                                                                                 cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATTCTTATTGGCTATGTACGCGTATCAACAAATGAc                                          >  1:220323/1‑149 (MQ=255)
                                                                                                   atggagccgttgctgcattgagtgttgctccttgtccaaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCa                                        <  2:364518‑M2/104‑1 (MQ=255)
                                                                                                         ggACACCAACAGGTAATGCAAAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTTCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAacac                                  <  1:400054/149‑1 (MQ=255)
                                                                                                                tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGGTTTTACCAGGGGGATTTAAc                                                                              <  2:125533‑M2/66‑1 (MQ=255)
                                                                                                                tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                                                              >  1:125533‑M2/33‑98 (MQ=255)
                                                                                                                         ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACACCGTAAt                  >  1:393224‑M2/24‑149 (MQ=255)
                                                                                                                                 cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                        >  1:577163‑M2/16‑71 (MQ=255)
                                                                                                                                 cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                        <  2:577163‑M2/56‑1 (MQ=255)
                                                                                                                                 cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACCACCTAAGAcgctaca          >  2:267756‑M2/16‑146 (MQ=255)
                                                                                                                                       ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc    <  1:82990‑M2/140‑1 (MQ=255)
                                                                                                                                       ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAACTACAACGTAATGCGCTGAACGGTGc    >  1:152307‑M2/10‑149 (MQ=255)
                                                                                                                                        caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                             <  2:483870‑M2/67‑1 (MQ=255)
                                                                                                                                        caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                             >  1:483870‑M2/9‑75 (MQ=255)
                                                                                                                                         aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAg  >  2:483978‑M2/8‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                             
TTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913_3_hsa_tpiA/1209475‑1209760

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.