New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA 1207790 =10 (0.260)27 (0.790) 24/248 0.5 76.2% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA 1209619 = 8 (0.230)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGC  >  NC_000913_3_hsa_tpiA/1209477‑1209758
                                                                                                                                              |                                                                                                                                           
gcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                                                                                       <  2:448404‑M2/7‑1 (MQ=255)
gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                                                                       >  2:430100/1‑149 (MQ=255)
              ggcattccacccacactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                         <  2:298129‑M2/21‑1 (MQ=255)
              ggCCAGTGGCGGGTATTCGGCAGCAGAAAGAGCGGCACCATTGCATTTGAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt                                                                                                                         <  2:267188/148‑1 (MQ=255)
                 attccacccccactgaatgattactgattgctcaataccaccactcctaaccggtatagctatgaaaccggtgagtcctaatgtagccgttgctgcatttggttttgctccttctccgcaaccaaCGTCTATGAAAGTGGAGAAATAAc                                                                                                                      >  1:46473‑M2/126‑149 (MQ=255)
                              tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                         >  1:153892‑M2/113‑149 (MQ=255)
                                       caaaaaGCCGCACCACTGCATTTCAGCCAGCCTGTTGGCGGAGTTGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                <  1:147361/148‑1 (MQ=255)
                                             tgctcaataccaccaatcataacctgtatagctatgaaaccggtgcgtcctaatgtagccgtttctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                                                          <  2:521716‑M2/52‑1 (MQ=255)
                                               ggCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                        >  2:153578/1‑149 (MQ=255)
                                                tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                       >  2:156140‑M2/95‑149 (MQ=255)
                                                tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAACAAATAACAAGCAAATGGCATCATTTCTGCTTTTACCAg                                                                                       >  2:293566‑M2/95‑149 (MQ=255)
                                                  aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                     >  1:416616‑M2/93‑149 (MQ=255)
                                                   ataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                    >  2:329215‑M2/92‑149 (MQ=255)
                                                                   ccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGGAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                                                                    >  1:240724‑M2/76‑149 (MQ=255)
                                                                      gtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                 <  1:485427‑M2/77‑1 (MQ=255)
                                                                      cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                 <  1:153578/149‑1 (MQ=255)
                                                                        gTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                               <  2:397432/149‑1 (MQ=255)
                                                                              ggAGTGGCTGAAGGCCACGGAACCGGTACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                         <  2:299925/149‑1 (MQ=255)
                                                                              atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                         <  1:429661‑M2/85‑1 (MQ=255)
                                                                                  aaccggtgcgtcctaatgtagccgttggtgcattgagttttgctccttctcccaaaccaaCGGTTATGAAAATGAAGAAATAACAAGCAAGTGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCAATGTACGCGt                                                     >  1:416609‑M2/61‑149 (MQ=255)
                                                                                    tcggtgcggcgtaaggtagccgttggtgcattaaggtttgatgcttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGGAt                                                   <  2:46473‑M2/91‑1 (MQ=255)
                                                                                       gtgcgtcctaatgtcgccgttggtgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAGATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                <  2:668766‑M2/94‑1 (MQ=255)
                                                                                       gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                >  1:429730‑M2/56‑149 (MQ=255)
                                                                                        tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccccCATTTATTTATAGGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                               <  1:329215‑M2/95‑1 (MQ=255)
                                                                                        tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGGAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                               <  2:240724‑M2/95‑1 (MQ=255)
                                                                                           gtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                                                      <  2:121839‑M2/40‑1 (MQ=255)
                                                                                           gtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                                                      >  1:121839‑M2/52‑91 (MQ=255)
                                                                                            ccACCGAACCGCGACACCAAAAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACGAGCAAATGGCATCATTCCTGCTTTTACCATGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                           <  2:510681/149‑1 (MQ=255)
                                                                                                    tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                   <  1:247874‑M2/107‑1 (MQ=255)
                                                                                                        cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAAcaca                               <  2:537057‑M2/111‑1 (MQ=255)
                                                                                                                  ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATGATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACATATTACCAGAACACAGATTTACAAc                     >  2:296744‑M2/29‑149 (MQ=255)
                                                                                                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa                 <  1:273867‑M2/125‑1 (MQ=255)
                                                                                                                       ttttgctccttctcccaaaccaaCGTTTATGAATAGAAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt                >  1:335808‑M2/24‑149 (MQ=255)
                                                                                                                       ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGCACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt                >  1:466990‑M2/24‑149 (MQ=255)
                                                                                                                            ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt           <  1:536223‑M2/131‑1 (MQ=255)
                                                                                                                                ttctcccaaaccaaCGTTTATGAAACTGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTGTACCAGGGAGATTTAACATGCTTATTGGCTATGTACGAGTATCAACAAATGACCAGAACACAGAACTACAACCTAATGCGCTgaag       >  2:656681‑M2/15‑148 (MQ=255)
                                                                                                                                     ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAAAACAGGTCGAGATGACAGATCTACAACGTAATGCGCTGAactgtcc  >  1:358866‑M2/10‑147 (MQ=255)
                                                                                                                                     ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCATTTACCAGGGGTATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc  >  1:588720‑M2/10‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                           
GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGC  >  NC_000913_3_hsa_tpiA/1209477‑1209758

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.