New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 3782597	21 (0.880)	4 (0.180)	4/270	2.6	15.7%	intergenic (‑35/+45)	cysE/gpsA	serine acetyltransferase/glycerol‑3‑phosphate dehydrogenase (NAD+)
	?	NC_000913_3_hsa_tpiA	= 3782604	23 (1.020)					intergenic (‑42/+38)	cysE/gpsA	serine acetyltransferase/glycerol‑3‑phosphate dehydrogenase (NAD+)
Rejected: Frequency below/above cutoff threshold.

AACTGGCCAGCATTGGCTCACAGTCCGCCAGCGTTCTGGCTTCGGCTTTAATATTGTTCCAGACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCA  >  NC_000913_3_hsa_tpiA/3782475‑3782751                                                                                                                                  |                                                                                                                                                    aaCTGGCCAGCATTGGCTCACAGTCCGCCAGCGTTCTGGCTTCGGCTTTAATATTGTTCCAGACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGt                                                                                                                                  <  2:407460/149‑1 (MQ=255)  aCTGGCCAGCATTGGCTCACAGTCCGCCAGCGTTCTGGCTTCGGCTTTAATATTGTTCCAGACAATTTCCAGTTCTTCACACGACATTGCTTACGCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTa                                                                                                                                 >  1:220177/1‑149 (MQ=255)        cAGCATTGGCTCACAGTCCGCCAGCGTTCTGGCTTCGGCTTTAATATTGTTCCAGACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTTCGCGGGCCGGGTCATGGCGGGACCAAAg                                                                                                                           >  1:293961/1‑149 (MQ=255)               ggCTCACAGTCCGCCAGCGTTCTGGCTTCGGCTTTAACATTGTTCCAGACAATTTCCAATTCCTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCt                                                                                                                    <  2:151869/149‑1 (MQ=255)                            ccAGCGTTCTGGCTTCGGCTTTAATATTGTTCCAGACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGgct                                                                                                       >  1:251289/1‑149 (MQ=255)                                  ttCTGGCTTCGGCTGTAATATGGTTCCAGACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTATGCGACGGCCGGGTTATATCCGTAACAAAGCTTCCCTGGGGCTAGTGGCTGCTgcg                                                                                                 >  2:87756/1‑149 (MQ=255)                                             gCTTTAATATTGTTCCAGACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCCCGTCCTTGc                                                                                      >  2:168759/1‑149 (MQ=255)                                               tttAATATTGTTCCAGACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGCCCCGTTCTACGCGGGCCGGGTCATAGCGGTCACACAGGTTCCCTAGGGTCAGTGGCTGCTGCGCTCGCCCTTGctg                                                                                    >  2:84380/1‑147 (MQ=255)                                                    tatTGTTCCAGACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCAct                                                                               >  2:236390/1‑148 (MQ=255)                                                             aGACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTTCACGACCTAGTat                                                                      >  2:215634/1‑148 (MQ=255)                                                             aGACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGCCCTAGTaa                                                                      >  1:145322/1‑149 (MQ=255)                                                              gaaaaTTTCCAGTTCTTCACACGACATTGTTTACTCCACACGATGAGATAATGACCGGCCCGGTCTGCGCGGCCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGGTATGCGCTGCTGCGCTCTTCCTGCCGTGCACGCCCTCGTaaa                                                                     >  2:178929/4‑149 (MQ=255)                                                              gACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTaaa                                                                     <  1:342135/149‑1 (MQ=255)                                                                    ttCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAAt                                                               <  1:385490/149‑1 (MQ=255)                                                                               cacaCGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGGCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAATGCTGCCTcccg                                                    >  2:224613/1‑149 (MQ=255)                                                                                  acGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAATGCTGCCTcgcgcgc                                                 >  2:33266/1‑149 (MQ=255)                                                                                           cTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGcgggccggtcattatctcatcgtgtggagtaagcaatgtcgtgtgaagaactggaaattgtctggaacaatattaaagccgaagccagaacgctggcggactgggagcc                                        <  1:352099‑M2/149‑110 (MQ=255)                                                                                             tACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGc                                      >  1:196738/1‑149 (MQ=255)                                                                                                  cacaCGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAatat                                 >  1:48123/1‑149 (MQ=255)                                                                                                                  gACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAAttt                 >  1:135406/1‑149 (MQ=255)                                                                                                                   aCCGGCCCGTTCTGCGcgggccggtcattatctcatcgtgtggagtaagcaatgtcgtgtgaa                                                                                                      >  2:180269‑M2/1‑16 (MQ=255)                                                                                                                   aCCGGCCCGTTCTGCGcgggccggtcattatctcatcgtgtggagtaagcaatgtcgtgtgaa                                                                                                      <  1:180269‑M2/63‑48 (MQ=255)                                                                                                                    ccGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAAtgct                                                            >  2:41474/1‑104 (MQ=255)                                                                                                                    ccGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAAtgct                                                            <  1:41474/104‑1 (MQ=255)                                                                                                                            tggcgcgcGGGGCGTGTCAGAGGCGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTAt       <  2:396503/145‑1 (MQ=255)                                                                                                                                 ggcgGGCCGGGCCATAGCGGTAACAAAGGATCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCa  <  2:330573/148‑1 (MQ=255)                                                                                                                                 cgcgggccggtcattatctcatcgtgtggagtaagcaatgtcgtgtgaagaactggaaattgtctggaacaatattaaagccgaagccagaacgctggcggactgtgagccaatgctggccagtttttaccacgcgacgctactcaagc  >  2:352099‑M2/1‑2 (MQ=255)                                                                                                                                  |                                                                                                                                                    AACTGGCCAGCATTGGCTCACAGTCCGCCAGCGTTCTGGCTTCGGCTTTAATATTGTTCCAGACAATTTCCAGTTCTTCACACGACATTGCTTACTCCACACGATGAGATAATGACCGGCCCGTTCTGCGCGGGCCGGGTCATAGCGGTAACAAAGGTTCCCTGGGGTTAGTGGCTGCTGCGCTCGTCCTTGCGTGCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCA  >  NC_000913_3_hsa_tpiA/3782475‑3782751

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