New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 454267216 (0.650)5 (0.220) 5/266 2.4 20.4% intergenic (+57/‑425) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913_3_hsa_tpiA = 4542968 24 (1.040)intergenic (+353/‑129) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

TTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGT  >  NC_000913_3_hsa_tpiA/4542528‑4542819
                                                                                                                                                |                                                                                                                                                   
ttATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACATTTGGGGCCAAACt                                                                                                                                                 >  1:253573/1‑149 (MQ=255)
   taCCGCCAGTAAAGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTc                                                                                                                                              <  1:108893/149‑1 (MQ=255)
                   gctCGTTTTGGCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                                                                                                              <  2:246583/149‑1 (MQ=255)
                   gctCGTTGTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                                                                                                              <  2:246580/149‑1 (MQ=255)
                             ccGGATTATGGAAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTAGGGCAAAACTGTCCATATCATACAAAAGTCACGCAtctt                                                                                                                    >  2:49090/1‑149 (MQ=255)
                                gagtcgggggAAAAAACAATCTCAGAAGAGTAACATTAAAAAGAGAAGGAGATCGATTTAACGTAAGGATAAAAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtc                                                                                                                 <  2:381937/143‑1 (MQ=255)
                                      gggAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCaaac                                                                                                           >  2:195879/1‑147 (MQ=255)
                                                       aTAAACAAAAAATTAAAAAGAGAAGAGTTTTTATTTAACTGATTGATAATAAAGTTAAAAAAACAAATAAAGACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                                                                                          <  2:286687‑M1/149‑60 (MQ=255)
                                                             gAAAAATTAAAAAGAGAAGAGGATTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATGATAAATAAGTTACGTATTTTTTATCAATCATAAAAATATTAAAAAACGaaaa                                                                                    >  2:161472/1‑149 (MQ=255)
                                                               aaaaataaaaatAAAATAGGATAGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctat                                                                                  <  2:253429‑M1/149‑68 (MQ=255)
                                                                   ttAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCa                                                                              >  2:407997/1‑149 (MQ=255)
                                                                      aaaaGAGAAGAGGTTTTATTTAAGTCATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACCATTGGGGCCAttttgatgcatagaggaaagcatcgcgaacaaactttttcagtttatttgttggcttaatattctattgttatc                                                                           <  1:387810‑M1/149‑75 (MQ=255)
                                                                        aaaagaagaggttTTATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGTGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctt                                                                         <  2:97038‑M1/146‑77 (MQ=255)
                                                                                                   gATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAtt                                              <  1:407997/149‑1 (MQ=255)
                                                                                                    ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgagg                                             >  1:186362‑M1/1‑45 (MQ=255)
                                                                                                                            aataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTCTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCg                     >  2:11966/1‑149 (MQ=255)
                                                                                                                             ataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaat                                        >  1:423574/1‑129 (MQ=255)
                                                                                                                             ataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaat                                        <  2:423574/129‑1 (MQ=255)
                                                                                                                                  aGACAACTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTatat               <  1:40859/149‑1 (MQ=255)
                                                                                                                                              ccAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACg   >  2:332817/1‑149 (MQ=255)
                                                                                                                                               cAAACTGTCCATATCATAAATAAGTTACGTATTTTTTTTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGt  <  1:54839/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGT  >  NC_000913_3_hsa_tpiA/4542528‑4542819

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.