New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_tpiA	3071319 =	22 (0.880)	4 (0.170)	4/274	2.9	17.0%	noncoding (3071319/4641682 nt)	repeat_region	repeat region
	?	NC_000913_3_vch_tpiA	3071351 =	18 (0.750)					noncoding (3071351/4641682 nt)	repeat_region	repeat region
Rejected: Frequency below/above cutoff threshold.

GCCAGGTTTTACGAAATCAAAAATCTTAGACATGTCTCTTGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTCCCTGCCCTGTGATTTTTTACTTCTTAGCGCGCTCTTCGAGCATCGCTACTGCAGG  >  NC_000913_3_vch_tpiA/3071212‑3071497                                                                                                                                            |                                                                                                                                                   gCCAGGTTTTACGAAGTCAAAAATCTTAGACATGTCTCTTGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCACTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCActctc                                                                                                                                           >  1:161218/1‑149 (MQ=255) gCCAGGTTTTACGAAATCAAAAATCTTAGACATGTCTCTTGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGTATTCActctc                                                                                                                                           >  1:318580/1‑149 (MQ=255)      acagggcagggaaacctgcccttgtttcagcgcgcttttagagcacgcaatttttcatactctaaataattcgagttgcaggaaggcgacaagagagtgaatCCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTc                                                                                                                                      >  1:64319‑M1/103‑149 (MQ=255)           gcagggaaacctgcccttgtttcagcgcgcttttagagcacgcaatttttcatactctaaataattcgagttgcaggaaggcgacaagagagtgaatCCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCa                                                                                                                                                >  2:395444‑M1/98‑134 (MQ=255)           gcagggaaacctgcccttgtttcagcgcgcttttagagcacgcaatttttcatactctaaataattcgagttgcaggaaggcgacaagagagtgaatCCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCa                                                                                                                                                <  1:395444‑M1/37‑1 (MQ=255)                     aaaTCTTAGACATGTCTCTTGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTc                                                                                                                       <  2:129725/149‑1 (MQ=255)                     aaaTCTTAGACATGTCTCTTGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTc                                                                                                                       >  2:120384/1‑149 (MQ=255)                                                         ccTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTaaa                                                                                   >  2:235154/1‑149 (MQ=255)                                                           ttATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTACTCaaaatt                                                                                 >  2:309745/1‑147 (MQ=255)                                                              taCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATCTAAAGTATGAAAAATTGCATACCCAAAAAgcac                                                                              >  1:201496/1‑147 (MQ=14)                                                                   gTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTAGTTCTGTAAACTCCTGGGGATTCACTGTCTTGTCGCCTTCCTGCAACTCGAATTATTTATAGTATGAAAAATTGCGTGCTATAAAAGCGAGCTGa                                                                         >  1:14479/1‑149 (MQ=17)                                                                          ttCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGCATGAAAAATTGCGTGCCCTCACAGCGCCCTGACACACgg                                                                  >  2:268576/1‑149 (MQ=14)                                                                            gaaGTTGTATGTGCTGCGTCTTCGTTCGCTTCACCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTTTCTTGTCGCCTTCCTGCGACTCGTATTATTTAGATTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGc                                                                <  2:85751/148‑1 (MQ=18)                                                                                ttGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGt                                                            >  1:121381/1‑149 (MQ=255)                                                                                    aTGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTcc                                                        >  1:61301/1‑149 (MQ=255)                                                                                        gctgcGTCTGCGTTCGCCCACCCCAGTCACTTACTAACGTAAGCTCCAGGGGATTCACTATCTTGTCGCCTTCCTGCAACCCGAATTATTTAGAGAATGAAAAATGACGGGACCTACCAACACCCTGAAACAAGTGCAGGATTccgtgc                                                    >  2:275837/1‑149 (MQ=255)                                                                                          tgcgtctgcgtTCGCTCACCCCAGTCACTTTTCTAAGCTACCTCGTGGGTATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTccctgccc                                                  <  2:201496/149‑1 (MQ=21)                                                                                             gacaagagagtgaatccccaggagcttacataagtaagtgactggggATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTCCCTGCCCtgt                                               >  1:64319‑M2/48‑149 (MQ=255)                                                                                                       tgaatccccaggagcttacataagtaagtgactggggATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTccctgc                                                    <  1:395444‑M2/97‑1 (MQ=255)                                                                                                       tgaatccccaggagcttacataagtaagtgactggggATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTccctgc                                                    >  2:395444‑M2/38‑134 (MQ=255)                                                                                                             cccaggagcttacataagtaagtgactggggATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGATTATGAAAAATTGCGTGTACTAAAAGCGCGCTGAAACAAGGGAAGGTTTTCCTGCCCCGTGATTTTTTACTACTTa                               >  2:99968‑M2/32‑149 (MQ=255)                                                                                                               CAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTCCCTGCCCTGTGATTTTTTACTTCTTAgc                             <  1:127119‑M1/149‑1 (MQ=255)                                                                                                                        acttaTGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGTATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTCCCTGCCCTGTGATTTTTTACTTCTTAGCGCGATCTTc                    >  2:174091/1‑149 (MQ=255)                                                                                                                          ttaTGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTCCCTGCCCTGTGATTTTTTACTTCTTAGCGCGCTCTTCGa                  <  1:235154/149‑1 (MQ=255)                                                                                                                               tAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTCCCTGCCCTGTGATTTTTTACTTCTTAGCGCGCTCTTCGAGCATc             >  1:419445/1‑149 (MQ=255)                                                                                                                               tAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTCCCTGCCCTGTGATTTTTTACTTCTTAGCGCGCTCTTCGAGCATc             <  2:61301/149‑1 (MQ=255)                                                                                                                                          ggATTCACTCTCTTGTCGCCTTCCTGCAACTCGGCTTATGTGTTTTATGAAAATTTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTCCCTGCCCTGTGATTTTTTACTTCTTAGCGCGCTCTTCGAGCATCGCTACTGCAgg  <  1:120384‑M2/147‑1 (MQ=255)                                                                                                                                          ggATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTCCCTGCCCTGTGATTTTTTACTTCTTAGCGCGCTCTTCGAGCATCGCTACTGCAgg  >  2:326727‑M2/3‑149 (MQ=255)                                                                                                                                          ggATTCACTCTCGTGTCGNCTCGCTGGAACCCGAAGTATTTAGAGTATGAAAAATTGCGTGCTCTAGAAGCGCGCTGAAACAAGGGCAGGTTTCCCTGCCCTGTGATTTTTTACTTCTTAGCGCGCTCTTCGAGCATCGCTACTGCAgg  <  1:268576‑M2/147‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                   GCCAGGTTTTACGAAATCAAAAATCTTAGACATGTCTCTTGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGAAAAATTGCGTGCTCTAAAAGCGCGCTGAAACAAGGGCAGGTTTCCCTGCCCTGTGATTTTTTACTTCTTAGCGCGCTCTTCGAGCATCGCTACTGCAGG  >  NC_000913_3_vch_tpiA/3071212‑3071497

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.