New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 3387786	31 (0.810)	6 (0.160)	6/306	NT	14.8%	coding (1832/3801 nt)	yhdP	DUF3971‑AsmA2 domains protein
	?	CP009273	= 3388496	39 (1.040)					coding (1122/3801 nt)	yhdP	DUF3971‑AsmA2 domains protein

TTACCCGGACCTTTAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/3387720‑3387786 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggCAAGATTAACACTCTGGCGCTGGATATCCCGCTTCAGGCGGCAGACAAGACCCGTTTTCAGGCATCGTGGAGCGATCTGGCCTGGAAGCAATGGAAATTATTACCGGGTGCGGAACACT  <  CP009273/3388496‑3388378                                                                                                                                                                                            TTACCCGGACCTTTAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGCAAGATTAACACTCTGGCGCTGGATATCCCGCTTCAG                                                                                    >  1:320074/1‑104     CCGGACCTTTAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGCAAGATTAACACTCTGGCGCTGGATATCCCGCTTCAGGCGG                                                                                <  2:320068/104‑1                     AGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGCAAGATTAACACTCTGGCGCTGGATATCCCGCTTCA                                                                                     >  2:147808/1‑83                         TCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGCAAGATTAACACTCTGGCGCTGGATATCCCGCTTCAGGCG                                                                                 <  1:147809/83‑1                                 CAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGCAAGATTAACACTCTGGCGCTGGATATCCCGCTTCAGGCGGCAGACAAGACCCGTTTTCAGGCATCGTGGAGCGATCTGGCCTGGAAGCAATGGAAATTATTACCGGGTGCGGAA      >  2:493987/1‑150                                     TTTTCTTTTGAGTAGTCAGGGATCACTGCGGCAAGATTAACACTCTGGCGCTGGATATCCCGCTTCAGGCGGCAGACAAGACCCGTTTTCAGGCATCGTGGAGCGATCTGGCCTGGAAGCAATGGAAATTATTACCGGGTGCGGAACACT  <  1:493994/150‑1                                                                   GCAAGATTAACACTCTGGCGCTGGATATCCCGCTTCAGGCGGCAGACAAGACCCGTTTTCAGGCA                                                         <  1:415802/65‑1                                                                                                                                                                                            TTACCCGGACCTTTAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/3387720‑3387786 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggCAAGATTAACACTCTGGCGCTGGATATCCCGCTTCAGGCGGCAGACAAGACCCGTTTTCAGGCATCGTGGAGCGATCTGGCCTGGAAGCAATGGAAATTATTACCGGGTGCGGAACACT  <  CP009273/3388496‑3388378

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