New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 4040683	42 (1.710)	7 (0.290)	4/280	2.9	18.3%	intergenic (+47/+223)	rrfA/mobB	5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
	?	NC_000913_3_hsa_tpiA	= 4213144	21 (0.870)					intergenic (+3/‑72)	rrfE/yjaA	5S ribosomal RNA of rrnE operon/stress‑induced protein
Rejected: Frequency below/above cutoff threshold.

TCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAG  >  NC_000913_3_hsa_tpiA/4213005‑4213291                                                                                                                                            |                                                                                                                                                    tCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacga                                                                                                                                            <  1:117826‑M2/149‑10 (MQ=255)   tGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaa                                                                                                                                          <  1:112083‑M2/149‑12 (MQ=255)     aaaaaacaGAATCTGCCTGGTGGCAGTAGCGCGGTGGCCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaagg                                                                                                                                        >  2:355743‑M2/3‑136 (MQ=255)        aaaCAGTACTTCCCTGGCTGCAGTAGCGCGGGGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCgg                                                                                                                                     <  2:299172/149‑1 (MQ=14)                   gCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAGataaaacgaaaggcccagtcttccgac                                                                                                                          >  1:6841‑M2/1‑122 (MQ=255)                   gCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGcc                                                                                                                          <  1:421958/149‑1 (MQ=255)                     cTGGCGGCAGTAGCGCGGTGGTACCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCgg                                                                                                                        <  1:91226/149‑1 (MQ=255)                                 gcgcGGTGGTCCCACCTGACCCCATGCCGGAATCAGAAGTGGAACGTCGTAGCGCCGGTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGCACTGAGAGGCATCAATTTAGAAAACGCGAGGTGCATAAGGGGGGGGGTTTTTgggg                                                                                                            >  2:99301/1‑146 (MQ=11)                                   gcgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAAGGACGTGTCGGTCATCATATTTGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCatat                                                                                                          <  2:265179/149‑1 (MQ=18)                                      gtggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCaa                                                                                                       >  2:402915/1‑149 (MQ=255)                                           cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAattatt                                                                                                  >  1:130927/1‑149 (MQ=255)                                               ccTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCa                                                                                              <  2:278985/149‑1 (MQ=255)                                                     cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtcttccgactgagccttt                                                                                                                 >  2:154765‑M2/1‑88 (MQ=255)                                                     cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtcttccgactgagccttt                                                                                                                 <  1:154765‑M2/124‑37 (MQ=255)                                                               aCTCAGAAGTGAGAAGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCAGATTAATTATTTGCATGCTGAAGGGAAtctc                                                                              >  2:331876/1‑149 (MQ=25)                                                                cTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCa                                                                             <  1:324309/149‑1 (MQ=255)                                                                                     cgcCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCaaa                                                        <  1:402915/149‑1 (MQ=255)                                                                                         gATGGTAGTGTGGGGTCTACATATGCGAGAGTAGGGAACTGCGAGGCATCAAATTAGAAAAACCACGGTACATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTTATTTGTGAATATGCAAATtcg                                                    >  2:378447/1‑149 (MQ=255)                                                                                          aTGGTTGTGTGGGGTCTCCCAAGGCGAGAGTAGTGACCTGTCAAGCAGCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgt                                                   <  1:295981/149‑1 (MQ=255)                                                                                              tAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTa                                               >  1:337155/1‑149 (MQ=255)                                                                                                       gTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGATCTGTATATCCAAATTCGTCGTCAGCAAAtga                                      >  2:179015/1‑147 (MQ=255)                                                                                                         ctcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCCACTTACt                                    >  2:61182/1‑149 (MQ=255)                                                                                                                 gCGAGCGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGa                            <  1:281902/149‑1 (MQ=255)                                                                                                                     gagTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCtt                        >  2:166090/1‑149 (MQ=255)                                                                                                                        tAGGGAACTGCCAGGCATCAAATTAGAAAAACCGCGGTCCATAAGGACGGGGTTTTTTGCATATCAAGTCCTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCAAAATTAGACGTGATAAAATTTCTATTCGGGATCTTGaa                     >  2:304399/1‑149 (MQ=255)                                                                                                                         aGGTAACTGCCAGGCAACACATCAGAAAACCCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATTTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaaa                    <  1:210750/149‑1 (MQ=255)                                                                                                                                     aacagataataaaacgaaaggcccagtcttccgactgagcctttcgttttatttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgg        >  2:366998‑M2/8‑8 (MQ=255)                                                                                                                                           aaaTTAGAAAAACCCCCGTCCATGAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCAGTATACCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAg  <  1:191925/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                    TCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAG  >  NC_000913_3_hsa_tpiA/4213005‑4213291

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ 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.