New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 4040683	92 (1.890)	11 (0.230)	6/276	5.7	13.7%	intergenic (+47/+223)	rrfA/mobB	5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
	?	NC_000913_3_hsa_tpiA	= 4213144	48 (1.000)					intergenic (+3/‑72)	rrfE/yjaA	5S ribosomal RNA of rrnE operon/stress‑induced protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGG  >  NC_000913_3_hsa_tpiA/4213000‑4213292                                                                                                                                                 |                                                                                                                                                      aGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGGGGGGTCTCCCCAAGCGAGAGTAGGGAACTGCCAGGCATCAAATTa                                                                                                                                                   <  1:685796/149‑1 (MQ=14)       atGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaaaaa                                                                                                                                             >  2:452644/2‑149 (MQ=25)        tGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATTCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAc                                                                                                                                            >  2:131645/1‑149 (MQ=37)          aTAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAccc                                                                                                                                          >  1:392061/1‑149 (MQ=38)          aTAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAccc                                                                                                                                          >  1:833578/1‑149 (MQ=38)           tAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAcccc                                                                                                                                         >  1:505887/1‑149 (MQ=39)            aaaaCAGAATGGGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCg                                                                                                                                        <  1:612825/149‑1 (MQ=255)              aaCAGCATTTGGCTGGAGCGGTTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGt                                                                                                                                      <  2:569915/149‑1 (MQ=14)               aCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTc                                                                                                                                     >  2:120113/1‑149 (MQ=255)                        gCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑cc                                                                                                                           >  1:554752/1‑149 (MQ=255)                          cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCAGCAAATTAGAAAAACCCCGGTCCATAAGG‑CCgg                                                                                                                         >  1:160816/1‑149 (MQ=255)                           tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAGAACCCGGGTCCATAAGG‑CCggg                                                                                                                        >  1:501519/1‑149 (MQ=255)                           tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACCCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                                                        <  2:769393/149‑1 (MQ=255)                                  aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGG‑CCTTATGGACCGGGGtttttc                                                                                                                 <  2:258504/149‑2 (MQ=14)                                   gTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTg                                                                                                                >  1:538880/1‑149 (MQ=255)                                   gTAGCGCGGTGGGCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTg                                                                                                                <  2:377235/149‑1 (MQ=255)                                                  cACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTg                                                                                                 <  1:120113/149‑1 (MQ=255)                                                   aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGc                                                                                                <  2:600035/149‑1 (MQ=255)                                                   aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGc                                                                                                >  1:523133/1‑149 (MQ=255)                                                           ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGaa                                                                                        >  1:613088/1‑149 (MQ=255)                                                              tGCCGAACGCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCCCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAggg                                                                                     <  1:500151/149‑1 (MQ=255)                                                                       cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATg                                                                            <  1:515793/149‑1 (MQ=255)                                                                       cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTGTTTTGCATAGTAAGTATTTGCATGATGAAGGGAACCTCATg                                                                            <  2:392061/149‑1 (MQ=21)                                                                       cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGAGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATg                                                                            <  1:24889/149‑1 (MQ=255)                                                                          aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCCTGATAAATAGTTGCATGATGAAGGGAATCTCATGTCa                                                                         <  2:833578/149‑1 (MQ=21)                                                                           aGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAg                                                                        <  1:616156/149‑1 (MQ=255)                                                                            gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCGGGCATCAAATTAGAAAAACCTCGGTCCATAAGG‑CCGGGGCTTTTTGCATATCAATTCTTTACATGATGAAGGGAATCTCATGttcgt                                                                       >  2:443529/1‑145 (MQ=25)                                                                                    ccGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGtata                                                               >  2:216528/1‑149 (MQ=255)                                                                                     cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTatat                                                              <  1:443247/149‑1 (MQ=255)                                                                                           ggcgggggTAGGGTGGGGTTTCCCCATGCGAGAGTAGGGAACCGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAAt                                                        <  2:498185/143‑1 (MQ=38)                                                                                             cGATGGACGTGGGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTTATGTCAGTTCTGTATATCCAAATtc                                                      <  2:185089/149‑1 (MQ=255)                                                                                               aTGGTCGTGTGGGGTCTCACCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgt                                                    >  2:750241/1‑149 (MQ=255)                                                                                                ttgTAGGGTGGGGGCACCCCATGGAAGAGTAGTGAAATGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                                                   <  2:162038/147‑1 (MQ=255)                                                                                                tGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCATTTCTGTATATCCAAATtcgtc                                                   <  2:762192/149‑1 (MQ=255)                                                                                                      tgtgGGGTCTCCCCAGGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATc                                             <  1:536099/149‑1 (MQ=255)                                                                                                          ggggCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAAt                                         <  1:131645/149‑1 (MQ=255)                                                                                                          gggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTATGTATATCCAAATTCGTCGTAATCAAAt                                         >  2:844042/1‑149 (MQ=255)                                                                                                               tcCCCATGCGAGAGTAGGGAACTGCCAAGCATCAAATTAGCAAAACCACGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCACATTCGTCGTAATCAAATTCCTg                                    >  2:105002/1‑149 (MQ=39)                                                                                                                 cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGtt                                  >  1:253270/1‑149 (MQ=255)                                                                                                                          gagTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCtt                         <  2:523133/149‑1 (MQ=255)                                                                                                                                  aaCTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGt                                                                  <  2:96223/100‑1 (MQ=255)                                                                                                                                  aaCTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGt                                                                  >  1:96223/1‑100 (MQ=255)                                                                                                                                   aCTGCCACGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCaaa                <  1:212838/149‑1 (MQ=255)                                                                                                                                         aGGCATCAAACTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa          <  2:711798/149‑1 (MQ=255)                                                                                                                                           gCACCAAATGAGAAACACCGCCGTACATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAgt        <  1:270990/149‑1 (MQ=255)                                                                                                                                            cATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAgtg       >  1:425159/1‑149 (MQ=255)                                                                                                                                               aaaaTTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTGTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGGCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCa    <  2:263653/148‑1 (MQ=255)                                                                                                                                                 aaTTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTCTTCGCGATCTTGAAAGCCAAGTAGAAGTGTCAgg  >  1:112022/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      AGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGG  >  NC_000913_3_hsa_tpiA/4213000‑4213292

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.