New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 4542672	16 (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)

GCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGG  >  NC_000913_3_hsa_tpiA/4542834‑4543115                                                                                                                                       |                                                                                                                                                    gCATCCGCCTCATGCCTCTCTCCAATAACCCCCTCATGCAATATAAACATCTATAAATAAAGATTACAATAGAATATTAAGGTAACAAATAAACTGAAAAAGATTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGtt                                                                                                                                       <  1:388140/149‑1 (MQ=255)   aTCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTc                                                                                                                                     >  1:235024/1‑149 (MQ=255)   aTCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTc                                                                                                                                     >  1:428214/1‑149 (MQ=255)             ttcttcTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                                                                                           <  1:275482/149‑1 (MQ=255)                       cAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCACCAAATGTTTCATCTTTTGGGGGAAAAAtgt                                                                                                                 >  1:427653/1‑149 (MQ=255)                       cAAAAACCACCTCATGCAATATAAACATCTAATGATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                                                                                                 >  2:309352/1‑149 (MQ=255)                               tccTCATGCAATATAAACATCTATAAATAAATATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                                                                                                         <  1:161472/148‑1 (MQ=255)                               accTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAGGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAACTGTTTCATCTTTTGGGGCAAAACTGGTCCGTGtt                                                                                                         >  2:175186/1‑149 (MQ=255)                                ccTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattat                                                                                                        >  1:186362‑M2/1‑104 (MQ=255)                                     tGCAATATAAACATCTATACATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTACGCGATGCTTTCATCTATGAGTCAAAATGGCCCCAAATGTTTAATCTTTTGGGGGAACAATTTGAAGTGTTGGCGGt                                                                                                   >  2:207576/1‑149 (MQ=255)                                               aCATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGtt                                                                                         >  2:74728/1‑149 (MQ=255)                                               aCATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTACGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGtt                                                                                         >  2:76438/1‑149 (MQ=255)                                                    tataAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTccaa                                                                                    >  2:419653/1‑149 (MQ=255)                                                     ataaataaAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaa                                                                                   >  2:389400/1‑149 (MQ=255)                                                            aaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccacaattgtcttgtatttatttgtttttttaactttattatcaataagttaaataaaacctcttctttt                                                                            <  2:97038‑M2/149‑74 (MQ=255)                                                              gATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTTCGCGATGCTTTCCTCTATGCATCAAAAtggccccaatggtcttgtatttatttgtttttttaactttattatcaatgacttaaataaaacctcttctctttt                                                                          <  1:387810‑M2/149‑76 (MQ=255)                                                                     aTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatctatcctattttatttttattttt                                                                   <  2:253429‑M2/149‑83 (MQ=255)                                                                     aTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAAc                                                                   >  2:184704/1‑149 (MQ=255)                                                                            atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                                                                                >  2:34789/1‑97 (MQ=255)                                                                            atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                                                                                <  1:34789/97‑1 (MQ=255)                                                                             ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtctttatttgtttttttaactttattatcaatcagttaaataaaaactcttctctttttaattttttgtttat                                                           <  2:286687‑M2/149‑91 (MQ=255)                                                                             ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGccc                                                           >  1:250205/1‑149 (MQ=255)                                                                                     acaaATAACCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAACCGAGTCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAt                                                   >  1:207751/1‑149 (MQ=255)                                                                                     acaaAAAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAt                                                   <  2:183449/149‑1 (MQ=255)                                                                                           aaaCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGGTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGaa                                             <  1:229756/149‑1 (MQ=255)                                                                                                              gcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAc                                                                        <  2:116219/103‑1 (MQ=255)                                                                                                              gcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAc                                                                        >  1:116219/1‑103 (MQ=255)                                                                                                               cgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaaggaaag                         <  2:250205/149‑1 (MQ=255)                                                                                                                                      aaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTgg  <  1:389400/149‑1 (MQ=255)                                                                                                                                       |                                                                                                                                                    GCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGG  >  NC_000913_3_hsa_tpiA/4542834‑4543115

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.