New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 114505	28 (0.670)	5 (0.130)	5/258	4.8	16.8%	intergenic (+18/+17)	guaC/hofC	GMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein
	?	NC_000913_3_hsa_tpiA	= 114511	24 (0.620)					intergenic (+24/+11)	guaC/hofC	GMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTAT  >  NC_000913_3_hsa_tpiA/114357‑114646                                                                                                                                                     |                                                                                                                                              cgcgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTg                                                                                                                                               <  2:508775/149‑1 (MQ=255)  gcgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTgg                                                                                                                                              <  1:734216/149‑1 (MQ=255)        tatTTTGGGCGGCCTGAGTTCAGCTTGTACATACGTTCGGGCTTCACGCCTGAAAGAGCAGACCAAGCGCACCACGTATATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAac                                                                                                                                        <  2:681818/149‑1 (MQ=255)          ttttGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGCAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAacac                                                                                                                                      >  1:653991/1‑149 (MQ=255)              gggCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCa                                                                                                                                  >  1:466056/1‑149 (MQ=255)                   gcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGtt                                                                                                                             >  1:413402/1‑149 (MQ=255)                       gCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCCACGCCAGCGttggga                                                                                                                         <  1:348144/149‑5 (MQ=255)                           tCACCTTGTACATACGTGGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccatc                                                                                                                     <  2:135156/149‑1 (MQ=255)                           gcAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGGACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTCCCACAACCTGTAATCCCCCAACGCTGGCGTGGAGAAACCCGCCAGCGTTatgccatc                                                                                                                     <  2:550465/148‑1 (MQ=255)                               cTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCAc                                                                                                                 <  1:481144/149‑1 (MQ=255)                                 tGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTc                                                                                                               <  1:498203/149‑1 (MQ=255)                                                         cTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTACATGGAAAAtt                                                                                       >  2:661877/1‑149 (MQ=255)                                                                 gCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATa                                                                               <  2:464105/149‑1 (MQ=255)                                                                              accacGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGccaccacc                                                                  >  1:350700/1‑149 (MQ=255)                                                                                      tATTCGTGCGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCCTCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTAcc                                                          <  1:92606/149‑1 (MQ=255)                                                                                           gtgtGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaataa                                                     >  1:393673/1‑149 (MQ=255)                                                                                                cAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATtcct                                                <  2:653991/149‑1 (MQ=255)                                                                                                    aaCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaataat                                                    <  1:170866/141‑1 (MQ=255)                                                                                                    aaCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaataat                                                    >  2:170866/1‑141 (MQ=255)                                                                                                    aaCAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCg                                            <  1:393503/149‑1 (MQ=255)                                                                                                        aGAAACGCGTATCTTCAACAACCTGTAATCTCCCAGCGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCTGCGTACCAATAATTCCTCCCGTTAt                                        <  2:372826/149‑1 (MQ=255)                                                                                                        aGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACAACAGACAGGTCGGGCAGCGCACGGTTCGTATCTCCTAAATGACAGATTGCAAAGGACGTTGCCACCACCAGCGTACCAATAATTCGTCCCGtaag                                        >  1:547092/1‑146 (MQ=255)                                                                                                              ccGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCTTCCCGTTATGATCAg                                  >  1:630242/1‑149 (MQ=255)                                                                                                               cGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGAGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGc                                 >  1:212990/1‑149 (MQ=255)                                                                                                                gCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCa                                >  1:641092/1‑149 (MQ=255)                                                                                                                         acaacCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTACCACCACCAGCGTACCACTCATGCTGCCCACTAGGATCAGCAACGCCgggt                       >  2:672536/1‑147 (MQ=255)                                                                                                                           aacCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTcc                     <  1:566861/149‑1 (MQ=255)                                                                                                                               tGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTa                 >  2:736414/1‑149 (MQ=255)                                                                                                                                     ctcCCAACGCTGGCGTGGAGCAACACGCCCGCGTTATCCCATCCCACTCTTCGCATCCCCTAAATGCAAAAGGGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCATCGCCGGATCCAGTAAGCCTg           >  2:598601/1‑149 (MQ=255)                                                                                                                                              cTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTAt  >  1:714705/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                              CGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTAT  >  NC_000913_3_hsa_tpiA/114357‑114646

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.